Global SNP analysis of 11,183 SARS-CoV-2 strains reveals high genetic diversity

Immunoinformatics design of a structural proteins driven multi-epitope candidate vaccine against different SARS-CoV-2 variants based on fynomer

The ideal vaccines for combating diseases that may emerge in the future require more than simply inactivating a few pathogenic strains. This study aims to provide a peptide-based multi-epitope vaccine effective against various severe acute respiratory syndrome coronavirus 2 strains. To design the vaccine, a library of peptides from the spike, nucleocapsid, membrane, and envelope structural proteins of various strains was prepared. Then, the final vaccine structure was optimized using the fully protected epitopes and the fynomer scaffold. Using bioinformatics tools, the antigenicity, allergenicity, toxicity, physicochemical properties, population coverage, and secondary and three-dimensional structures of the vaccine candidate were evaluated. The bioinformatic analyses confirmed the high quality of the vaccine. According to further investigations, this structure is similar to native protein and there is a stable and strong interaction between vaccine and receptors. Based on molecular dynamics simulation, structural compactness and stability in binding were also observed. In addition, the immune simulation showed that the vaccine can stimulate immune responses similar to real conditions. Finally, codon optimization and in silico cloning confirmed efficient expression in Escherichia coli. In conclusion, the fynomer-based vaccine can be considered as a new style in designing and updating vaccines to protect against coronavirus disease.

Design and Development of an Antigen Test for SARS-CoV-2 Nucleocapsid Protein to Validate the Viral Quality Assurance Panels

The continuing mutability of the SARS-CoV-2 virus can result in failures of diagnostic assays. To address this, we describe a generalizable bioinformatics-to-biology pipeline developed for the calibration and quality assurance of inactivated SARS-CoV-2 variant panels provided to Radical Acceleration of Diagnostics programs (RADx)-radical program awardees. A heuristic genetic analysis based on variant-defining mutations demonstrated the lowest genetic variance in the Nucleocapsid protein (Np)-C-terminal domain (CTD) across all SARS-CoV-2 variants. We then employed the Shannon entropy method on (Np) sequences collected from the major variants, verifying the CTD with lower entropy (less prone to mutations) than other Np regions. Polyclonal and monoclonal antibodies were raised against this target CTD antigen and used to develop an Enzyme-linked immunoassay (ELISA) test for SARS-CoV-2. Blinded Viral Quality Assurance (VQA) panels comprised of UV-inactivated SARS-CoV-2 variants (XBB.1.5, BF.7, BA.1, B.1.617.2, and WA1) and distractor respiratory viruses (CoV 229E, CoV OC43, RSV A2, RSV B, IAV H1N1, and IBV) were assembled by the RADx-rad Diagnostics core and tested using the ELISA described here. The assay tested positive for all variants with high sensitivity (limit of detection: 1.72-8.78 ng/mL) and negative for the distractor virus panel. Epitope mapping for the monoclonal antibodies identified a 20 amino acid antigenic peptide on the Np-CTD that an in-silico program also predicted for the highest antigenicity. This work provides a template for a bioinformatics pipeline to select genetic regions with a low propensity for mutation (low Shannon entropy) to develop robust 'pan-variant' antigen-based assays for viruses prone to high mutational rates.

Sequence analysis of the Spike, RNA-dependent RNA polymerase, and protease genes reveals a distinct evolutionary pattern of SARS-CoV-2 variants circulating in Yogyakarta and Central Java provinces, Indonesia

During the Covid-19 pandemic, the resurgence of SARS-CoV-2 was due to the development of novel variants of concern (VOC). Thus, genomic surveillance is essential to monitor continuing evolution of SARS-CoV-2 and to track the emergence of novel variants. In this study, we performed phylogenetic, mutation, and selection pressure analyses of the Spike, nsp12, nsp3, and nsp5 genes of SARS-CoV-2 isolates circulating in Yogyakarta and Central Java provinces, Indonesia from May 2021 to February 2022. Various bioinformatics tools were employed to investigate the evolutionary dynamics of distinct SARS-CoV-2 isolates. During the study period, 213 and 139 isolates of Omicron and Delta variants were identified, respectively. Particularly in the Spike gene, mutations were significantly more abundant in Omicron than in Delta variants. Consistently, in all of four genes studied, the substitution rates of Omicron were higher than that of Delta variants, especially in the Spike and nsp12 genes. In addition, selective pressure analysis revealed several sites that were positively selected in particular genes, implying that these sites were functionally essential for virus evolution. In conclusion, our study demonstrated a distinct evolutionary pattern of SARS-CoV-2 variants circulating in Yogyakarta and Central Java provinces, Indonesia.

Significance of understanding the genomics of host-pathogen interaction in limiting antibiotic resistance development: lessons from COVID-19 pandemic

The entire world is facing the stiff challenge of COVID-19 pandemic. To overcome the spread of this highly infectious disease, several short-sighted strategies were adopted such as the use of broad-spectrum antibiotics and antifungals. However, the misuse and/or overuse of antibiotics have accentuated the emergence of the next pandemic: antimicrobial resistance (AMR). It is believed that pathogens while transferring between humans and the environment carry virulence and antibiotic-resistant factors from varied species. It is presumed that all such genetic factors are quantifiable and predictable, a better understanding of which could be a limiting step for the progression of AMR. Herein, we have reviewed how genomics-based understanding of host-pathogen interactions during COVID-19 could reduce the non-judicial use of antibiotics and prevent the eruption of an AMR-based pandemic in future.

Some common deleterious mutations are shared in SARS-CoV-2 genomes from deceased COVID-19 patients across continents

