Phylo-geo-network and haplogroup analysis of 611 novel coronavirus (SARS-CoV-2) genomes from India

Epidemiological Characterization and Genetic Variation of the SARS-CoV-2 Delta Variant in Palestine

The emergence of new SARS-CoV-2 variants in Palestine highlights the need for continuous genetic surveillance and accurate screening strategies. This case series study aimed to investigate the geographic distribution and genetic variation of the SARS-CoV-2 Delta Variant in Palestine in August 2021. Samples were collected at random in August 2021 (n = 571) from eight districts in the West Bank, Palestine. All samples were confirmed as positive for COVID-19 by RT-PCR. The samples passed the quality control test and were successfully sequenced using the ARTIC protocol. The Delta Variant was revealed to have four dominant lineages: B.1.617 (19%), AY.122 (18%), AY.106 (17%), and AY.121 (13%). The study revealed eight significant purely spatial clusters (p < 0.005) distributed in the northern and southern parts of Palestine. Phylogenetic analysis of SARS-CoV-2 genomes (n = 552) showed no geographically specific clades. The haplotype network revealed three haplogroups without any geographic distribution. Chronologically, the Delta Variant peak in Palestine was shortly preceded by the one in the neighboring Israeli community and shortly followed by the peak in Jordan. In addition, the study revealed an extremely intense transmission network of the Delta Variant circulating between the Palestinian districts as hubs (SHR ≈ 0.5), with Al-Khalil, the district with the highest prevalence of COVID-19, witnessing the highest frequency of transitions. Genetic diversity analysis indicated closely related haplogroups, as haplotype diversity (Hd) is high but has low nucleotide diversity (π). However, nucleotide diversity (π) in Palestine is still higher than the global figures. Neutrality tests were significantly (p < 0.05) low, including Tajima's D, Fu-Li's F, and Fu-Li's D, suggesting one or more of the following: population expansion, selective sweep, and natural negative selection. Wright's F-statistic (Fst) showed genetic differentiation (Fst > 0.25) with low to medium gene flow (Nm). Recombination events were minimal between clusters (Rm) and between adjacent sites (Rs). The study confirms the utility of the whole genome sequence as a surveillance system to track the emergence of new SARS-CoV-2 variants for any possible geographical association and the use of genetic variation analysis and haplotype networking to delineate any minimal change or slight deviation in the viral genome from a reference strain.

Phylo-geo haplotype network-based characterization of SARS-CoV-2 strains circulating in India (2020-2022)

Background & objectives Genetic analysis of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) strains circulating in India during 2020-2022 was carried out to understand the evolution of potentially expanding and divergent clades. Methods SARS-CoV-2 sequences (n=612) randomly selected from among the sequences of samples collected through a nationwide network of Virus Research Diagnostic Laboratories during 2020 (n=1532) and Indian sequences available in Global Initiative on Sharing All Influenza Data during March 2020-March 2022 (n=53077), were analyzed using the phylo-geo haplotype network approach with reference to the Wuhan prototype sequence. Results On haplotype analysis, 420 haplotypes were revealed from 643 segregating sites among the sequences. Haplotype sharing was noted among the strains from different geographical regions. Nevertheless, the genetic distance among the viral haplotypes from different clades could differentiate the strains into distinct haplo groups regarding variant emergence. Interpretation & conclusions The haplotype analysis revealed that the G and GR clades were co-evolved and an epicentrefor the evolution of the GH, GK and GRA clades. GH was more frequently identified in northern parts of India than in other parts, whereas GK was detected less in north India than in other parts. Thus, the network analysis facilitated a detailed illustration of the pathways of evolution and circulation of SARS-CoV-2 variants.

Elucidating the evolution of monkeypox virus genomes through phylo-geo-network and haplogroup analysis

BACKGROUND

As the world settles down from the COVID-19 pandemic, many countries are faced with an unexpected outbreak of monkeypox infection. Monkeypox is a zoonotic disease caused by monkeypox virus (MPXV), which is an enveloped, double stranded DNA virus belonging to the Poxviridae family. Presently, we construct and analyze the phylo-geo-network and the corresponding haplogroups. Presently, we performed the haplogroup analysis with their defining mutations and phylogenetic lineage study along with geographical distributions with the aim to understand the evolutionary path of the MPXV across the world.

RESULTS

