Coronavirus genomes carry the signatures of their habitats

SARS-CoV-2 and innate immunity: the good, the bad, and the "goldilocks"

An ancient conflict between hosts and pathogens has driven the innate and adaptive arms of immunity. Knowledge about this interplay can not only help us identify biological mechanisms but also reveal pathogen vulnerabilities that can be leveraged therapeutically. The humoral response to SARS-CoV-2 infection has been the focus of intense research, and the role of the innate immune system has received significantly less attention. Here, we review current knowledge of the innate immune response to SARS-CoV-2 infection and the various means SARS-CoV-2 employs to evade innate defense systems. We also consider the role of innate immunity in SARS-CoV-2 vaccines and in the phenomenon of long COVID.

Genetic diversity and evolutionary patterns of SARS-CoV-2 among the Bhutanese population during the pandemic

BACKGROUND

The coronavirus disease 2019 (COVID-19) pandemic, caused by a dynamic virus, has had a profound global impact. Despite declining global COVID-19 cases and mortality rates, the emergence of new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants remains a major concern. This study provides a comprehensive analysis of the genomic sequences of SARS-CoV-2 within the Bhutanese population during the pandemic. The primary aim was to elucidate the molecular epidemiology and evolutionary patterns of SARS-CoV-2 in Bhutan, with a particular focus on genetic variations and lineage dynamics.

METHODS

Whole-genome sequences of SARS-CoV-2 collected from Bhutan between May 2020 and February 2023 (n=135) were retrieved from the Global Initiative on Sharing All Influenza Database.

RESULTS

The SARS-CoV-2 variants in Bhutan were predominantly classified within the Nextstrain clade 20A (31.1%), followed by clade 21L (20%) and clade 22D (15.6%). We identified 26 Pangolin lineages with variations in their spatial and temporal distribution. Bayesian time-scaled phylogenetic analysis estimated the time to the most recent common ancestor as February 15, 2020, with a substitution rate of 0.97×10-3 substitutions per site per year. Notably, the spike glycoprotein displayed the highest mutation frequency among major viral proteins, with 116 distinct mutations, including D614G. The Bhutanese isolates also featured mutations such as E484K, K417N, and S477N in the spike protein, which have implications for altered viral properties.

CONCLUSION

This is the first study to describe the genetic diversity of SARS-CoV-2 circulating in Bhutan during the pandemic, and this data can inform public health policies and strategies for preventing future outbreaks in Bhutan.

Similar deamination activities but different phenotypic outcomes induced by APOBEC3 enzymes in breast epithelial cells

APOBEC3 (A3) enzymes deaminate cytosine to uracil in viral single-stranded DNA as a mutagenic barrier for some viruses. A3-induced deaminations can also occur in human genomes resulting in an endogenous source of somatic mutations in multiple cancers. However, the roles of each A3 are unclear since few studies have assessed these enzymes in parallel. Thus, we developed stable cell lines expressing A3A, A3B, or A3H Hap I using non-tumorigenic MCF10A and tumorigenic MCF7 breast epithelial cells to assess their mutagenic potential and cancer phenotypes in breast cells. The activity of these enzymes was characterized by γH2AX foci formation and in vitro deamination. Cell migration and soft agar colony formation assays assessed cellular transformation potential. We found that all three A3 enzymes had similar γH2AX foci formation, despite different deamination activities in vitro. Notably, in nuclear lysates, the in vitro deaminase activity of A3A, A3B, and A3H did not require digestion of cellular RNA, in contrast to that of A3B and A3H in whole-cell lysates. Their similar activities in cells, nonetheless, resulted in distinct phenotypes where A3A decreased colony formation in soft agar, A3B decreased colony formation in soft agar after hydroxyurea treatment, and A3H Hap I promoted cell migration. Overall, we show that in vitro deamination data do not always reflect cell DNA damage, all three A3s induce DNA damage, and the impact of each is different.

