Genomic Surveillance and Phylodynamic Analyses Reveal the Emergence of Novel Mutations and Co-mutation Patterns Within SARS-CoV-2 Variants Prevalent in India

Optimizing next-generation sequencing efficiency in clinical settings: analysis of read length impact on cost and performance

BACKGROUND

The expansion of sequencing technologies as a result of the response to the COVID-19 pandemic enabled pathogen (meta)genomics to be deployed as a routine component of surveillance in many countries. Scaling genomic surveillance, however, comes with associated costs in both equipment and sequencing reagents, which should be optimized. Here, we evaluate the cost efficiency and performance of different read lengths in identifying pathogens in metagenomic samples. We carefully evaluated performance metrics, costs, and time requirements relative to choices of 75, 150 and 300 base pairs (bp) read lengths in pathogen identification.

RESULTS

Our findings revealed that moving from 75 bp to 150 bp read length approximately doubles both the cost and sequencing time. Opting for 300 bp reads leads to approximately two- and three-fold increases, respectively, in cost and sequencing time compared to 75 bp reads. For viral pathogen detection, the sensitivity median ranged from 99% with 75 bp reads to 100% with 150-300 bp reads. However, bacterial pathogens detection was less effective with shorter reads: 87% with 75 bp, 95% with 150 bp, and 97% with 300 bp reads. These findings were consistent across different levels of taxa abundance. The precision of pathogen detection using shorter reads was comparable to that of longer reads across most viral and bacterial taxa.

CONCLUSIONS

During disease outbreak situations, when swift responses are required for pathogen identification, we suggest prioritizing 75 bp read lengths, especially if detection of viral pathogens is aimed. This practical approach allows better use of resources, enabling the sequencing of more samples using streamlined workflows, while maintaining a reliable response capability.

The TRAF3-DYRK1A-RAD54L2 complex maintains ACE2 expression to promote SARS-CoV-2 infection

Angiotensin converting enzyme 2 (ACE2), the host receptor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, is differentially expressed in a wide variety of tissues and cell types. The expression of ACE2 is under tight regulation, but the mechanisms regulating ACE2 expression have not yet been well defined. Through a genome-wide CRISPR knockout screen, we discovered that host factors TRAF3, DYRK1A, and RAD54L2 (TDR) form a complex to regulate the expression of ACE2. Knockout of TRAF3, DYRK1A, or RAD54L2 reduces the mRNA levels of ACE2 and inhibits the cellular entry of SARS-CoV-2. On the other hand, SARS-CoV-2 continuously evolves by genetic mutations for the adaption to the host. We have identified mutations in spike (S) (P1079T) and nucleocapsid (N) (S194L) that enhance the replication of SARS-CoV-2 in cells that express ACE2 at a low level. Our results have revealed the mechanisms for the transcriptional regulation of ACE2 and the adaption of SARS-CoV-2.

IMPORTANCE

The expression of ACE2 is essential for the entry of SARS-CoV-2 into host cells. We identify a new complex-the TDR complex-that acts to maintain the abundance of ACE2 in host cells. The identification and characterization of the TDR complex provide new targets for the development of therapeutics against SARS-CoV-2 infection. By analysis of SARS-CoV-2 virus replicating in cells expressing low levels of ACE2, we identified mutations in spike (P1079T) and nucleocapsid (S194L) that overcome the restriction of limited ACE2. Functional analysis of these key amino acids in S and N extends our knowledge of the impact of SARS-CoV-2 variants on virus infection and transmission.

Emergence of a unique SARS-CoV-2 Delta sub-cluster harboring a constellation of co-appearing non-Spike mutations

Accumulation of diverse mutations across the structural and nonstructural genes is leading to rapid evolution of SARS-CoV-2, altering its pathogenicity. We performed whole genome sequencing of 239 SARS-CoV-2 RNA samples collected from both adult and pediatric patients across eastern India (West Bengal), during the second pandemic wave in India (April-May 2021). In addition to several common spike mutations within the Delta variant, a unique constellation of eight co-appearing non-Spike mutations was identified, which revealed a high degree of positive mutual correlation. Our results also demonstrated the dynamics of SARS-CoV-2 variants among unvaccinated pediatric patients. 41.4% of our studied Delta strains harbored this signature set of eight co-appearing non-Spike mutations and phylogenetically out-clustered other Delta sub-lineages like 21J, 21A, or 21I. This is the first report from eastern India that portrayed a landscape of co-appearing mutations in the non-Spike proteins, which might have led to the evolution of a distinct Delta subcluster. Accumulation of such mutations in SARS-CoV-2 may lead to the emergence of "vaccine-evading variants." Hence, monitoring of such non-Spike mutations will be significant in the formulation of any future vaccines against those SARS-CoV-2 variants that might evade the current vaccine-induced immunity, among both the pediatric and adult populations.

