Missense mutations in SARS-CoV2 genomes from Indian patients

Non-uniform aspects of the SARS-CoV-2 intraspecies evolution reopen question of its origin

Several hypotheses have been presented on the origin of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) from its identification as the agent causing the current coronavirus disease 19 (COVID-19) pandemic. So far, no solid evidence has been found to support any hypothesis on the origin of this virus, and the issue continue to resurface over and over again. Here we have unfolded a pattern of distribution of several mutations in the SARS-CoV-2 proteins in 24 geo-locations across different continents. The results showed an evenly uneven distribution of the unique protein variants, distinct mutations, unique frequency of common conserved residues, and mutational residues across these 24 geo-locations. Furthermore, ample mutations were identified in the evolutionarily conserved invariant regions in the SARS-CoV-2 proteins across almost all geo-locations studied. This pattern of mutations potentially breaches the law of evolutionary conserved functional units of the beta-coronavirus genus. These mutations may lead to several novel SARS-CoV-2 variants with a high degree of transmissibility and virulence. A thorough investigation on the origin and characteristics of SARS-CoV-2 needs to be conducted in the interest of science and for the preparation of meeting the challenges of potential future pandemics.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

BepiTBR: T-B reciprocity enhances B cell epitope prediction

The ability to predict B cell epitopes is critical for biomedical research and many clinical applications. Investigators have observed the phenomenon of T-B reciprocity, in which candidate B cell epitopes with nearby CD4+ T cell epitopes have higher chances of being immunogenic. To our knowledge, existing B cell epitope prediction algorithms have not considered this interesting observation. We developed a linear B cell epitope prediction model, BepiTBR, based on T-B reciprocity. We showed that explicitly including the enrichment of putative CD4+ T cell epitopes (predicted HLA class II epitopes) in the model leads to significant enhancement in the prediction of linear B cell epitopes. Curiously, the positive impact on B cell epitope generation is specific to the enrichment of DQ allele binders. Overall, our work provides interesting mechanistic insights into the generation of B cell epitopes and points to a new avenue to improve B cell epitope prediction for the field.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

Emergence of unique SARS-CoV-2 ORF10 variants and their impact on protein structure and function

The devastating impact of the ongoing coronavirus disease 2019 (COVID-19) on public health, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has made targeting the COVID-19 pandemic a top priority in medical research and pharmaceutical development. Surveillance of SARS-CoV-2 mutations is essential for the comprehension of SARS-CoV-2 variant diversity and their impact on virulence and pathogenicity. The SARS-CoV-2 open reading frame 10 (ORF10) protein interacts with multiple human proteins CUL2, ELOB, ELOC, MAP7D1, PPT1, RBX1, THTPA, TIMM8B, and ZYG11B expressed in lung tissue. Mutations and co-occurring mutations in the emerging SARS-CoV-2 ORF10 variants are expected to impact the severity of the virus and its associated consequences. In this article, we highlight 128 single mutations and 35 co-occurring mutations in the unique SARS-CoV-2 ORF10 variants. The possible predicted effects of these mutations and co-occurring mutations on the secondary structure of ORF10 variants and host protein interactomes are presented. The findings highlight the possible effects of mutations and co-occurring mutations on the emerging 140 ORF10 unique variants from secondary structure and intrinsic protein disorder perspectives.

Amantadine has potential for the treatment of COVID-19 because it inhibits known and novel ion channels encoded by SARS-CoV-2

