Movement of the A-strain maize streak virus in and out of Madagascar.

The maize streak virus belongs in the genus Mastrevirus, in the family Geminiviridae. The A-strain of the virus (MSV-A) is recognised as the principal causative agent of the most severe manifestation of maize streak disease (MSD). This disease continues to be a persistent limitation on maize output across sub-Saharan Africa and the nearby Indian Ocean islands. Irrespective of the causes behind the spread of MSV-A, we can determine the paths and speeds with which MSV-A spreads by analysing MSV genome sequence data along with information on when and where samples were taken. This information is valuable for identifying the geographical origins of viral strains that cause sporadic MSD epidemics in specific places and the geographical regions where viruses remain in reservoirs and contribute to prolonged epidemics during outbreaks. Our aim is to utilise these analyses to estimate the timing and origin of the MSV-A that arrived on the island of Madagascar in the Indian Ocean. Specifically, we employ model-based phylogeographic analyses on 524 complete MSV-A genome sequences, which consist of 56 newly obtained genomes from infected maize plants collected in Madagascar. These studies allow us to reconstruct the most likely paths of MSV-A to Madagascar. We found strong evidence for the existence of at least four separate movements of MSV-A variants from East and southern Africa to Madagascar. These movements took place between roughly 1979 (with a 95% highest probability density interval [HPD] ranging from 1976 to 1982) and 2003 (with a 95% HPD ranging from 2002 to 2003). While we inferred that MSV-A variants are spreading at an average rate of 38.9 km/year (with a 95% highest posterior density interval of 34.0-44.4) across their geographical range. Since their arrival in Madagascar, MSV-A variants have been migrating at an average rate of 47.6 km/year (with a 95% highest posterior density interval of 36.05-61.70). Human influences are likely significant contributors to both sporadic long-range movements of MSV-A between mainland Africa and Madagascar, as well as shorter to medium range movements within the island.

Network-based integrative multi-omics approach reveals biosignatures specific to COVID-19 disease phases.

Background: COVID-19 disease is characterized by a spectrum of disease phases (mild, moderate, and severe). Each disease phase is marked by changes in omics profiles with corresponding changes in the expression of features (biosignatures). However, integrative analysis of multiple omics data from different experiments across studies to investigate biosignatures at various disease phases is limited. Exploring an integrative multi-omics profile analysis through a network approach could be used to determine biosignatures associated with specific disease phases and enable the examination of the relationships between the biosignatures. Aim: To identify and characterize biosignatures underlying various COVID-19 disease phases in an integrative multi-omics data analysis. Method: We leveraged a multi-omics network-based approach to integrate transcriptomics, metabolomics, proteomics, and lipidomics data. The World Health Organization Ordinal Scale WHO Ordinal Scale was used as a disease severity reference to harmonize COVID-19 patient metadata across two studies with independent data. A unified COVID-19 knowledge graph was constructed by assembling a disease-specific interactome from the literature and databases. Disease-state specific omics-graphs were constructed by integrating multi-omics data with the unified COVID-19 knowledge graph. We expanded on the network layers of multiXrank, a random walk with restart on multilayer network algorithm, to explore disease state omics-specific graphs and perform enrichment analysis. Results: Network analysis revealed the biosignatures involved in inducing chemokines and inflammatory responses as hubs in the severe and moderate disease phases. We observed distinct biosignatures between severe and moderate disease phases as compared to mild-moderate and mild-severe disease phases. Mild COVID-19 cases were characterized by a unique biosignature comprising C-C Motif Chemokine Ligand 4 (CCL4), and Interferon Regulatory Factor 1 (IRF1). Hepatocyte Growth Factor (HGF), Matrix Metallopeptidase 12 (MMP12), Interleukin 10 (IL10), Nuclear Factor Kappa B Subunit 1 (NFKB1), and suberoylcarnitine form hubs in the omics network that characterizes the moderate disease state. The severe cases were marked by biosignatures such as Signal Transducer and Activator of Transcription 1 (STAT1), Superoxide Dismutase 2 (SOD2), HGF, taurine, lysophosphatidylcholine, diacylglycerol, triglycerides, and sphingomyelin that characterize the disease state. Conclusion: This study identified both biosignatures of different omics types enriched in disease-related pathways and their associated interactions (such as protein-protein, protein-transcript, protein-metabolite, transcript-metabolite, and lipid-lipid interactions) that are unique to mild, moderate, and severe COVID-19 disease states. These biosignatures include molecular features that underlie the observed clinical heterogeneity of COVID-19 and emphasize the need for disease-phase-specific treatment strategies. The approach implemented here can be used to find associations between transcripts, proteins, lipids, and metabolites in other diseases.