The identification of deleterious mutations in different variants of SARS-CoV-2 and their roles in the morbidity of COVID-19 patients has yet to be thoroughly investigated. To unravel the spectrum of mutations and their effects within SARS-CoV-2 genomes, we analyzed 5,724 complete genomes from deceased COVID-19 patients sourced from the GISAID database. This analysis was conducted using the Nextstrain platform, applying a generalized time-reversible model for evolutionary phylogeny. These genomes were compared to the reference strain (hCoV-19/Wuhan/WIV04/2019) using MAFFT v7.470. Our findings revealed that SARS-CoV-2 genomes from deceased individuals belonged to 21 Nextstrain clades, with clade 20I (Alpha variant) being the most predominant, followed by clade 20H (Beta variant) and clade 20J (Gamma variant). The majority of SARS-CoV-2 genomes from deceased patients (33.4%) were sequenced in North America, while the lowest percentage (0.98%) came from Africa. The 'G' clade was dominant in the SARS-CoV-2 genomes of Asian, African, and North American regions, while the 'GRY' clade prevailed in Europe. In our analysis, we identified 35,799 nucleotide (NT) mutations throughout the genome, with the highest frequency (11,402 occurrences) found in the spike protein. Notably, we observed 4150 point-specific amino acid (AA) mutations in SARS-CoV-2 genomes, with D614G (20%) and N501Y (14%) identified as the top two deleterious mutations in the spike protein on a global scale. Furthermore, we detected five common deleterious AA mutations, including G18V, W45S, I33T, P30L, and Q418H, which play a key role in defining each clade of SARS-CoV-2. Our novel findings hold potential value for genomic surveillance, enabling the monitoring of the evolving pattern of SARS-CoV-2 infection, its emerging variants, and their impact on the development of effective vaccination and control strategies.

Detection of SARS-CoV-2 Δ426 ORF8 Deletion Mutant Cluster in NGS Screening

Next-generation sequencing (NGS) from SARS-CoV-2-positive swabs collected during the last months of 2022 revealed a large deletion spanning ORF7b and ORF8 (426 nt) in six patients infected with the BA.5.1 Omicron variant. This extensive genome loss removed a large part of these two genes, maintaining in frame the first 22 aminoacids of ORF7b and the last three aminoacids of ORF8. Interestingly, the deleted region was flanked by two small repeats, which were likely involved in the formation of a hairpin structure. Similar rearrangements, comparable in size and location to the deletion, were also identified in 15 sequences in the NCBI database. In this group, seven out of 15 cases from the USA and Switzerland presented both the BA.5.1 variant and the same 426 nucleotides deletion. It is noteworthy that three out of six cases were detected in patients with immunodeficiency, and it is conceivable that this clinical condition could promote the replication and selection of these mutations.

Transient SARS-CoV-2 RNA-Dependent RNA Polymerase Mutations after Remdesivir Treatment for Chronic COVID-19 in Two Transplant Recipients: Case Report and Intra-Host Viral Genomic Investigation

Remdesivir is the first FDA-approved drug for treating severe SARS-CoV-2 infection and targets RNA-dependent RNA polymerase (RdRp) that is required for viral replication. To monitor for the development of mutations that may result in remdesivir resistance during prolonged treatment, we sequenced SARS-CoV-2 specimens collected at different treatment time points in two transplant patients with severe COVID-19. In the first patient, an allogeneic hematopoietic stem cell transplant recipient, a transient RdRp catalytic subunit mutation (nsp12:A449V) was observed that has not previously been associated with remdesivir resistance. As no in vitro study had been conducted to elucidate the phenotypic effect of nsp12:A449V, its clinical significance is unclear. In the second patient, two other transient RdRp mutations were detected: one in the catalytic subunit (nsp12:V166A) and the other in an accessory subunit important for processivity (nsp7:D67N). This is the first case report for a potential link between the nsp12:V166A mutation and remdesivir resistance in vivo, which had only been previously described by in vitro studies. The nsp7:D67N mutation has not previously been associated with remdesivir resistance, and whether it has a phenotypic effect is unknown. Our study revealed SARS-CoV-2 genetic dynamics during remdesivir treatment in transplant recipients that involved mutations in the RdRp complex (nsp7 and nsp12), which may be the result of selective pressure. These results suggest that close monitoring for potential resistance during the course of remdesivir treatment in highly vulnerable patient populations may be beneficial. Development and utilization of diagnostic RdRp genotyping tests may be a future direction for improving the management of chronic COVID-19.

Genomic Epidemiology of SARS-CoV-2 in Urban Settings in Senegal

We used whole genome sequencing to identify and analyze mutations in SARS-CoV-2 in urban settings during the deadliest wave of the COVID-19 epidemic-from March to April 2021-in Senegal. Nasopharyngeal samples testing positive for SARS-CoV-2 were sequenced on the Illumina NovaSeq 6000 sequencing system using the COVIDSeq protocol. A total of 291 genotypable consensus genome sequences were obtained. Phylogenetic analyses grouped the genomes into 16 distinct PANGOLIN lineages. The major lineage was B.1.1.420, despite circulation of the Alpha variant of concern (VOC). A total of 1125 different SNPs, identified relative to the Wuhan reference genome, were detected. These included 13 SNPs in non-coding regions. An average density of 37.2 SNPs per 1000 nucleotides was found, with the highest density occurring in ORF10. This analysis allowed, for the first time, the detection of a Senegalese SARS-CoV-2 strain belonging to the P.1.14 (GR/20J, Gamma V3) sublineage of the Brazilian P.1 lineage (or Gamma VOC). Overall, our results highlight substantial SARS-CoV-2 diversification in Senegal during the study period.

Antitarget, Anti-SARS-CoV-2 Leads, Drugs, and the Drug Discovery-Genetics Alliance Perspective

