Clinical and whole genome characterization of SARS-CoV-2 in India

Genomic epidemiology of SARS-CoV-2 from Uttar Pradesh, India

The various strains and mutations of SARS-CoV-2 have been tracked using several forms of genomic classification systems. The present study reports high-throughput sequencing and analysis of 99 SARS-CoV-2 specimens from Western Uttar Pradesh using sequences obtained from the GISAID database, followed by phylogeny and clade classification. Phylogenetic analysis revealed that Omicron lineages BA-2-like (55.55%) followed by Delta lineage-B.1.617.2 (45.5%) were predominantly circulating in this area Signature substitution at positions S: N501Y, S: D614G, S: T478K, S: K417N, S: E484A, S: P681H, and S: S477N were commonly detected in the Omicron variant-BA-2-like, however S: D614G, S: L452R, S: P681R and S: D950N were confined to Delta variant-B.1.617.2. We have also identified three escape variants in the S gene at codon position 19 (T19I/R), 484 (E484A/Q), and 681 (P681R/H) during the fourth and fifth waves in India. Based on the phylogenetic diversification studies and similar changes in other lineages, our analysis revealed indications of convergent evolution as the virus adjusts to the shifting immunological profile of its human host. To the best of our knowledge, this study is an approach to comprehensively map the circulating SARS-CoV-2 strains from Western Uttar Pradesh using an integrated approach of whole genome sequencing and phylogenetic analysis. These findings will be extremely valuable in developing a structured approach toward pandemic preparedness and evidence-based intervention plans in the future.

Whole-Genome Comparison of Representatives of All Variants of SARS-CoV-2, Including Subvariant BA.2 and the GKA Clade

Since its discovery at the end of 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly evolved into many variants, including the subvariant BA.2 and the GKA clade. Genomic clarification is needed for better management of the current pandemic as well as the possible reemergence of novel variants. The sequence of the reference genome Wuhan-Hu-1 and approximately 20 representatives of each variant were downloaded from GenBank and GISAID. Two representatives with no track of in-definitive nucleotides were selected. The sequences were aligned using muscle. The location of insertion/deletion (indel) in the genome was mapped following the open reading frame (ORF) of Wuhan-Hu-1. The phylogeny of the spike protein coding region was constructed using the maximum likelihood method. Amino acid substitutions in all ORFs were analyzed separately. There are two indel sites in ORF1AB, eight in spike, and one each in ORF3A, matrix (MA), nucleoprotein (NP), and the 3'-untranslated regions (3'UTR). Some indel sites and residues/substitutions are not unique, and some are variant-specific. The phylogeny shows that Omicron, Deltacron, and BA2 are clustered together and separated from other variants with 100% bootstrap support. In conclusion, whole-genome comparison of representatives of all variants revealed indel patterns that are specific to SARS-CoV-2 variants or subvariants. Polymorphic amino acid comparison across all coding regions also showed amino acid residues shared by specific groups of variants. Finally, the higher transmissibility of BA.2 might be due at least in part to the 48 nucleotide deletions in the 3'UTR, while the seem-to-be extinction of GKA clade is due to the lack of genetic advantages as a consequence of amino acid substitutions in various genes.

Proficiency testing for SARS-CoV-2 whole genome sequencing

Extensive studies and analyses into the molecular features of severe acute respiratory syndrome related coronavirus 2 (SARS-CoV-2) have enhanced the surveillance and investigation of its clusters and transmission worldwide. The whole genome sequencing (WGS) approach is crucial in identifying the source of infection and transmission routes by monitoring the emergence of variants over time and through communities. Varying SARS-CoV-2 genomics capacity and capability levels have been established in public health laboratories across different Australian states and territories. Therefore, laboratories performing SARS-CoV-2 WGS for public health purposes are recommended to participate in an external proficiency testing program (PTP). This study describes the development of a SARS-CoV-2 WGS PTP. The PTP assessed the performance of laboratories while providing valuable insight into the current state of SARS-CoV-2 genomics in public health across Australia. Part 1 of the PTP contained eight simulated SARS-CoV-2 positive and negative specimens to assess laboratories' wet and dry laboratory capacity. Part 2 involved the analysis of a genomic dataset that consisted of a multi-FASTA file of 70 consensus genomes of SARS-CoV-2. Participating laboratories were required to (1) submit raw data for independent bioinformatics analysis, (2) analyse the data with their processes, and (3) answer relevant questions about the data. The performance of the laboratories was commendable, despite some variation in the reported results due to the different sequencing and bioinformatics approaches used by laboratories. The overall outcome is positive and demonstrates the critical role of the PTP in supporting the implementation and validation of SARS-CoV-2 WGS processes. The data derived from this PTP will contribute to the development of SARS-CoV-2 bioinformatic quality control (QC) and performance benchmarking for accreditation.

Multifaceted Assessment of Wastewater-Based Epidemiology for SARS-CoV-2 in Selected Urban Communities in Davao City, Philippines: A Pilot Study