Evolutionary Pattern Comparisons of the SARS-CoV-2 Delta Variant in Countries/Regions with High and Low Vaccine Coverage

It has been argued that vaccine-breakthrough infections of SARS-CoV-2 would likely accelerate the emergence of novel variants with immune evasion. This study explored the evolutionary patterns of the Delta variant in countries/regions with relatively high and low vaccine coverage based on large-scale sequences. Our results showed that (i) the sequences were grouped into two clusters (L and R); the R cluster was dominant, its proportion increased over time and was higher in the high-vaccine-coverage areas; (ii) genetic diversities in the countries/regions with low vaccine coverage were higher than those in the ones with high vaccine coverage; (iii) unique mutations and co-mutations were detected in different countries/regions; in particular, common co-mutations were exhibited in highly occurring frequencies in the areas with high vaccine coverage and presented in increasing frequencies over time in the areas with low vaccine coverage; (iv) five sites on the S protein were under strong positive selection in different countries/regions, with three in non-C to U sites (I95T, G142D and T950N), and the occurring frequencies of I95T in high vaccine coverage areas were higher, while G142D and T950N were potentially immune-pressure-selected sites; and (v) mutation at the N⁶-methyladenosine site 4 on ORF7a (C27527T, P45L) was detected and might be caused by immune pressure. Our study suggested that certain variation differences existed between countries/regions with high and low vaccine coverage, but they were not likely caused by host immune pressure. We inferred that no extra immune pressures on SARS-CoV-2 were generated with high vaccine coverage, and we suggest promoting and strengthening the uptake of the COVID-19 vaccine worldwide, especially in less developed areas.

Analysis of co-occurring and mutually exclusive amino acid changes and detection of convergent and divergent evolution events in SARS-CoV-2

The inflation of SARS-CoV-2 lineages with a high number of accumulated mutations (such as the recent case of Omicron) has risen concerns about the evolutionary capacity of this virus. Here, we propose a computational study to examine non-synonymous mutations gathered within genomes of SARS-CoV-2 from the beginning of the pandemic until February 2022. We provide both qualitative and quantitative descriptions of such corpus, focusing on statistically significant co-occurring and mutually exclusive mutations within single genomes. Then, we examine in depth the distributions of mutations over defined lineages and compare those of frequently co-occurring mutation pairs. Based on this comparison, we study mutations' convergence/divergence on the phylogenetic tree. As a result, we identify 1,818 co-occurring pairs of non-synonymous mutations showing at least one event of convergent evolution and 6,625 co-occurring pairs with at least one event of divergent evolution. Notable examples of both types are shown by means of a tree-based representation of lineages, visually capturing mutations' behaviors. Our method confirms several well-known cases; moreover, the provided evidence suggests that our workflow can explain aspects of the future mutational evolution of SARS-CoV-2.

SARS-CoV-2: Emergence of New Variants and Effectiveness of Vaccines

The emergence of SARS-CoV-2 variants may cause resistance at the immunity level against current vaccines. Some emergent new variants have increased transmissibility, infectivity, hospitalization, and mortality. Since the administration of the first SARS-CoV-2 vaccine to a human in March 2020, there is an ongoing global race against SARS-CoV-2 to control the current pandemic situation. Spike (S) glycoprotein of SARS-CoV-2 is the main target for current vaccine development, which can neutralize the infection. Companies and academic institutions have developed vaccines based on the S glycoprotein, as well as its antigenic domains and epitopes, which have been proven effective in generating neutralizing antibodies. The effectiveness of SARS-CoV-2 vaccines and other therapeutics developments are limited by the new emergent variants at the global level. We have discussed the emergent variants of SARS-CoV-2 on the efficacy of developed vaccines. Presently, most of the vaccines have been tremendously effective in severe diseases. However, there are still noteworthy challenges in certifying impartial vaccines; the stories of re-infections are generating more stressful conditions, and this needs further clinical evaluation.