The dire need for COVID-19 treatments has inspired strategies of repurposing approved drugs. Amantadine has been suggested as a candidate, and cellular as well as clinical studies have indicated beneficial effects of this drug. We demonstrate that amantadine and hexamethylene-amiloride (HMA), but not rimantadine, block the ion channel activity of Protein E from SARS-CoV-2, a conserved viroporin among coronaviruses. These findings agree with their binding to Protein E as evaluated by solution NMR and molecular dynamics simulations. Moreover, we identify two novel viroporins of SARS-CoV-2; ORF7b and ORF10, by showing ion channel activity in a X. laevis oocyte expression system. Notably, amantadine also blocks the ion channel activity of ORF10, thereby providing two ion channel targets in SARS-CoV-2 for amantadine treatment in COVID-19 patients. A screen of known viroporin inhibitors on Protein E, ORF7b, ORF10 and Protein 3a from SARS-CoV-2 revealed inhibition of Protein E and ORF7b by emodin and xanthene, the latter also blocking Protein 3a. This illustrates a general potential of well-known ion channel blockers against SARS-CoV-2 and specifically a dual molecular basis for the promising effects of amantadine in COVID-19 treatment.

Toft-Bertelsen et al. describe repurposing of anti-influenza drug amantadine and its derivatives for the treatment of SARS-CoV-2. They show that Amantadine, Emodin and Xanthene show significant blockage of ionchannels formed by SARS-CoV-2 which are crucial for its assembly and pathophysiology.

Implications derived from S-protein variants of SARS-CoV-2 from six continents

The spike (S) protein is a critical determinant of the infectivity and antigenicity of SARS-CoV-2. Several mutations in the S protein of SARS-CoV-2 have already been detected, and their effect in immune system evasion and enhanced transmission as a cause of increased morbidity and mortality are being investigated. From pathogenic and epidemiological perspectives, S proteins are of prime interest to researchers. This study focused on the unique variants of S proteins from six continents: Asia, Africa, Europe, Oceania, South America, and North America. In comparison to the other five continents, Africa had the highest percentage of unique S proteins (29.1%). The phylogenetic relationship implies that unique S proteins from North America are significantly different from those of the other five continents. They are most likely to spread to the other geographic locations through international travel or naturally by emerging mutations. It is suggested that restriction of international travel should be considered, and massive vaccination as an utmost measure to combat the spread of the COVID-19 pandemic. It is also further suggested that the efficacy of existing vaccines and future vaccine development must be reviewed with careful scrutiny, and if needed, further re-engineered based on requirements dictated by new emerging S protein variants.

A next generation sequencing (NGS) analysis to reveal genomic and proteomic mutation landscapes of SARS-CoV-2 in South Asia

•

48 SNPs were identified from the genome-wide analysis of 410 South Asian SARS-CoV-2 sequences.

•

About 85% SNPs are packaged in ORF1ab, spike protein, and nucleocapsid.

•

South Asian strains are highly related to the South American and European strains according to the phylogenetic analysis.

•

Unlike other countries, frequency of 1163A>T missense mutation is very high (78.80%) in Bangladeshi samples.

Counts for SARS-CoV-2 associated infections and fatalities are on the rise globally even in regions which contained the spread momentarily. The pattern of infections has been found to be controlled by the distinctive selection pressures exerted by fluctuating environmental nature and hosts. A total of 410 whole-genome sequences submitted by the South Asian countries were retrieved from the GISAID database and analyzed to assess the impact and pattern of mutations in this region. Most common and frequent mutations in the South Asian countries are 241C > T, 3037C > T, 14408C > T, and 23403A > G and about 85% SNPs are localized in ORF1ab, spike protein, and nucleocapsid. Among the identified mutations, the proportion of missense type (54.17%) was highest, followed by the synonymous (41.66%) and the non-coding types (4.17%). While analyzing transmission source in terms of geolocation, the largest clustered group from the South Asian countries was based on the G-clade (D614G) (81.7%; 335/410 samples), tracing the inception and transmission of SARS-CoV-2 infections in the South Asian countries from European regions. Phylogenetic analysis also revealed that the South Asian strains are highly related to the South American and European strains. We found that G-clade mutations are more prevalent (96.19%) in the samples of Bangladesh which were also prevalent in the European isolates. Surprisingly, one missense mutation (1163A > T) in ORF1ab gene became dominant only in Bangladesh (78.8%), which led to debates regarding effects on the pathogenicity and transmissibility of the virus. Overall, the findings of this study highlight the frequently mutated SARS-CoV-2 variants among the COVID-19 patients in the South Asian countries which might ease containment of the disease in this region through investigating the virulence reducing factors as the identified mutations are strongly correlated with low infection and mortality rate.

