Genomic surveillance reveals multiple introductions of SARS-CoV-2 into Northern California

Investigating the cause of a 2021 winter wave of COVID-19 in a border region in eastern Germany: a mixed-methods study, August to November 2021

It is so far unclear how the COVID-19 winter waves started and what should be done to prevent possible future waves. In this study, we deciphered the dynamic course of a winter wave in 2021 in Saxony, a state in Eastern Germany neighbouring the Czech Republic and Poland. The study was carried out through the integration of multiple virus genomic epidemiology approaches to track transmission chains, identify emerging variants and investigate dynamic changes in transmission clusters. For identified local variants of interest, functional evaluations were performed. Multiple long-lasting community transmission clusters have been identified acting as driving force for the winter wave 2021. Analysis of the dynamic courses of two representative clusters indicated a similar transmission pattern. However, the transmission cluster caused by a locally occurring new Delta variant AY.36.1 showed a distinct transmission pattern, and functional analyses revealed a replication advantage of it. This study indicated that long-lasting community transmission clusters starting since early autumn caused by imported or locally occurring variants all contributed to the development of the 2021 winter wave. The information we achieved might help future pandemic prevention.

Sequencing analysis of SARS-CoV-2 cases in Slovenian long-term care facilities to support outbreak control

Introduction

Residents of long-term care facilities (LTCFs) are at high risk of morbidity and mortality due to COVID-19, especially when new variants of concern (VOC) emerge. To provide intradisciplinary data in order to tailor public health interventions during future epidemics, available epidemiologic and genomic data from Slovenian LTCFs during the initial phases of the COVID-19 pandemic was analyzed.

Methods

The first part of the study included SARS-CoV-2 reverse-transcription Real-Time PCR (rtRT-PCR) positive LTCF residents, from 21 facilities with COVID-19 outbreaks occurring in October 2020. The second part of the study included SARS-CoV-2 rtRT-PCR positive LTCF residents and staff between January and April 2021, when VOC Alpha emerged in Slovenia. Next-generation sequencing (NGS) was used to acquire SARS-CoV-2 genomes, and lineage determination. In-depth phylogenetic and mutational profile analysis were performed and coupled with available field epidemiological data to assess the dynamics of SARS-CoV-2 introduction and transmission.

Results

370/498 SARS-CoV-2 positive residents as well as 558/699 SARS-CoV-2 positive residents and 301/358 staff were successfully sequenced in the first and second part of the study, respectively. In October 2020, COVID-19 outbreaks in the 21 LTCFs were caused by intra-facility transmission as well as multiple independent SARS-CoV-2 introductions. The Alpha variant was confirmed in the first LTCF resident approximately 1.5 months after the first Alpha case was identified in Slovenia. The data also showed a slower replacement of existing variants by Alpha in residents compared to staff and the general population.

Discussion

Multiple SARS CoV-2 introductions as well as intra-facility spreading impacted disease transmission in Slovenian LTCFs. Timely implementation of control measures aimed at limiting new introductions while controlling in-facility transmission are of paramount importance, especially as new VOCs emerge. Sequencing, in conjunction with epidemiological data, can facilitate the determination of the need for future improvements in control measures to protect LTCF residents from COVID-19 or other respiratory infections.

Deciphering male influence in gynogenetic Pengze crucian carp (Carassius auratus var. pengsenensis): insights from Nanopore sequencing of structural variations

In this study, we investigate gynogenetic reproduction in Pengze Crucian Carp (Carassius auratus var. pengsenensis) using third-generation Nanopore sequencing to uncover structural variations (SVs) in offspring. Our objective was to understand the role of male genetic material in gynogenesis by examining the genomes of both parents and their offspring. We discovered a notable number of male-specific structural variations (MSSVs): 1,195 to 1,709 MSSVs in homologous offspring, accounting for approximately 0.52%-0.60% of their detected SVs, and 236 to 350 MSSVs in heterologous offspring, making up about 0.10%-0.13%. These results highlight the significant influence of male genetic material on the genetic composition of offspring, particularly in homologous pairs, challenging the traditional view of asexual reproduction. The gene annotation of MSSVs revealed their presence in critical gene regions, indicating potential functional impacts. Specifically, we found 5 MSSVs in the exonic regions of protein-coding genes in homologous offspring, suggesting possible direct effects on protein structure and function. Validation of an MSSV in the exonic region of the polyunsaturated fatty acid 5-lipoxygenase gene confirmed male genetic material transmission in some offspring. This study underscores the importance of further research on the genetic diversity and gynogenesis mechanisms, providing valuable insights for reproductive biology, aquaculture, and fostering innovation in biological research and aquaculture practices.

Viral genome sequencing to decipher in-hospital SARS-CoV-2 transmission events

The SARS-CoV-2 pandemic has highlighted the need to better define in-hospital transmissions, a need that extends to all other common infectious diseases encountered in clinical settings. To evaluate how whole viral genome sequencing can contribute to deciphering nosocomial SARS-CoV-2 transmission 926 SARS-CoV-2 viral genomes from 622 staff members and patients were collected between February 2020 and January 2021 at a university hospital in Munich, Germany, and analysed along with the place of work, duration of hospital stay, and ward transfers. Bioinformatically defined transmission clusters inferred from viral genome sequencing were compared to those inferred from interview-based contact tracing. An additional dataset collected at the same time at another university hospital in the same city was used to account for multiple independent introductions. Clustering analysis of 619 viral genomes generated 19 clusters ranging from 3 to 31 individuals. Sequencing-based transmission clusters showed little overlap with those based on contact tracing data. The viral genomes were significantly more closely related to each other than comparable genomes collected simultaneously at other hospitals in the same city (n = 829), suggesting nosocomial transmission. Longitudinal sampling from individual patients suggested possible cross-infection events during the hospital stay in 19.2% of individuals (14 of 73 individuals). Clustering analysis of SARS-CoV-2 whole genome sequences can reveal cryptic transmission events missed by classical, interview-based contact tracing, helping to decipher in-hospital transmissions. These results, in line with other studies, advocate for viral genome sequencing as a pathogen transmission surveillance tool in hospitals.

Integration of individualized and population-level molecular epidemiology data to model COVID-19 outcomes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with enhanced transmissibility and immune escape have emerged periodically throughout the coronavirus disease 2019 (COVID-19) pandemic, but the impact of these variants on disease severity has remained unclear. In this single-center, retrospective cohort study, we examined the association between SARS-CoV-2 clade and patient outcome over a two-year period in Chicago, Illinois. Between March 2020 and March 2022, 14,252 residual diagnostic specimens were collected from SARS-CoV-2-positive inpatients and outpatients alongside linked clinical and demographic metadata, of which 2,114 were processed for viral whole-genome sequencing. When controlling for patient demographics and vaccination status, several viral clades were associated with risk for hospitalization, but this association was negated by the inclusion of population-level confounders, including case count, sampling bias, and shifting standards of care. These data highlight the importance of integrating non-virological factors into disease severity and outcome models for the accurate assessment of patient risk.

Genomic epidemiology reveals the dominance of Hennepin County in the transmission of SARS-CoV-2 in Minnesota from 2020 to 2022

IMPORTANCE

We analyzed over 22,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes of patient samples tested at Mayo Clinic Laboratories during a 2-year period in the COVID-19 pandemic, which included Alpha, Delta, and Omicron variants of concern to examine the roles and relationships of Minnesota virus transmission. We found that Hennepin County, the most populous county, drove the transmission of SARS-CoV-2 viruses in the state after including the formation of earlier clades including 20A, 20C, and 20G, as well as variants of concern Alpha and Delta. We also found that Hennepin County was the source for most of the county-to-county introductions after an initial predicted introduction with the virus in early 2020 from an international source, while other counties acted as transmission "sinks." In addition, major policies, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, did not appear to have an impact on virus diversity across individual counties.

Motifs in SARS-CoV-2 evolution