Virome release of an invasive exotic plant species in southern France.

The increase in human-mediated introduction of plant species to new regions has resulted in a rise of invasive exotic plant species (IEPS) that has had significant effects on biodiversity and ecosystem processes. One commonly accepted mechanism of invasions is that proposed by the enemy release hypothesis (ERH), which states that IEPS free from their native herbivores and natural enemies in new environments can outcompete indigenous species and become invasive. We here propose the virome release hypothesis (VRH) as a virus-centered variant of the conventional ERH that is only focused on enemies. The VRH predicts that vertically transmitted plant-associated viruses (PAV, encompassing phytoviruses and mycoviruses) should be co-introduced during the dissemination of the IEPS, while horizontally transmitted PAV of IEPS should be left behind or should not be locally transmitted in the introduced area due to a maladaptation of local vectors. To document the VRH, virome richness and composition as well as PAV prevalence, co-infection, host range, and transmission modes were compared between indigenous plant species and an invasive grass, cane bluestem (Bothriochloa barbinodis), in both its introduced range (southern France) and one area of its native range (Sonoran Desert, Arizona, USA). Contrary to the VRH, we show that invasive populations of B. barbinodis in France were not associated with a lower PAV prevalence or richness than native populations of B. barbinodis from the USA. However, comparison of virome compositions and network analyses further revealed more diverse and complex plant-virus interactions in the French ecosystem, with a significant richness of mycoviruses. Setting mycoviruses apart, only one putatively vertically transmitted phytovirus (belonging to the Amalgaviridae family) and one putatively horizontally transmitted phytovirus (belonging to the Geminiviridae family) were identified from B. barbinodis plants in the introduced area. Collectively, these characteristics of the B. barbinodis-associated PAV community in southern France suggest that a virome release phase may have immediately followed the introduction of B. barbinodis to France in the 1960s or 1970s, and that, since then, the invasive populations of this IEPS have already transitioned out of this virome release phase, and have started interacting with several local mycoviruses and a few local plant viruses.

An Oxford Nanopore Technology-Based Hepatitis B Virus Sequencing Protocol Suitable For Genomic Surveillance Within Clinical Diagnostic Settings.

Chronic hepatitis B virus (HBV) infection remains a significant public health concern, particularly in Africa, where there is a substantial burden. HBV is an enveloped virus, with isolates being classified into ten phylogenetically distinct genotypes (A - J) determined based on full-genome sequence data or reverse hybridization-based diagnostic tests. In practice, limitations are noted in that diagnostic sequencing, generally using Sanger sequencing, tends to focus only on the S-gene, yielding little or no information on intra-patient HBV genetic diversity with very low-frequency variants and reverse hybridization detects only known genotype-specific mutations. To resolve these limitations, we developed an Oxford Nanopore Technology (ONT)-based HBV genotyping protocol suitable for clinical virology, yielding complete HBV genome sequences and extensive data on intra-patient HBV diversity. Specifically, the protocol involves tiling-based PCR amplification of HBV sequences, library preparation using the ONT Rapid Barcoding Kit, ONT GridION sequencing, genotyping using Genome Detective software, recombination analysis using jpHMM and RDP5 software, and drug resistance profiling using Geno2pheno software. We prove the utility of our protocol by efficiently generating and characterizing high-quality near full-length HBV genomes from 148 left-over diagnostic Hepatitis B patient samples obtained in the Western Cape province of South Africa, providing valuable insights into the genetic diversity and epidemiology of HBV in this region of the world.

