Dominant clade-featured SARS-CoV-2 co-occurring mutations reveal plausible epistasis: An in silico based hypothetical model

Alternative genome sequencing approaches of SARS-CoV-2 using Ion AmpliSeq Technology

A thorough understanding of SARS-CoV-2 genetic features is compulsory to track the ongoing pandemic across multiple geographical locations of the world. Thermo Fisher Scientific USA has developed the Ion AmpliSeq SARS-CoV-2 Research Panel for the targeted sequencing of SARS-CoV-2 complete genome with high coverage and lower error rate. In this study an alternative approach of complete genome sequencing has been validated using different commercial sequencing kits to sequence the SARS-CoV-2. Amplification of cDNA with the SARS-CoV-2 primer pool was performed separately using two different master mixes: 2X environmental master mix (EM) and Platinum™ PCR SuperMix High Fidelity master mix (PM) instead of 5X Ion AmpliSeq™ HiFi Mix whereas NEBNext® Fast DNA Library Prep Set for Ion Torrent™ kit was used as an alternative to Ion AmpliSeq Library Kit Plus for other reagents. This study demonstrated a successful procedure to sequence the SARS-CoV-2 whole genome with average ∼2351 depth and 98.1% of total the reads aligned against the reference sequence (SARS-CoV-2, isolate Wuhan-Hu-1, complete genome). Although genome coverage varied, complete genomes were retrieved for both reagent sets with a reduced cost. This study proposed an alternative approach of high throughput sequencing using Ion torrent technology for the sequencing of SARS-CoV-2 in developing countries where sequencing facilities are low. This blended sequencing technique also offers a low cost protocol in developing countries like Bangladesh.

Analysis of the SARS-CoV-2 nsp12 P323L/A529V mutations: coeffect in the transiently peaking lineage C.36.3 on protein structure and response to treatment in Egyptian records

SARS-CoV-2 nsp12, the RNA-dependent RNA-polymerase plays a crucial role in virus replication. Monitoring the effect of its emerging mutants on viral replication and response to antiviral drugs is important. Nsp12 of two Egyptian isolates circulating in 2020 and 2021 were sequenced. Both isolates included P323L, one included the A529V. Tracking A529V mutant frequency, it relates to the transience peaked C.36.3 variant and its parent C.36, both peaked worldwide on February-August 2021, enlisted as high transmissible variants under investigation (VUI) on May 2021. Both Mutants were reported to originate from Egypt and showed an abrupt low frequency upon screening, we analyzed all 1104 nsp12 Egyptian sequences. A529V mutation was in 36 records with an abrupt low frequency on June 2021. As its possible reappearance might obligate actions for a candidate VUI, we analyzed the predicted co-effect of P323L and A529V mutations on protein stability and dynamics through protein structure simulations. Three available structures for drug-nsp12 interaction were used representing remdesivir, suramin and favipiravir drugs. Remdesivir and suramin showed an increase in structure stability and considerable change in flexibility while favipiravir showed an extreme interaction. Results predict a favored efficiency of the drugs except for favipiravir in case of the reported mutations.

Characterization of SARS-CoV-2 Variants in Military and Civilian Personnel of an Air Force Airport during Three Pandemic Waves in Italy

We investigated SARS-CoV-2 variants circulating, from November 2020 to March 2022, among military and civilian personnel at an Air Force airport in Italy in order to classify viral isolates in a potential hotspot for virus spread. Positive samples were subjected to Next-Generation Sequencing (NGS) of the whole viral genome and Sanger sequencing of the spike coding region. Phylogenetic analysis classified viral isolates and traced their evolutionary relationships. Clusters were identified using 70% cut-off. Sequencing methods yielded comparable results in terms of variant classification. In 2020 and 2021, we identified several variants, including B.1.258 (4/67), B.1.177 (9/67), Alpha (B.1.1.7, 9/67), Gamma (P.1.1, 4/67), and Delta (4/67). In 2022, only Omicron and its sub-lineage variants were observed (37/67). SARS-CoV-2 isolates were screened to detect naturally occurring resistance in genomic regions, the target of new therapies, comparing them to the Wuhan Hu-1 reference strain. Interestingly, 2/30 non-Omicron isolates carried the G15S 3CLpro substitution responsible for reduced susceptibility to protease inhibitors. On the other hand, Omicron isolates carried unusual substitutions A1803V, D1809N, and A949T on PLpro, and the D216N on 3CLpro. Finally, the P323L substitution on RdRp coding regions was not associated with the mutational pattern related to polymerase inhibitor resistance. This study highlights the importance of continuous genomic surveillance to monitor SARS-CoV-2 evolution in the general population, as well as in restricted communities.