Comparative susceptibility of SARS-CoV-2, SARS-CoV, and MERS-CoV across mammals

Exploring wild reservoirs of pathogenic viruses is critical for their long-term control and for predicting future pandemic scenarios. Here, a comparative in vitro infection analysis was first performed on 83 cell cultures derived from 55 mammalian species using pseudotyped viruses bearing S proteins from SARS-CoV-2, SARS-CoV, and MERS-CoV. Cell cultures from Thomas's horseshoe bats, king horseshoe bats, green monkeys, and ferrets were found to be highly susceptible to SARS-CoV-2, SARS-CoV, and MERS-CoV pseudotyped viruses. Moreover, five variants (del69-70, D80Y, S98F, T572I, and Q675H), that beside spike receptor-binding domain can significantly alter the host tropism of SARS-CoV-2. An examination of phylogenetic signals of transduction rates revealed that closely related taxa generally have similar susceptibility to MERS-CoV but not to SARS-CoV and SARS-CoV-2 pseudotyped viruses. Additionally, we discovered that the expression of 95 genes, e.g., PZDK1 and APOBEC3, were commonly associated with the transduction rates of SARS-CoV, MERS-CoV, and SARS-CoV-2 pseudotyped viruses. This study provides basic documentation of the susceptibility, variants, and molecules that underlie the cross-species transmission of these coronaviruses.

The mutational spectrum of SARS-CoV-2 genomic and antigenomic RNA

The raw material for viral evolution is provided by intra-host mutations occurring during replication, transcription or post-transcription. Replication and transcription of Coronaviridae proceed through the synthesis of negative-sense 'antigenomes' acting as templates for positive-sense genomic and subgenomic RNA. Hence, mutations in the genomes of SARS-CoV-2 and other coronaviruses can occur during (and after) the synthesis of either negative-sense or positive-sense RNA, with potentially distinct patterns and consequences. We explored for the first time the mutational spectrum of SARS-CoV-2 (sub)genomic and anti(sub)genomic RNA. We use a high-quality deep sequencing dataset produced using a quantitative strand-aware sequencing method, controlled for artefacts and sequencing errors, and scrutinized for accurate detection of within-host diversity. The nucleotide differences between negative- and positive-sense strand consensus vary between patients and do not show dependence on age or sex. Similarities and differences in mutational patterns between within-host minor variants on the two RNA strands suggested strand-specific mutations or editing by host deaminases and oxidative damage. We observe generally neutral and slight negative selection on the negative strand, contrasting with purifying selection in ORF1a, ORF1b and S genes of the positive strand of the genome.

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

The low abundance of CpG in the SARS-CoV-2 genome is not an evolutionarily signature of ZAP

The zinc finger antiviral protein (ZAP) is known to restrict viral replication by binding to the CpG rich regions of viral RNA, and subsequently inducing viral RNA degradation. This enzyme has recently been shown to be capable of restricting SARS-CoV-2. These data have led to the hypothesis that the low abundance of CpG in the SARS-CoV-2 genome is due to an evolutionary pressure exerted by the host ZAP. To investigate this hypothesis, we performed a detailed analysis of many coronavirus sequences and ZAP RNA binding preference data. Our analyses showed neither evidence for an evolutionary pressure acting specifically on CpG dinucleotides, nor a link between the activity of ZAP and the low CpG abundance of the SARS-CoV-2 genome.