The most advanced antiviral molecules addressing major SARS-CoV-2 targets (Main protease, Spike protein, and RNA polymerase), compared with proteins of other human pathogenic coronaviruses, may have a short-lasting clinical efficacy. Accumulating knowledge on the mechanisms underlying the target structural basis, its mutational progression, and the related biological significance to virus replication allows envisaging the development of better-targeted therapies in the context of COVID-19 epidemic and future coronavirus outbreaks. The identification of evolutionary patterns based solely on sequence information analysis for those targets can provide meaningful insights into the molecular basis of host-pathogen interactions and adaptation, leading to drug resistance phenomena. Herein, we will explore how the study of observed and predicted mutations may offer valuable suggestions for the application of the so-called "synthetic lethal" strategy to SARS-CoV-2 Main protease and Spike protein. The synergy between genetics evidence and drug discovery may prioritize the development of novel long-lasting antiviral agents.

SARS-CoV-2 mutations on diagnostic gene targets in the second wave in Zimbabwe: A retrospective genomic analysis

BACKGROUND

SARS-CoV-2 continues to be a major issue in resource-limited settings, particularly owing to the limited supply of vaccinescaused by inequitable distribution.

OBJECTIVE

To monitor diagnostic gene targets to identify potential test failures caused by mutations, which is important for public health.

METHODS

Here we analysed the genome sequence of SARS-CoV-2 from the second wave in Zimbabwe. A total of 377 samples weresequenced at Quadram Institute Bioscience. After quality control, 192 sequences passed and were analysed.

RESULTS

The Beta variant was dominant during this period, contributing 77.6% (149) of the genomes sequenced and having a total of 2994mutations in diagnostic polymerase chain reaction target genes. Many single nucleotide polymorphism mutations resulted in amino acidsubstitution that had the potential to impact viral fitness by increasing the rate of transmission or evading the immune response to previousinfection or vaccination.

CONCLUSION

There were nine lineages circulating in Zimbabwe during the second wave. The B.1.351 was dominant, accounting for >75%.There were over 3 000 mutations on the diagnostic genes and lineage B.1.351, contributing almost two-thirds of the mutations. The S-genehad the most mutations and the E-gene was the least mutated.

Computational Analysis Predicts Correlations among Amino Acids in SARS-CoV-2 Proteomes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a serious global challenge requiring urgent and permanent therapeutic solutions. These solutions can only be engineered if the patterns and rate of mutations of the virus can be elucidated. Predicting mutations and the structure of proteins based on these mutations have become necessary for early drug and vaccine design purposes in anticipation of future viral mutations. The amino acid composition (AAC) of proteomes and individual viral proteins provide avenues for exploitation since AACs have been previously used to predict structure, shape and evolutionary rates. Herein, the frequency of amino acid residues found in 1637 complete proteomes belonging to 11 SARS-CoV-2 variants/lineages were analyzed. Leucine is the most abundant amino acid residue in the SARS-CoV-2 with an average AAC of 9.658% while tryptophan had the least abundance of 1.11%. The AAC and ranking of lysine and glycine varied in the proteome. For some variants, glycine had higher frequency and AAC than lysine and vice versa in other variants. Tryptophan was also observed to be the most intolerant to mutation in the various proteomes for the variants used. A correlogram revealed a very strong correlation of 0.999992 between B.1.525 (Eta) and B.1.526 (Iota) variants. Furthermore, isoleucine and threonine were observed to have a very strong negative correlation of -0.912, while cysteine and isoleucine had a very strong positive correlation of 0.835 at p < 0.001. Shapiro-Wilk normality test revealed that AAC values for all the amino acid residues except methionine showed no evidence of non-normality at p < 0.05. Thus, AACs of SARS-CoV-2 variants can be predicted using probability and z-scores. AACs may be beneficial in classifying viral strains, predicting viral disease types, members of protein families, protein interactions and for diagnostic purposes. They may also be used as a feature along with other crucial factors in machine-learning based algorithms to predict viral mutations. These mutation-predicting algorithms may help in developing effective therapeutics and vaccines for SARS-CoV-2.

Genomic Evidence of Multiple Introductions of SARS-CoV-2 in Mauritania