Information about 719 full length genomes of MPXV were collected from GISAID repository and the sequences extracted from NCBI. The alignment of 719 MPXV genomes and their subsequent analysis revealed a total of 1530 segregating sites of which 330 were parsimony informative (PI) sites. The variations had a positive value of Tajima's D statistic indicating some mutations being prevalent and hence balancing selection. A total of 39 haplogroups were observed in the phylo-geo-network and their defining mutations along with the evolutionary path has been discussed. The phylo-geo-network revealed the nodal haplogroup is represented by GISAID ID 13889450, haplogroup A1, an isolate from Germany, having a total of 296 identical sequences in the study incident across 22 countries. The localized evolution is highlighted by country specific sequences and haplogroups. USA had a total of 58 genomes and 13 haplogroups as compared to Peru (89 genomes, 7 haplogroups) and Germany (26 genomes, 6 haplogroups).

CONCLUSIONS

The evolution of MPXV can be happening in a localized manner and hence accumulation of variations in the MPXV genomes needs to be monitored in order to be prepared for any possible threats.

BA.1, BA.2 and BA.2.75 variants show comparable replication kinetics, reduced impact on epithelial barrier and elicit cross-neutralizing antibodies

The Omicron variant of SARS-CoV-2 is capable of infecting unvaccinated, vaccinated and previously-infected individuals due to its ability to evade neutralization by antibodies. With multiple sub-lineages of Omicron emerging in the last 12 months, there is inadequate information on the quantitative antibody response generated upon natural infection with Omicron variant and whether these antibodies offer cross-protection against other sub-lineages of Omicron variant. In this study, we characterized the growth kinetics of Kappa, Delta and Omicron variants of SARS-CoV-2 in Calu-3 cells. Relatively higher amounts infectious virus titers, cytopathic effect and disruption of epithelial barrier functions was observed with Delta variant whereas infection with Omicron sub-lineages led to a more robust induction of interferon pathway, lower level of virus replication and mild effect on epithelial barrier. The replication kinetics of BA.1, BA.2 and BA.2.75 sub-lineages of the Omicron variant were comparable in cell culture and natural infection in a subset of individuals led to a significant increase in binding and neutralizing antibodies to the Delta variant and all the three sub-lineages of Omicron but the level of neutralizing antibodies were lowest against the BA.2.75 variant. Finally, we show that Cu2+, Zn2+ and Fe2+ salts inhibited in vitro RdRp activity but only Cu2+ and Fe2+ inhibited both the Delta and Omicron variants in cell culture. Thus, our results suggest that high levels of interferons induced upon infection with Omicron variant may counter virus replication and spread. Waning neutralizing antibody titers rendered subjects susceptible to infection by Omicron variants and natural Omicron infection elicits neutralizing antibodies that can cross-react with other sub-lineages of Omicron and other variants of concern.

Evaluating the transmission feasibility of SARS-CoV-2 Omicron (B.1.1.529) variant to 143 mammalian hosts: insights from S protein RBD and host ACE2 interaction studies

In comparison to previously known severe respiratory syndrome coronavirus 2 (SARS-CoV-2) variants, the newly emerged Omicron (B.1.1.529) variant shows higher infectivity in humans. Exceptionally high infectivity of this variant raises concern of its possible transmission via other intermediate hosts. The SARS-CoV-2 infectivity is established via the association of spike (S) protein receptor binding domain (RBD) with host angiotensin I converting enzyme 2 (hACE2) receptor. In the course of this study, we investigated the interaction between Omicron S protein RBD with the ACE2 receptor of 143 mammalian hosts including human by protein-protein interaction analysis. The goal of this study was to forecast the likelihood that the virus may infect other mammalian species that coexist with or are close to humans in the household, rural, agricultural, or zoological environments. The Omicron RBD was found to interact with higher binding affinity with the ACE2 receptor of 122 mammalian hosts via different amino acid residues from the human ACE2 (hACE2). The rat (Rattus rattus) ACE2 was found to show the strongest interaction with Omicron RBD with a binding affinity of -1393.6 kcal/mol. These distinct strong binding affinity of RBD of Omicron with host ACE2 indicates a greater potential of new host transmissibility and infection via intermediate hosts. Though expected but the phylogenetic position of the mammalian species may not dictate the Omicron RBD binding to the host ACE2 receptor suggesting an involvement of multiple factors in guiding host divergence of the variant.

Analyzing the interaction of human ACE2 and RBD of spike protein of SARS-CoV-2 in perspective of Omicron variant