The substitution spectra of coronavirus genomes

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has triggered an unprecedented international effort to sequence complete viral genomes. We leveraged this wealth of information to characterize the substitution spectrum of SARS-CoV-2 and to compare it with those of other human and animal coronaviruses. We show that, once nucleotide composition is taken into account, human and most animal coronaviruses display a mutation spectrum dominated by C to U and G to U substitutions, a feature that is not shared by other positive-sense RNA viruses. However, the proportions of C to U and G to U substitutions tend to decrease as divergence increases, suggesting that, whatever their origin, a proportion of these changes is subsequently eliminated by purifying selection. Analysis of the sequence context of C to U substitutions showed little evidence of apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC)-mediated editing and such contexts were similar for SARS-CoV-2 and Middle East respiratory syndrome coronavirus sampled from different hosts, despite different repertoires of APOBEC3 proteins in distinct species. Conversely, we found evidence that C to U and G to U changes affect CpG dinucleotides at a frequency higher than expected. Whereas this suggests ongoing selective reduction of CpGs, this effect alone cannot account for the substitution spectra. Finally, we show that, during the first months of SARS-CoV-2 pandemic spread, the frequency of both G to U and C to U substitutions increased. Our data suggest that the substitution spectrum of SARS-CoV-2 is determined by an interplay of factors, including intrinsic biases of the replication process, avoidance of CpG dinucleotides and other constraints exerted by the new host.

Striking lineage diversity of severe acute respiratory syndrome coronavirus 2 from non-human sources

Due to the necessity to control human-to-human spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the overwhelming majority of the generated data on this virus was solely related to the genomic characteristics of strains infecting humans; conversely, this work aimed to recover and analyze the diversity of viral genomes from non-human sources. From a set of 3595 publicly available SARS-CoV-2 genome sequences, 128 lineages were identified in non-human hosts, the majority represented by the variants of concern Delta (n = 1105, 30.7%) and Alpha (n = 466, 12.9%), followed by B.1.1.298 lineage (n = 458, 12.7%). Environment, Neovison vison, Odocoileus virginianus and Felis catus were the non-human sources with the highest number of lineages (14, 12 and 10, respectively). Phylogenomic analyses showed viral clusters from environmental sources, N. vison, O. virginianus, Panthera tigris, and Panthera leo. These clusters were collectively related to human viruses as well as all other non-human sources that were heterogeneously distributed in the phylogenetic tree. Further, the genetic details of viral genomes from bats and pangolins were independently investigated owing to their high divergence, revealing five distinct clusters. Cluster 4 exclusively included bat-sourced genomes and the SARS-CoV-2 reference strain Wuhan-01. In summary, this study identified new genetic landmarks of SARS-CoV-2 evolution. We propose potential interspecies transmission routes of SARS-CoV-2 between animals and humans, which should be considered in order to establish better pathogen surveillance and containment strategies.

The Roles of APOBEC-mediated RNA Editing in SARS-CoV-2 Mutations, Replication and Fitness.. [PMID: 34981048 PMCID: PMC8722585 DOI: 10.1101/2021.12.18.473309]

During COVID-19 pandemic, mutations of SARS-CoV-2 produce new strains that can be more infectious or evade vaccines. Viral RNA mutations can arise from misincorporation by RNA-polymerases and modification by host factors. Analysis of SARS-CoV-2 sequence from patients showed a strong bias toward C-to-U mutation, suggesting a potential mutational role by host APOBEC cytosine deaminases that possess broad anti-viral activity. We report the first experimental evidence demonstrating that APOBEC3A, APOBEC1, and APOBEC3G can edit on specific sites of SARS-CoV-2 RNA to produce C-to-U mutations. However, SARS-CoV-2 replication and viral progeny production in Caco-2 cells are not inhibited by the expression of these APOBECs. Instead, expression of wild-type APOBEC3 greatly promotes viral replication/propagation, suggesting that SARS-CoV-2 utilizes the APOBEC-mediated mutations for fitness and evolution. Unlike the random mutations, this study suggests the predictability of all possible viral genome mutations by these APOBECs based on the UC/AC motifs and the viral genomic RNA structure.

Efficient Editing of SARS-CoV-2 genomic RNA by Host APOBEC deaminases and Its Potential Impacts on the Viral Replication and Emergence of New Strains in COVID-19 Pandemic

Intra-host evolution during SARS-CoV-2 prolonged infection