The rapid and global spread of the novel coronavirus severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has raised serious public health concerns, including in Mauritania. We sequenced and analyzed the entire genome of 13 SARS-CoV-2 virus strains isolated from polymerase chain reaction (PCR)-positive symptomatic patients sampled from March 3 to May 31, 2021 to better understand SARS-CoV-2 introduction, propagation, and evolution in Mauritania. A phylogenetic tree using available data from the EpiCoV GISAID database and a variant network with non-Mauritanian sequences were constructed. Variant analysis of the 13 Mauritanian SARS-CoV-2 genome sequences indicated an average mutational percentage of 0.39, which is similar to that in other countries. Phylogenetic analysis revealed multiple spatiotemporal introductions, mainly from Europe (France, Belgium) and Africa (Senegal, Côte d'Ivoire), which also provided evidence of early community transmission. A total of 2 unique mutations, namely, NSP6_Q208K and NSP15_S273T, were detected in the NSP6 and NSP15 genes, respectively, confirming the aforementioned introduction of SARS-CoV-2 in Mauritania. These findings highlight the relevance of continuous genomic monitoring strategies for understanding virus transmission dynamics and acquiring knowledge to address forthcoming sources of infection in Africa.

Bioinformatics pipeline unveils genetic variability to synthetic vaccine design for Indian SARS-CoV-2 genomes

In the worrisome scenarios of various waves of SARS-CoV-2 pandemic, a comprehensive bioinformatics pipeline is essential to analyse the virus genomes in order to understand its evolution, thereby identifying mutations as signature SNPs, conserved regions and subsequently to design epitope based synthetic vaccine. We have thus performed multiple sequence alignment of 4996 Indian SARS-CoV-2 genomes as a case study using MAFFT followed by phylogenetic analysis using Nextstrain to identify virus clades. Furthermore, based on the entropy of each genomic coordinate of the aligned sequences, conserved regions are identified. After refinement of the conserved regions, based on its length, one conserved region is identified for which the primers and probes are reported for virus detection. The refined conserved regions are also used to identify T-cell and B-cell epitopes along with their immunogenic and antigenic scores. Such scores are used for selecting the most immunogenic and antigenic epitopes. By executing this pipeline, 40 unique signature SNPs are identified resulting in 23 non-synonymous signature SNPs which provide 28 amino acid changes in protein. On the other hand, 12 conserved regions are selected based on refinement criteria out of which one is selected as the potential target for virus detection. Additionally, 22 MHC-I and 21 MHC-II restricted T-cell epitopes with 10 unique HLA alleles each and 17 B-cell epitopes are obtained for 12 conserved regions. All the results are validated both quantitatively and qualitatively which show that from genetic variability to synthetic vaccine design, the proposed pipeline can be used effectively to combat SARS-CoV-2.

Cellulosic copper nanoparticles and a hydrogen peroxide-based disinfectant trigger rapid inactivation of pseudoviral particles expressing the Spike protein of SARS-CoV-2, SARS-CoV, and MERS-CoV

Severe acute respiratory syndrome (SARS) is a viral respiratory infection caused by human coronaviruses that include SARS-CoV-2, SARS-CoV, and Middle East respiratory syndrome coronavirus (MERS-CoV). Although their primary mode of transmission is through contaminated respiratory droplets from infected carriers, the deposition of expelled virus particles onto surfaces and fomites could contribute to viral transmission. Here, we use replication-deficient murine leukemia virus (MLV) pseudoviral particles expressing SARS-CoV-2, SARS-CoV, or MERS-CoV Spike (S) protein on their surface. These surrogates of native coronavirus counterparts serve as a model to analyze the S-mediated entry into target cells. Carboxymethyl cellulose (CMC) nanofibers that are combined with copper (Cu) exhibit strong antimicrobial properties. S-pseudovirions that are exposed to CMC-Cu nanoparticles (30 s) display a dramatic reduction in their ability to infect target Vero E6 cells, with ∼97% less infectivity as compared to untreated pseudovirions. In contrast, addition of the Cu chelator tetrathiomolybdate protects S-pseudovirions from CMC-Cu-mediated inactivation. When S-pseudovirions were treated with a hydrogen peroxide-based disinfectant (denoted SaberTM) used at 1:250 dilution, their infectivity was dramatically reduced by ∼98%. However, the combined use of SaberTM and CMC-Cu is the most effective approach to restrict infectivity of SARS-CoV-2-S, SARS-CoV-S, and MERS-CoV-S pseudovirions in Vero E6 cell assays. Together, these results show that cellulosic Cu nanoparticles enhance the effectiveness of diluted SaberTM sanitizer, setting up an improved strategy to lower the risk of surface- and fomite-mediated transmission of enveloped respiratory viruses.

Shielding reduced incidence of COVID-19 in patients with inflammatory arthritis but vulnerability is associated with increased mortality

OBJECTIVES

Investigate whether individuals with inflammatory arthritis (IA), their treatments and shielding status affect the risk of adverse outcomes from COVID-19 for the entire population of Wales, UK.

METHODS

Retrospective, population-based cohort study using linked, anonymized electronic health data from SAIL Databank, including primary/secondary care, rheumatology, Office for National Statistics Mortality and COVID-19 laboratory data. Individuals aged 18 years and over testing positive for COVID-19 between March 2020 and May 2021 with READ Codes present for rheumatoid arthritis, psoriatic arthritis and ankylosing spondylitis formed the study cases.

RESULTS

A total of 1966 people with IA and 166 602 without tested positive for COVID-19. The incidence rate was 3.5% (1966/56 914) in IA, vs 6% in the general population (166 602/2 760 442), (difference: 2.5%, 95% CI: 2.4%, 2.7%, P ≤0.001). In an adjusted Cox proportional hazard model, IA was not associated with higher mortality (HR: 0.56, 95% CI: 0.18, 1.64, P=0.286). Significant risk factors included shielding (HR: 1.52, 95% CI: 1.40, 1.64, P ≤0.001), hospitalization for previous infections (HR: 1.20, 95% CI: 1.12, 1.28, P ≤0.001), hospitalizations one year pre-pandemic (HR: 1.34, 95% CI: 1.25, 1.44, P ≤0.001) and glucocorticoid use (HR: 1.17, 95% CI: 1.09, 1.25, P ≤0.001).

CONCLUSIONS

Individuals with IA had a lower incidence of COVID-19, probably due to shielding. IA was not associated with increased mortality following COVID-19 infection; being vulnerable (shielded), comorbidities and other factors were associated with increased risk. These key risk factors can identify individuals with IA at greater risk from COVID-19 and advised to shield during high community prevalence.

Prediction of the effects of the top 10 nonsynonymous variants from 30229 SARS-CoV-2 strains on their proteins

Background: SARS-CoV-2 virus is a highly transmissible pathogen that causes COVID-19. The outbreak originated in Wuhan, China in December 2019. A number of nonsynonymous mutations located at different SARS-CoV-2 proteins have been reported by multiple studies. However, there are limited computational studies on the biological impacts of these mutations on the structure and function of the proteins.  

Methods: In our study nonsynonymous mutations of the SARS-CoV-2 genome and their frequencies were identified from 30,229 sequences. Subsequently, the effects of the top 10 highest frequency nonsynonymous mutations of different SARS-CoV-2 proteins were analyzed using bioinformatics tools including co-mutation analysis, prediction of the protein structure stability and flexibility analysis, and prediction of the protein functions.  

Results: A total of 231 nonsynonymous mutations were identified from 30,229 SARS-CoV-2 genome sequences. The top 10 nonsynonymous mutations affecting nine amino acid residues were ORF1a nsp5 P108S, ORF1b nsp12 P323L and A423V, S protein N501Y and D614G, ORF3a Q57H, N protein P151L, R203K and G204R. Many nonsynonymous mutations showed a high concurrence ratio, suggesting these mutations may evolve together and interact functionally. Our result showed that ORF1a nsp5 P108S, ORF3a Q57H and N protein P151L mutations may be deleterious to the function of SARS-CoV-2 proteins. In addition, ORF1a nsp5 P108S and S protein D614G may destabilize the protein structures while S protein D614G may have a more open conformation compared to the wild type.  

Conclusion: The biological consequences of these nonsynonymous mutations of SARS-CoV-2 proteins should be further validated by in vivo and in vitro experimental studies in the future.

Whole genome sequencing of SARS-CoV2 strains circulating in Iran during five waves of pandemic

Purpose

Whole genome sequencing of SARS-CoV2 is important to find useful information about the viral lineages, variants of interests and variants of concern. As there are not enough data about the circulating SARS-CoV2 variants in Iran, we sequenced 54 SARS-CoV2 genomes during the 5 waves of pandemic in Iran.

Methods

After viral RNA extraction from clinical samples collected during the COVID-19 pandemic, next generation sequencing was performed using the Nextseq platform. The sequencing data were analyzed and compared with reference sequences.

Results

During the 1st wave, V and L clades were detected. The second wave was recognized by G, GH and GR clades. Circulating clades during the 3rd wave were GH and GR. In the fourth wave GRY (alpha variant), GK (delta variant) and one GH clade (beta variant) were detected. All viruses in the fifth wave were in clade GK (delta variant). There were different mutations in all parts of the genomes but Spike-D614G, NSP12-P323L, N-R203K and N-G204R were the most frequent mutants in these studied viruses.

Conclusions

These findings display the significance of SARS-CoV2 monitoring to help on time detection of possible variants for pandemic control and vaccination plans.

Microbial ecology and evolution is key to pandemics: using the coronavirus model to mitigate future public health challenges

Pandemics are global challenges that lead to total disruption of human activities. From the inception of human existence, all pandemics have resulted in loss of human lives. The coronavirus disease caused by SAR-CoV-2 began in China and is now at the global scale with an increase in mortality and morbidity. Numerous anthropogenic activities have been implicated in the emergence and severity of pandemics, including COVID-19. These activities cause changes in microbial ecology, leading to evolution due to mutation and recombination. This review hypothesized that an understanding of these anthropogenic activities would explain the dynamics of pandemics. The recent coronavirus model was used to study issues leading to microbial evolution, towards preventing future pandemics. Our review highlighted anthropogenic activities, including deforestation, mining activities, waste treatment, burning of fossil fuel, as well as international travels as drivers of microbial evolution leading to pandemics. Furthermore, human-animal interaction has also been implicated in pandemic incidents. Our study recommends substantial control of such anthropogenic activities as having been highlighted as ways to reduce the frequency of mutation, reduce pathogenic reservoirs, and the emergence of infectious diseases.

Whole-genome analysis of SARS-CoV-2 in a 2020 infection cluster in a nursing home of Southern Italy

Background

Nursing homes have represented important hotspots of viral spread during the initial wave of COVID-19 pandemics. The proximity of patients inside nursing homes allows investigate the dynamics of viral transmission, which may help understand SARS-Cov2 biology and spread.

Methods

SARS-CoV-2 viral genomes obtained from 46 patients infected in an outbreak inside a nursing home in Calabria region (South Italy) were analyzed by Next Generation Sequencing. We also investigated the evolution of viral genomes in 8 patients for which multiple swabs were available. Phylogenetic analysis and haplotype reconstruction were carried out with IQ-TREE software and RegressHaplo tool, respectively.

Results

All viral strains isolated from patients infected in the nursing home were classified as B.1 lineage, clade G. Overall, 14 major single nucleotide variations (SNVs) (frequency > 80%) and 12 minor SNVs (frequency comprised between 20% and 80%) were identified with reference to the Wuhan-H-1 sequence (NC_045512.2).

All patients presented the same 6 major SNVs: D614G in the S gene; P4715L, ntC3037T (F924F) and S5398P in Orf1ab gene; ntC26681T (F53F) in the M gene; and ntC241T in the non-coding UTR region. However, haplotype reconstruction identified a founder haplotype (Hap A) in 36 patients carrying only the 6 common SNVs indicated above, and 10 other haplotypes (Hap B—K) derived from Hap A in the remaining 10 patients. Notably, no significant association between a specific viral haplotype and clinical parameters was found.

Conclusion

The predominant viral strain responsible for the infection in a nursing home in Calabria was the B.1 lineage (clade G). Viral genomes were classified into 11 haplotypes (Hap A in 36 patients, Hap B—K in the remaining patients).

Identifying SARS-CoV-2 Lineage Mutation Hallmarks and Correlating Them With Clinical Outcomes in Egypt: A Pilot Study

The SARS-CoV-2 pandemic has led to over 4.9 million deaths as of October 2021. One of the main challenges of creating vaccines, treatment, or diagnostic tools for the virus is its mutations and emerging variants. A couple of variants were declared as more virulent and infectious than others. Some approaches were used as nomenclature for SARS-CoV-2 variants and lineages. One of the most used is the Pangolin nomenclature. In our study, we enrolled 35 confirmed SARS-CoV-2 patients and sequenced the viral RNA in their samples. We also aimed to highlight the hallmark mutations in the most frequent lineage. We identified a seven-mutation signature for the SARS-CoV-2 C36 lineage, detected in 56 countries and an emerging lineage in Egypt. In addition, we identified one mutation which was highly negatively correlated with the lineage. On the other hand, we found no significant correlation between our clinical outcomes and the C36 lineage. In conclusion, the C36 lineage is an emerging SARS-CoV-2 variant that needs more investigation regarding its clinical outcomes compared to other strains. Our study paves the way for easier diagnosis of variants of concern using mutation signatures.

Genomic Surveillance of SARS-CoV-2 Lineages Indicates Early Circulation of P.1 (Gamma) Variant of Concern in Southern Brazil

The SARS-CoV-2 P.1 lineage emerged in Amazonas (AM), North Brazil and its evolution has been dynamically reported associated with increased transmissibility and/or immune evasion. Here, we evaluated the lineages circulating in 29 cities in Rio Grande do Sul (RS), Southern Brazil between March 2020 and May 2021 and investigated the genetic events associated with the emergence of the P.1. A total of 202 oro/nasopharyngeal SARS-CoV-2 specimens from patients during routine hospital care were submitted to whole-genome sequencing. Phylogenetic and Bayesian Evolutionary Analyses of the P.1 lineage were carried out to determine the relationship between sequences from RS and AM and dated their common ancestor and origin. One hundred six (53%) sequences were assigned as P.1 and most carried the 22 lineage-defining mutations. All the P.1 sequences included other important mutations, such as P314L and R203K/G204R, and revealed a high genetic diversity in the phylogenetic tree. The time-scaled inference suggests that the oldest P.1 sequences from different Brazilian states share a ancestor with those from AM, but the origin of some sequences from RS is unknown. Further, the common ancestor of sequences from RS is dated to mid-June/July 2020, earlier than those previously reported from AM. Our results demonstrate that there is a high degree of genetic diversity among P.1 sequences, which suggests a continuous evolution and community spread of the virus. Although the first P.1 outbreak was reported in AM, the lineage was associated with multiple introductory events and had already been circulating in Southern Brazil prior to November 2020.

IMPORTANCE The SARS-CoV-2 P.1 lineage is associated with increased transmissibility and/or immune evasion and presents a dynamic evolution in Brazil. The significance of our research relies in the fact that we evaluated the SARS-CoV-2 lineages circulating in Southern Brazil between March 2020 and May 2021. This evaluation allowed us to detect the genetic events associated with the emergence of the P.1 and its sublineages. This study is important because we were able to establish that the common ancestor of P.1 sequences from Rio Grande do Sul, Southern Brazil, is dated of mid-June/July 2020, earlier than the P.1 sequences previously reported from Amazonas (AM) state. Noteworthy, the high degree of genetic diversity among P.1 sequences found in this study suggests a continuous evolution and community spread of the virus. Moreover, the oldest P.1 sequences from different Brazilian states share a ancestor with those from AM.

Highlight of potential impact of new viral genotypes of SARS-CoV-2 on vaccines and anti-viral therapeutics

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causal agent of the coronavirus disease (COVID-19) pandemic, has infected millions of people globally. Genetic variation and selective pressures lead to the accumulation of single nucleotide polymorphism (SNP) within the viral genome that may affect virulence, transmission rate, viral recognition and the efficacy of prophylactic and interventional measures. To address these concerns at the genomic level, we assessed the phylogeny and SNPs of the SARS-CoV-2 mutant population collected to date in Iran in relation to globally reported variants. Phylogenetic analysis of mutant strains revealed the occurrence of the variants known as B.1.1.7 (Alpha), B.1.525 (Eta), and B.1.617 (Delta) that appear to have delineated independently in Iran. SNP analysis of the Iranian sequences revealed that the mutations were predominantly positioned within the S protein-coding region, with most SNPs localizing to the S1 subunit. Seventeen S1-localizing SNPs occurred in the RNA binding domain that interacts with ACE2 of the host cell. Importantly, many of these SNPs are predicted to influence the binding of antibodies and anti-viral therapeutics, indicating that the adaptive host response appears to be imposing a selective pressure that is driving the evolution of the virus in this closed population through enhancing virulence. The SNPs detected within these mutant cohorts are addressed with respect to current prophylactic measures and therapeutic interventions.

Insights Into the Changing Landscape of Coronavirus Disease 2019

Coronavirus disease 2019 (COVID-19) is a highly contagious, infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which emerged in late 2019 in Wuhan China. A year after the World Health Organization declared COVID-19 a global pandemic, over 215 million confirmed cases and approximately 5 million deaths have been reported worldwide. In this multidisciplinary review, we summarize important insights for COVID-19, ranging from its origin, pathology, epidemiology, to clinical manifestations and treatment. More importantly, we also highlight the foundational connection between genetics and the development of personalized medicine and how these aspects have an impact on disease treatment and management in the dynamic landscape of this pandemic.

Modeling SARS-CoV-2 spike/ACE2 protein-protein interactions for predicting the binding affinity of new spike variants for ACE2, and novel ACE2 structurally related human protein targets, for COVID-19 handling in the 3PM context

Aims

The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM.

Methods

In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein–protein interface.

Results

Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike “boat-shaped” receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the “stern” or “bow” regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the “hull” region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein–protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities.

Conclusion

In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-021-00267-w.

Detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its first variants in fourplex real-time quantitative reverse transcription-PCR assays

The early diagnosis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is required to identify and isolate contagious patients to prevent further transmission of SARS-CoV-2. In this study, we present a multitarget real-time TaqMan reverse transcription PCR (rRT-PCR) assay for the quantitative detection of SARS-CoV-2 and some of its circulating variants harboring mutations that give the virus a selective advantage. Seven different primer-probe sets that included probes containing locked nucleic acid (LNA) nucleotides were designed to amplify specific wild-type and mutant sequences in Orf1ab, Envelope (E), Spike (S), and Nucleocapsid (N) genes. Furthermore, a newly developed primer-probe set targeted human β₂-microglobulin (B2M) as a highly sensitive internal control for RT efficacy. All singleplex and fourplex assays detected ≤ 14 copies/reaction of quantified synthetic RNA transcripts, with a linear amplification range of nine logarithmic orders. Primer-probe sets for detection of SARS-CoV-2 exhibited no false-positive amplifications with other common respiratory pathogens, including human coronaviruses NL63, 229E, OC43, and HKU-1. Fourplex assays were evaluated using 160 clinical samples positive for SARS-CoV-2. Results showed that SARS-CoV-2 viral RNA was detected in all samples, including viral strains harboring mutations in the Spike coding sequence that became dominant in the pandemic. Given the emergence of SARS-CoV-2 variants and their rapid spread in some populations, fourplex rRT-PCR assay containing four primer-probe sets represents a reliable approach to allow quicker detection of circulating relevant variants in a single reaction.

Hotspot Mutations in SARS-CoV-2

Since its emergence in Wuhan, China, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has spread very rapidly around the world, resulting in a global pandemic. Though the vaccination process has started, the number of COVID-affected patients is still quite large. Hence, an analysis of hotspot mutations of the different evolving virus strains needs to be carried out. In this regard, multiple sequence alignment of 71,038 SARS-CoV-2 genomes of 98 countries over the period from January 2020 to June 2021 is performed using MAFFT followed by phylogenetic analysis in order to visualize the virus evolution. These steps resulted in the identification of hotspot mutations as deletions and substitutions in the coding regions based on entropy greater than or equal to 0.3, leading to a total of 45 unique hotspot mutations. Moreover, 10,286 Indian sequences are considered from 71,038 global SARS-CoV-2 sequences as a demonstrative example that gives 52 unique hotspot mutations. Furthermore, the evolution of the hotspot mutations along with the mutations in variants of concern is visualized, and their characteristics are discussed as well. Also, for all the non-synonymous substitutions (missense mutations), the functional consequences of amino acid changes in the respective protein structures are calculated using PolyPhen-2 and I-Mutant 2.0. In addition to this, SSIPe is used to report the binding affinity between the receptor-binding domain of Spike protein and human ACE2 protein by considering L452R, T478K, E484Q, and N501Y hotspot mutations in that region.

A Global Mutational Profile of SARS-CoV-2: A Systematic Review and Meta-Analysis of 368,316 COVID-19 Patients

Since its first detection in December 2019, more than 232 million cases of COVID-19, including 4.7 million deaths, have been reported by the WHO. The SARS-CoV-2 viral genomes have evolved rapidly worldwide, causing the emergence of new variants. This systematic review and meta-analysis was conducted to provide a global mutational profile of SARS-CoV-2 from December 2019 to October 2020. The review was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA), and a study protocol was lodged with PROSPERO. Data from 62 eligible studies involving 368,316 SARS-CoV-2 genomes were analyzed. The mutational data analyzed showed most studies detected mutations in the Spike protein (n = 50), Nucleocapsid phosphoprotein (n = 34), ORF1ab gene (n = 29), 5′-UTR (n = 28) and ORF3a (n = 25). Under the random-effects model, pooled prevalence of SARS-CoV-2 variants was estimated at 95.1% (95% CI; 93.3–96.4%; I² = 98.952%; p = 0.000) while subgroup meta-analysis by country showed majority of the studies were conducted ‘Worldwide’ (n = 10), followed by ‘Multiple countries’ (n = 6) and the USA (n = 5). The estimated prevalence indicated a need to continuously monitor the prevalence of new mutations due to their potential influence on disease severity, transmissibility and vaccine effectiveness.

