Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes

Impact of vaccination on SARS-CoV-2 evolution and immune escape variants

Vaccines and host genetic factors can influence the SARS-CoV-2 evolution and emergence of new variants. Even vaccinated cases get affected as virus continues to evolve, raising concerns about vaccine efficacy and the emergence of immune escape variants. Here, we have analyzed 2295 whole-genome sequences of SARS-CoV-2 collected from vaccinated and unvaccinated cases to evaluate the impact of vaccines on virus diversity within hosts. Our comparative analysis revealed a significant higher incidence of intra-host single nucleotides variants (iSNVs) in vaccinated cases compared to unvaccinated ones (p value<0.0001). Furthermore, we have found that specific mutational processes, including APOBEC (C > T) mediated and ADAR1 (A > G) mediated mutations, were found more prevalent in vaccinated cases. Vaccinated cases exhibited higher accumulation of nonsynonymous mutation than unvaccinated cases. Fixed iSNVs were predominantly located in the ORF1ab and spike genes, several key omicron defining immune escape variants S477N, Q493R, Q498R, Y505H, L452R, and N501Y were identified in the RBD domain of spike gene in vaccinated cases. Our findings suggest that vaccine plays an important role in the evolution of the virus genome. The virus genome acquires random mutations due to error-prone replication of the virus, host modification through APOBEC and ADAR1 mediated editing mechanism, and oxidative stress. These mutations become fixed in the viral population due to the selective pressure imposed by vaccination.

PHDtools: A platform for pathogen detection and multi-dimensional genetic signatures decoding to realize pathogen genomics data analyses online

OBJECTIVES

Facing the emerging diseases, rapid identification of the pathogen and multi-dimensional characterization of the genomic features at both isolate-level and population-level through high-throughput sequencing data can provide invaluable information to guide the development of antiviral agents and strategies. However, a user-friendly program is in urgent need for clinical laboratories without bioinformatics background to decode the complex big genomics data.

METHODS

In this study, we developed an interactive online platform named PHDtools with a total of 15 functions to analyze metagenomics data to identify the potential pathogen and decode multi-dimensional genetic signatures including intra-/inter-host variations and lineage-level variations. The platform was applied to analyze the meta-genomic data of the samples collected from the 172 imported COVID-19 cases.

RESULTS

According to the analytical results of mNGS, 27 patients were found to have the co-infections of SARS-CoV-2 with either influenza virus (n = 9) or human picobirnavirus (n = 19). Enough coverages of all the assembled SARS-CoV-2 genomes provided the sub-lineages of Omicron variant, and the number of mutations in the non-structural genes and M gene was increased, as well as the intra-host variations occurred in E and M gene were under positive selection (Ka/Ks > 1). These findings of increased or changed mutations in the SARS-CoV-2 genome characterized the current adaptive evolution patterns of Omicron sub-lineages, and revealed the evolution speed of these sub-lineages might increase.

CONCLUSIONS

Consequently, the application of PHDtools has proved that this platform is accurate, user-friendly and convenient for clinical users who are deficient in bioinformatics, and the clinical monitor of SARS-CoV-2 genomes by PHDtools also highlighted the potential evolution features of current SARS-CoV-2 and indicated that the development of anti-SARS-CoV-2 agents and new-designed vaccines should incorporate the gene variations other than S gene.

Rapid identification of full-length genome and tracing variations of monkeypox virus in clinical specimens based on mNGS and amplicon sequencing

The monkeypox virus (MPXV) has triggered a current outbreak globally. Genome sequencing of MPXV and rapid tracing of genetic variants will benefit disease diagnosis and control. It is a significant challenge but necessary to optimize the strategy and application of rapid full-length genome identification and to track variations of MPXV in clinical specimens with low viral loads, as it is one of the DNA viruses with the largest genome and the most AT-biased, and has a significant number of tandem repeats. Here we evaluated the performance of metagenomic and amplicon sequencing techniques, and three sequencing platforms in MPXV genome sequencing based on multiple clinical specimens of five mpox cases in Chinese mainland. We rapidly identified the full-length genome of MPXV with the assembly of accurate tandem repeats in multiple clinical specimens. Amplicon sequencing enables cost-effective and rapid sequencing of clinical specimens to obtain high-quality MPXV genomes. Third-generation sequencing facilitates the assembly of the terminal tandem repeat regions in the monkeypox virus genome and corrects a common misassembly in published sequences. Besides, several intra-host single nucleotide variations were identified in the first imported mpox case. This study offers an evaluation of various strategies aimed at identifying the complete genome of MPXV in clinical specimens. The findings of this study will significantly enhance the surveillance of MPXV.