Modular characterization of SARS-CoV-2 nucleocapsid protein domain functions in nucleocapsid-like assembly

SARS-CoV-2 and its variants, with the Omicron subvariant XBB currently prevailing the global infections, continue to pose threats on public health worldwide. This non-segmented positive-stranded RNA virus encodes the multi-functional nucleocapsid protein (N) that plays key roles in viral infection, replication, genome packaging and budding. N protein consists of two structural domains, NTD and CTD, and three intrinsically disordered regions (IDRs) including the N_IDR, the serine/arginine rich motif (SR_IDR), and the C_IDR. Previous studies revealed functions of N protein in RNA binding, oligomerization, and liquid-liquid phase separation (LLPS), however, characterizations of individual domains and their dissected contributions to N protein functions remain incomplete. In particular, little is known about N protein assembly that may play essential roles in viral replication and genome packing. Here, we present a modular approach to dissect functional roles of individual domains in SARS-CoV-2 N protein that reveals inhibitory or augmented modulations of protein assembly and LLPS in the presence of viral RNAs. Intriguingly, full-length N protein (N_FL) assembles into ring-like architecture whereas the truncated SR_IDR-CTD-C_IDR (N_182-419) promotes filamentous assembly. Moreover, LLPS droplets of N_FL and N_182-419 are significantly enlarged in the presence of viral RNAs, and we observed filamentous structures in the N_182-419 droplets using correlative light and electron microscopy (CLEM), suggesting that the formation of LLPS droplets may promote higher-order assembly of N protein for transcription, replication and packaging. Together this study expands our understanding of the multiple functions of N protein in SARS-CoV-2.

Evaluation of RT-PCR assays for detection of SARS-CoV-2 variants of concern

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) pandemic has been considered with great importance on correct screening procedure. The detection efficiency of recent variants of concern were observed by comparing 5 commercial RT-PCR kits and a SYBR-green method developed and validated in our laboratory. The RNA was extracted from nasopharyngeal samples from suspected COVID-19 patients and RT-PCR assay was performed according to the instruction of the respective manufacturers. The specificity and sensitivity of Maccura kit was 81.8% and 82.5%, A*Star kit was 100% and 75.4%, Da An Gene kit was 100% and 68.4%, Sansure kit was 54.5% and 91.2% and TaqPath kit was 100% and 70.2% respectively. Our in house SYBR-Green method showed a consistent detection result with 90.9% specificity and 91.2% sensitivity. We also found that detection kits targeting more genes showed better accuracy which facilitates less false positive results (< 20%). Our study found a significant difference (p < 0.005) in Ct value reported for common target genes shared by the RT-PCR kits in relation with different variants of COVID-19 infection. Recent variants of concerns contain more than 30 mutations in the spike proteins including 2 deletion and a unique insertion mutation by which makes detection of these variants difficult and these facilitates the variants to escape from being detected.

Genomic Analysis of SARS-CoV-2 Alpha, Beta and Delta Variants of Concern Uncovers Signatures of Neutral and Non-Neutral Evolution

Due to the emergence of new variants of the SARS-CoV-2 coronavirus, the question of how the viral genomes evolved, leading to the formation of highly infectious strains, becomes particularly important. Three major emergent strains, Alpha, Beta and Delta, characterized by a significant number of missense mutations, provide a natural test field. We accumulated and aligned 4.7 million SARS-CoV-2 genomes from the GISAID database and carried out a comprehensive set of analyses. This collection covers the period until the end of October 2021, i.e., the beginnings of the Omicron variant. First, we explored combinatorial complexity of the genomic variants emerging and their timing, indicating very strong, albeit hidden, selection forces. Our analyses show that the mutations that define variants of concern did not arise gradually but rather co-evolved rapidly, leading to the emergence of the full variant strain. To explore in more detail the evolutionary forces at work, we developed time trajectories of mutations at all 29,903 sites of the SARS-CoV-2 genome, week by week, and stratified them into trends related to (i) point substitutions, (ii) deletions and (iii) non-sequenceable regions. We focused on classifying the genetic forces active at different ranges of the mutational spectrum. We observed the agreement of the lowest-frequency mutation spectrum with the Griffiths-Tavaré theory, under the Infinite Sites Model and neutrality. If we widen the frequency range, we observe the site frequency spectra much more consistently with the Tung-Durrett model assuming clone competition and selection. The coefficients of the fitting model indicate the possibility of selection acting to promote gradual growth slowdown, as observed in the history of the variants of concern. These results add up to a model of genomic evolution, which partly fits into the classical drift barrier ideas. Certain observations, such as mutation "bands" persistent over the epidemic history, suggest contribution of genetic forces different from mutation, drift and selection, including recombination or other genome transformations. In addition, we show that a "toy" mathematical model can qualitatively reproduce how new variants (clones) stem from rare advantageous driver mutations, and then acquire neutral or disadvantageous passenger mutations which gradually reduce their fitness so they can be then outcompeted by new variants due to other driver mutations.