A Global Mutational Profile of SARS-CoV-2: A Systematic Review and Meta-Analysis of 368,316 COVID-19 Patients

Since its first detection in December 2019, more than 232 million cases of COVID-19, including 4.7 million deaths, have been reported by the WHO. The SARS-CoV-2 viral genomes have evolved rapidly worldwide, causing the emergence of new variants. This systematic review and meta-analysis was conducted to provide a global mutational profile of SARS-CoV-2 from December 2019 to October 2020. The review was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA), and a study protocol was lodged with PROSPERO. Data from 62 eligible studies involving 368,316 SARS-CoV-2 genomes were analyzed. The mutational data analyzed showed most studies detected mutations in the Spike protein (n = 50), Nucleocapsid phosphoprotein (n = 34), ORF1ab gene (n = 29), 5′-UTR (n = 28) and ORF3a (n = 25). Under the random-effects model, pooled prevalence of SARS-CoV-2 variants was estimated at 95.1% (95% CI; 93.3–96.4%; I² = 98.952%; p = 0.000) while subgroup meta-analysis by country showed majority of the studies were conducted ‘Worldwide’ (n = 10), followed by ‘Multiple countries’ (n = 6) and the USA (n = 5). The estimated prevalence indicated a need to continuously monitor the prevalence of new mutations due to their potential influence on disease severity, transmissibility and vaccine effectiveness.

Comparative study of predicted miRNA between Indonesia and China (Wuhan) SARS-CoV-2: a bioinformatics analysis

BACKGROUND

Several reports on the discovery of SARS-CoV-2 mutations and variations in Indonesia COVID-19 cases led to genomic dysregulation with the first pandemic cases in Wuhan, China. MicroRNA (miRNA) plays an important role in this genetic regulation and contributes to the enhancement of viral RNA binding through the host mRNA.

OBJECTIVE

This research is aimed to detect miRNA targets of SARS-CoV-2 and examines their role in Indonesia cases against Wuhan cases.

METHODS

