Rate variation and recurrent sequence errors in pandemic-scale phylogenetics.

Prediction analysis of porcine AXL protein as a potential receptor for SARS-CoV-2

MotivationBuilding reliable phylogenies from very large collections of sequences with a limited number of phylogenetically informative sites is challenging because sequencing errors and recurrent/backward mutations interfere with the phylogenetic signal, confounding true evolutionary relationships. Massive global efforts of sequencing genomes and reconstructing the phylogeny of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains exemplify these difficulties since there are only hundreds of phylogenetically informative sites but millions of genomes. For such datasets, we set out to develop a method for building the phylogenetic tree of genomic haplotypes consisting of positions harboring common variants to improve the signal-to-noise ratio for more accurate and fast phylogenetic inference of resolvable phylogenetic features.ResultsWe present the TopHap approach that determines spatiotemporally common haplotypes of common variants and builds their phylogeny at a fraction of the computational time of traditional methods. We develop a bootstrap strategy that resamples genomes spatiotemporally to assess topological robustness. The application of TopHap to build a phylogeny of 68 057 SARS-CoV-2 genomes (68KG) from the first year of the pandemic produced an evolutionary tree of major SARS-CoV-2 haplotypes. This phylogeny is concordant with the mutation tree inferred using the co-occurrence pattern of mutations and recovers key phylogenetic relationships from more traditional analyses. We also evaluated alternative roots of the SARS-CoV-2 phylogeny and found that the earliest sampled genomes in 2019 likely evolved by four mutations of the most recent common ancestor of all SARS-CoV-2 genomes. An application of TopHap to more than 1 million SARS-CoV-2 genomes reconstructed the most comprehensive evolutionary relationships of major variants, which confirmed the 68KG phylogeny and provided evolutionary origins of major and recent variants of concern.Availability and implementationTopHap is available at https://github.com/SayakaMiura/TopHap.Supplementary informationSupplementary data are available at Bioinformatics online.

Evidence from the evolutionary analysis of nucleotide sequences for a recombinant history of SARS-CoV

BackgroundThe emergence and rapid spread of new severe acute respiratory syndrome coronavirus 2 (SARS-COV-2) variants have challenged the control of the COVID-19 pandemic globally. Burundi was not spared by that pandemic, but the genetic diversity, evolution, and epidemiology of those variants in the country remained poorly understood. The present study sought to investigate the role of different SARS-COV-2 variants in the successive COVID-19 waves experienced in Burundi and the impact of their evolution on the course of that pandemic.We conducted a cross-sectional descriptive study using positive SARS-COV-2 samples for genomic sequencing. Subsequently, we performed statistical and bioinformatics analyses of the genome sequences in light of available metadata.ResultsIn total, we documented 27 PANGO lineages of which BA.1, B.1.617.2, AY.46, AY.122, and BA.1.1, all VOCs, accounted for 83.15% of all the genomes isolated in Burundi from May 2021 to January 2022. Delta (B.1.617.2) and its descendants predominated the peak observed in July–October 2021. It replaced the previously predominant B.1.351 lineage. It was itself subsequently replaced by Omicron (B.1.1.529, BA.1, and BA.1.1). Furthermore, we identified amino acid mutations including E484K, D614G, and L452R known to increase infectivity and immune escape in the spike proteins of Delta and Omicron variants isolated in Burundi. The SARS-COV-2 genomes from imported and community-detected cases were genetically closely related.ConclusionThe global emergence of SARS-COV-2 VOCs and their subsequent introductions in Burundi was accompanied by new peaks (waves) of COVID-19. The relaxation of travel restrictions and the mutations occurring in the virus genome played an important role in the introduction and the spread of new SARS-COV-2 variants in the country. It is of utmost importance to strengthen the genomic surveillance of SARS-COV-2, enhance the protection by increasing the SARS-COV-2 vaccine coverage, and adjust the public health and social measures ahead of the emergence or introduction of new SARS-COV-2 VOCs in the country.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) evolved slowly over the first year of the Coronavirus Disease 19 (COVID-19) pandemic with differential mutation rates across lineages. Here, we explore how this variation arose. Whether evolutionary change accumulated gradually within lineages or during viral lineage branching is unclear. Using phylogenetic regression models, we show that ~13% of SARS-CoV-2 genomic divergence up to May 2020 is attributable to lineage branching events (punctuated evolution). The net number of branching events along lineages predicts ~5% of the deviation from the strict molecular clock. We did not detect punctuated evolution in SARS-CoV-1, possibly due to the small sample size, and in sarbecovirus broadly, likely due to a different evolutionary process altogether. Punctuation in SARS-CoV-2 is probably neutral because most mutations were not positively selected and because the strength of the punctuational effect remained constant over time, at least until May 2020, and across continents. However, the small punctuational contribution to SARS-CoV-2 diversity is consistent with the founder effect arising from narrow transmission bottlenecks. Therefore, punctuation in SARS-CoV-2 may represent the macroevolutionary consequence (rate variation) of a microevolutionary process (transmission bottleneck).

