Sequential viral introductions and spread of BA.1 across Pakistan provinces during the Omicron wave

BACKGROUND

COVID-19 waves caused by specific SARS-CoV-2 variants have occurred globally at different times. We focused on Omicron variants to understand the genomic diversity and phylogenetic relatedness of SARS-CoV-2 strains in various regions of Pakistan.

METHODS

We studied 276,525 COVID-19 cases and 1,031 genomes sequenced from December 2021 to August 2022. Sequences were analyzed and visualized using phylogenetic trees.

RESULTS

The highest case numbers and deaths were recorded in Sindh and Punjab, the most populous provinces in Pakistan. Omicron variants comprised 93% of all genomes, with BA.2 (32.6%) and BA.5 (38.4%) predominating. The first Omicron wave was associated with the sequential identification of BA.1 in Sindh, then Islamabad Capital Territory, Punjab, Khyber Pakhtunkhwa (KP), Azad Jammu Kashmir (AJK), Gilgit-Baltistan (GB) and Balochistan. Phylogenetic analysis revealed Sindh to be the source of BA.1 and BA.2 introductions into Punjab and Balochistan during early 2022. BA.4 was first introduced in AJK and BA.5 in Punjab. Most recent common ancestor (MRCA) analysis revealed relatedness between the earliest BA.1 genome from Sindh with Balochistan, AJK, Punjab and ICT, and that of first BA.1 from Punjab with strains from KPK and GB.

CONCLUSIONS

Phylogenetic analysis provides insights into the introduction and transmission dynamics of the Omicron variant in Pakistan, identifying Sindh as a hotspot for viral dissemination. Such data linked with public health efforts can help limit surges of new infections.

Evolutionary History and Introduction of SARS-CoV-2 Alpha VOC/B.1.1.7 in Pakistan Through International Travelers