SARS-CoV-2 sequences were obtained from GISAID ( https://www.gisaid.org/ ), NCBI ( https://ncbi.nlm.nih.gov ), and National Genomics Data Center ( https://bigd.big.ac.cn/gwh/ ) databases. MiRDB ( https://github.com/gbnegrini/mirdb-custom-target-search ) was used to annotate and predict target human mature miRNAs. For statistical analysis, we utilized a series chi-square test to obtain significant miRNA. DIANA-miRPath v3.0 ( http://www.microrna.gr/miRPathv3 ) analyzed the Gene Ontology of mature miRNAs.

RESULT

The statistical results detected five significant miRNAs. Two miRNAs: hsa-miR-4778-5p and hsa-miR-4531 were consistently found in the majority of Wuhan samples, while they were only found in less than half of the Indonesia samples. The other three miRNA, hsa-miR-6844, hsa-miR-627-5p, and hsa-miR-3674, were discovered in most samples in both groups but with a significant difference ratio. Among these five significant miRNA targets, hsa-miR-6844 is the only miRNA that has an association with the ORF1ab gene of SARS-CoV-2.

CONCLUSION

The Gene Ontology analysis of five significant miRNA targets indicates a significant role in inflammation and the immune system. The specific detection of host miRNAs in this study shows that there are differences in the characteristics of SARS-CoV-2 between Indonesia and Wuhan.

Determination of Ligand Binding Modes in Hydrated Viral Ion Channels to Foster Drug Design and Repositioning

Target-based design and repositioning are mainstream strategies of drug discovery. Numerous drug design and repositioning projects have been launched to fight the ongoing COVID-19 pandemic. The resulting drug candidates have often failed due to the misprediction of their target-bound structures. The determination of water positions of such structures is particularly challenging due to the large number of possible drugs and the diversity of their hydration patterns. To answer this challenge and help correct predictions, we introduce a new protocol HydroDock, which can build hydrated drug-target complexes from scratch. HydroDock requires only the dry target and drug structures and produces their complexes with appropriately positioned water molecules. As a test application of the protocol, we built the structures of amantadine derivatives in complex with the influenza M2 transmembrane ion channel. The repositioning of amantadine derivatives from this influenza target to the SARS-CoV-2 envelope protein was also investigated. Excellent agreement was observed between experiments and the structures determined by HydroDock. The atomic resolution complex structures showed that water plays a similar role in the binding of amphipathic amantadine derivatives to transmembrane ion channels of both influenza A and SARS-CoV-2. While the hydrophobic regions of the channels capture the bulky hydrocarbon group of the ligand, the surrounding waters direct its orientation parallel with the axes of the channels via bridging interactions with the ionic ligand head. As HydroDock supplied otherwise undetermined structural details, it can be recommended to improve the reliability of future design and repositioning of antiviral drug candidates and many other ligands with an influence of water structure on their mechanism of action.

SARS-CoV-2 genomics: An Indian perspective on sequencing viral variants

Since its emergence as a pneumonia-like outbreak in the Chinese city of Wuhan in late 2019, the novel coronavirus disease COVID-19 has spread widely to become a global pandemic. The first case of COVID-19 in India was reported on 30 January 2020 and since then it has affected more than ten million people and resulted in around 150,000 deaths in the country. Over time, the viral genome has accumulated mutations as it passes through its human hosts, a common evolutionary mechanism found in all microorganisms. This has implications for disease surveillance and management, vaccines and therapeutics, and the emergence of reinfections. Sequencing the viral genome can help monitor these changes and provides an extraordinary opportunity to understand the genetic epidemiology and evolution of the virus as well as tracking its spread in a population. Here we review the past year in the context of the phylogenetic analysis of variants isolated over the course of the pandemic in India and highlight the importance of continued sequencing-based surveillance in the country.

Amantadine Inhibits SARS-CoV-2 In Vitro

Since the SARS-CoV-2 pandemic started in late 2019, the search for protective vaccines and for drug treatments has become mandatory to fight the global health emergency. Travel restrictions, social distancing, and face masks are suitable counter measures, but may not bring the pandemic under control because people will inadvertently or at a certain degree of restriction severity or duration become incompliant with the regulations. Even if vaccines are approved, the need for antiviral agents against SARS-CoV-2 will persist. However, unequivocal evidence for efficacy against SARS-CoV-2 has not been demonstrated for any of the repurposed antiviral drugs so far. Amantadine was approved as an antiviral drug against influenza A, and antiviral activity against SARS-CoV-2 has been reasoned by analogy but without data. We tested the efficacy of amantadine in vitro in Vero E6 cells infected with SARS-CoV-2. Indeed, amantadine inhibited SARS-CoV-2 replication in two separate experiments with IC50 concentrations between 83 and 119 µM. Although these IC50 concentrations are above therapeutic amantadine levels after systemic administration, topical administration by inhalation or intranasal instillation may result in sufficient amantadine concentration in the airway epithelium without high systemic exposure. However, further studies in other models are needed to prove this hypothesis.

Genomic Variations in SARS-CoV-2 Genomes From Gujarat: Underlying Role of Variants in Disease Epidemiology

Humanity has seen numerous pandemics during its course of evolution. The list includes several incidents from the past, such as measles, Ebola, severe acute respiratory syndrome (SARS), and Middle East respiratory syndrome (MERS), etc. The latest edition to this is coronavirus disease 2019 (COVID-19), caused by the novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). As of August 18, 2020, COVID-19 has affected over 21 million people from 180 + countries with 0.7 million deaths across the globe. Genomic technologies have enabled us to understand the genomic constitution of pathogens, their virulence, evolution, and rate of mutation, etc. To date, more than 83,000 viral genomes have been deposited in public repositories, such as GISAID and NCBI. While we are writing this, India is the third most affected country by COVID-19, with 2.7 million cases and > 53,000 deaths. Gujarat is the 11th highest affected state with a 3.48% death rate compared to the national average of 1.91%. In this study, a total of 502 SARS-CoV-2 genomes from Gujarat were sequenced and analyzed to understand its phylogenetic distribution and variants against global and national sequences. Further variants were analyzed from diseased and recovered patients from Gujarat and the world to understand its role in pathogenesis. Among the missense mutations present in the Gujarat SARS-CoV-2 genomes, C28854T (Ser194Leu) had an allele frequency of 47.62 and 7.25% in deceased patients from the Gujarat and global datasets, respectively. In contrast, the allele frequency of 35.16 and 3.20% was observed in recovered patients from the Gujarat and global datasets, respectively. It is a deleterious mutation present in the nucleocapsid (N) gene and is significantly associated with mortality in Gujarat patients with a p-value of 0.067 and in the global dataset with a p-value of 0.000924. The other deleterious variant identified in deceased patients from Gujarat (p-value of 0.355) and the world (p-value of 2.43E-06) is G25563T, which is located in Orf3a and plays a potential role in viral pathogenesis. SARS-CoV-2 genomes from Gujarat are forming distinct clusters under the GH clade of GISAID. This study will shed light on the viral haplotype in SARS-CoV-2 samples from Gujarat, India.