Rapid dynamic changes of FL.2 variant: A case report of COVID-19 breakthrough infection.

We investigated intra-host genetic evolution using two SARS-CoV-2 isolates from a fully vaccinated (primary schedule x2 doses of AstraZeneca plus a booster of Pfizer), >70-year-old woman with a history of lymphoma and hypertension who presented a SARS-CoV-2 infection for 3 weeks prior to death due to COVID-19. Two full genome sequences were determined from samples taken 13 days apart with both belonging to Pango lineage FL.2: the first detection of this Omicron sub-variant in Botswana. FL.2 is a sub-lineage of XBB.1.9.1. The repertoire of mutations and minority variants in the Spike protein differed between the two time points. Notably, we also observed deletions within the ORF1a and Membrane proteins; both regions are associated with high T-cell epitope density. The internal milieu of immune-suppressed individuals may accelerate SARS-CoV-2 evolution; hence, close monitoring is warranted.

Complete genome reconstruction of the global and European regional dispersal history of the lumpy skin disease virus.

IMPORTANCE: Lumpy skin disease virus (LSDV) has a complex epidemiology involving multiple strains, recombination, and vaccination. Its DNA genome provides limited genetic variation to trace outbreaks in space and time. Sequencing of LSDV whole genomes has also been patchy at global and regional scales. Here, we provide the first fine-grained whole genome sequence sampling of a constrained LSDV outbreak (southeastern Europe, 2015-2017), which we analyze along with global publicly available genomes. We formally evaluate the past occurrence of recombination events as well as the temporal signal that is required for calibrating molecular clock models and subsequently conduct a time-calibrated spatially explicit phylogeographic reconstruction. Our study further illustrates the importance of accounting for recombination events before reconstructing global and regional dynamics of DNA viruses. More LSDV whole genomes from endemic areas are needed to obtain a comprehensive understanding of global LSDV dispersal dynamics.

Metagenomics reveals the structure of Mastrevirus-host interaction network within an agro-ecosystem.

As highly pervasive parasites that sometimes cause disease, viruses are likely major components of all natural ecosystems. An important step towards both understanding the precise ecological roles of viruses and determining how natural communities of viral species are assembled and evolve is obtaining full descriptions of viral diversity and distributions at ecosystem scales. Here, we focused on obtaining such 'community-scale' data for viruses in a single genus. We chose the genus Mastrevirus (family Geminiviridae), members of which have predominantly been found infecting uncultivated grasses (family Poaceae) throughout the tropical and sub-tropical regions of the world. We sampled over 3 years, 2,884 individual Poaceae plants belonging to thirty different species within a 2-ha plot which included cultivated and uncultivated areas on the island of Reunion. Mastreviruses were found in ∼8 per cent of the samples, of which 96 per cent did not have any discernible disease symptoms. The multitude of host-virus associations that we uncovered reveals both the plant species that most commonly host mastreviruses and the mastrevirus species (such as maize streak virus and maize streak Reunion virus) that have especially large host ranges. Our findings are consistent with the hypothesis that perennial plant species capable of hosting years-long mixed mastrevirus infections likely play a disproportionately important role in the generation of inter-species and inter-strain mastrevirus recombinants.

Evolution of anelloviruses from a circovirus-like ancestor through gradual augmentation of the jelly-roll capsid protein.