Characterisation of SARS-CoV-2 clades based on signature SNPs unveils continuous evolution

Since the emergence of SARS-CoV-2 in Wuhan, China more than a year ago, it has spread across the world in a very short span of time. Although, different forms of vaccines are being rolled out for vaccination programs around the globe, the mutation of the virus is still a cause of concern among the research communities. Hence, it is important to study the constantly evolving virus and its strains in order to provide a much more stable form of cure. This fact motivated us to conduct this research where we have initially carried out multiple sequence alignment of 15359 and 3033 global dataset without Indian and the dataset of exclusive Indian SARS-CoV-2 genomes respectively, using MAFFT. Subsequently, phylogenetic analyses are performed using Nextstrain to identify virus clades. Consequently, the virus strains are found to be distributed among 5 major clades or clusters viz. 19A, 19B, 20A, 20B and 20C. Thereafter, mutation points as SNPs are identified in each clade. Henceforth, from each clade top 10 signature SNPs are identified based on their frequency i.e. number of occurrences in the virus genome. As a result, 50 such signature SNPs are individually identified for global dataset without Indian and dataset of exclusive Indian SARS-CoV-2 genomes respectively. Out of each 50 signature SNPs, 39 and 41 unique SNPs are identified among which 25 non-synonymous signature SNPs (out of 39) resulted in 30 amino acid changes in protein while 27 changes in amino acid are identified from 22 non-synonymous signature SNPs (out of 41). These 30 and 27 amino acid changes for the non-synonymous signature SNPs are visualised in their respective protein structure as well. Finally, in order to judge the characteristics of the identified clades, the non-synonymous signature SNPs are considered to evaluate the changes in proteins as biological functions with the sequences using PROVEAN and PolyPhen-2 while I-Mutant 2.0 is used to evaluate their structural stability. As a consequence, for global dataset without Indian sequences, G251V in ORF3a in clade 19A, F308Y and G196V in NSP4 and ORF3a in 19B are the unique amino acid changes which are responsible for defining each clade as they are all deleterious and unstable. Such changes which are common for both global dataset without Indian and dataset of exclusive Indian sequences are R203M in Nucleocapsid for 20B, T85I and Q57H in NSP2 and ORF3a respectively for 20C while for exclusive Indian sequences such unique changes are A97V in RdRp, G339S and G339C in NSP2 in 19A and Q57H in ORF3a in 20A.