SARS-CoV-2 Genome from the Khyber Pakhtunkhwa Province of Pakistan

Among viral outbreaks, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is one of the deadliest ones, and it has triggered the global COVID-19 pandemic. In Pakistan, until 5th September 2020, a total of 6342 deaths have been reported, of which 1255 were from the Khyber Pakhtunkhwa (KPK) province. To understand the disease progression and control and also to produce vaccines and therapeutic efforts, whole genome sequence analysis is important. In the current investigation, we sequenced a single sample of SARS-CoV-2 genomes (accession no. MT879619) from a male suspect from Peshawar, the KPK capital city, during the first wave of infection. The local SARS-CoV-2 strain shows some unique characteristics compared to neighboring Iranian and Chinese isolates in phylogenetic tree and mutations. The circulating strains of SARS-CoV-2 represent an intermediate evolution from China and Iran. Furthermore, eight complete whole genome sequences, including the current Pakistani isolates which have been submitted to Global Initiative on Sharing All Influenza Data (GSAID), were also investigated for specific mutations and characters. Some novel mutations [NSP2 (D268del), NSP5 (N228K), and NS3 (F105S)] and specific characters have been detected in the coding regions, which may affect viral transmission, epidemiology, and disease severity. The computational modeling revealed that a majority of these mutations may have a stabilizing effect on the viral protein structure. In conclusion, the genome sequencing of local strains is important for better understanding the pathogenicity, immunogenicity, and epidemiology of causative agents.

Sequence analysis of Indian SARS-CoV-2 isolates shows a stronger interaction of mutant receptor-binding domain with ACE2

OBJECTIVE

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has affected the whole world, including Odisha, a state in eastern India. Many people have migrated to the state from different countries as well as other states during this SARS-CoV-2 pandemic. The aim of this study was to analyse the receptor-binding domain (RBD) sequence of the spike protein from isolates collected from throat swab samples of COVID-19-positive patients and further to assess the RBD affinity for angiotensin-converting enzyme 2 (ACE2) of different species, including humans.

METHODS

Whole-genome sequencing for 35 clinical SARS-CoV-2 isolates from COVID-19-positive patients was performed by ARTIC amplicon-based sequencing. Sequence analysis and phylogenetic analysis were performed for the spike region and the RBD region of all isolates. The interaction between the RBD and ACE2 of five different species was also analysed.

RESULTS