SARS-CoV-2 variants continue to emerge, and their identification is important for the public health response to COVID-19. Genomic sequencing provides robust information but may not always be accessible and therefore rapid mutation-based PCR approaches can be used to identify known variants. International travelers arriving in Karachi between December 2020 and February 2021 were tested for SARS-CoV-2 by PCR. A subset of positive samples was tested for S-Gene Target Failure (SGTF) on TaqPathTM COVID-19 (Thermo Fisher Scientific) and for mutations using the GSD NovaType SARS-CoV-2 (Eurofins Technologies) assays. Sequencing was conducted on the MinION platform (Oxford Nanopore Technologies (ONT). Bayesian phylogeographic inference was performed integrating the patients’ travel history information. Of the thirty-five COVID-19 cases screened, thirteen had isolates with SGTF. The travelers transmitted infection to sixty-eight contact cases. The B.1.1.7 lineage was confirmed through sequencing and PCR. Phylogenetic analysis of sequence data available for six cases included four B.1.1.7 strains and one B.1.36 and B.1.1.212 lineage isolate, respectively. Phylogeographic modeling estimated at least three independent B.1.1.7 introductions into Karachi, Pakistan, originating from the UK. B.1.1.212 and B.1.36 were inferred to be introduced either from the UK or the travelers’ layover location. We report the introduction of SARS-CoV-2 B.1.1.7 and other lineages in Pakistan by international travelers arriving via different flight routes. This highlights SARS-CoV-2 transmission through travel, importance of testing and quarantine post-travel to prevent transmission of new strains, as well as recording detailed patients’ metadata. Such results help inform policies on restricting travel from destinations where new highly transmissible variants have emerged.

Genomic diversity of SARS-CoV-2 during early introduction into the Baltimore-Washington metropolitan area

The early COVID-19 pandemic was characterized by rapid global spread. In Maryland and Washington, DC, United States, more than 2500 cases were reported within 3 weeks of the first COVID-19 detection in March 2020. We aimed to use genomic sequencing to understand the initial spread of SARS-CoV-2 — the virus that causes COVID-19 — in the region. We analyzed 620 samples collected from the Johns Hopkins Health System during March 11–31, 2020, comprising 28.6% of the total cases in Maryland and Washington, DC. From these samples, we generated 114 complete viral genomes. Analysis of these genomes alongside a subsampling of over 1000 previously published sequences showed that the diversity in this region rivaled global SARS-CoV-2 genetic diversity at that time and that the sequences belong to all of the major globally circulating lineages, suggesting multiple introductions into the region. We also analyzed these regional SARS-CoV-2 genomes alongside detailed clinical metadata and found that clinically severe cases had viral genomes belonging to all major viral lineages. We conclude that efforts to control local spread of the virus were likely confounded by the number of introductions into the region early in the epidemic and the interconnectedness of the region as a whole.

Reference Genome for the Highly Transformable Setaria viridis ME034V

Setaria viridis (green foxtail) is an important model system for improving cereal crops due to its diploid genome, ease of cultivation, and use of C4 photosynthesis. The S. viridis accession ME034V is exceptionally transformable, but the lack of a sequenced genome for this accession has limited its utility. We present a 397 Mb highly contiguous de novo assembly of ME034V using ultra-long nanopore sequencing technology (read N50 = 41kb). We estimate that this genome is largely complete based on our updated k-mer based genome size estimate of 401 Mb for S. viridis Genome annotation identified 37,908 protein-coding genes and >300k repetitive elements comprising 46% of the genome. We compared the ME034V assembly with two other previously sequenced Setaria genomes as well as to a diversity panel of 235 S. viridis accessions. We found the genome assemblies to be largely syntenic, but numerous unique polymorphic structural variants were discovered. Several ME034V deletions may be associated with recent retrotransposition of copia and gypsy LTR repeat families, as evidenced by their low genotype frequencies in the sampled population. Lastly, we performed a phylogenomic analysis to identify gene families that have expanded in Setaria, including those involved in specialized metabolism and plant defense response. The high continuity of the ME034V genome assembly validates the utility of ultra-long DNA sequencing to improve genetic resources for emerging model organisms. Structural variation present in Setaria illustrates the importance of obtaining the proper genome reference for genetic experiments. Thus, we anticipate that the ME034V genome will be of significant utility for the Setaria research community.

Genomic Diversity of SARS-CoV-2 During Early Introduction into the United States National Capital Region

BACKGROUND

The early COVID-19 pandemic has been characterized by rapid global spread. In the United States National Capital Region, over 2,000 cases were reported within three weeks of its first detection in March 2020. We aimed to use genomic sequencing to understand the initial spread of SARS-CoV-2, the virus that causes COVID-19, in the region. By correlating genetic information to disease phenotype, we also aimed to gain insight into any correlation between viral genotype and case severity or transmissibility.

METHODS

We performed whole genome sequencing of clinical SARS-CoV-2 samples collected in March 2020 by the Johns Hopkins Health System. We analyzed these regional SARS-CoV-2 genomes alongside detailed clinical metadata and the global phylogeny to understand early establishment of the virus within the region.

RESULTS

We analyzed 620 samples from the Johns Hopkins Health System collected between March 11-31, 2020, comprising 37.3% of the total cases in Maryland during this period. We selected 143 of these samples for sequencing, generating 114 complete viral genomes. These genomes belong to all five major Nextstrain-defined clades, suggesting multiple introductions into the region and underscoring the diversity of the regional epidemic. We also found that clinically severe cases had genomes belonging to all of these clades.

CONCLUSIONS

We established a pipeline for SARS-CoV-2 sequencing within the Johns Hopkins Health system, which enabled us to capture the significant viral diversity present in the region as early as March 2020. Efforts to control local spread of the virus were likely confounded by the number of introductions into the region early in the epidemic and interconnectedness of the region as a whole.

Global marine biodiversity in the context of achieving the Aichi Targets: ways forward and addressing data gaps

In 2010, the Conference of the Parties of the Convention on Biological Diversity agreed on the Strategic Plan for Biodiversity 2011–2020 in Aichi Prefecture, Japan. As this plan approaches its end, we discussed whether marine biodiversity and prediction studies were nearing the Aichi Targets during the 4th World Conference on Marine Biodiversity held in Montreal, Canada in June 2018. This article summarises the outcome of a five-day group discussion on how global marine biodiversity studies should be focused further to better understand the patterns of biodiversity. We discussed and reviewed seven fundamental biodiversity priorities related to nine Aichi Targets focusing on global biodiversity discovery and predictions to improve and enhance biodiversity data standards (quantity and quality), tools and techniques, spatial and temporal scale framing, and stewardship and dissemination. We discuss how identifying biodiversity knowledge gaps and promoting efforts have and will reduce such gaps, including via the use of new databases, tools and technology, and how these resources could be improved in the future. The group recognised significant progress toward Target 19 in relation to scientific knowledge, but negligible progress with regard to Targets 6 to 13 which aimed to safeguard and reduce human impacts on biodiversity.

Evolution on the Biophysical Fitness Landscape of an RNA Virus

Viral evolutionary pathways are determined by the fitness landscape, which maps viral genotype to fitness. However, a quantitative description of the landscape and the evolutionary forces on it remain elusive. Here, we apply a biophysical fitness model based on capsid folding stability and antibody binding affinity to predict the evolutionary pathway of norovirus escaping a neutralizing antibody. The model is validated by experimental evolution in bulk culture and in a drop-based microfluidics that propagates millions of independent small viral subpopulations. We demonstrate that along the axis of binding affinity, selection for escape variants and drift due to random mutations have the same direction, an atypical case in evolution. However, along folding stability, selection and drift are opposing forces whose balance is tuned by viral population size. Our results demonstrate that predictable epistatic tradeoffs between molecular traits of viral proteins shape viral evolution.

A high-throughput drop microfluidic system for virus culture and analysis

High mutation rates and short replication times lead to rapid evolution in RNA viruses. New tools for high-throughput culture and analysis of viral phenotypes will enable more effective studies of viral evolutionary processes. A water-in-oil drop microfluidic system to study virus-cell interactions at the single event level on a massively parallel scale is described here. Murine norovirus (MNV-1) particles were co-encapsulated with individual RAW 264.7 cells in 65 pL aqueous drops formed by flow focusing in 50 μm microchannels. At low multiplicity of infection (MOI), viral titers increased greatly, reaching a maximum 18 h post-encapsulation. This system was employed to evaluate MNV-1 escape from a neutralizing monoclonal antibody (clone A6.2). Further, the system was validated as a means for testing escape from antibody neutralization using a series of viral point mutants. Finally, the replicative capacity of single viral particles in drops under antibody stress was tested. Under standard conditions, many RNA virus stocks harbor minority populations of genotypic and phenotypic variants, resulting in quasispecies. These data show that when single cells are encapsulated with single viral particles under antibody stress without competition from other virions, the number of resulting infectious particles is nearly equivalent to the number of viral genomes present. These findings suggest that lower fitness virions can infect cells successfully and replicate, indicating that the microfluidics system may serve as an effective tool for isolating mutants that escape evolutionary stressors.