Mutational Asymmetries in the SARS-CoV-2 Genome May Lead to Increased Hydrophobicity of Virus Proteins

The genomic diversity of SARS-CoV-2 has been a focus during the ongoing COVID-19 pandemic. Here, we analyzed the distribution and character of emerging mutations in a data set comprising more than 95,000 virus genomes covering eight major SARS-CoV-2 lineages in the GISAID database, including genotypes arising during COVID-19 therapy. Globally, the C>U transitions and G>U transversions were the most represented mutations, accounting for the majority of single-nucleotide variations. Mutational spectra were not influenced by the time the virus had been circulating in its host or medical treatment. At the amino acid level, we observed about a 2-fold excess of substitutions in favor of hydrophobic amino acids over the reverse. However, most mutations constituting variants of interests of the S-protein (spike) lead to hydrophilic amino acids, counteracting the global trend. The C>U and G>U substitutions altered codons towards increased amino acid hydrophobicity values in more than 80% of cases. The bias is explained by the existing differences in the codon composition for amino acids bearing contrasting biochemical properties. Mutation asymmetries apparently influence the biochemical features of SARS CoV-2 proteins, which may impact protein-protein interactions, fusion of viral and cellular membranes, and virion assembly.

Analysis of COVID-19 pandemics in Kazakhstan