Prolonged shedding of infectious viruses with haplotype switches of SARS-CoV-2 in an immunocompromised patient

A concern has been raised that the persistent COVID-19 infection in an immunocompromised host can be the source of the SARS-CoV-2 variants. This is the case of a 61-year-old man in complete remission of a follicular lymphoma after six cycles of rituximab and bendamustine with additional two cycles of rituximab completed eight months prior to the episode of COVID-19 pneumonia. The patient's respiratory failure was long-lasting, and required mechanical ventilation until day 75. Acquired immunity tested negative throughout the observational period. The viral RNA was detectable until day 100 while the infectious virus was isolated until day 79. Seven haplotypes were identified and the non-synonymous mutations accumulated in the spike gene which included E484Q and S494P. In the management of COVID-19 cases with suppressed immune statuses, initial evaluation of existing immunity and monitoring for infectiousness throughout the clinical course including the convalescent stage may be necessary.

Plasticity in structure and assembly of SARS-CoV-2 nucleocapsid protein

Worldwide SARS-CoV-2 sequencing efforts track emerging mutations in its spike protein, as well as characteristic mutations in other viral proteins. Besides their epidemiological importance, the observed SARS-CoV-2 sequences present an ensemble of viable protein variants, and thereby a source of information on viral protein structure and function. Charting the mutational landscape of the nucleocapsid (N) protein that facilitates viral assembly, we observe variability exceeding that of the spike protein, with more than 86% of residues that can be substituted, on average by three to four different amino acids. However, mutations exhibit an uneven distribution that tracks known structural features but also reveals highly protected stretches of unknown function. One of these conserved regions is in the central disordered linker proximal to the N-G215C mutation that has become dominant in the Delta variant, outcompeting G215 variants without further spike or N-protein substitutions. Structural models suggest that the G215C mutation stabilizes conserved transient helices in the disordered linker serving as protein-protein interaction interfaces. Comparing Delta variant N-protein to its ancestral version in biophysical experiments, we find a significantly more compact and less disordered structure. N-G215C exhibits substantially stronger self-association, shifting the unliganded protein from a dimeric to a tetrameric oligomeric state, which leads to enhanced coassembly with nucleic acids. This suggests that the sequence variability of N-protein is mirrored by high plasticity of N-protein biophysical properties, which we hypothesize can be exploited by SARS-CoV-2 to achieve greater efficiency of viral assembly, and thereby enhanced infectivity.

Plasticity in structure and assembly of SARS-CoV-2 nucleocapsid protein

Worldwide SARS-CoV-2 sequencing efforts track emerging mutations in its spike protein, as well as characteristic mutations in other viral proteins. Besides their epidemiological importance, the observed SARS-CoV-2 sequences present an ensemble of viable protein variants, and thereby a source of information on viral protein structure and function. Charting the mutational landscape of the nucleocapsid (N) protein that facilitates viral assembly, we observe variability exceeding that of the spike protein, with more than 86% of residues that can be substituted, on average by 3-4 different amino acids. However, mutations exhibit an uneven distribution that tracks known structural features but also reveals highly protected stretches of unknown function. One of these conserved regions is in the central disordered linker proximal to the N-G215C mutation that has become dominant in the Delta variant, outcompeting G215 variants without further spike or N-protein substitutions. Structural models suggest that the G215C mutation stabilizes conserved transient helices in the disordered linker serving as protein-protein interaction interfaces. Comparing Delta variant N-protein to its ancestral version in biophysical experiments, we find a significantly more compact and less disordered structure. N-G215C exhibits substantially stronger self-association, shifting the unliganded protein from a dimeric to a tetrameric oligomeric state, which leads to enhanced co-assembly with nucleic acids. This suggests that the sequence variability of N-protein is mirrored by high plasticity of N-protein biophysical properties, which we hypothesize can be exploited by SARS-CoV-2 to achieve greater efficiency of viral assembly, and thereby enhanced infectivity.