We present a novel framework enhancing the prediction of whether novel lineage poses the threat of eventually dominating the viral population. The framework is based purely on genomic sequence data, without requiring prior established biological analysis. Its building blocks are sets of coevolving sites in the alignment (motifs), identified via coevolutionary signals. The collection of such motifs forms a relational structure over the polymorphic sites. Motifs are constructed using distances quantifying the coevolutionary coupling of pairs and manifest as coevolving clusters of sites. We present an approach to genomic surveillance based on this notion of relational structure. Our system will issue an alert regarding a lineage, based on its contribution to drastic changes in the relational structure. We then conduct a comprehensive retrospective analysis of the COVID-19 pandemic based on SARS-CoV-2 genomic sequence data in GISAID from October 2020 to September 2022, across 21 lineages and 27 countries with weekly resolution. We investigate the performance of this surveillance system in terms of its accuracy, timeliness, and robustness. Lastly, we study how well each lineage is classified by such a system.

Online Phylogenetics with matOptimize Produces Equivalent Trees and is Dramatically More Efficient for Large SARS-CoV-2 Phylogenies than de novo and Maximum-Likelihood Implementations

Phylogenetics has been foundational to SARS-CoV-2 research and public health policy, assisting in genomic surveillance, contact tracing, and assessing emergence and spread of new variants. However, phylogenetic analyses of SARS-CoV-2 have often relied on tools designed for de novo phylogenetic inference, in which all data are collected before any analysis is performed and the phylogeny is inferred once from scratch. SARS-CoV-2 data sets do not fit this mold. There are currently over 14 million sequenced SARS-CoV-2 genomes in online databases, with tens of thousands of new genomes added every day. Continuous data collection, combined with the public health relevance of SARS-CoV-2, invites an "online" approach to phylogenetics, in which new samples are added to existing phylogenetic trees every day. The extremely dense sampling of SARS-CoV-2 genomes also invites a comparison between likelihood and parsimony approaches to phylogenetic inference. Maximum likelihood (ML) and pseudo-ML methods may be more accurate when there are multiple changes at a single site on a single branch, but this accuracy comes at a large computational cost, and the dense sampling of SARS-CoV-2 genomes means that these instances will be extremely rare because each internal branch is expected to be extremely short. Therefore, it may be that approaches based on maximum parsimony (MP) are sufficiently accurate for reconstructing phylogenies of SARS-CoV-2, and their simplicity means that they can be applied to much larger data sets. Here, we evaluate the performance of de novo and online phylogenetic approaches, as well as ML, pseudo-ML, and MP frameworks for inferring large and dense SARS-CoV-2 phylogenies. Overall, we find that online phylogenetics produces similar phylogenetic trees to de novo analyses for SARS-CoV-2, and that MP optimization with UShER and matOptimize produces equivalent SARS-CoV-2 phylogenies to some of the most popular ML and pseudo-ML inference tools. MP optimization with UShER and matOptimize is thousands of times faster than presently available implementations of ML and online phylogenetics is faster than de novo inference. Our results therefore suggest that parsimony-based methods like UShER and matOptimize represent an accurate and more practical alternative to established ML implementations for large SARS-CoV-2 phylogenies and could be successfully applied to other similar data sets with particularly dense sampling and short branch lengths.

Implementation of California COVIDNet - a multi-sector collaboration for statewide SARS-CoV-2 genomic surveillance

Introduction

The SARS-CoV-2 pandemic represented a formidable scientific and technological challenge to public health due to its rapid spread and evolution. To meet these challenges and to characterize the virus over time, the State of California established the California SARS-CoV-2 Whole Genome Sequencing (WGS) Initiative, or "California COVIDNet". This initiative constituted an unprecedented multi-sector collaborative effort to achieve large-scale genomic surveillance of SARS-CoV-2 across California to monitor the spread of variants within the state, to detect new and emerging variants, and to characterize outbreaks in congregate, workplace, and other settings.

Methods

California COVIDNet consists of 50 laboratory partners that include public health laboratories, private clinical diagnostic laboratories, and academic sequencing facilities as well as expert advisors, scientists, consultants, and contractors. Data management, sample sourcing and processing, and computational infrastructure were major challenges that had to be resolved in the midst of the pandemic chaos in order to conduct SARS-CoV-2 genomic surveillance. Data management, storage, and analytics needs were addressed with both conventional database applications and newer cloud-based data solutions, which also fulfilled computational requirements.

Results

Representative and randomly selected samples were sourced from state-sponsored community testing sites. Since March of 2021, California COVIDNet partners have contributed more than 450,000 SARS-CoV-2 genomes sequenced from remnant samples from both molecular and antigen tests. Combined with genomes from CDC-contracted WGS labs, there are currently nearly 800,000 genomes from all 61 local health jurisdictions (LHJs) in California in the COVIDNet sequence database. More than 5% of all reported positive tests in the state have been sequenced, with similar rates of sequencing across 5 major geographic regions in the state.

Discussion

Implementation of California COVIDNet revealed challenges and limitations in the public health system. These were overcome by engaging in novel partnerships that established a successful genomic surveillance program which provided valuable data to inform the COVID-19 public health response in California. Significantly, California COVIDNet has provided a foundational data framework and computational infrastructure needed to respond to future public health crises.

A fast and accurate method for SARS-CoV-2 genomic tracing

To contain infectious diseases, it is crucial to determine the origin and transmission routes of the pathogen, as well as how the virus evolves. With the development of genome sequencing technology, genome epidemiology has emerged as a powerful approach for investigating the source and transmission of pathogens. In this study, we first presented the rationale for genomic tracing of SARS-CoV-2 and the challenges we currently face. Identifying the most genetically similar reference sequence to the query sequence is a critical step in genome tracing, typically achieved using either a phylogenetic tree or a sequence similarity search. However, these methods become inefficient or computationally prohibitive when dealing with tens of millions of sequences in the reference database, as we encountered during the COVID-19 pandemic. To address this challenge, we developed a novel genomic tracing algorithm capable of processing 6 million SARS-CoV-2 sequences in less than a minute. Instead of constructing a giant phylogenetic tree, we devised a weighted scoring system based on mutation characteristics to quantify sequences similarity. The developed method demonstrated superior performance compared to previous methods. Additionally, an online platform was developed to facilitate genomic tracing and visualization of the spatiotemporal distribution of sequences. The method will be a valuable addition to standard epidemiological investigations, enabling more efficient genomic tracing. Furthermore, the computational framework can be easily adapted to other pathogens, paving the way for routine genomic tracing of infectious diseases.

Genomic epidemiology of SARS-CoV-2 from Uttar Pradesh, India

The various strains and mutations of SARS-CoV-2 have been tracked using several forms of genomic classification systems. The present study reports high-throughput sequencing and analysis of 99 SARS-CoV-2 specimens from Western Uttar Pradesh using sequences obtained from the GISAID database, followed by phylogeny and clade classification. Phylogenetic analysis revealed that Omicron lineages BA-2-like (55.55%) followed by Delta lineage-B.1.617.2 (45.5%) were predominantly circulating in this area Signature substitution at positions S: N501Y, S: D614G, S: T478K, S: K417N, S: E484A, S: P681H, and S: S477N were commonly detected in the Omicron variant-BA-2-like, however S: D614G, S: L452R, S: P681R and S: D950N were confined to Delta variant-B.1.617.2. We have also identified three escape variants in the S gene at codon position 19 (T19I/R), 484 (E484A/Q), and 681 (P681R/H) during the fourth and fifth waves in India. Based on the phylogenetic diversification studies and similar changes in other lineages, our analysis revealed indications of convergent evolution as the virus adjusts to the shifting immunological profile of its human host. To the best of our knowledge, this study is an approach to comprehensively map the circulating SARS-CoV-2 strains from Western Uttar Pradesh using an integrated approach of whole genome sequencing and phylogenetic analysis. These findings will be extremely valuable in developing a structured approach toward pandemic preparedness and evidence-based intervention plans in the future.

Population genomics of diarrheagenic Escherichia coli uncovers high connectivity between urban and rural communities in Ecuador

Human movement may be an important driver of transmission dynamics for enteric pathogens but has largely been underappreciated except for international 'travelers' diarrhea or cholera. Phylodynamic methods, which combine genomic and epidemiological data, are used to examine rates and dynamics of disease matching underlying evolutionary history and biogeographic distributions, but these methods often are not applied to enteric bacterial pathogens. We used phylodynamics to explore the phylogeographic and evolutionary patterns of diarrheagenic E. coli in northern Ecuador to investigate the role of human travel in the geographic distribution of strains across the country. Using whole genome sequences of diarrheagenic E. coli isolates, we built a core genome phylogeny, reconstructed discrete ancestral states across urban and rural sites, and estimated migration rates between E. coli populations. We found minimal structuring based on site locations, urban vs. rural locality, pathotype, or clinical status. Ancestral states of phylogenomic nodes and tips were inferred to have 51% urban ancestry and 49% rural ancestry. Lack of structuring by location or pathotype E. coli isolates imply highly connected communities and extensive sharing of genomic characteristics across isolates. Using an approximate structured coalescent model, we estimated rates of migration among circulating isolates were 6.7 times larger for urban towards rural populations compared to rural towards urban populations. This suggests increased inferred migration rates of diarrheagenic E. coli from urban populations towards rural populations. Our results indicate that investments in water and sanitation prevention in urban areas could limit the spread of enteric bacterial pathogens among rural populations.

