Different selection dynamics of S and RdRp between SARS-CoV-2 genomes with and without the dominant mutations

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.

Genome analysis of SARS-CoV-2 isolates from a population reveals the rapid selective sweep of a haplotype carrying many pre-existing and new mutations

To understand the mechanism underlying the evolution of SARS-CoV-2 in a population, we sequenced 92 viral genomes from Assam, India. Analysis of these and database sequences revealed a complete selective sweep of a haplotype in Assam carrying 13 pre-existing variants, including a high leap in frequency of a variant on ORF8, which is involved in immune evasion. A comparative study between sequences of same lineage and similar time frames in and outside Assam showed that 10 of the 13 pre-existing variants had a frequency ranging from 96 to 99%, and the remaining 3 had a low frequency outside Assam. Using a phylogenetic approach to infer sequential occurrences of variants we found that the variant Phe120del on ORF8, which had a low frequency (1.75%) outside Assam, is at the base of the phylogenetic tree of variants and became totally fixed (100%) in Assam population. Based on this observation, we inferred that the variant on ORF8 had a selective advantage, so it carried the haplotype to reach the100% frequency. The haplotype also carried 32 pre-existing variants at a frequency from 1.00 to 80.00% outside Assam. Those of these variants that are more closely linked to the S-protein locus, which often carries advantageous mutations and is tightly linked to the ORF8 locus, retained higher frequencies, while the less tightly linked variants showed lower frequencies, likely due to recombination among co- circulating variants in Assam. The ratios of non-synonymous substitutions to synonymous substitutions suggested that some genes such as those coding for the S-protein and non-structural proteins underwent positive selection while others were subject to purifying selection during their evolution in Assam. Furthermore, we observed negative correlation of the Ct value of qRT-PCR of the patients with abundant ORF6 transcripts, suggesting that ORF6 can be used as a marker for estimating viral titer. In conclusion, our in-depth analysis of SARS-CoV-2 genomes in a regional population reveals the mechanism and dynamics of viral evolution.

Molecular Evolution of SARS-CoV-2 during the COVID-19 Pandemic

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) produced diverse molecular variants during its recent expansion in humans that caused different transmissibility and severity of the associated disease as well as resistance to monoclonal antibodies and polyclonal sera, among other treatments. In order to understand the causes and consequences of the observed SARS-CoV-2 molecular diversity, a variety of recent studies investigated the molecular evolution of this virus during its expansion in humans. In general, this virus evolves with a moderate rate of evolution, in the order of 10-3-10-4 substitutions per site and per year, which presents continuous fluctuations over time. Despite its origin being frequently associated with recombination events between related coronaviruses, little evidence of recombination was detected, and it was mostly located in the spike coding region. Molecular adaptation is heterogeneous among SARS-CoV-2 genes. Although most of the genes evolved under purifying selection, several genes showed genetic signatures of diversifying selection, including a number of positively selected sites that affect proteins relevant for the virus replication. Here, we review current knowledge about the molecular evolution of SARS-CoV-2 in humans, including the emergence and establishment of variants of concern. We also clarify relationships between the nomenclatures of SARS-CoV-2 lineages. We conclude that the molecular evolution of this virus should be monitored over time for predicting relevant phenotypic consequences and designing future efficient treatments.

Sequence analysis of SARS-CoV-2 Delta variant isolated from Makassar, South Sulawesi, Indonesia

Introduction

This study aimed to perform mutation and phylogenetic analyses of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Delta variants and analyze the characteristic signs and symptoms of patients infected with SARS-CoV-2 Delta variant originated from Makassar during the Delta outbreak.Methods: We collected samples from patients who were infected with coronavirus disease 2019 (COVID-19) between June and October 2021. We selected the Quantitative Reverse Transcription-Polymerase Chain Reaction (PCR)-positive samples with a cycle threshold value of <30 for whole genome sequencing. Total viral ribonucleic acid (RNA) was isolated from 34 PCR-positive nasopharyngeal swab samples, and whole genome sequencing was performed using the Oxford Nanopore GridlON sequencer. Phylogenetic and maximum clade credibility analyses were performed using the Bayesian Markov chain Monte Carlo method.