The spike region of 32 isolates showed one or multiple alterations in nucleotide bases in comparison with the Wuhan reference strain. One of the identified mutations, at position 1204 (Ref A, RMRC 22 C), in the RBD coding region of the spike protein showed stronger binding affinity for human ACE2. Furthermore, RBDs of all the Indian isolates showed binding affinity for ACE2 of different species.

CONCLUSION

As mutant RBD showed stronger interaction with human ACE2, it could potentially result in higher infectivity. The binding affinity of the RBDs for ACE2 of all five species studied suggests that the virus can infect a wide variety of animals, which could also act as natural reservoir for SARS-CoV-2.

Rare mutations in the accessory proteins ORF6, ORF7b, and ORF10 of the SARS-CoV-2 genomes

A total number of 3080 SARS-CoV-2 genomes from all continents are considered from the NCBI database. Every accessory protein ORF6, ORF7b, and ORF10 of SARS-CoV-2 possess a single missense mutation in less than 1.5% of the 3080 genomes. It has now been observed that different non-synonymous mutations occurred in these three accessory proteins. Most of these rare mutations are changing the amino acids such as hydrophilic to hydrophobic, acidic or basic to hydrophobic, and vice versa etc. So these highly conserved proteins might play an essential role in virus pathogenicity. This study opens a question whether it carries some messages about the virus rapid replications, and virulence.

Mutational analysis and assessment of its impact on proteins of SARS-CoV-2 genomes from India

Mutational status of SARS-CoV-2 genomes from India along with their impact on proteins was ascertained through multiple tools including MEGA, Genome Detective, SIFT, PROVEAN and ws-SNPs&GO. Excluding gaps and ambiguous sequences, 493 variable sites (152 parsimony informative and 341 singleton) were observed. NSP3 had the highest incidence of 101 sites followed by S protein (74), NSP12b (43) and ORF3a (31). Average mutations per sample for males and females was 2.56 and 2.88 respectively. Non-uniform geographical distribution of mutations suggests that sequences in some regions are mutating faster than others. There were 281 mutations (198 Neutral and 83 Disease) affecting amino acid sequence. NSP13 has a maximum of 14 Disease variants followed by S protein and ORF3a with 13 each. Disease mutations in genomes from asymptomatic people was mere 11% but those from deceased patients was at 38% indicating contribution of these mutations to the pathophysiology of the SARS-CoV-2.

SARS-CoV-2 genomics: An Indian perspective on sequencing viral variants

Since its emergence as a pneumonia-like outbreak in the Chinese city of Wuhan in late 2019, the novel coronavirus disease COVID-19 has spread widely to become a global pandemic. The first case of COVID-19 in India was reported on 30 January 2020 and since then it has affected more than ten million people and resulted in around 150,000 deaths in the country. Over time, the viral genome has accumulated mutations as it passes through its human hosts, a common evolutionary mechanism found in all microorganisms. This has implications for disease surveillance and management, vaccines and therapeutics, and the emergence of reinfections. Sequencing the viral genome can help monitor these changes and provides an extraordinary opportunity to understand the genetic epidemiology and evolution of the virus as well as tracking its spread in a population. Here we review the past year in the context of the phylogenetic analysis of variants isolated over the course of the pandemic in India and highlight the importance of continued sequencing-based surveillance in the country.

Comprehensive analysis of genomic diversity of SARS-CoV-2 in different geographic regions of India: an endeavour to classify Indian SARS-CoV-2 strains on the basis of co-existing mutations