Tracking SARS-CoV-2 introductions in Mozambique using pandemic-scale phylogenies: a retrospective observational study

BACKGROUND

From the start of the SARS-CoV-2 outbreak, global sequencing efforts have generated an unprecedented amount of genomic data. Nonetheless, unequal sampling between high-income and low-income countries hinders the implementation of genomic surveillance systems at the global and local level. Filling the knowledge gaps of genomic information and understanding pandemic dynamics in low-income countries is essential for public health decision making and to prepare for future pandemics. In this context, we aimed to discover the timing and origin of SARS-CoV-2 variant introductions in Mozambique, taking advantage of pandemic-scale phylogenies.

METHODS

We did a retrospective, observational study in southern Mozambique. Patients from Manhiça presenting with respiratory symptoms were recruited, and those enrolled in clinical trials were excluded. Data were included from three sources: (1) a prospective hospital-based surveillance study (MozCOVID), recruiting patients living in Manhiça, attending the Manhiça district hospital, and fulfilling the criteria of suspected COVID-19 case according to WHO; (2) symptomatic and asymptomatic individuals with SARS-CoV-2 infection recruited by the National Surveillance system; and (3) sequences from SARS-CoV-2-infected Mozambican cases deposited on the Global Initiative on Sharing Avian Influenza Data database. Positive samples amenable for sequencing were analysed. We used Ultrafast Sample placement on Existing tRees to understand the dynamics of beta and delta waves, using available genomic data. This tool can reconstruct a phylogeny with millions of sequences by efficient sample placement in a tree. We reconstructed a phylogeny (~7·6 million sequences) adding new and publicly available beta and delta sequences.

FINDINGS

A total of 5793 patients were recruited between Nov 1, 2020, and Aug 31, 2021. During this time, 133 328 COVID-19 cases were reported in Mozambique. 280 good quality new SARS-CoV-2 sequences were obtained after the inclusion criteria were applied and an additional 652 beta (B.1.351) and delta (B.1.617.2) public sequences were included from Mozambique. We evaluated 373 beta and 559 delta sequences. We identified 187 beta introductions (including 295 sequences), divided in 42 transmission groups and 145 unique introductions, mostly from South Africa, between August, 2020 and July, 2021. For delta, we identified 220 introductions (including 494 sequences), with 49 transmission groups and 171 unique introductions, mostly from the UK, India, and South Africa, between April and November, 2021.

INTERPRETATION

The timing and origin of introductions suggests that movement restrictions effectively avoided introductions from non-African countries, but not from surrounding countries. Our results raise questions about the imbalance between the consequences of restrictions and health benefits. This new understanding of pandemic dynamics in Mozambique can be used to inform public health interventions to control the spread of new variants.

FUNDING

European and Developing Countries Clinical Trials, European Research Council, Bill & Melinda Gates Foundation, and Agència de Gestió d'Ajuts Universitaris i de Recerca.

Genetic determination of regional connectivity in modelling the spread of COVID-19 outbreak for more efficient mitigation strategies

For the COVID-19 pandemic, viral transmission has been documented in many historical and geographical contexts. Nevertheless, few studies have explicitly modeled the spatiotemporal flow based on genetic sequences, to develop mitigation strategies. Additionally, thousands of SARS-CoV-2 genomes have been sequenced with associated records, potentially providing a rich source for such spatiotemporal analysis, an unprecedented amount during a single outbreak. Here, in a case study of seven states, we model the first wave of the outbreak by determining regional connectivity from phylogenetic sequence information (i.e. "genetic connectivity"), in addition to traditional epidemiologic and demographic parameters. Our study shows nearly all of the initial outbreak can be traced to a few lineages, rather than disconnected outbreaks, indicative of a mostly continuous initial viral flow. While the geographic distance from hotspots is initially important in the modeling, genetic connectivity becomes increasingly significant later in the first wave. Moreover, our model predicts that isolated local strategies (e.g. relying on herd immunity) can negatively impact neighboring regions, suggesting more efficient mitigation is possible with unified, cross-border interventions. Finally, our results suggest that a few targeted interventions based on connectivity can have an effect similar to that of an overall lockdown. They also suggest that while successful lockdowns are very effective in mitigating an outbreak, less disciplined lockdowns quickly decrease in effectiveness. Our study provides a framework for combining phylodynamic and computational methods to identify targeted interventions.

Genomic epidemiology reveals the dominance of Hennepin County in transmission of SARS-CoV-2 in Minnesota from 2020-2022

SARS-CoV-2 has had an unprecedented impact on human health and highlights the need for genomic epidemiology studies to increase our understanding of virus evolution and spread, and to inform policy decisions. We sequenced viral genomes from over 22,000 patient samples tested at Mayo Clinic Laboratories between 2020-2022 and use Bayesian phylodynamics to describe county and regional spread in Minnesota. The earliest introduction into Minnesota was to Hennepin County from a domestic source around January 22, 2020; six weeks before the first confirmed case in the state. This led to the virus spreading to Northern Minnesota, and eventually, the rest of the state. International introductions were most abundant in Hennepin (home to the Minneapolis/St. Paul International (MSP) airport) totaling 45 (out of 107) over the two-year period. Southern Minnesota counties were most common for domestic introductions with 19 (out of 64), potentially driven by bordering states such as Iowa and Wisconsin as well as Illinois which is nearby. Hennepin also was, by far, the most dominant source of in-state transmissions to other Minnesota locations (n=772) over the two-year period. We also analyzed the diversity of the location source of SARS-CoV-2 viruses in each county and noted the timing of state-wide policies as well as trends in clinical cases. Neither the number of clinical cases or major policy decisions, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, appeared to have impact on virus diversity across each individual county.

Evaluation of miniaturized Illumina DNA preparation protocols for SARS-CoV-2 whole genome sequencing

The global pandemic caused by SARS-CoV-2 has increased the demand for scalable sequencing and diagnostic methods, especially for genomic surveillance. Although next-generation sequencing has enabled large-scale genomic surveillance, the ability to sequence SARS-CoV-2 in some settings has been limited by the cost of sequencing kits and the time-consuming preparations of sequencing libraries. We compared the sequencing outcomes, cost and turn-around times obtained using the standard Illumina DNA Prep kit protocol to three modified protocols with fewer clean-up steps and different reagent volumes (full volume, half volume, one-tenth volume). We processed a single run of 47 samples under each protocol and compared the yield and mean sequence coverage. The sequencing success rate and quality for the four different reactions were as follows: the full reaction was 98.2%, the one-tenth reaction was 98.0%, the full rapid reaction was 97.5% and the half-reaction, was 97.1%. As a result, uniformity of sequence quality indicated that libraries were not affected by the change in protocol. The cost of sequencing was reduced approximately seven-fold and the time taken to prepare the library was reduced from 6.5 hours to 3 hours. The sequencing results obtained using the miniaturised volumes showed comparability to the results obtained using full volumes as described by the manufacturer. The adaptation of the protocol represents a lower-cost, streamlined approach for SARS-CoV-2 sequencing, which can be used to produce genomic data quickly and more affordably, especially in resource-constrained settings.

Maximum likelihood pandemic-scale phylogenetics

Phylogenetics has a crucial role in genomic epidemiology. Enabled by unparalleled volumes of genome sequence data generated to study and help contain the COVID-19 pandemic, phylogenetic analyses of SARS-CoV-2 genomes have shed light on the virus's origins, spread, and the emergence and reproductive success of new variants. However, most phylogenetic approaches, including maximum likelihood and Bayesian methods, cannot scale to the size of the datasets from the current pandemic. We present 'MAximum Parsimonious Likelihood Estimation' (MAPLE), an approach for likelihood-based phylogenetic analysis of epidemiological genomic datasets at unprecedented scales. MAPLE infers SARS-CoV-2 phylogenies more accurately than existing maximum likelihood approaches while running up to thousands of times faster, and requiring at least 100 times less memory on large datasets. This extends the reach of genomic epidemiology, allowing the continued use of accurate phylogenetic, phylogeographic and phylodynamic analyses on datasets of millions of genomes.