Results

It was found that 33 patients were infected with the SARS-CoV-2 Delta variant in this cohort study, among whom 63.6% (21) patients were female. According to the clinical data, 24 (72.7%), 7 (21.2%), and 2 (6.1%) patients had mild, moderate, and severe COVID-19 infections. Phylogenetic analysis based on the spike and RNA-dependent RNA polymerase (RdRp) genes showed that the collected samples were clustered in the main lineage of B.1.617.2 (Delta variant). The Delta variants had a high frequency of distinct mutations in the spike protein region, including T19R (94.12%), L452R (88.23%), T478K (91.17%), D614G (97%), P681R (97%), and D950 N (97%). Other unique mutations found in a smaller frequency in our samples were present in the N-terminal domain, including A27T (2.94%) and A222V (14.70%), and in the receptor-binding domain, including Q414K (5.88%), G446V (2.94%), and T470 N (2.94%).

Conclusion

This study revealed the unique mutations in the S protein region of Delta variants. T19R, L452R, T478K/T478R, D614G, P681R, and D950 N were the most common substitutions in Makassar's Delta variant.

Evolutionary history of the SARS-CoV-2 Gamma variant of concern (P.1): a perfect storm

Our goal was to describe in more detail the evolutionary history of Gamma and two derived lineages (P.1.1 and P.1.2), which are part of the arms race that SARS-CoV-2 wages with its host. A total of 4,977 sequences of the Gamma strain of SARS-CoV-2 from Brazil were analyzed. We detected 194 sites under positive selection in 12 genes/ORFs: Spike, N, M, E, ORF1a, ORF1b, ORF3, ORF6, ORF7a, ORF7b, ORF8, and ORF10. Some diagnostic sites for Gamma lacked a signature of positive selection in our study, but these were not fixed, apparently escaping the action of purifying selection. Our network analyses revealed branches leading to expanding haplotypes with sites under selection only detected when P.1.1 and P.1.2 were considered. The P.1.2 exclusive haplotype H_5 originated from a non-synonymous mutational step (H3509Y) in H_1 of ORF1a. The selected allele, 3509Y, represents an adaptive novelty involving ORF1a of P.1. Finally, we discuss how phenomena such as epistasis and antagonistic pleiotropy could limit the emergence of new alleles (and combinations thereof) in SARS-COV-2 lineages, maintaining infectivity in humans, while providing rapid response capabilities to face the arms race triggered by host immuneresponses.

Genomic Bootstrap Barcodes and Their Application to Study the Evolution of Sarbecoviruses

Recombination creates mosaic genomes containing regions with mixed ancestry, and the accumulation of such events over time can complicate greatly many aspects of evolutionary inference. Here, we developed a sliding window bootstrap (SWB) method to generate genomic bootstrap (GB) barcodes to highlight the regions supporting phylogenetic relationships. The method was applied to an alignment of 56 sarbecoviruses, including SARS-CoV and SARS-CoV-2, responsible for the SARS epidemic and COVID-19 pandemic, respectively. The SWB analyses were also used to construct a consensus tree showing the most reliable relationships and better interpret hidden phylogenetic signals. Our results revealed that most relationships were supported by just a few genomic regions and confirmed that three divergent lineages could be found in bats from Yunnan: SCoVrC, which groups SARS-CoV related coronaviruses from China; SCoV2rC, which includes SARS-CoV-2 related coronaviruses from Southeast Asia and Yunnan; and YunSar, which contains a few highly divergent viruses recently described in Yunnan. The GB barcodes showed evidence for ancient recombination between SCoV2rC and YunSar genomes, as well as more recent recombination events between SCoVrC and SCoV2rC genomes. The recombination and phylogeographic patterns suggest a strong host-dependent selection of the viral RNA-dependent RNA polymerase. In addition, SARS-CoV-2 appears as a mosaic genome composed of regions sharing recent ancestry with three bat SCoV2rCs from Yunnan (RmYN02, RpYN06, and RaTG13) or related to more ancient ancestors in bats from Yunnan and Southeast Asia. Finally, our results suggest that viral circular RNAs may be key molecules for the mechanism of recombination.