In 2019, the novel highly infectious severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak rapidly led to a global pandemic with more than 346 million confirmed cases worldwide, resulting in 5.5 million associated deaths (January 2022). Entry of all SARS-CoV-2 variants is mediated by the cellular angisin-converting enzyme 2 (ACE2). The virus abundantly replicates in the epithelia of the upper respiratory tract. Beyond vaccines for immunization, there is an imminent need for novel treatment options in COVID-19 patients. So far, only a few drugs have found their way into the clinics, often with modest success. Specific gene silencing based on small interfering RNA (siRNA) has emerged as a promising strategy for therapeutic intervention, preventing/limiting SARS-CoV-2 entry into host cells or interfering with viral replication. Here, we pursued both strategies. We designed and screened nine siRNAs (siA1-9) targeting the viral entry receptor ACE2. SiA1, (siRNA against exon1 of ACE2 mRNA) was most efficient, with up to 90% knockdown of the ACE2 mRNA and protein for at least six days. In vitro, siA1 application was found to protect Vero E6 and Huh-7 cells from infection with SARS-CoV-2 with an up to ∼92% reduction of the viral burden indicating that the treatment targets both the endosomal and the viral entry at the cytoplasmic membrane. Since the RNA-encoded genome makes SARS-CoV-2 vulnerable to RNA interference (RNAi), we designed and analysed eight siRNAs (siV1-8) directly targeting the Orf1a/b region of the SARS-CoV-2 RNA genome, encoding for non-structural proteins (nsp). As a significant hallmark of this study, we identified siV1 (siRNA against leader protein of SARS-CoV-2), which targets the nsp1-encoding sequence (a.k.a. ‘host shutoff factor’) as particularly efficient. SiV1 inhibited SARS-CoV-2 replication in Vero E6 or Huh-7 cells by more than 99% or 97%, respectively. It neither led to toxic effects nor induced type I or III interferon production. Of note, sequence analyses revealed the target sequence of siV1 to be highly conserved in SARS-CoV-2 variants. Thus, our results identify the direct targeting of the viral RNA genome (ORF1a/b) by siRNAs as highly efficient and introduce siV1 as a particularly promising drug candidate for therapeutic intervention.

SARS-CoV-2 Outbreak Dynamics in an Isolated US Military Recruit Training Center With Rigorous Prevention Measures

This study conduct viral genome sequencing among severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected pregnant and postpartum individuals, and investigates disease severity and maternal and perinatal outcomes considering variant of concern (VOC) and non-VOC groups. The study was designed as a prospective cohort study of unvaccinated pregnant and postpartum women with SARS-CoV-2 infection diagnosed by quantitative reverse transcription polymerase chain reaction between May 2020 and April 2021 at a maternal referral center. Initially, 111 participants were identified, and 50 presented with conditions for genome sequencing and further analysis. The SARS-CoV-2 genome was subjected to next-generation sequencing in eligible cases. Sociodemographic characteristics, background, and clinical and obstetric outcomes were compared according to the SARS-CoV-2 genomic data (VOC vs. non-VOC). Among the sequenced cases, 18 were classified as VOC, among which 14 were Gamma cases and four Alpha. Oxygen desaturation, hospitalization, intensive care unit admission, invasive ventilation, and maternal death were significantly higher in SARS-CoV-2 VOC-infected than in non-VOC-infected cases. All three cases of maternal death were classified as VOC, including two Gamma and one Alpha. Pregnant and postpartum people presented with an increased morbidity and mortality due to COVID-19 caused by SARS-CoV-2 VOC, especially Gamma. The increased risk of clinical severity in the unvaccinated obstetric population underscores the importance of addressing the repercussions of the SARS-CoV-2 genomic evolution.