Anelloviruses are highly prevalent in diverse mammals, including humans, but so far have not been linked to any disease and are considered to be part of the 'healthy virome'. These viruses have small circular single-stranded DNA (ssDNA) genomes and encode several proteins with no detectable sequence similarity to proteins of other known viruses. Thus, anelloviruses are the only family of eukaryotic ssDNA viruses currently not included in the realm Monodnaviria. To gain insights into the provenance of these enigmatic viruses, we sequenced more than 250 complete genomes of anelloviruses from nasal and vaginal swab samples of Weddell seal (Leptonychotes weddellii) from Antarctica and a fecal sample of grizzly bear (Ursus arctos horribilis) from the USA and performed a comprehensive family-wide analysis of the signature anellovirus protein ORF1. Using state-of-the-art remote sequence similarity detection approaches and structural modeling with AlphaFold2, we show that ORF1 orthologs from all Anelloviridae genera adopt a jelly-roll fold typical of viral capsid proteins (CPs), establishing an evolutionary link to other eukaryotic ssDNA viruses, specifically, circoviruses. However, unlike CPs of other ssDNA viruses, ORF1 encoded by anelloviruses from different genera display remarkable variation in size, due to insertions into the jelly-roll domain. In particular, the insertion between ß-strands H and I forms a projection domain predicted to face away from the capsid surface and function at the interface of virus-host interactions. Consistent with this prediction and supported by recent experimental evidence, the outermost region of the projection domain is a mutational hotspot, where rapid evolution was likely precipitated by the host immune system. Collectively, our findings further expand the known diversity of anelloviruses and explain how anellovirus ORF1 proteins likely diverged from canonical jelly-roll CPs through gradual augmentation of the projection domain. We suggest assigning Anelloviridae to a new phylum, 'Commensaviricota', and including it into the kingdom Shotokuvirae (realm Monodnaviria), alongside Cressdnaviricota and Cossaviricota.

RASCL: Rapid Assessment of Selection in CLades through molecular sequence analysis.

An important unmet need revealed by the COVID-19 pandemic is the near-real-time identification of potentially fitness-altering mutations within rapidly growing SARS-CoV-2 lineages. Although powerful molecular sequence analysis methods are available to detect and characterize patterns of natural selection within modestly sized gene-sequence datasets, the computational complexity of these methods and their sensitivity to sequencing errors render them effectively inapplicable in large-scale genomic surveillance contexts. Motivated by the need to analyze new lineage evolution in near-real time using large numbers of genomes, we developed the Rapid Assessment of Selection within CLades (RASCL) pipeline. RASCL applies state of the art phylogenetic comparative methods to evaluate selective processes acting at individual codon sites and across whole genes. RASCL is scalable and produces automatically updated regular lineage-specific selection analysis reports: even for lineages that include tens or hundreds of thousands of sampled genome sequences. Key to this performance is (i) generation of automatically subsampled high quality datasets of gene/ORF sequences drawn from a selected "query" viral lineage; (ii) contextualization of these query sequences in codon alignments that include high-quality "background" sequences representative of global SARS-CoV-2 diversity; and (iii) the extensive parallelization of a suite of computationally intensive selection analysis tests. Within hours of being deployed to analyze a novel rapidly growing lineage of interest, RASCL will begin yielding JavaScript Object Notation (JSON)-formatted reports that can be either imported into third-party analysis software or explored in standard web-browsers using the premade RASCL interactive data visualization dashboard. By enabling the rapid detection of genome sites evolving under different selective regimes, RASCL is well-suited for near-real-time monitoring of the population-level selective processes that will likely underlie the emergence of future variants of concern in measurably evolving pathogens with extensive genomic surveillance.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance.

Investment in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequencing in Africa over the past year has led to a major increase in the number of sequences that have been generated and used to track the pandemic on the continent, a number that now exceeds 100,000 genomes. Our results show an increase in the number of African countries that are able to sequence domestically and highlight that local sequencing enables faster turnaround times and more-regular routine surveillance. Despite limitations of low testing proportions, findings from this genomic surveillance study underscore the heterogeneous nature of the pandemic and illuminate the distinct dispersal dynamics of variants of concern-particularly Alpha, Beta, Delta, and Omicron-on the continent. Sustained investment for diagnostics and genomic surveillance in Africa is needed as the virus continues to evolve while the continent faces many emerging and reemerging infectious disease threats. These investments are crucial for pandemic preparedness and response and will serve the health of the continent well into the 21st century.

Conserved recombination patterns across coronavirus subgenera.