Tracking the Genomic Evolution of SARS-CoV-2 for 29 Months in South Korea

The pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued, with the persistent emergence of variants of concern (VOCs). Therefore, this study aimed to track the genomic evolution of SARS-CoV-2 strains by sequencing the spike protein for 29 months, which accounted for the majority of the COVID-19 pandemic period. A total of 109 swabs from patients with confirmed coronavirus disease 2019 (COVID-19) infection were randomly collected between March 2020 and July 2022. After genomic sequencing, we analyzed the naming systems and phylogenetic trees. Five surge peaks of COVID-19 cases have been reported in South Korea, resulting in 14,000,000 cumulative confirmed cases and 17,000 deaths. Among the sequenced samples, 34 wild-type strains and 75 VOCs, including 4 Alpha, 33 Delta, 2 Epsilon, and 36 Omicron VOCs, were identified. Omicron strains were comprised of 8 BA.1.1 (21 K), 27 BA.2 (21 L), and 1 BA.2.12.1 (22C). Phylogenetic analysis of the identified isolates and representative sequences of SARS-CoV-2 strains revealed clusters that presented the WHO VOCs. Specific or unique mutations for each VOC waxed and waned according to the variant waves. Our findings allowed recognition of the overall trends of SARS-CoV-2 isolates, which implicated replication advantage, immune evasion, and disease management.

Faster SARS-CoV-2 sequence validation and annotation for GenBank using VADR

In 2020 and 2021, >1.5 million SARS-CoV-2 sequences were submitted to GenBank. The initial version (v1.0) of the VADR (Viral Annotation DefineR) software package that GenBank uses to automatically validate and annotate incoming viral sequences is too slow and memory intensive to process many thousands of SARS-CoV-2 sequences in a reasonable amount of time. Additionally, long stretches of ambiguous N nucleotides, which are common in many SARS-CoV-2 sequences, prevent VADR from accurate validation and annotation. VADR has been updated to more accurately and rapidly annotate SARS-CoV-2 sequences. Stretches of consecutive Ns are now identified and temporarily replaced with expected nucleotides to facilitate processing, and the slowest steps have been overhauled using blastn and glsearch, increasing speed, reducing the memory requirement from 64Gb to 2Gb per thread, and allowing simple, coarse-grained parallelization on multiple processors per host. VADR is now nearly 1000 times faster than it was in early 2020 SARS-CoV-2 sequence processing. It has been used to screen and annotate more than 1.5 million SARS-CoV-2 sequences since June 2020, and it is now efficient enough to cope with the current rate of hundreds of thousands of submitted sequences per month.

Dissecting Phenotype from Genotype with Clinical Isolates of SARS-CoV-2 First Wave Variants

The emergence and availability of closely related clinical isolates of SARS-CoV-2 offers a unique opportunity to identify novel nonsynonymous mutations that may impact phenotype. Global sequencing efforts show that SARS-CoV-2 variants have emerged and then been replaced since the beginning of the pandemic, yet we have limited information regarding the breadth of variant-specific host responses. Using primary cell cultures and the K18-hACE2 mouse, we investigated the replication, innate immune response, and pathology of closely related, clinical variants circulating during the first wave of the pandemic. Mathematical modeling of the lung viral replication of four clinical isolates showed a dichotomy between two B.1. isolates with significantly faster and slower infected cell clearance rates, respectively. While isolates induced several common immune host responses to infection, one B.1 isolate was unique in the promotion of eosinophil-associated proteins IL-5 and CCL11. Moreover, its mortality rate was significantly slower. Lung microscopic histopathology suggested further phenotypic divergence among the five isolates showing three distinct sets of phenotypes: (i) consolidation, alveolar hemorrhage, and inflammation, (ii) interstitial inflammation/septal thickening and peribronchiolar/perivascular lymphoid cells, and (iii) consolidation, alveolar involvement, and endothelial hypertrophy/margination. Together these findings show divergence in the phenotypic outcomes of these clinical isolates and reveal the potential importance of nonsynonymous mutations in nsp2 and ORF8.

Rapid System to Detect Variants of SARS-CoV-2 in Nasopharyngeal Swabs

Currently, the reference method for identifying the presence of variants of SARS-CoV-2 is whole genome sequencing. Although it is less expensive than in the past, it is still time-consuming, and interpreting the results is difficult, requiring staff with specific skills who are not always available in diagnostic laboratories. The test presented in this study aimed to detect, using traditional real-time PCR, the presence of the main variants described for the spike protein of the SARS-CoV-2 genome. The primers and probes were designed to detect the main deletions that characterize the different variants. The amplification targets were deletions in the S gene: 25-27, 69-70, 241-243, and 157-158. In the ORF1a gene, the deletion 3675-3677 was chosen. Some of these mutations can be considered specific variants, while others can be identified by the simultaneous presence of one or more deletions. We avoided using point mutations in order to improve the speed of the test. Our test can help clinical and medical microbiologists quickly recognize the presence of variants in biological samples (particularly nasopharyngeal swabs). The test can also be used to identify variants of the virus that could potentially be more diffusive as well as not responsive to the vaccine.

Retrospective Insights of the COVID-19 Epidemic in the Major Latin American City, São Paulo, Southeastern Brazil

São Paulo is the financial center of Brazil, with a population of over 12 million, that receives travelers from all over the world for business and tourism. It was the first city in Brazil to report a case of COVID-19 that rapidly spread across the city despite the implementation of the restriction measures. Despite many reports, much is still unknown regarding the genomic diversity and transmission dynamics of this virus in the city of São Paulo. Thus, in this study, we provide a retrospective overview of the COVID-19 epidemic in São Paulo City, Southeastern, Brazil, by generating a total of 9995 near-complete genome sequences from all the city's different macro-regions (North, West, Central, East, South, and Southeast). Our analysis revealed that multiple independent introduction events of different variants (mainly Gamma, Delta, and Omicron) occurred throughout time. Additionally, our estimates of viral movement within the different macro-regions further suggested that the East and the Southeast regions were the largest contributors to the Gamma and Delta viral exchanges to other regions. Meanwhile, the North region had a higher contribution to the dispersion of the Omicron variant. Together, our results reinforce the importance of increasing SARS-CoV-2 genomic monitoring within the city and the country to track the real-time evolution of the virus and to detect earlier any eventual emergency of new variants of concern that could undermine the fight against COVID-19 in Brazil and worldwide.

Origin, Genetic Variation and Molecular Epidemiology of SARS-CoV-2 Strains Circulating in Sardinia (Italy) during the First and Second COVID-19 Epidemic Waves

Understanding how geography and human mobility shape the patterns and spread of infectious diseases such as COVID-19 is key to control future epidemics. An interesting example is provided by the second wave of the COVID-19 epidemic in Europe, which was facilitated by the intense movement of tourists around the Mediterranean coast in summer 2020. The Italian island of Sardinia is a major tourist destination and is widely believed to be the origin of the second Italian wave. In this study, we characterize the genetic variation among SARS-CoV-2 strains circulating in northern Sardinia during the first and second Italian waves using both Illumina and Oxford Nanopore Technologies Next Generation Sequencing methods. Most viruses were placed into a single clade, implying that despite substantial virus inflow, most outbreaks did not spread widely. The second epidemic wave on the island was actually driven by local transmission of a single B.1.177 subclade. Phylogeographic analyses further suggest that those viral strains circulating on the island were not a relevant source for the second epidemic wave in Italy. This result, however, does not rule out the possibility of intense mixing and transmission of the virus among tourists as a major contributor to the second Italian wave.

The chloroplast genome of Cuphea hookeriana Walp. (Lythraceae), a Mexico ornamental plant