Puzzles, challenges, and information reservoir of SARS-CoV-2 quasispecies

Upon the emergence of SARS-CoV-2 in the human population, it was conjectured that for this coronavirus the dynamic intra-host heterogeneity typical of RNA viruses would be toned down. Nothing of this sort is observed. Here we review the main observations on the complexity and diverse composition of SARS-CoV-2 mutant spectra sampled from infected patients, within the framework of quasispecies dynamics. The analyses suggest that the information provided by myriads of genomic sequences within infected individuals may have a predictive value of the genomic sequences that acquire epidemiological relevance. Possibilities to reconcile the presence of broad mutant spectra in the large RNA coronavirus genome with its encoding a 3' to 5' exonuclease proofreading-repair activity are considered. Indeterminations in the behavior of individual viral genomes provide a benefit for the survival of the ensemble. We propose that this concept falls in the domain of "stochastic thinking," a notion that applies also to cellular processes, as a means for biological systems to face unexpected needs.

SARS-CoV-2 within-host diversity of human hosts and its implications for viral immune evasion

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is continuously evolving, bringing great challenges to the control of the virus. In the present study, we investigated the characteristics of SARS-CoV-2 within-host diversity of human hosts and its implications for immune evasion using about 2,00,000 high-depth next-generation genome sequencing data of SARS-CoV-2. A total of 44% of the samples showed within-host variations (iSNVs), and the average number of iSNVs in the samples with iSNV was 1.90. C-to-U is the dominant substitution pattern for iSNVs. C-to-U/G-to-A and A-to-G/U-to-C preferentially occur in 5'-CG-3' and 5'-AU-3' motifs, respectively. In addition, we found that SARS-CoV-2 within-host variations are under negative selection. About 15.6% iSNVs had an impact on the content of the CpG dinucleotide (CpG) in SARS-CoV-2 genomes. We detected signatures of faster loss of CpG-gaining iSNVs, possibly resulting from zinc-finger antiviral protein-mediated antiviral activities targeting CpG, which could be the major reason for CpG depletion in SARS-CoV-2 consensus genomes. The non-synonymous iSNVs in the S gene can largely alter the S protein's antigenic features, and many of these iSNVs are distributed in the amino-terminal domain (NTD) and receptor-binding domain (RBD). These results suggest that SARS-CoV-2 interacts actively with human hosts and attempts to take different evolutionary strategies to escape human innate and adaptive immunity. These new findings further deepen and widen our understanding of the within-host evolutionary features of SARS-CoV-2. IMPORTANCE Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative pathogen of the coronavirus disease 2019, has evolved rapidly since it was discovered. Recent studies have pointed out that some mutations in the SARS-CoV-2 S protein could confer SARS-CoV-2 the ability to evade the human adaptive immune system. In addition, it is observed that the content of the CpG dinucleotide in SARS-CoV-2 genome sequences has decreased over time, reflecting the adaptation to the human host. The significance of our research is revealing the characteristics of SARS-CoV-2 within-host diversity of human hosts, identifying the causes of CpG depletion in SARS-CoV-2 consensus genomes, and exploring the potential impacts of non-synonymous within-host variations in the S gene on immune escape, which could further deepen and widen our understanding of the evolutionary features of SARS-CoV-2.

Nanopore sequencing technology and its applications

Since the development of Sanger sequencing in 1977, sequencing technology has played a pivotal role in molecular biology research by enabling the interpretation of biological genetic codes. Today, nanopore sequencing is one of the leading third-generation sequencing technologies. With its long reads, portability, and low cost, nanopore sequencing is widely used in various scientific fields including epidemic prevention and control, disease diagnosis, and animal and plant breeding. Despite initial concerns about high error rates, continuous innovation in sequencing platforms and algorithm analysis technology has effectively addressed its accuracy. During the coronavirus disease (COVID-19) pandemic, nanopore sequencing played a critical role in detecting the severe acute respiratory syndrome coronavirus-2 virus genome and containing the pandemic. However, a lack of understanding of this technology may limit its popularization and application. Nanopore sequencing is poised to become the mainstream choice for preventing and controlling COVID-19 and future epidemics while creating value in other fields such as oncology and botany. This work introduces the contributions of nanopore sequencing during the COVID-19 pandemic to promote public understanding and its use in emerging outbreaks worldwide. We discuss its application in microbial detection, cancer genomes, and plant genomes and summarize strategies to improve its accuracy.