Long-term infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represents a challenge to virus dispersion and the control of coronavirus disease 2019 (COVID-19) pandemic. The reason why some people have prolonged infection and how the virus persists for so long are still not fully understood. Recent studies suggested that the accumulation of intra-host single nucleotide variants (iSNVs) over the course of the infection might play an important role in persistence as well as emergence of mutations of concern. For this reason, we aimed to investigate the intra-host evolution of SARS-CoV-2 during prolonged infection. Thirty-three patients who remained reverse transcription polymerase chain reaction (RT-PCR) positive in the nasopharynx for on average 18 days from the symptoms onset were included in this study. Whole-genome sequences were obtained for each patient at two different time points. Phylogenetic, populational, and computational analyses of viral sequences were consistent with prolonged infection without evidence of coinfection in our cohort. We observed an elevated within-host genomic diversity at the second time point samples positively correlated with cycle threshold (Ct) values (lower viral load). Direct transmission was also confirmed in a small cluster of healthcare professionals that shared the same workplace by the presence of common iSNVs. A differential accumulation of missense variants between the time points was detected targeting crucial structural and non-structural proteins such as Spike and helicase. Interestingly, longitudinal acquisition of iSNVs in Spike protein coincided in many cases with SARS-CoV-2 reactive and predicted T cell epitopes. We observed a distinguishing pattern of mutations over the course of the infection mainly driven by increasing A→U and decreasing G→A signatures. G→A mutations may be associated with RNA-editing enzyme activities; therefore, the mutational profiles observed in our analysis were suggestive of innate immune mechanisms of the host cell defense. Therefore, we unveiled a dynamic and complex landscape of host and pathogen interaction during prolonged infection of SARS-CoV-2, suggesting that the host’s innate immunity shapes the increase of intra-host diversity. Our findings may also shed light on possible mechanisms underlying the emergence and spread of new variants resistant to the host immune response as recently observed in COVID-19 pandemic.

Molecular characterization of SARS-CoV-2 from Bangladesh: implications in genetic diversity, possible origin of the virus, and functional significance of the mutations

In a try to understand the pathogenesis, evolution and epidemiology of the SARS-CoV-2 virus, scientists from all over the world are tracking its genomic changes in real-time. Genomic studies can be helpful in understanding the disease dynamics. We have downloaded 324 complete and near complete SARS-CoV-2 genomes submitted in GISAID database from Bangladesh which were isolated between 30 March to 7 September, 2020. We then compared these genomes with Wuhan reference sequence and found 4160 mutation events including 2253 missense single nucleotide variations, 38 deletions and 10 insertions. The C>T nucleotide change was most prevalent (41% of all mutations) possibly due to selective mutation pressure to reduce CpG sites to evade CpG targeted host immune response. The most frequent mutation that occurred in 98% isolates was 3037C>T which is a synonymous change that usually accompanied 3 other mutations that include 241C>T, 14408C>T (P323L in RdRp) and 23403A>G (D614G in spike protein). The P323L was reported to increase mutation rate and D614G is associated with increased viral replication and currently most prevalent variant circulating all over the world. We identified multiple missense mutations in B-cell and T-cell predicted epitope regions and/or PCR target regions (including R203K and G204R that occurred in 86% of the isolates) that may impact immunogenicity and/or RT-PCR based diagnosis. Our analysis revealed 5 large deletion events in ORF7a and ORF8 gene products that may be associated with less severity of the disease and increased viral clearance. Our phylogeny analysis identified most of the isolates belonged to the Nextstrain clade 20B (86%) and GISAID clade GR (88%). Most of our isolates shared common ancestors either directly with European countries or jointly with middle eastern countries as well as Australia and India. Interestingly, the 19B clade (GISAID S clade) was unique to Chittagong, which was originally prevalent in China. This reveals possible multiple introductions of the virus in Bangladesh via different routes. Hence, more genome sequencing and analysis with related clinical data is needed to interpret functional significance and better predict the disease dynamics that may be helpful for policy makers to control the COVID-19 pandemic.