Cuphea hookeriana Walp. is an ornamental plant belonging to the Lythraceae. In this study, we reported the complete chloroplast (cp) genome sequence here and analyzed the phylogenetic relationship among Lythraceae plants. The length of the cp genome was 158,999 bp, including a large single-copy (LSC, 89,311 bp) region and a small single-copy (SSC, 18,436 bp) region separated by a pair of inverted repeats (IRs, 25,626 bp). There were 72 unique protein-coding genes (PCGs), 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes in the cp genome of C. hookeriana. A total of 223 simple sequence repeats (SSRs) and 34 long repeat sequences were identified. Phylogenetic analyses using maximum-likelihood (ML) revealed that C. hookeriana was close to C. hyssopifolia. In addition, the two Cuphea species were the sister group of Woodfordia fruticosa.

Local-Scale phylodynamics reveal differential community impact of SARS-CoV-2 in metropolitan US county

SARS-CoV-2 transmission is largely driven by heterogeneous dynamics at a local scale, leaving local health departments to design interventions with limited information. We analyzed SARS-CoV-2 genomes sampled between February 2020 and March 2022 jointly with epidemiological and cell phone mobility data to investigate fine scale spatiotemporal SARS-CoV-2 transmission dynamics in King County, Washington, a diverse, metropolitan US county. We applied an approximate structured coalescent approach to model transmission within and between North King County and South King County alongside the rate of outside introductions into the county. Our phylodynamic analyses reveal that following stay-at-home orders, the epidemic trajectories of North and South King County began to diverge. We find that South King County consistently had more reported and estimated cases, COVID-19 hospitalizations, and longer persistence of local viral transmission when compared to North King County, where viral importations from outside drove a larger proportion of new cases. Using mobility and demographic data, we also find that South King County experienced a more modest and less sustained reduction in mobility following stay-at-home orders than North King County, while also bearing more socioeconomic inequities that might contribute to a disproportionate burden of SARS-CoV-2 transmission. Overall, our findings suggest a role for local-scale phylodynamics in understanding the heterogeneous transmission landscape.

ViralVar: A Web Tool for Multilevel Visualization of SARS-CoV-2 Genomes

The unprecedented growth of publicly available SARS-CoV-2 genome sequence data has increased the demand for effective and accessible SARS-CoV-2 data analysis and visualization tools. The majority of the currently available tools either require computational expertise to deploy them or limit user input to preselected subsets of SARS-CoV-2 genomes. To address these limitations, we developed ViralVar, a publicly available, point-and-click webtool that gives users the freedom to investigate and visualize user-selected subsets of SARS-CoV-2 genomes obtained from the GISAID public database. ViralVar has two primary features that enable: (1) the visualization of the spatiotemporal dynamics of SARS-CoV-2 lineages and (2) a structural/functional analysis of genomic mutations. As proof-of-principle, ViralVar was used to explore the evolution of the SARS-CoV-2 pandemic in the USA in pediatric, adult, and elderly populations (n > 1.7 million genomes). Whereas the spatiotemporal dynamics of the variants did not differ between these age groups, several USA-specific sublineages arose relative to the rest of the world. Our development and utilization of ViralVar to provide insights on the evolution of SARS-CoV-2 in the USA demonstrates the importance of developing accessible tools to facilitate and accelerate the large-scale surveillance of circulating pathogens.

Evaluation of co-circulating pathogens and microbiome from COVID-19 infections

Co-infections or secondary infections with SARS-CoV-2 have the potential to affect disease severity and morbidity. Additionally, the potential influence of the nasal microbiome on COVID-19 illness is not well understood. In this study, we analyzed 203 residual samples, originally submitted for SARS-CoV-2 testing, for the presence of viral, bacterial, and fungal pathogens and non-pathogens using a comprehensive microarray technology, the Lawrence Livermore Microbial Detection Array (LLMDA). Eighty-seven percent of the samples were nasopharyngeal samples, and 23% of the samples were oral, nasal and oral pharyngeal swabs. We conducted bioinformatics analyses to examine differences in microbial populations of these samples, as a proxy for the nasal and oral microbiome, from SARS-CoV-2 positive and negative specimens. We found 91% concordance with the LLMDA relative to a diagnostic RT-qPCR assay for detection of SARS-CoV-2. Sixteen percent of all the samples (32/203) revealed the presence of an opportunistic bacterial or frank viral pathogen with the potential to cause co-infections. The two most detected bacteria, Streptococcus pyogenes and Streptococcus pneumoniae, were present in both SARS-CoV-2 positive and negative samples. Human metapneumovirus was the most prevalent viral pathogen in the SARS-CoV-2 negative samples. Sequence analysis of 16S rRNA was also conducted to evaluate bacterial diversity and confirm LLMDA results.

An Overview of Anthropogenic Actions as Drivers for Emerging and Re-Emerging Zoonotic Diseases

Population growth and industrialization have led to a race for greater food and supply productivity. As a result, the occupation and population of forest areas, contact with wildlife and their respective parasites and vectors, the trafficking and consumption of wildlife, the pollution of water sources, and the accumulation of waste occur more frequently. Concurrently, the agricultural and livestock production for human consumption has accelerated, often in a disorderly way, leading to the deforestation of areas that are essential for the planet's climatic and ecological balance. The effects of human actions on other ecosystems such as the marine ecosystem cause equally serious damage, such as the pollution of this habitat, and the reduction of the supply of fish and other animals, causing the coastal population to move to the continent. The sum of these factors leads to an increase in the demands such as housing, basic sanitation, and medical assistance, making these populations underserved and vulnerable to the effects of global warming and to the emergence of emerging and re-emerging diseases. In this article, we discuss the anthropic actions such as climate changes, urbanization, deforestation, the trafficking and eating of wild animals, as well as unsustainable agricultural intensification which are drivers for emerging and re-emerging of zoonotic pathogens such as viral (Ebola virus, hantaviruses, Hendravirus, Nipah virus, rabies, and severe acute respiratory syndrome coronavirus disease-2), bacterial (leptospirosis, Lyme borreliosis, and tuberculosis), parasitic (leishmaniasis) and fungal pathogens, which pose a substantial threat to the global community. Finally, we shed light on the urgent demand for the implementation of the One Health concept as a collaborative global approach to raise awareness and educate people about the science behind and the battle against zoonotic pathogens to mitigate the threat for both humans and animals.

Neutralizing Immunity Induced Against the Omicron BA.1 and BA.2 Variants in Vaccine Breakthrough Infections

BACKGROUND

As of early 2022, the Omicron variants are the predominant circulating lineages globally. Understanding neutralizing antibody responses against Omicron BA.1 and BA.2 after vaccine breakthrough infections will provide insights into BA.2 infectivity and susceptibility to subsequent reinfection.

METHODS

Live virus neutralization assays were used to study immunity against Delta and Omicron BA.1 and BA.2 variants in samples from 86 individuals, 24 unvaccinated (27.9%) and 62 vaccinated (72.1%), who were infected with Delta (n = 42, 48.8%) or BA.1 (n = 44, 51.2%). Among the 62 vaccinated individuals, 39 were unboosted (62.9%), whereas 23 were boosted (37.1%).

RESULTS

In unvaccinated infections, neutralizing antibodies (nAbs) against the three variants were weak or undetectable, except against Delta for Delta-infected individuals. Both Delta and BA.1 breakthrough infections resulted in strong nAb responses against ancestral wild-type and Delta lineages, but moderate nAb responses against BA.1 and BA.2, with similar titers between unboosted and boosted individuals. Antibody titers against BA.2 were generally higher than those against BA.1 in breakthrough infections.

CONCLUSIONS

These results underscore the decreased immunogenicity of BA.1 compared to BA.2, insufficient neutralizing immunity against BA.2 in unvaccinated individuals, and moderate to strong neutralizing immunity induced against BA.2 in Delta and BA.1 breakthrough infections.