Within-host genetic diversity of SARS-CoV-2 lineages in unvaccinated and vaccinated individuals

Viral and host factors can shape SARS-CoV-2 evolution. However, little is known about lineage-specific and vaccination-specific mutations that occur within individuals. Here, we analysed deep sequencing data from 2,820 SARS-CoV-2 respiratory samples with different viral lineages to describe the patterns of within-host diversity under different conditions, including vaccine-breakthrough infections. In unvaccinated individuals, variant of Concern (VOC) Alpha, Delta, and Omicron respiratory samples were found to have higher within-host diversity and were under neutral to purifying selection at the full genome level compared to non-VOC SARS-CoV-2. Breakthrough infections in 2-dose or 3-dose Comirnaty and CoronaVac vaccinated individuals did not increase levels of non-synonymous mutations and did not change the direction of selection pressure. Vaccine-induced antibody or T cell responses did not appear to have significant impact on within-host SARS-CoV-2 sequence diversification. Our findings suggest that vaccination does not increase exploration of SARS-CoV-2 protein sequence space and may not facilitate emergence of viral variants.

The race to understand immunopathology in COVID-19: Perspectives on the impact of quantitative approaches to understand within-host interactions

The COVID-19 pandemic has revealed the need for the increased integration of modelling and data analysis to public health, experimental, and clinical studies. Throughout the first two years of the pandemic, there has been a concerted effort to improve our understanding of the within-host immune response to the SARS-CoV-2 virus to provide better predictions of COVID-19 severity, treatment and vaccine development questions, and insights into viral evolution and the impacts of variants on immunopathology. Here we provide perspectives on what has been accomplished using quantitative methods, including predictive modelling, population genetics, machine learning, and dimensionality reduction techniques, in the first 26 months of the COVID-19 pandemic approaches, and where we go from here to improve our responses to this and future pandemics.

Distribution of intra-host variations and mutations in the genomes of SARS-CoV-2 and their implications on detection and therapeutics

The ongoing circulation of SARS-CoV-2 variants of concern (VOCs) has caused global concerns, because VOCs could escape current vaccines, antiviral drugs, and diagnosis. Analyzing mutations and intra-host diversities in different and widespread VOCs can provide important insights to virus adaptive evolution and validity of vaccines, antiviral drugs, and diagnosis. In this study, by analyzing 1744 high-throughput sequencing data for intra-host single-nucleotide variations (iSNVs) and 3,668,205 genome sequences for mutations in different VOCs, it was found that Omicron variant is still evolving at high speed, especially having high iSNVs frequency in its S and N genes. The efficacies of antibodies or detection primers targeting these two genes are at high risks to be invalid. Instead, highly conserved regions such as NSP8 gene could be better therapeutic and detection targets. Furthermore, mutations in later VOCs could be traced to the minor alleles in the previous variant samples such as Alpha and Delta in different countries. Finally, it was found that mutations C14408T in RdRp and A18163G in NSP14 gene might be associated with the higher genetic diversity in Omicron. Our findings not only contribute to understanding the adaptive evolution of SARS-CoV-2 VOCs, but also provide useful information for both drugs and diagnostic kits development.

Within-host genetic diversity of SARS-CoV-2 in the context of large-scale hospital-associated genomic surveillance

The COVID-19 pandemic has resulted in extensive surveillance of the genomic diversity of SARS-CoV-2. Sequencing data generated as part of these efforts can also capture the diversity of the SARS-CoV-2 virus populations replicating within infected individuals. To assess this within-host diversity of SARS-CoV-2 we quantified low frequency (minor) variants from deep sequence data of thousands of clinical samples collected by a large urban hospital system over the course of a year. Using a robust analytical pipeline to control for technical artefacts, we observe that at comparable viral loads, specimens from patients hospitalized due to COVID-19 had a greater number of minor variants than samples from outpatients. Since individuals with highly diverse viral populations could be disproportionate drivers of new viral lineages in the patient population, these results suggest that transmission control should pay special attention to patients with severe or protracted disease to prevent the spread of novel variants.

Heterologous immunity induced by 1st generation COVID-19 vaccines and its role in developing a pan-coronavirus vaccine