The newly identified Omicron (B.1.1.529) variant of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has steered concerns across the world due to the possession of a large number of mutations leading to high infectivity and vaccine escape potential. The Omicron variant houses 32 mutations in spike (S) protein alone. The viral infectivity is determined mainly by the ability of S protein Receptor Binding Domain (RBD) to bind to the human Angiotensin I Converting Enzyme 2 (hACE2) receptor. In this paper, the interaction of the RBDs of SARS-CoV-2 variants with hACE2 was analyzed by using protein-protein docking and compared with the novel Omicron variant. Our findings reveal that the Omicron RBD interacts strongly with hACE2 receptor via unique amino acid residues as compared to the Wuhan and many other variants. However, the interacting residues of RBD are found to be the same in Lamda (C.37) variant. This unique binding of Omicron RBD with hACE2 suggests an increased potential of infectivity and vaccine evasion potential of the new variant. The evolutionary drive of the SARS-CoV-2 may not be exclusively driven by RBD variants but surely provides for the platform for emergence of new variants.

Development of an efficient Sanger sequencing-based assay for detecting SARS-CoV-2 spike mutations

Novel strains of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) harboring nucleotide changes (mutations) in the spike gene have emerged and are spreading rapidly. These mutations are associated with SARS-CoV-2 transmissibility, virulence, or resistance to some neutralizing antibodies. Thus, the accurate detection of spike mutants is crucial for controlling SARS-CoV-2 transmission and identifying neutralizing antibody-resistance caused by amino acid changes in the receptor-binding domain. Here, we developed five SARS-CoV-2 spike gene primer pairs (5-SSG primer assay; 69S, 144S, 417S, 484S, and 570S) and verified their ability to detect nine key spike mutations (ΔH69/V70, T95I, G142D, ΔY144, K417T/N, L452R, E484K/Q, N501Y, and H655Y) using a Sanger sequencing-based assay. The 5-SSG primer assay showed 100% specificity and a conservative limit of detection with a median tissue culture infective dose (TCID₅₀) values of 1.4 × 10² TCID₅₀/mL. The accuracy of the 5-SSG primer assay was confirmed by next generation sequencing. The results of these two approaches showed 100% consistency. Taken together, the ability of the 5-SSG primer assay to accurately detect key SARS-CoV-2 spike mutants is reliable. Thus, it is a useful tool for detecting SARS-CoV-2 spike gene mutants in a clinical setting, thereby helping to improve the management of patients with COVID-19.

Differential mutation profile of SARS-CoV-2 proteins across deceased and asymptomatic patients

BACKGROUND

The SARS-CoV-2 infection has spread at an alarming rate with many places showing multiple peaks in incidence. Present study analyzes a total of 332 SARS-CoV-2 genome sequences from 114 asymptomatic and 218 deceased patients from twenty-one different countries to assess the mutation profile therein in order to establish the correlation between the clinical status and the observed mutations.

METHODS

The mining of mutations was carried out using the GISAID CoVSurver (www.gisaid.org/epiflu-applications/covsurver-mutations-app) with the reference sequence 'hCoV-19/Wuhan/WIV04/2019' present in NCBI with Accession number NC-045512.2. The impact of the mutations on SARS-CoV-2 proteins mutation was predicted using PredictSNP1(loschmidt.chemi.muni.cz/predictsnp1) which is a meta-server integrating six predictor tools: SIFT, PhD-SNP, PolyPhen-1, PolyPhen-2, MAPP and SNAP. The iStable integrated server (predictor.nchu.edu.tw/iStable) was used to predict shifts in the protein stability due to mutations.

RESULTS

A total of 372 variants were observed in the 332 SARS-CoV-2 sequences with several variants present in multiple patients accounting for a total of 1596 incidences. Asymptomatic and deceased specific mutants constituted 32% and 62% of the repertoire respectively indicating their partial exclusivity. However, the most prevalent mutations were those present in both. Though some parts of the genome are more variable than others but there was clear difference between incidence and prevalence. Non-structural protein 3 (NSP3) with 68 variants had a total of only 105 incidences whereas Spike protein had 346 incidences with just 66 variants. Amongst the Deleterious variants, NSP3 had the highest incidence of 25 followed by NSP2 (16), ORF3a (14) and N (14). Spike protein had just 7 Deleterious variants out of 66.

CONCLUSION

Deceased patients have more Deleterious than Neutral variants as compared to the asymptomatic ones. Further, it appears that the Deleterious variants which decrease protein stability are more significant in pathogenicity of SARS-CoV-2.