HIV-1 and SARS-CoV-2: Patterns in the evolution of two pandemic pathogens

Humanity is currently facing the challenge of two devastating pandemics caused by two very different RNA viruses: HIV-1, which has been with us for decades, and SARS-CoV-2, which has swept the world in the course of a single year. The same evolutionary strategies that drive HIV-1 evolution are at play in SARS-CoV-2. Single nucleotide mutations, multi-base insertions and deletions, recombination, and variation in surface glycans all generate the variability that, guided by natural selection, enables both HIV-1's extraordinary diversity and SARS-CoV-2's slower pace of mutation accumulation. Even though SARS-CoV-2 diversity is more limited, recently emergent SARS-CoV-2 variants carry Spike mutations that have important phenotypic consequences in terms of both antibody resistance and enhanced infectivity. We review and compare how these mutational patterns manifest in these two distinct viruses to provide the variability that fuels their evolution by natural selection.

AID and APOBECs as Multifaceted Intrinsic Virus-Restricting Factors: Emerging Concepts in the Light of COVID-19

The AID (activation-induced cytidine deaminase)/APOBEC (apolipoprotein B mRNA editing enzyme catalytic subunit) family with its multifaceted mode of action emerges as potent intrinsic host antiviral system that acts against a variety of DNA and RNA viruses including coronaviruses. All family members are cytosine-to-uracil deaminases that either have a profound role in driving a strong and specific humoral immune response (AID) or restricting the virus itself by a plethora of mechanisms (APOBECs). In this article, we highlight some of the key aspects apparently linking the AID/APOBECs and SARS-CoV-2. Among those is our discovery that APOBEC4 shows high expression in cell types and anatomical parts targeted by SARS-CoV-2. Additional focus is given by us to the lymphoid structures and AID as the master regulator of germinal center reactions, which result in antibody production by plasma and memory B cells. We propose the dissection of the AID/APOBECs gene signature towards decisive determinants of the patient-specific and/or the patient group-specific antiviral response. Finally, the patient-specific mapping of the AID/APOBEC polymorphisms should be considered in the light of COVID-19.

The Heterogeneous Landscape and Early Evolution of Pathogen-Associated CpG Dinucleotides in SARS-CoV-2

COVID-19 can lead to acute respiratory syndrome, which can be due to dysregulated immune signaling. We analyze the distribution of CpG dinucleotides, a pathogen-associated molecular pattern, in the SARS-CoV-2 genome. We characterize CpG content by a CpG force that accounts for statistical constraints acting on the genome at the nucleotidic and amino acid levels. The CpG force, as the CpG content, is overall low compared with other pathogenic betacoronaviruses; however, it widely fluctuates along the genome, with a particularly low value, comparable with the circulating seasonal HKU1, in the spike coding region and a greater value, comparable with SARS and MERS, in the highly expressed nucleocapside coding region (N ORF), whose transcripts are relatively abundant in the cytoplasm of infected cells and present in the 3'UTRs of all subgenomic RNA. This dual nature of CpG content could confer to SARS-CoV-2 the ability to avoid triggering pattern recognition receptors upon entry, while eliciting a stronger response during replication. We then investigate the evolution of synonymous mutations since the outbreak of the COVID-19 pandemic, finding a signature of CpG loss in regions with a greater CpG force. Sequence motifs preceding the CpG-loss-associated loci in the N ORF match recently identified binding patterns of the zinc finger antiviral protein. Using a model of the viral gene evolution under human host pressure, we find that synonymous mutations seem driven in the SARS-CoV-2 genome, and particularly in the N ORF, by the viral codon bias, the transition-transversion bias, and the pressure to lower CpG content.

Glimpses into evolutionary trajectories of SARS-CoV-2: emerging variants and potential immune evasion routes

Tweetable abstract An opinion on the coronaviruses' evolution paradoxes, the continuing adaptation of the SARS-CoV-2 in humans following the zoonotic transmission, and clues into escape routes from host immune responses.