A Robust, Highly Multiplexed Mass Spectrometry Assay to Identify SARS-CoV-2 Variants

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants are characterized by differences in transmissibility and response to therapeutics. Therefore, discriminating among them is vital for surveillance, infection prevention, and patient care. While whole-genome sequencing (WGS) is the "gold standard" for variant identification, molecular variant panels have become increasingly available. Most, however, are based on limited targets and have not undergone comprehensive evaluation. We assessed the diagnostic performance of the highly multiplexed Agena MassARRAY SARS-CoV-2 Variant Panel v3 to identify variants in a diverse set of 391 SARS-CoV-2 clinical RNA specimens collected across our health systems in New York City, USA and Bogotá, Colombia (September 2, 2020 to March 2, 2022). We demonstrated almost perfect levels of interrater agreement between this assay and WGS for 9 of 11 variant calls (κ ≥ 0.856) and 25 of 30 targets (κ ≥ 0.820) tested on the panel. The assay had a high diagnostic sensitivity (≥93.67%) for contemporary variants (e.g., Iota, Alpha, Delta, and Omicron [BA.1 sublineage]) and a high diagnostic specificity for all 11 variants (≥96.15%) and all 30 targets (≥94.34%) tested. Moreover, we highlighted distinct target patterns that could be utilized to identify variants not yet defined on the panel, including the Omicron BA.2 and other sublineages. These findings exemplified the power of highly multiplexed diagnostic panels to accurately call variants and the potential for target result signatures to elucidate new ones. IMPORTANCE The continued circulation of SARS-CoV-2 amid limited surveillance efforts and inconsistent vaccination of populations has resulted in the emergence of variants that uniquely impact public health systems. Thus, in conjunction with functional and clinical studies, continuous detection and identification are quintessential to informing diagnostic and public health measures. Furthermore, until WGS becomes more accessible in the clinical microbiology laboratory, the ideal assay for identifying variants must be robust, provide high resolution, and be adaptable to the evolving nature of viruses like SARS-CoV-2. Here, we highlighted the diagnostic capabilities of a highly multiplexed commercial assay to identify diverse SARS-CoV-2 lineages that circulated from September 2, 2020 to March 2, 2022 among patients seeking care in our health systems. This assay demonstrated variant-specific signatures of nucleotide/amino acid polymorphisms and underscored its utility for the detection of contemporary and emerging SARS-CoV-2 variants of concern.

[Physical inactivity in young people during home confinement due to COVID-19]

Background

Physical inactivity is a highly prevalent condition in the world and has been associated with increased susceptibility to develop comorbidities and present with severe respiratory distress syndrome due to COVID-19.

Objective

To identify the factors present in the family environment and the personal reasons associated with physical inactivity in young people during confinement at home due to COVID-19.

Material and methods

A cross-sectional study analyzed data from 1,326 young people, ages 15 - 18. To collect information, a questionnaire was constructed using the Google Forms tool and distributed through the WhatsApp application and email to collect the information.

Results

The prevalence of physical inactivity was 43.4%. Approximately 24.4% were overweight, and 8.8% were obese. Near 43.0% of young people reported living in an environment with a dysfunctional family. The multivariate logistic regression analysis showed that suffering from obesity, does not have space at home, or devices to exercise and present a change in emotions, are related to the physical inactivity of young people during confinement at home due to the COVID-19 pandemic.

Conclusions

It is important to promote a harmonious environment within the family and the personal development of a healthy lifestyle, during the period of application of the contingency plan due to the presence of a pandemic, in order to maintain a better healthy physical and mental state.

[Epidemiological and clinical characteristics of the COVID-19 epidemic in Mexico: Quintana Roo case]

Objective

Identify risk factors for severe outcome in Mexican patients with COVID-19 in the population of Quintana Roo.

Material and methods

Study of 5,916 who met the criteria for suspected cases of COVID-19, 2,531 confirmed by qrTPCR-Sars-CoV-2 tests, of which 1,486 were positive, among which they were classified as hospitalized (severe COVID-19) and outpatients. Multivariate logistic regression analysis was performed to explore the factors associated with the severity of COVID-19 and death as clinical outcomes. The basic reproduction number (R0) was calculated Statistical analysis) Endorsement of the ethics committee 2301.

Results

SARS-CoV-2 positive patients presented a high prevalence of hypertension 29.1%, diabetes 23.5%, obesity 24%, and 48.5% have at least one chronic disease. There is a high risk of severity for COVID-19 in patients with diabetes OR=3.14, hypertension OR=1.88, obesity OR=1.68, kidney disease OR=3.2, older than 65 years OR=13.6 and men OR=1.7. These factors also increase the risk of death up to 7.7 times. The maximum R0 during the epidemic was 2.4.

Conclusion

Liver and kidney disease, diabetes, hypertension, and obesity are significantly associated with severe COVID-19 and death.

Contact tracing reveals community transmission of COVID-19 in New York City

Understanding SARS-CoV-2 transmission within and among communities is critical for tailoring public health policies to local context. However, analysis of community transmission is challenging due to a lack of high-resolution surveillance and testing data. Here, using contact tracing records for 644,029 cases and their contacts in New York City during the second pandemic wave, we provide a detailed characterization of the operational performance of contact tracing and reconstruct exposure and transmission networks at individual and ZIP code scales. We find considerable heterogeneity in reported close contacts and secondary infections and evidence of extensive transmission across ZIP code areas. Our analysis reveals the spatial pattern of SARS-CoV-2 spread and communities that are tightly interconnected by exposure and transmission. We find that locations with higher vaccination coverage and lower numbers of visitors to points-of-interest had reduced within- and cross-ZIP code transmission events, highlighting potential measures for curtailing SARS-CoV-2 spread in urban settings.

Estimation of local time-varying reproduction numbers in noisy surveillance data

A valuable metric in understanding local infectious disease dynamics is the local time-varying reproduction number, i.e. the expected number of secondary local cases caused by each infected individual. Accurate estimation of this quantity requires distinguishing cases arising from local transmission from those imported from elsewhere. Realistically, we can expect identification of cases as local or imported to be imperfect. We study the propagation of such errors in estimation of the local time-varying reproduction number. In addition, we propose a Bayesian framework for estimation of the true local time-varying reproduction number when identification errors exist. And we illustrate the practical performance of our estimator through simulation studies and with outbreaks of COVID-19 in Hong Kong and Victoria, Australia. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.

Ecological and evolutionary dynamics of multi-strain RNA viruses

Potential interactions among co-circulating viral strains in host populations are often overlooked in the study of virus transmission. However, these interactions probably shape transmission dynamics by influencing host immune responses or altering the relative fitness among co-circulating strains. In this Review, we describe multi-strain dynamics from ecological and evolutionary perspectives, outline scales in which multi-strain dynamics occur and summarize important immunological, phylogenetic and mathematical modelling approaches used to quantify interactions among strains. We also discuss how host-pathogen interactions influence the co-circulation of pathogens. Finally, we highlight outstanding questions and knowledge gaps in the current theory and study of ecological and evolutionary dynamics of multi-strain viruses.

SARS-COV-2/COVID-19: scenario, epidemiology, adaptive mutations, and environmental factors

The coronavirus pandemic of 2019 has already exerted an enormous impact. For over a year, the worldwide pandemic has ravaged the whole globe, with approximately 250 million verified human infection cases and a mortality rate surpassing 4 million. While the genetic makeup of the related pathogen (SARS-CoV-2) was identified, many unknown facets remain a mystery, comprising the virus's origin and evolutionary trend. There were many rumors that SARS-CoV-2 was human-borne and its evolution was predicted many years ago, but scientific investigation proved them wrong and concluded that bats might be the origin of SARS-CoV-2 and pangolins act as intermediary species to transmit the virus from bats to humans. Airborne droplets were found to be the leading cause of human-to-human transmission of this virus, but later studies showed that contaminated surfaces and other environmental factors are also involved in its transmission. The evolution of different SARS-CoV-2 variants worsens the condition and has become a challenge to overcome this pandemic. The emergence of COVID-19 is still a mystery, and scientists are unable to explain the exact origin of SARS-CoV-2. This review sheds light on the possible origin of SARS-CoV-2, its transmission, and the key factors that worsen the situation.

SARS-CoV-2 chronological genomic evolution and epidemiology in the Middle East and North Africa (MENA) region as affected by vaccination, conflict and socioeconomical disparities: a population-based cohort study

OBJECTIVE

To describe the chronological genomic evolution of SARS-CoV-2 and its impact on public health in the Middle East and North Africa (MENA) region.

METHODS

This study analysed all available SARS-CoV-2 genomic sequences, metadata and rates of COVID-19 infection from the MENA region retrieved from the Global Initiative on Sharing All Influenza Data database from January 2020 to August 2021. Inferential and ‎descriptive statistics were conducted to describe the epidemiology of SARS-CoV-2.

RESULTS