Recombination contributes to the genetic diversity found in coronaviruses and is known to be a prominent mechanism whereby they evolve. It is apparent, both from controlled experiments and in genome sequences sampled from nature, that patterns of recombination in coronaviruses are non-random and that this is likely attributable to a combination of sequence features that favour the occurrence of recombination break points at specific genomic sites, and selection disfavouring the survival of recombinants within which favourable intra-genome interactions have been disrupted. Here we leverage available whole-genome sequence data for six coronavirus subgenera to identify specific patterns of recombination that are conserved between multiple subgenera and then identify the likely factors that underlie these conserved patterns. Specifically, we confirm the non-randomness of recombination break points across all six tested coronavirus subgenera, locate conserved recombination hot- and cold-spots, and determine that the locations of transcriptional regulatory sequences are likely major determinants of conserved recombination break-point hotspot locations. We find that while the locations of recombination break points are not uniformly associated with degrees of nucleotide sequence conservation, they display significant tendencies in multiple coronavirus subgenera to occur in low guanine-cytosine content genome regions, in non-coding regions, at the edges of genes, and at sites within the Spike gene that are predicted to be minimally disruptive of Spike protein folding. While it is apparent that sequence features such as transcriptional regulatory sequences are likely major determinants of where the template-switching events that yield recombination break points most commonly occur, it is evident that selection against misfolded recombinant proteins also strongly impacts observable recombination break-point distributions in coronavirus genomes sampled from nature.

Virion-Associated Nucleic Acid-Based Metagenomics: A Decade of Advances in Molecular Characterization of Plant Viruses.

Over the last decade, viral metagenomic studies have resulted in the discovery of thousands of previously unknown viruses. These studies are likely to play a pivotal role in obtaining an accurate and robust understanding of how viruses affect the stability and productivity of ecosystems. Among the metagenomics-based approaches that have been developed since the beginning of the 21st century, shotgun metagenomics applied specifically to virion-associated nucleic acids (VANA) has been used to disentangle the diversity of the viral world. We summarize herein the results of 24 VANA-based studies, focusing on plant and insect samples conducted over the last decade (2010 to 2020). Collectively, viruses from 85 different families were reliably detected in these studies, including capsidless RNA viruses that replicate in fungi, oomycetes, and plants. Finally, strengths and weaknesses of the VANA approach are summarized and perspectives of applications in detection, epidemiological surveillance, environmental monitoring, and ecology of plant viruses are provided. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Detection of SARS-CoV-2 intra-host recombination during superinfection with Alpha and Epsilon variants in New York City.

Recombination is an evolutionary process by which many pathogens generate diversity and acquire novel functions. Although a common occurrence during coronavirus replication, detection of recombination is only feasible when genetically distinct viruses contemporaneously infect the same host. Here, we identify an instance of SARS-CoV-2 superinfection, whereby an individual was infected with two distinct viral variants: Alpha (B.1.1.7) and Epsilon (B.1.429). This superinfection was first noted when an Alpha genome sequence failed to exhibit the classic S gene target failure behavior used to track this variant. Full genome sequencing from four independent extracts reveals that Alpha variant alleles comprise around 75% of the genomes, whereas the Epsilon variant alleles comprise around 20% of the sample. Further investigation reveals the presence of numerous recombinant haplotypes spanning the genome, specifically in the spike, nucleocapsid, and ORF 8 coding regions. These findings support the potential for recombination to reshape SARS-CoV-2 genetic diversity.

Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa.