Genomic characterization and evolution of SARS-CoV-2 of a Canadian population

COVID-19 has greatly affected public health and world economy. In this study, we analyzed 129 full-length genomes of SARS-CoV-2 viruses of a Canadian population during early phase of the pandemic. Phylogenetic analysis revealed three major paths of transmission of SARS-CoV-2 viruses into Canada. Twenty-one substitutions that have frequencies greater than 3% of viral population were identified. Analysis of these substitutions indicated that P1427I (ORF1b), Y1464C (ORF1b), and Q57H (ORF3a) might affect functions of the corresponding SARS-CoV-2 encoded proteins. Additionally, we found the evidence of positive selection on the ORF3a and codon 614 of Spike protein, suggesting the viral components responsible for host entry and activation of inflammation response were targeted by host immune responses. The study showed genomic variation and evolution of SARS-CoV-2 in a Canadian population. These information may help develop preventive strategies and be used for further study of SARS-CoV-2 pathogenesis and therapeutics development.

Host-directed editing of the SARS-CoV-2 genome

The extensive sequence data generated from SARS-CoV-2 during the 2020 pandemic has facilitated the study of viral genome evolution over a brief period of time. This has highlighted instances of directional mutation pressures exerted on the SARS-CoV-2 genome from host antiviral defense systems. In this brief review we describe three such human defense mechanisms, the apolipoprotein B mRNA editing catalytic polypeptide-like proteins (APOBEC), adenosine deaminase acting on RNA proteins (ADAR), and reactive oxygen species (ROS), and discuss their potential implications on SARS-CoV-2 evolution.

Domains and Functions of Spike Protein in Sars-Cov-2 in the Context of Vaccine Design

The spike protein in SARS-CoV-2 (SARS-2-S) interacts with the human ACE2 receptor to gain entry into a cell to initiate infection. Both Pfizer/BioNTech's BNT162b2 and Moderna's mRNA-1273 vaccine candidates are based on stabilized mRNA encoding prefusion SARS-2-S that can be produced after the mRNA is delivered into the human cell and translated. SARS-2-S is cleaved into S1 and S2 subunits, with S1 serving the function of receptor-binding and S2 serving the function of membrane fusion. Here, I dissect in detail the various domains of SARS-2-S and their functions discovered through a variety of different experimental and theoretical approaches to build a foundation for a comprehensive mechanistic understanding of how SARS-2-S works to achieve its function of mediating cell entry and subsequent cell-to-cell transmission. The integration of structure and function of SARS-2-S in this review should enhance our understanding of the dynamic processes involving receptor binding, multiple cleavage events, membrane fusion, viral entry, as well as the emergence of new viral variants. I highlighted the relevance of structural domains and dynamics to vaccine development, and discussed reasons for the spike protein to be frequently featured in the conspiracy theory claiming that SARS-CoV-2 is artificially created.

Cytidine deamination-induced perpetual immunity to SAR-CoV-2 infection is a potential new therapeutic target

As the world is racing to develop perpetual immunity to the SARS-CoV-2 virus. The emergence of new viral strains, together with vaccination and reinfections, are all contributing to a long-term immunity against the deadly virus that has taken over the world since its introduction to humans in late December 2019. The discovery that more than 95 percent of people who recovered from COVID-19 had long-lasting immunity and that asymptomatic people have a different immune response to SARS-CoV-2 than symptomatic people has shifted attention to how our immune system initiates such diverse responses. These findings have provided reason to believe that SARS-CoV-2 days are numbered. Hundreds of research papers have been published on the causes of long-lasting immune responses and variations in the numbers of different immune cell types in COVID 19 survivors, but the main reason of these differences has still not been adequately identified. In this article, we focus on the activation-induced cytidine deaminase (AID), which initiates molecular processes that allow our immune system to generate antibodies against SARS-CoV-2. To establish lasting immunity to SARS-CoV-2, we suggest that AID could be the key to unlocking it.