Background: This study aimed to analyze the demographic and epidemiological features of identified COVID-19 cases in Kazakhstan.

Study design: A cross-sectional study.

Methods: This cross-sectional study aimed to analyze COVID-19 cases (n=5116) collected from March 13 to June 6, 2020, in Kazakhstan. The data were obtained from a state official medical electronic database. The study investigated the geographic and demographic data of patients as well as the association of COVID-19 cases with gender and age. The prevalence of symptoms, the presence of comorbidities, complications, and COVID-19 mortality were determined for all patients.

Results: The mean ±SD age of the patients in this study was 34.8 ±17.6 years, and the majority (55.7%) of COVID-19 cases were male and residents of cities (79.6%). In total, 80% of the cases had the asymptomatic/mild form of the disease. Cough (20.8 %) and sore throat (17.1%) were the most common symptoms among patients, and pneumonia was diagnosed in 1 out of 5 cases. Acute respiratory distress syndrome (ARDS) was recorded in 1.2% of the patients. The fatality rate was 1% in the study population and lethality was 2.6 times higher in males compared to females. Each additional year in age increased the probability of COVID-19 infection by 1.06 times. The presence of cardiovascular, diabetes, respiratory, and kidney diseases affected the rate of mortality (P<0.05).

Conclusions: The results demonstrated a high proportion (40%) of the asymptomatic type of coronavirus infection in the Kazakhstan population. The severity of COVID-19 symptoms and lethality were directly related to the age of patients and the presence of comorbidities.