Genomic surveillance of SARS-CoV-2 in the MENA region indicated that the variants in January 2020 predominately belonged to the G, GR, GH or O clades and that the most common variant of concern was Alpha. By August 2021, however, the GK clade dominated (57.4% of all sequenced genomes), followed by the G clade (18.7%) and the GR clade (11.6%). In August, the most commonly sequenced variants of concern were Delta in the Middle East region (91%); Alpha (44.3%) followed by Delta (29.7%) and Beta (25.3%) in the North Africa region; and Alpha (88.9%), followed by Delta (10%) in the fragile and conflict-affected regions of MENA. The mean proportion of the variants of concern among the total sequenced samples differed significantly by country (F=1.93, P=0.0112) but not by major MENA region (F=0.14, P=0.27) or by vaccination coverage (F=1.84, P=0.176).

CONCLUSION

This analysis of the genomic surveillance of SARS-CoV-2 provides an essential description the virus evolution and its impact on public health safety in the MENA region. As of August 2021, the Delta variant showed a genomic advantage in the MENA region. The MENA region includes several fragile and conflict-affected countries with extremely low levels of vaccination coverage and little genomic surveillance, which may soon exacerbate the existing health crisis within those countries and globally.

Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic

Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.

Genomic epidemiology of the first two waves of SARS-CoV-2 in Canada

Tracking the emergence and spread of SARS-CoV-2 lineages using phylogenetics has proven critical to inform the timing and stringency of COVID-19 public health interventions. We investigated the effectiveness of international travel restrictions at reducing SARS-CoV-2 importations and transmission in Canada in the first two waves of 2020 and early 2021. Maximum likelihood phylogenetic trees were used to infer viruses' geographic origins, enabling identification of 2263 (95% confidence interval: 2159-2366) introductions, including 680 (658-703) Canadian sublineages, which are international introductions resulting in sampled Canadian descendants, and 1582 (1501-1663) singletons, introductions with no sampled descendants. Of the sublineages seeded during the first wave, 49% (46-52%) originated from the USA and were primarily introduced into Quebec (39%) and Ontario (36%), while in the second wave, the USA was still the predominant source (43%), alongside a larger contribution from India (16%) and the UK (7%). Following implementation of restrictions on the entry of foreign nationals on 21 March 2020, importations declined from 58.5 (50.4-66.5) sublineages per week to 10.3-fold (8.3-15.0) lower within 4 weeks. Despite the drastic reduction in viral importations following travel restrictions, newly seeded sublineages in summer and fall 2020 contributed to the persistence of COVID-19 cases in the second wave, highlighting the importance of sustained interventions to reduce transmission. Importations rebounded further in November, bringing newly emergent variants of concern (VOCs). By the end of February 2021, there had been an estimated 30 (19-41) B.1.1.7 sublineages imported into Canada, which increasingly displaced previously circulating sublineages by the end of the second wave.Although viral importations are nearly inevitable when global prevalence is high, with fewer importations there are fewer opportunities for novel variants to spark outbreaks or outcompete previously circulating lineages.

matOptimize: a parallel tree optimization method enables online phylogenetics for SARS-CoV-2

MOTIVATION

Phylogenetic tree optimization is necessary for precise analysis of evolutionary and transmission dynamics, but existing tools are inadequate for handling the scale and pace of data produced during the coronavirus disease 2019 (COVID-19) pandemic. One transformative approach, online phylogenetics, aims to incrementally add samples to an ever-growing phylogeny, but there are no previously existing approaches that can efficiently optimize this vast phylogeny under the time constraints of the pandemic.

RESULTS

Here, we present matOptimize, a fast and memory-efficient phylogenetic tree optimization tool based on parsimony that can be parallelized across multiple CPU threads and nodes, and provides orders of magnitude improvement in runtime and peak memory usage compared to existing state-of-the-art methods. We have developed this method particularly to address the pressing need during the COVID-19 pandemic for daily maintenance and optimization of a comprehensive SARS-CoV-2 phylogeny. matOptimize is currently helping refine on a daily basis possibly the largest-ever phylogenetic tree, containing millions of SARS-CoV-2 sequences.

AVAILABILITY AND IMPLEMENTATION

The matOptimize code is freely available as part of the UShER package (https://github.com/yatisht/usher) and can also be installed via bioconda (https://bioconda.github.io/recipes/usher/README.html). All scripts we used to perform the experiments in this manuscript are available at https://github.com/yceh/matOptimize-experiments.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Contact tracing reveals community transmission of COVID-19 in New York City

Understanding SARS-CoV-2 transmission within and among communities is critical for tailoring public health policies to local context. However, analysis of community transmission is challenging due to a lack of high-resolution surveillance and testing data. Here, using contact tracing records for 644,029 cases and their contacts in New York City during the second pandemic wave, we provide a detailed characterization of the operational performance of contact tracing and reconstruct exposure and transmission networks at individual and ZIP code scales. We find considerable heterogeneity in reported close contacts and secondary infections and evidence of extensive transmission across ZIP code areas. Our analysis reveals the spatial pattern of SARS-CoV-2 spread and communities that are tightly interconnected by exposure and transmission. We find that higher vaccination coverage and reduced numbers of visitors to points-of-interest are associated with fewer within- and cross-ZIP code transmission events, highlighting potential measures for curtailing SARS-CoV-2 spread in urban settings.

COVID-19 Variant Detection with a High-Fidelity CRISPR-Cas12 Enzyme

Laboratory tests for the accurate and rapid identification of SARS-CoV-2 variants can potentially guide the treatment of COVID-19 patients and inform infection control and public health surveillance efforts. Here, we present the development and validation of a rapid COVID-19 variant DETECTR assay incorporating loop-mediated isothermal amplification (LAMP) followed by CRISPR-Cas12 based identification of single nucleotide polymorphism (SNP) mutations in the SARS-CoV-2 spike (S) gene. This assay targets the L452R, E484K/Q/A, and N501Y mutations, at least one of which is found in nearly all major variants. In a comparison of three different Cas12 enzymes, only the newly identified enzyme CasDx1 was able to accurately identify all targeted SNP mutations. An analysis pipeline for CRISPR-based SNP identification from 261 clinical samples yielded a SNP concordance of 97.3% and agreement of 98.9% (258 of 261) for SARS-CoV-2 lineage classification, using SARS-CoV-2 whole-genome sequencing and/or real-time RT-PCR as test comparators. We also showed that detection of the single E484A mutation was necessary and sufficient to accurately identify Omicron from other major circulating variants in patient samples. These findings demonstrate the utility of CRISPR-based DETECTR as a faster and simpler diagnostic method compared with sequencing for SARS-CoV-2 variant identification in clinical and public health laboratories.

Clinical performance characteristics of the Swift Normalase Amplicon Panel for sensitive recovery of SARS-CoV-2 genomes

Amplicon-based sequencing methods are central in characterizing the diversity, transmission, and evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), but need to be rigorously assessed for clinical utility. Herein, we validated the Swift Biosciences' SARS-CoV-2 Swift Normalase Amplicon Panels using remnant clinical specimens. High-quality genomes meeting our established library and sequence quality criteria were recovered from positive specimens, with 95% limit of detection of 40.08 SARS-CoV-2 copies/PCR. Breadth of genome recovery was evaluated across a range of C_T values (11.3 to 36.7; median, 21.6). Of 428 positive samples, 413 (96.5%) generated genomes with <10% unknown bases, with a mean genome coverage of 13,545× ± SD 8382×. No genomes were recovered from PCR-negative specimens (n = 30) or from specimens positive for non–SARS-CoV-2 respiratory viruses (n = 20). Compared with whole-genome shotgun metagenomic sequencing (n = 14) or Sanger sequencing for the spike gene (n = 11), pairwise identity between consensus sequences was 100% in all cases, with highly concordant allele frequencies (R² = 0.99) between Swift and shotgun libraries. When samples from different clades were mixed at varying ratios, expected variants were detected even in 1:99 mixtures. When deployed as a clinical test, 268 tests were performed in the first 23 weeks, with a median turnaround time of 11 days, ordered primarily for outbreak investigations and infection control.

Comparative Dynamics of Delta and Omicron SARS-CoV-2 Variants across and between California and Mexico