Coronavirus disease 2019 (COVID-19) is an emerging life-threatening pulmonary disease caused by infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which originated in Wuhan, Hubei Province, China, in December 2019. COVID-19 develops after close contact via inhalation of respiratory droplets containing SARS-CoV-2 during talking, coughing, or sneezing by asymptomatic, presymptomatic, and symptomatic carriers. This virus evolved over time, and numerous genetic variants have been reported to have increased disease severity, mortality, and transmissibility. Variants have also developed resistance to antivirals and vaccination and can escape the immune response of humans. Reverse transcription polymerase chain reaction (RT–PCR) is the method of choice among diagnostic techniques, including nucleic acid amplification tests (NAATs), serological tests, and diagnostic imaging, such as computed tomography (CT). The limitation of RT–PCR is that it cannot distinguish fragmented RNA genomes from live transmissible viruses. Thus, SARS-CoV-2 isolation by using cell culture has been developed and makes important contributions in the field of diagnosis, development of antivirals, vaccines, and SARS-CoV-2 virology research. In this research, two SARS-CoV-2 strains were isolated from four RT–PCR-positive nasopharyngeal swabs using VERO E6 cell culture. One isolate was cultured successfully with a blind passage on day 3 post inoculation from a swab with a Ct > 35, while the cells did not develop cytopathic effects without a blind passage until day 14 post inoculation. Our results indicated that infectious SARS-CoV-2 virus particles existed, even with a Ct > 35. Cultivable viruses could provide additional consideration for releasing the patient from quarantine. The results of the whole genome sequencing and bioinformatic analysis suggested that these two isolates contain a spike 68-76del+spike 675-679del double-deletion variation. The double deletion was confirmed by amplification of the regions spanning the spike gene deletion using Sanger sequencing. Phylogenetic analysis revealed that this double-deletion variant was rare (one per million in public databases, including GenBank and GISAID). The impact of this double deletion in the spike gene on the SARS-CoV-2 virus itself as well as on cultured cells and/or humans remains to be further elucidated.

Countries around the world are gearing for the transition of the coronavirus disease 2019 (COVID-19) from pandemic to endemic phase but the emergence of new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants could lead to a prolonged pandemic. SARS-CoV-2 has continued to evolve as it optimizes its adaptation to the human host and the successive waves of COVID-19 have been linked to the explosion of particular variant of concern. As the genetic diversity and epidemiological landscape of SARS-CoV-2 differ from country to country, this study aims to provide insights into the variants that are circulating in Malaysia. Whole genome sequencing was performed for 204 SARS-CoV-2 from COVID-19 cases and an additional 18,667 SARS-CoV-2 genome sequences were retrieved from the GISAID EpiCoV database for clade, lineage and genetic variation analyses. Complete genome sequences with high coverage were then used for phylogeny investigation and the resulting phylogenetic tree was constructed from 8,716 sequences. We found that the different waves of COVID-19 in Malaysia were dominated by different clades with the L and O clade for first and second wave, respectively, whereas the progressive replacement by G, GH, and GK of the GRA clade were observed in the subsequence waves. Continuous monitoring of the genetic diversity of SARS-CoV-2 is important to identify the emergence and dominance of new variant in different locality so that the appropriate countermeasures can be taken to effectively contain the spread of SARS-CoV-2.