SARS-CoV-2 main protease suppresses type I interferon production by preventing nuclear translocation of phosphorylated IRF3

Suppression of type I interferon (IFN) response is one pathological outcome of the infection of highly pathogenic human coronaviruses. To effect this, severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 encode multiple IFN antagonists. In this study, we reported on the IFN antagonism of SARS-CoV-2 main protease NSP5. NSP5 proteins of both SARS-CoV and SARS-CoV-2 counteracted Sendai virus-induced IFN production. NSP5 variants G15S and K90R commonly seen in circulating strains of SARS-CoV-2 retained the IFN-antagonizing property. The suppressive effect of NSP5 on IFN-β gene transcription induced by RIG-I, MAVS, TBK1 and IKKϵ suggested that NSP5 likely acts at a step downstream of IRF3 phosphorylation in the cytoplasm. NSP5 did not influence steady-state expression or phosphorylation of IRF3, suggesting that IRF3, regardless of its phosphorylation state, might not be the substrate of NSP5 protease. However, nuclear translocation of phosphorylated IRF3 was severely compromised in NSP5-expressing cells. Taken together, our work revealed a new mechanism by which NSP5 proteins encoded by SARS-CoV and SARS-CoV-2 antagonize IFN production by retaining phosphorylated IRF3 in the cytoplasm. Our findings have implications in rational design and development of antiviral agents against SARS-CoV-2.

SARS-CoV-2 genomic diversity and the implications for qRT-PCR diagnostics and transmission

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 lineages that have spread throughout the world. In this study, we investigated 129 RNA-seq data sets and 6928 consensus genomes to contrast the intra-host and inter-host diversity of SARS-CoV-2. Our analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights intra-host single nucleotide variant (iSNV) and SNP similarity, albeit with differences in C > U changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of insertions and deletions contribute to the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.

Intra-host variation and evolutionary dynamics of SARS-CoV-2 populations in COVID-19 patients

BACKGROUND

Since early February 2021, the causative agent of COVID-19, SARS-CoV-2, has infected over 104 million people with more than 2 million deaths according to official reports. The key to understanding the biology and virus-host interactions of SARS-CoV-2 requires the knowledge of mutation and evolution of this virus at both inter- and intra-host levels. However, despite quite a few polymorphic sites identified among SARS-CoV-2 populations, intra-host variant spectra and their evolutionary dynamics remain mostly unknown.

METHODS

Using high-throughput sequencing of metatranscriptomic and hybrid captured libraries, we characterized consensus genomes and intra-host single nucleotide variations (iSNVs) of serial samples collected from eight patients with COVID-19. The distribution of iSNVs along the SARS-CoV-2 genome was analyzed and co-occurring iSNVs among COVID-19 patients were identified. We also compared the evolutionary dynamics of SARS-CoV-2 population in the respiratory tract (RT) and gastrointestinal tract (GIT).

RESULTS

The 32 consensus genomes revealed the co-existence of different genotypes within the same patient. We further identified 40 intra-host single nucleotide variants (iSNVs). Most (30/40) iSNVs presented in a single patient, while ten iSNVs were found in at least two patients or identical to consensus variants. Comparing allele frequencies of the iSNVs revealed a clear genetic differentiation between intra-host populations from the respiratory tract (RT) and gastrointestinal tract (GIT), mostly driven by bottleneck events during intra-host migrations. Compared to RT populations, the GIT populations showed a better maintenance and rapid development of viral genetic diversity following the suspected intra-host bottlenecks.

CONCLUSIONS

Our findings here illustrate the intra-host bottlenecks and evolutionary dynamics of SARS-CoV-2 in different anatomic sites and may provide new insights to understand the virus-host interactions of coronaviruses and other RNA viruses.