Over 60 countries have integrated wastewater-based epidemiology (WBE) in their COVID-19 surveillance programs, focusing on wastewater treatment plants (WWTP). In this paper, we piloted the assessment of SARS-CoV-2 WBE as a complementary public health surveillance method in susceptible communities in a highly urbanized city without WWTP in the Philippines by exploring the extraction and detection methods, evaluating the contribution of physico-chemical-anthropogenic factors, and attempting whole-genome sequencing (WGS). Weekly wastewater samples were collected from sewer pipes or creeks in six communities with moderate-to-high risk of COVID-19 transmission, as categorized by the City Government of Davao from November to December 2020. Physico-chemical properties of the wastewater and anthropogenic conditions of the sites were noted. Samples were concentrated using a PEG-NaCl precipitation method and analyzed by RT-PCR to detect the SARS-CoV-2 N, RdRP, and E genes. A subset of nine samples were subjected to WGS using the Minion sequencing platform. SARS-CoV-2 RNA was detected in twenty-two samples (91.7%) regardless of the presence of new cases. Cycle threshold values correlated with RNA concentration and attack rate. The lack of a sewershed map in the sampled areas highlights the need to integrate this in the WBE planning. A combined analysis of wastewater physico-chemical parameters such as flow rate, surface water temperature, salinity, dissolved oxygen, and total dissolved solids provided insights on the ideal sampling location, time, and method for WBE, and their impact on RNA recovery. The contribution of fecal matter in the wastewater may also be assessed through the coliform count and in the context of anthropogenic conditions in the area. Finally, our attempt on WGS detected single-nucleotide polymorphisms (SNPs) in wastewater which included clinically reported and newly identified mutations in the Philippines. This exploratory report provides a contextualized framework for applying WBE surveillance in low-sanitation areas.

Molecular Analysis of SARS-CoV-2 Lineages in Armenia

The sequencing of SARS-CoV-2 provides essential information on viral evolution, transmission, and epidemiology. In this paper, we performed the whole-genome sequencing of SARS-CoV-2 using nanopore and Illumina sequencing to describe the circulation of the virus lineages in Armenia. The analysis of 145 full genomes identified six clades (19A, 20A, 20B, 20I, 21J, and 21K) and considerable intra-clade PANGO lineage diversity. Phylodynamic and transmission analysis allowed to attribute specific clades as well as infer their importation routes. Thus, the first two waves of positive case increase were caused by the 20B clade, the third peak caused by the 20I (Alpha), while the last two peaks were caused by the 21J (Delta) and 21K (Omicron) variants. The functional analyses of mutations in sequences largely affected epitopes associated with protective HLA loci and did not cause the loss of the signal in PCR tests targeting ORF1ab and N genes as confirmed by RT-PCR. We also compared the performance of nanopore and Illumina short-read sequencing and showed the utility of nanopore sequencing as an efficient and affordable alternative for large-scale molecular epidemiology research. Thus, our paper describes new data on the genomic diversity of SARS-CoV-2 variants in Armenia in the global context of the virus molecular genomic surveillance.

Genome Sequencing Reveals a Mixed Picture of SARS-CoV-2 Variant of Concern Circulation in Eastern Uttar Pradesh, India

Uttar Pradesh is the densely populated state of India and is the sixth highest COVID-19 affected state with 22,904 deaths recorded on November 12, 2021. Whole-genome sequencing (WGS) is being used as a potential approach to investigate genomic evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. In this study, a total of 87 SARS-CoV-2 genomes-49 genomes from the first wave (March 2020 to February 2021) and 38 genomes from the second wave (March 2021 to July 2021) from Eastern Uttar Pradesh (E-UP) were sequenced and analyzed to understand its evolutionary pattern and variants against publicaly available sequences. The complete genome analysis of SARS-CoV-2 during the first wave in E-UP largely reported transmission of G, GR, and GH clades with specific mutations. In contrast, variants of concerns (VOCs) such as Delta (71.0%) followed by Delta AY.1 (21.05%) and Kappa (7.9%) lineages belong to G clade with prominent signature amino acids were introduced in the second wave. Signature substitution at positions S:L452R, S:P681R, and S:D614G were commonly detected in the Delta, Delta AY.1, and Kappa variants whereas S:T19R and S:T478K were confined to Delta and Delta AY.1 variants only. Vaccine breakthrough infections showed unique mutational changes at position S:D574Y in the case of the Delta variant, whereas position S:T95 was conserved among Kappa variants compared to the Wuhan isolate. During the transition from the first to second waves, a shift in the predominant clade from GH to G clade was observed. The identified spike protein mutations in the SARS-CoV-2 genome could be used as the potential target for vaccine and drug development to combat the effects of the COVID-19 disease.

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

Identification of the genomic diversity and the phylodynamic profiles of prevalent variants is critical to understand the evolution and spread of SARS-CoV-2 variants. We performed whole-genome sequencing of 54 SARS-CoV-2 variants collected from COVID-19 patients in Kolkata, West Bengal during August-October 2020. Phylogeographic and phylodynamic analyses were performed using these 54 and other sequences from India and abroad that are available in the GISAID database. We estimated the clade dynamics of the Indian variants and compared the clade-specific mutations and the co-mutation patterns across states and union territories of India over the time course. Frequent mutations and co-mutations observed within the major clades across time periods do not show much overlap, indicating the emergence of newer mutations in the viral population prevailing in the country. Furthermore, we explored the possible association of specific mutations and co-mutations with the infection outcomes manifested in Indian patients.