The outbreak of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is spreading fast worldwide. There is a pressing need to understand how the virus counteracts host innate immune responses. Deleterious clinical manifestations of coronaviruses have been associated with virus-induced direct dysregulation of innate immune responses occurring via viral macrodomains located within nonstructural protein-3 (Nsp3). However, no substantial information is available concerning the relationship of macrodomains to the unusually high pathogenicity of SARS-CoV-2. Here, we show that structural evolution of macrodomains may impart a critical role to the unique pathogenicity of SARS-CoV-2. Using sequence, structural, and phylogenetic analysis, we identify a specific set of historical substitutions that recapitulate the evolution of the macrodomains that counteract host immune response. These evolutionary substitutions may alter and reposition the secondary structural elements to create new intra-protein contacts and, thereby, may enhance the ability of SARS-CoV-2 to inhibit host immunity. Further, we find that the unusual virulence of this virus is potentially the consequence of Darwinian selection‐driven epistasis in protein evolution. Our findings warrant further characterization of macrodomain-specific evolutionary substitutions in in vitro and in vivo models to determine their inhibitory effects on the host immune system.

Novel variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continue to emerge as the coronavirus disease 2019 (COVID-19) pandemic extends into its fourth year. Understanding SARS-CoV-2 circulation in university populations is vital for effective interventions in higher education settings and will inform public health policy during pandemics. In this study, we performed whole-genome sequencing of 537 of 1717 SARS-CoV-2-positive nasopharyngeal/nasal swab samples collected over a nearly 20-month period from two university populations in Wisconsin, USA. We observed that the viral sequences were distributed into 57 lineages/sub-lineages belonging to 15 clades, of which the majority were from 21K (omicron, 36.13 %) and 21J (delta, 30.91 %). Nearly 40 % (213) of the sequences were omicron, of which BA.1 and its eight descendent lineages accounted for 91 %, while the remaining belonged to BA.2 and its six descendent lineages. Independent analysis of the sequences from these two universities revealed significant differences in the circulating SARS-CoV-2 variants. Phylogenetic analysis of university sequences with a global sub-dataset demonstrated that the sequences of the same lineages from the university populations were more closely related. Genome-based analysis of closely related strains, along with phylogenetic clusters and mutational differences, identified that potential virus transmission occurred within and between universities, as well as between the university and the local community. Although this study improves our understanding of the distinct transmission patterns of circulating variants in local universities, expanding genomic surveillance capacity will aid local jurisdictions not only in identifying emerging SARS-CoV-2 variants, but also in improving data-driven public health mitigation and policy efforts.

Unique SARS-CoV protein nsp1: bioinformatics, biochemistry and potential effects on virulence

In the traditional view of molecular evolution, the rate of point mutation is uniform over the genome of an organism and variation in the rate of nucleotide substitution among DNA regions reflects differential selective constraints. Here we provide evidence for significant variation in mutation rate among regions in the mammalian genome. We show first that substitutions at silent (degenerate) sites in protein-coding genes in mammals seem to be effectively neutral (or nearly so) as they do not occur significantly less frequently than substitutions in pseudogenes. We then show that the rate of silent substitution varies among genes and is correlated with the base composition of genes and their flanking DNA. This implies that the variation in both silent substitution rate and base composition can be attributed to systematic differences in the rate and pattern of mutation over regions of the genome. We propose that the differences arise because mutation patterns vary with the timing of replication of different chromosomal regions in the germline. This hypothesis can account for both the origin of isochores in mammalian genomes and the observation that silent nucleotide substitutions in different mammalian genes do not have the same molecular clock.

Along with real-time reverse transcriptase PCR (RT-PCR) diagnostic testing, whole-genomic sequencing (WGS) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been an irreplaceable tool for epidemiological evaluation, and genomic sequence analysis has provided essential information in the development of antiviral therapeutics and vaccines (1, 2). Good sequencing coverage needs to be obtained to get information: the default threshold for lineage call at Pangolin is at least 50% of non-N bases (3). Obtaining sufficient WGS coverage from low-viral-load samples is challenging and often not subject to WGS in many labs. However, WGS of low-viral-load patient samples will enable a more complete picture of viral transmission and viral evolution (4).