Accumulation of mutations within the genome is the primary driving force in viral evolution within an endemic setting. This inherent feature often leads to altered virulence, infectivity and transmissibility, and antigenic shifts to escape host immunity, which might compromise the efficacy of vaccines and antiviral drugs. Therefore, we carried out a genome-wide analysis of circulating SARS-CoV-2 strains to detect the emergence of novel co-existing mutations and trace their geographical distribution within India. Comprehensive analysis of whole genome sequences of 837 Indian SARS-CoV-2 strains revealed the occurrence of 33 different mutations, 18 of which were unique to India. Novel mutations were observed in the S glycoprotein (6/33), NSP3 (5/33), RdRp/NSP12 (4/33), NSP2 (2/33), and N (1/33). Non-synonymous mutations were found to be 3.07 times more prevalent than synonymous mutations. We classified the Indian isolates into 22 groups based on their co-existing mutations. Phylogenetic analysis revealed that the representative strains of each group were divided into various sub-clades within their respective clades, based on the presence of unique co-existing mutations. The A2a clade was found to be dominant in India (71.34%), followed by A3 (23.29%) and B (5.36%), but a heterogeneous distribution was observed among various geographical regions. The A2a clade was highly predominant in East India, Western India, and Central India, whereas the A2a and A3 clades were nearly equal in prevalence in South and North India. This study highlights the divergent evolution of SARS-CoV-2 strains and co-circulation of multiple clades in India. Monitoring of the emerging mutations will pave the way for vaccine formulation and the design of antiviral drugs.

Genetics and genomics of SARS-CoV-2: A review of the literature with the special focus on genetic diversity and SARS-CoV-2 genome detection

The outbreak of 2019-novel coronavirus disease (COVID-19), caused by SARS-CoV-2, started in late 2019; in a short time, it has spread rapidly all over the world. Although some possible antiviral and anti-inflammatory medications are available, thousands of people are dying daily. Well-understanding of the SARS-CoV-2 genome is not only essential for the development of new treatments/vaccines, but it also can be used for improving the sensitivity and specificity of current approaches for virus detection. Accordingly, we reviewed the most critical findings related to the genetics of the SARS-CoV-2, with a specific focus on genetic diversity and reported mutations, molecular-based diagnosis assays, using interfering RNA technology for the treatment of patients, and genetic-related vaccination strategies. Additionally, considering the unanswered questions or uncertainties in these regards, different topics were discussed.

SARS-CoV2 envelope protein: non-synonymous mutations and its consequences

In the NCBI database, as on June 6, 2020, total number of available complete genome sequences of SARS-CoV2 across the world is 3617. The envelope (E) protein of SARS-CoV2 possesses several non-synonymous mutations over the transmembrane and C-terminus domains in 15 (0.414%) genomes among 3617 SARS-CoV2 genomes, analyzed. More precisely, 10(0.386%) out of 2588 genomes from the USA, 3(0.806%) from Asia, 1 (0.348%) from Europe and 1 (0.274%) from Oceania contained the missense mutations over the E-protein of SARS-CoV2 genomes. The C-terminus motif DLLV has been to DFLV and YLLV in the proteins from QJR88103 (Australia: Victoria) and QKI36831 (China: Guangzhou) respectively, which might affect the binding of this motif with the host protein PALS1.

Molecular conservation and differential mutation on ORF3a gene in Indian SARS-CoV2 genomes

A global emergency due to the COVID-19 pandemic demands various studies related to genes and genomes of the SARS-CoV2. Among other important proteins, the role of accessory proteins are of immense importance in replication, regulation of infections of the coronavirus in the hosts. The largest accessory protein in the SARS-CoV2 genome is ORF3a which modulates the host response to the virus infection and consequently it plays an important role in pathogenesis. In this study, an attempt is made to decipher the conservation of nucleotides, dimers, codons and amino acids in the ORF3a genes across thirty-two genomes of Indian patients. ORF3a gene possesses single and double point mutations in Indian SARS-CoV2 genomes suggesting the change of SARS-CoV2's virulence property in Indian patients. We find that the parental origin of the ORF3a gene over the genomes of SARS-CoV2 and Pangolin-CoV is same from the phylogenetic analysis based on conservation of nucleotides and so on. This study highlights the accumulation of mutation on ORF3a in Indian SARS-CoV2 genomes which may provide the designing therapeutic approach against SARS-CoV2.