Severe acute respiratory syndrome virus-2 (SARS-CoV-2), the causative infectious agent of the COVID-19 pandemic, has led to multiple (4-6) waves of infections worldwide during the past two years. The development of vaccines against SARS-CoV-2 has led to successful mass immunizations worldwide, mitigating the worldwide mortality due the pandemic to a great extent. Yet the evolution of new variants highlights a need to develop a universal vaccine which can prevent infections from all virulent SARS-CoV-2. Most of the current first generation COVID-19 vaccines are based on the Spike protein from the original Wuhan-hu-1 virus strain. It is encouraging that they still protect from serious illnesses, hospitalizations and mortality against a number of mutated viral strains, to varying degrees. Understanding the mechanisms by which these vaccines provide heterologous protection against multiple highly mutated variants can reveal strategies to develop a universal vaccine. In addition, many unexposed individuals have been found to harbor T cells that are cross-reactive against SARS-CoV-2 antigens, with a possible protective role. In this review, we will discuss various aspects of natural or vaccine-induced heterologous (cross-reactive) adaptive immunity against SARS-CoV-2 and other coronaviruses, and their role in achieving the concept of a pan-coronavirus vaccine.

Within-host diversity of SARS-CoV-2 lineages and effect of vaccination

Viral and host factors can shape SARS-CoV-2 within-host viral diversity and virus evolution. However, little is known about lineage-specific and vaccination-specific mutations that occur within individuals. Here we analysed deep sequencing data from 2,146 SARS-CoV-2 samples with different viral lineages to describe the patterns of within-host diversity in different conditions, including vaccine-breakthrough infections. Variant of Concern (VOC) Alpha, Delta, and Omicron samples were found to have higher within-host nucleotide diversity while being under weaker purifying selection at full genome level compared to non-VOC SARS-CoV-2 viruses. Breakthrough Delta and Omicron infections in Comirnaty and CoronaVac vaccinated individuals appeared to have higher within-host purifying selection at the full-genome and/or Spike gene levels. Vaccine-induced antibody or T cell responses did not appear to have significant impact on within-host SARS-CoV-2 evolution. Our findings suggest that vaccination does not increase SARS-CoV-2 protein sequence space and may not facilitate emergence of more viral variants.

How SARS-CoV-2 dodges immune surveillance and facilitates infection: an analytical review

INTRODUCTION

Effective treatments for the ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic are limited. The virus has evolved strategies to evade the immune system or hijack immune responses to facilitate infection and escape immune surveillance. Mechanistically, SARS-CoV-2 takes advantage of TLR4 and cytokine-induced integrins to promote its entrance into the cell. Furthermore, the activation of pattern recognition receptors (PRR)-mediated signaling pathways is compromised by SARS-CoV-2 non-structural proteins (NSPs), accessory protein open reading frames (ORFs), and structural proteins upon infection, contributing to viral infection and replication. Host factors necessary for cellular protein synthesis, metabolism, and viral replication can also be inhibited by the SARS-CoV-2 proteins. Exploring specific mechanisms would optimize the therapy methods and benefit drug research and development.

AREAS COVERED

: We describe pathways and mechanisms by which SARS-CoV-2 evades immune system; these include the mechanisms that operate during virus entry, signaling pathways involved, and processes at RNA and protein levels.

EXPERT OPINION

: Increased understanding of how viruses interfere with immune responses would provide more evidence for drug development. Drugs targeting conserved viral proteins to inhibit their replication or host factors to enhance immune responses would minimize the impact of virus mutations and prepare for future coronavirus outbreaks.

Rescuing low frequency variants within intra-host viral populations directly from Oxford Nanopore sequencing data

Infectious disease monitoring on Oxford Nanopore Technologies (ONT) platforms offers rapid turnaround times and low cost. Tracking low frequency intra-host variants provides important insights with respect to elucidating within-host viral population dynamics and transmission. However, given the higher error rate of ONT, accurate identification of intra-host variants with low allele frequencies remains an open challenge with no viable computational solutions available. In response to this need, we present Variabel, a novel approach and first method designed for rescuing low frequency intra-host variants from ONT data alone. We evaluate Variabel on both synthetic data (SARS-CoV-2) and patient derived datasets (Ebola virus, norovirus, SARS-CoV-2); our results show that Variabel can accurately identify low frequency variants below 0.5 allele frequency, outperforming existing state-of-the-art ONT variant callers for this task. Variabel is open-source and available for download at: www.gitlab.com/treangenlab/variabel.

Tracking low frequency intra-host variants has helped understanding within-host viral population dynamics and transmission. Precise tracking, however, depends partially on the error rate of the sequencing platforms used. Here, Liu et al. present Variabel, a method to rescue low frequency intra-host variants from Oxford Nanopore Technologies (ONT) platforms and validate their approach on Ebola virus, norovirus, and SARS-CoV-2 datasets.