Three lineages (BA.1, BA.2 and BA.3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern predominantly drove South Africa's fourth Coronavirus Disease 2019 (COVID-19) wave. We have now identified two new lineages, BA.4 and BA.5, responsible for a fifth wave of infections. The spike proteins of BA.4 and BA.5 are identical, and similar to BA.2 except for the addition of 69-70 deletion (present in the Alpha variant and the BA.1 lineage), L452R (present in the Delta variant), F486V and the wild-type amino acid at Q493. The two lineages differ only outside of the spike region. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure, on the background of variants not possessing this feature. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa by the first week of April 2022. Using a multinomial logistic regression model, we estimated growth advantages for BA.4 and BA.5 of 0.08 (95% confidence interval (CI): 0.08-0.09) and 0.10 (95% CI: 0.09-0.11) per day, respectively, over BA.2 in South Africa. The continued discovery of genetically diverse Omicron lineages points to the hypothesis that a discrete reservoir, such as human chronic infections and/or animal hosts, is potentially contributing to further evolution and dispersal of the virus.

Recent Zoonotic Spillover and Tropism Shift of a Canine Coronavirus Is Associated with Relaxed Selection and Putative Loss of Function in NTD Subdomain of Spike Protein.

A canine coronavirus (CCoV) has now been reported from two independent human samples from Malaysia (respiratory, collected in 2017-2018; CCoV-HuPn-2018) and Haiti (urine, collected in 2017); these two viruses were nearly genetically identical. In an effort to identify any novel adaptations associated with this apparent shift in tropism we carried out detailed evolutionary analyses of the spike gene of this virus in the context of related Alphacoronavirus 1 species. The spike 0-domain retains homology to CCoV2b (enteric infections) and Transmissible Gastroenteritis Virus (TGEV; enteric and respiratory). This domain is subject to relaxed selection pressure and an increased rate of molecular evolution. It contains unique amino acid substitutions, including within a region important for sialic acid binding and pathogenesis in TGEV. Overall, the spike gene is extensively recombinant, with a feline coronavirus type II strain serving a prominent role in the recombinant history of the virus. Molecular divergence time for a segment of the gene where temporal signal could be determined, was estimated at around 60 years ago. We hypothesize that the virus had an enteric origin, but that it may be losing that particular tropism, possibly because of mutations in the sialic acid binding region of the spike 0-domain.

Emergence and phenotypic characterization of the global SARS-CoV-2 C.1.2 lineage.

Global genomic surveillance of SARS-CoV-2 has identified variants associated with increased transmissibility, neutralization resistance and disease severity. Here we report the emergence of the PANGO lineage C.1.2, detected at low prevalence in South Africa and eleven other countries. The initial C.1.2 detection is associated with a high substitution rate, and includes changes within the spike protein that have been associated with increased transmissibility or reduced neutralization sensitivity in SARS-CoV-2 variants of concern or variants of interest. Like Beta and Delta, C.1.2 shows significantly reduced neutralization sensitivity to plasma from vaccinees and individuals infected with the ancestral D614G virus. In contrast, convalescent donors infected with either Beta or Delta show high plasma neutralization against C.1.2. These functional data suggest that vaccine efficacy against C.1.2 will be equivalent to Beta and Delta, and that prior infection with either Beta or Delta will likely offer protection against C.1.2.

HIV-1 Evolutionary Dynamics under Nonsuppressive Antiretroviral Therapy.

Prolonged virologic failure on 2nd-line protease inhibitor (PI)-based antiretroviral therapy (ART) without emergence of major protease mutations is well recognized and provides an opportunity to study within-host evolution in long-term viremic individuals. Using next-generation sequencing and in silico haplotype reconstruction, we analyzed whole-genome sequences from longitudinal plasma samples of eight chronically infected HIV-1-positive individuals failing 2nd-line regimens from the French National Agency for AIDS and Viral Hepatitis Research (ANRS) 12249 Treatment as Prevention (TasP) trial. On nonsuppressive ART, there were large fluctuations in synonymous and nonsynonymous variant frequencies despite stable viremia. Reconstructed haplotypes provided evidence for selective sweeps during periods of partial adherence, and viral haplotype competition, during periods of low drug exposure. Drug resistance mutations in reverse transcriptase (RT) were used as markers of viral haplotypes in the reservoir, and their distribution over time indicated recombination. We independently observed linkage disequilibrium decay, indicative of recombination. These data highlight dramatic changes in virus population structure that occur during stable viremia under nonsuppressive ART. IMPORTANCE HIV-1 infections are most commonly initiated with a single founder virus and are characterized by extensive inter- and intraparticipant genetic diversity. However, existing literature on HIV-1 intrahost population dynamics is largely limited to untreated infections, predominantly in subtype B-infected individuals. The manuscript characterizes viral population dynamics in long-term viremic treatment-experienced individuals, which has not been previously characterized. These data are particularly relevant for understanding HIV dynamics but can also be applied to other RNA viruses. With this unique data set we propose that the virus is highly unstable, and we have found compelling evidence of HIV-1 within-host viral diversification, recombination, and haplotype competition during nonsuppressive ART.

Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function.

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Neutralisation sensitivity of the SARS-CoV-2 omicron (B.1.1.529) variant: a cross-sectional study.

BACKGROUND: The SARS-CoV-2 omicron (B.1.1.529) variant, which was first identified in November, 2021, spread rapidly in many countries, with a spike protein highly diverged from previously known variants, and raised concerns that this variant might evade neutralising antibody responses. We therefore aimed to characterise the sensitivity of the omicron variant to neutralisation. METHODS: For this cross-sectional study, we cloned the sequence encoding the omicron spike protein from a diagnostic sample to establish an omicron pseudotyped virus neutralisation assay. We quantified the neutralising antibody ID50 (the reciprocal dilution that produces 50% inhibition) against the omicron spike protein, and the fold-change in ID50 relative to the spike of wild-type SARS-CoV-2 (ie, the pandemic founder variant), for one convalescent reference plasma pool (WHO International Standard for anti-SARS-CoV-2 immunoglobulin [20/136]), three reference serum pools from vaccinated individuals, and two cohorts from Stockholm, Sweden: one comprising previously infected hospital workers (17 sampled in November, 2021, after vaccine rollout and nine in June or July, 2020, before vaccination) and one comprising serum from 40 randomly sampled blood donors donated during week 48 (Nov 29-Dec 5) of 2021. Furthermore, we assessed the neutralisation of omicron by five clinically relevant monoclonal antibodies (mAbs). FINDINGS: Neutralising antibody responses in reference sample pools sampled shortly after infection or vaccination were substantially less potent against the omicron variant than against wild-type SARS-CoV-2 (seven-fold to 42-fold reduction in ID50 titres). Similarly, for sera obtained before vaccination in 2020 from a cohort of convalescent hospital workers, neutralisation of the omicron variant was low to undetectable (all ID50 titres <20). However, in serum samples obtained in 2021 from two cohorts in Stockholm, substantial cross-neutralisation of the omicron variant was observed. Sera from 17 hospital workers after infection and subsequent vaccination had a reduction in average potency of only five-fold relative to wild-type SARS-CoV-2 (geometric mean ID50 titre 495 vs 105), and two donors had no reduction in potency. A similar pattern was observed in randomly sampled blood donors (n=40), who had an eight-fold reduction in average potency against the omicron variant compared with wild-type SARS-CoV-2 (geometric mean ID50 titre 369 vs 45). We found that the omicron variant was resistant to neutralisation (50% inhibitory concentration [IC50] >10 µg/mL) by mAbs casirivimab (REGN-10933), imdevimab (REGN-10987), etesevimab (Ly-CoV016), and bamlanivimab (Ly-CoV555), which form part of antibody combinations used in the clinic to treat COVID-19. However, S309, the parent of sotrovimab, retained most of its activity, with only an approximately two-fold reduction in potency against the omicron variant compared with ancestral D614G SARS-CoV-2 (IC50 0·1-0·2 µg/mL). INTERPRETATION: These data highlight the extensive, but incomplete, evasion of neutralising antibody responses by the omicron variant, and suggest that boosting with licensed vaccines might be sufficient to raise neutralising antibody titres to protective levels. FUNDING: European Union Horizon 2020 research and innovation programme, European and Developing Countries Clinical Trials Partnership, SciLifeLab, and the Erling-Persson Foundation.

Selection analysis identifies unusual clustered mutational changes in Omicron lineage BA.1 that likely impact Spike function.

Among the 30 non-synonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (i) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (ii) interactions of Spike with ACE2 receptors, and (iii) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any genomes within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron over all previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.