Identification and phylogenetic analysis of two canine coronavirus strains

Canine coronavirus (CCoV), a member of the genus Alphacoronavirus, is an enveloped, single-stranded positive-sense RNA virus that responsible for gastroenteritis in dogs. In this study, two CCoV isolates were successfully propagated from 53 CCoV-positive clinical specimens by serial passaging in A-72 cells. These two strains, CCoV JS1706 and CCoV JS1712, caused cytopathic effects in A-72 cells. The sizes of virus plaque formed by them differed in early passages. Electron microscopy revealed a large quantity of typical coronavirus particles with 80–120 nm in diameter in cell culture media and cytoplasm of infected cells, in which they appeared as inclusion bodies. RT-PCR analysis of S gene indicated that these two isolates were belonged to CCoV IIa subtype. Homology of RdRp, S, M and N proteins between the two strains were 100, 99.6, 99.2 and 100.0%, respectively, whereas they were 99.4–100%, 83.1–95.2%, 88.5–99.2% and 91.9–99.7% identity compared to CCoV II reference strains. Phylogenetic analysis of RdRp, S, M and N protein showed that they were closely related to CCoV II strains. These two subtype IIa isolates will be useful for evaluating the pathogenesis and evolution of CCoV and for developing diagnostic reagents and vaccines.

Neutralising antibody escape of SARS-CoV-2 spike protein: Risk assessment for antibody-based Covid-19 therapeutics and vaccines

The Spike protein is the target of both antibody-based therapeutics (convalescent plasma, polyclonal serum, monoclonal antibodies) and vaccines. Mutations in Spike could affect efficacy of those treatments. Hence, monitoring of mutations is necessary to forecast and readapt the inventory of therapeutics. Different phylogenetic nomenclatures have been used for the currently circulating SARS-CoV-2 clades. The Spike protein has different hotspots of mutation and deletion, the most dangerous for immune escape being the ones within the receptor binding domain (RBD), such as K417N/T, N439K, L452R, Y453F, S477N, E484K, and N501Y. Convergent evolution has led to different combinations of mutations among different clades. In this review we focus on the main variants of concern, that is, the so-called UK (B.1.1.7), South African (B.1.351) and Brazilian (P.1) strains.

Genomic characterization of SARS-CoV-2 isolates from patients in Turkey reveals the presence of novel mutations in spike and nsp12 proteins

Novel mutations have been emerging in the genome of severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2); consequently, the evolving of more virulent and treatment resistance strains have the potential to increase transmissibility and mortality rates. The characterization of full-length SARS-CoV-2 genomes is critical for understanding the origin and transmission pathways of the virus, as well as identifying mutations that affect the transmissibility and pathogenicity of the virus. We present an analysis of the mutation pattern and clade distribution of full-length SARS-CoV-2 genome sequences obtained from specimens tested at Gazi University Medical Virology Laboratory. Viral RNA was extracted from nasopharyngeal specimens. Next-generation sequencing libraries were prepared and sequenced on Illumina iSeq 100 platform. Raw sequencing data were processed to obtain full-length genome sequences and variant calling was performed to analyze amino acid changes. Clade distribution was determined to understand the phylogenetic background in relation to global data. A total of 293 distinct mutations were identified, of which 152 missense, 124 synonymous, 12 noncoding, and 5 deletions. The most frequent mutations were P323L (nsp12), D614G (ORF2/S), and 2421C>T (5'-untranslated region) found simultaneously in all sequences. Novel mutations were found in nsp12 (V111A, H133R, Y453C, M626K) and ORF2/S (R995G, V1068L). Nine different Pangolin lineages were detected. The most frequently assigned lineage was B.1.1 (17 sequences), followed by B.1 (7 sequences) and B.1.1.36 (3 sequences). Sequence information is essential for revealing genomic diversity. Mutations might have significant functional implications and analysis of these mutations provides valuable information for therapeutic and vaccine development studies. Our findings point to the introduction of the virus into Turkey through various sources and the subsequent spread of several key variants.