Evolutionary analysis using viral sequence data can elucidate the epidemiology of transmission. Using publicly available SARS-CoV-2 sequence and epidemiological data, we developed discrete phylogeographic models to interrogate the emergence and dispersal of the Delta and Omicron variants in 2021 between and across California and Mexico. External introductions of Delta and Omicron in the region peaked in early July (2021-07-10 [95% CI: 2021-04-20, 2021-11-01]) and mid-December (2021-12-15 [95% CI: 2021-11-14, 2022-01-09]), respectively, 3 months and 2 weeks after first detection. These repeated introductions coincided with domestic migration events with no evidence of a unique transmission hub. The spread of Omicron was most consistent with gravity centric patterns within Mexico. While cross-border events accounted for only 5.1% [95% CI: 4.3-6] of all Delta migration events, they accounted for 20.6% [95% CI: 12.4-29] of Omicron movements, paralleling the increase in international travel observed in late 2021. Our investigations of the Delta and Omicron epidemics in the California/Mexico region illustrate the complex interplay and the multiplicity of viral and structural factors that need to be considered to limit viral spread, even as vaccination is reducing disease burden. Understanding viral transmission patterns may help intra-governmental responses to viral epidemics.

Benchmarking and Assessment of Eight De Novo Genome Assemblers on Viral Next-Generation Sequencing Data, Including the SARS-CoV-2

Viral genomics has become crucial in clinical diagnostics and ecology, not to mention to stem the COVID-19 pandemic. Whole-genome sequencing (WGS) is pivotal in gaining an improved understanding of viral evolution, genomic epidemiology, infectious outbreaks, pathobiology, clinical management, and vaccine development. Genome assembly is one of the crucial steps in WGS data analyses. A series of different assemblers has been developed with the advent of high-throughput next-generation sequencing (NGS). Various studies have reported the evaluation of these assembly tools on distinct datasets; however, these lack data from viral origin. In this study, we performed a comparative evaluation and benchmarking of eight de novo assemblers: SOAPdenovo, Velvet, assembly by short sequences (ABySS), iterative De Bruijn graph assembler (IDBA), SPAdes, Edena, iterative virus assembler, and VICUNA on the viral NGS data from distinct Illumina (GAIIx, Hiseq, Miseq, and Nextseq) platforms. WGS data of diverse viruses, that is, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), dengue virus 3, human immunodeficiency virus 1, hepatitis B virus, human herpesvirus 8, human papillomavirus 16, rhinovirus A, and West Nile virus, were utilized to assess these assemblers. Performance metrics such as genome fraction recovery, assembly lengths, NG50, N50, contig length, contig numbers, mismatches, and misassemblies were analyzed. Overall, three assemblers, that is, SPAdes, IDBA, and ABySS, performed consistently well, including for genome assembly of SARS-CoV-2. These assembly methods should be considered and recommended for future studies of viruses. The study also suggests that implementing two or more assembly approaches should be considered in viral NGS studies, especially in clinical settings. Taken together, the benchmarking of eight de novo genome assemblers reported in this study can inform future public health and ecology research concerning the viruses, the COVID-19 pandemic, and viral outbreaks.

SARS-CoV-2 variants preferentially emerge at intrinsically disordered protein sites helping immune evasion

The SARS‐CoV‐2 pandemic is maintained by the emergence of successive variants, highlighting the flexibility of the protein sequences of the virus. We show that experimentally determined intrinsically disordered regions (IDRs) are abundant in the SARS‐CoV‐2 viral proteins, making up to 28% of disorder content for the S1 subunit of spike and up to 51% for the nucleoprotein, with the vast majority of mutations occurring in the 13 major variants mapped to these IDRs. Strikingly, antigenic sites are enriched in IDRs, in the receptor‐binding domain (RBD) and in the N‐terminal domain (NTD), suggesting a key role of structural flexibility in the antigenicity of the SARS‐CoV‐2 protein surface. Mutations occurring in the S1 subunit and nucleoprotein (N) IDRs are critical for immune evasion and antibody escape, suggesting potential additional implications for vaccines and monoclonal therapeutic strategies. Overall, this suggests the presence of variable regions on S1 and N protein surfaces, which confer sequence and antigenic flexibility to the virus without altering its protein functions.

A robust, highly multiplexed mass spectrometry assay to identify SARS-CoV-2 variants

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants are characterized by differences in transmissibility and response to therapeutics. Therefore, discriminating among them is vital for surveillance, infection prevention, and patient care. While whole viral genome sequencing (WGS) is the "gold standard" for variant identification, molecular variant panels have become increasingly available. Most, however, are based on limited targets and have not undergone comprehensive evaluation. We assessed the diagnostic performance of the highly multiplexed Agena MassARRAY ® SARS-CoV-2 Variant Panel v3 to identify variants in a diverse set of 391 SARS-CoV-2 clinical RNA specimens collected across our health systems in New York City, USA as well as in Bogotá, Colombia (September 2, 2020 - March 2, 2022). We demonstrate almost perfect levels of interrater agreement between this assay and WGS for 9 of 11 variant calls (κ ≥ 0.856) and 25 of 30 targets (κ ≥ 0.820) tested on the panel. The assay had a high diagnostic sensitivity (≥93.67%) for contemporary variants (e.g., Iota, Alpha, Delta, Omicron [BA.1 sublineage]) and a high diagnostic specificity for all 11 variants (≥96.15%) and all 30 targets (≥94.34%) tested. Moreover, we highlight distinct target patterns that can be utilized to identify variants not yet defined on the panel including the Omicron BA.2 and other sublineages. These findings exemplify the power of highly multiplexed diagnostic panels to accurately call variants and the potential for target result signatures to elucidate new ones.

Importance

The continued circulation of SARS-CoV-2 amidst limited surveillance efforts and inconsistent vaccination of populations has resulted in emergence of variants that uniquely impact public health systems. Thus, in conjunction with functional and clinical studies, continuous detection and identification are quintessential to inform diagnostic and public health measures. Furthermore, until WGS becomes more accessible in the clinical microbiology laboratory, the ideal assay for identifying variants must be robust, provide high resolution, and be adaptable to the evolving nature of viruses like SARS-CoV-2. Here, we highlight the diagnostic capabilities of a highly multiplexed commercial assay to identify diverse SARS-CoV-2 lineages that circulated at over September 2, 2020 - March 2, 2022 among patients seeking care at our health systems. This assay demonstrates variant-specific signatures of nucleotide/amino acid polymorphisms and underscores its utility for detection of contemporary and emerging SARS-CoV-2 variants of concern.

Role of genomics in combating COVID-19 pandemic

The coronavirus disease 2019 (COVID-19) quickly swept over the world, becoming one of the most devastating outbreaks in human history. Being the first pandemic in the post-genomic era, advancements in genomics contributed significantly to scientific understanding and public health response to COVID-19. Genomic technologies have been employed by researchers all over the world to better understand the biology of SARS-CoV-2 and its origin, genomic diversity, and evolution. Worldwide genomic resources have greatly aided in the investigation of the COVID-19 pandemic. The pandemic has ushered in a new era of genomic surveillance, wherein scientists are tracking the changes of the SARS-CoV-2 genome in real-time at the international and national levels. Availability of genomic and proteomic information enables the rapid development of molecular diagnostics and therapeutics. The advent of high-throughput sequencing and genome editing technologies led to the development of modern vaccines. We briefly discuss the impact of genomics in the ongoing COVID-19 pandemic in this review.

Epidemiologic and Genomic Analysis of the Severe Acute Respiratory Syndrome Coronavirus 2 Epidemic in the Nebraska Region of the United States, March 2020-2021

COVID-19 emerged at varying intervals in different regions of the United States in 2020. This report details the epidemiologic and genetic evolution of Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during the first year of the epidemic in the state of Nebraska using data collected from the Creighton Catholic Health Initiatives (CHI) health system. Statistical modelling identified age, gender, and previous history of diabetes and/or stroke as significant risk factors associated with mortality in COVID-19 patients. In parallel, the viral genomes of over 1,000 samples were sequenced. The overall rate of viral variation in the population was 0.07 mutations/day. Genetically, the first 9 months of the outbreak, which include the initial outbreak, a small surge in August and a major outbreak in November 2020 were primarily characterized by B.1. lineage viruses. In early 2021, the United Kingdom variant (B.1.1.7 or alpha) quickly became the dominant variant. Notably, surveillance of non-consensus variants detected B.1.1.7 defining mutations months earlier in Fall 2020. This work provides insights into the regional variance and evolution of SARS-CoV-2 in the Nebraska region during the first year of the pandemic.