Phylogenetic analysis of SARS-CoV-2 genomes in Turkey

Mutational landscape of SARS-CoV-2 genome in Turkey and impact of mutations on spike protein structure

The Coronavirus Disease 2019 (COVID-19) was declared a pandemic in March 2020 by the World Health Organization (WHO). As of May 25th, 2021 there were 2.059.941 SARS-COV2 genome sequences that have been submitted to the GISAID database, with numerous variations. Here, we aim to analyze the SARS-CoV-2 genome data submitted to the GISAID database from Turkey and to determine the variant and clade distributions by the end of May 2021, in accordance with their appearance timeline. We compared these findings to USA, Europe, and Asia data as well. We have also evaluated the effects of spike protein variations, detected in a group of genome sequences of 13 patients who applied to our clinic, by using 3D modeling algorithms. For this purpose, we analyzed 4607 SARS-CoV-2 genome sequences submitted by different lab centers from Turkey to the GISAID database between March 2020 and May 2021. Described mutations were also introduced in silico to the spike protein structure to analyze their isolated impacts on the protein structure. The most abundant clade was GR followed by G, GH, and GRY and we did not detect any V clade. The most common variant was B.1, followed by B.1.1, and the UK variant, B.1.1.7. Our results clearly show a concordance between the variant distributions, the number of cases, and the timelines of different variant accumulations in Turkey. The 3D simulations indicate an increase in the surface hydrophilicity of the reference spike protein and the detected mutations. There was less surface hydrophilicity increase in the Asp614Gly mutation, which exhibits a more compact conformation around the ACE-2 receptor binding domain region, rendering the structure in a “down” conformation. Our genomic findings can help to model vaccination programs and protein modeling may lead to different approaches for COVID-19 treatment strategies.

SARS-CoV-2 surveillance in Italy through phylogenomic inferences based on Hamming distances derived from pan-SNPs, -MNPs and -InDels

BACKGROUND

Faced with the ongoing global pandemic of coronavirus disease, the 'National Reference Centre for Whole Genome Sequencing of microbial pathogens: database and bioinformatic analysis' (GENPAT) formally established at the 'Istituto Zooprofilattico Sperimentale dell'Abruzzo e del Molise' (IZSAM) in Teramo (Italy) is in charge of the SARS-CoV-2 surveillance at the genomic scale. In a context of SARS-CoV-2 surveillance requiring correct and fast assessment of epidemiological clusters from substantial amount of samples, the present study proposes an analytical workflow for identifying accurately the PANGO lineages of SARS-CoV-2 samples and building of discriminant minimum spanning trees (MST) bypassing the usual time consuming phylogenomic inferences based on multiple sequence alignment (MSA) and substitution model.

RESULTS

GENPAT constituted two collections of SARS-CoV-2 samples. The first collection consisted of SARS-CoV-2 positive swabs collected by IZSAM from the Abruzzo region (Italy), then sequenced by next generation sequencing (NGS) and analyzed in GENPAT (n = 1592), while the second collection included samples from several Italian provinces and retrieved from the reference Global Initiative on Sharing All Influenza Data (GISAID) (n = 17,201). The main results of the present work showed that (i) GENPAT and GISAID detected the same PANGO lineages, (ii) the PANGO lineages B.1.177 (i.e. historical in Italy) and B.1.1.7 (i.e. 'UK variant') are major concerns today in several Italian provinces, and the new MST-based method (iii) clusters most of the PANGO lineages together, (iv) with a higher dicriminatory power than PANGO lineages, (v) and faster that the usual phylogenomic methods based on MSA and substitution model.

CONCLUSIONS

The genome sequencing efforts of Italian provinces, combined with a structured national system of NGS data management, provided support for surveillance SARS-CoV-2 in Italy. We propose to build phylogenomic trees of SARS-CoV-2 variants through an accurate, discriminant and fast MST-based method avoiding the typical time consuming steps related to MSA and substitution model-based phylogenomic inference.

SARS-CoV-2 Genome Analysis of Pediatric Patients in Konya Region, Turkey

Objective Genome sequencing is useful for following the change in mutation and variants in viral agent during pandemics. In this study, we performed next-generation sequencing of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) complete genomes on pediatric patients.

Methods Six pediatric patients aged 0 to 18 years who were positive for SARS-CoV-2 by reverse transcription polymerase chain reaction were included in this study. SARS-CoV-2 genome sequencing was performed using Oxford Nanopore Technologies MinION, following the ARTIC Network protocols. Sequencing data were obtained using the FASTQ program and quality assessment was evaluated. The sequence information of all samples was uploaded to the Global Initiative on Sharing All Influenza Data (GISAID) database. Genome, variant, clade, and phylogenetic tree analyses were performed with bioinformatic analysis.

Results Two of these six samples were at 20A, two were at 20B, and two were at 19A in the nextstrain clade. According to Pango lineages, B.1.36, B.1.218, B.1, and B.1.260 lineages were detected. A total of 84 mutations were observed in all samples. None of the variants were classified as variants of concern (VOC) nor variants of interest (VOI) according to the Pango database.

Conclusion This study is the first comprehensively sequence analysis registered in the GISAID database reported from the Konya region in Turkey. Similar studies will be informative to track changes in the virus genome, obtain epidemiological data, guide studies on diagnosis and treatment, and evaluate vaccine efficacy.

Genomic characterization of SARS-CoV-2 isolates from patients in Turkey reveals the presence of novel mutations in spike and nsp12 proteins

Novel mutations have been emerging in the genome of severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2); consequently, the evolving of more virulent and treatment resistance strains have the potential to increase transmissibility and mortality rates. The characterization of full-length SARS-CoV-2 genomes is critical for understanding the origin and transmission pathways of the virus, as well as identifying mutations that affect the transmissibility and pathogenicity of the virus. We present an analysis of the mutation pattern and clade distribution of full-length SARS-CoV-2 genome sequences obtained from specimens tested at Gazi University Medical Virology Laboratory. Viral RNA was extracted from nasopharyngeal specimens. Next-generation sequencing libraries were prepared and sequenced on Illumina iSeq 100 platform. Raw sequencing data were processed to obtain full-length genome sequences and variant calling was performed to analyze amino acid changes. Clade distribution was determined to understand the phylogenetic background in relation to global data. A total of 293 distinct mutations were identified, of which 152 missense, 124 synonymous, 12 noncoding, and 5 deletions. The most frequent mutations were P323L (nsp12), D614G (ORF2/S), and 2421C>T (5'-untranslated region) found simultaneously in all sequences. Novel mutations were found in nsp12 (V111A, H133R, Y453C, M626K) and ORF2/S (R995G, V1068L). Nine different Pangolin lineages were detected. The most frequently assigned lineage was B.1.1 (17 sequences), followed by B.1 (7 sequences) and B.1.1.36 (3 sequences). Sequence information is essential for revealing genomic diversity. Mutations might have significant functional implications and analysis of these mutations provides valuable information for therapeutic and vaccine development studies. Our findings point to the introduction of the virus into Turkey through various sources and the subsequent spread of several key variants.

Unlocking capacities of viral genomics for the COVID-19 pandemic response

More than any other infectious disease epidemic, the COVID-19 pandemic has been characterized by the generation of large volumes of viral genomic data at an incredible pace due to recent advances in high-throughput sequencing technologies, the rapid global spread of SARS-CoV-2, and its persistent threat to public health. However, distinguishing the most epidemiologically relevant information encoded in these vast amounts of data requires substantial effort across the research and public health communities. Studies of SARS-CoV-2 genomes have been critical in tracking the spread of variants and understanding its epidemic dynamics, and may prove crucial for controlling future epidemics and alleviating significant public health burdens. Together, genomic data and bioinformatics methods enable broad-scale investigations of the spread of SARS-CoV-2 at the local, national, and global scales and allow researchers the ability to efficiently track the emergence of novel variants, reconstruct epidemic dynamics, and provide important insights into drug and vaccine development and disease control. Here, we discuss the tremendous opportunities that genomics offers to unlock the effective use of SARS-CoV-2 genomic data for efficient public health surveillance and guiding timely responses to COVID-19.

The Mutation Profile of SARS-CoV-2 Is Primarily Shaped by the Host Antiviral Defense

Understanding SARS-CoV-2 evolution is a fundamental effort in coping with the COVID-19 pandemic. The virus genomes have been broadly evolving due to the high number of infected hosts world-wide. Mutagenesis and selection are two inter-dependent mechanisms of virus diversification. However, which mechanisms contribute to the mutation profiles of SARS-CoV-2 remain under-explored. Here, we delineate the contribution of mutagenesis and selection to the genome diversity of SARS-CoV-2 isolates. We generated a comprehensive phylogenetic tree with representative genomes. Instead of counting mutations relative to the reference genome, we identified each mutation event at the nodes of the phylogenetic tree. With this approach, we obtained the mutation events that are independent of each other and generated the mutation profile of SARS-CoV-2 genomes. The results suggest that the heterogeneous mutation patterns are mainly reflections of host (i) antiviral mechanisms that are achieved through APOBEC, ADAR, and ZAP proteins, and (ii) probable adaptation against reactive oxygen species.

Structure of SARS-CoV-2 Spike Glycoprotein for Therapeutic and Preventive Target

The global crisis caused by the coronavirus disease 2019 (COVID-19) led to the most significant economic loss and human deaths after World War II. The pathogen causing this disease is a novel virus called the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). As of December 2020, there have been 80.2 million confirmed patients, and the mortality rate is known as 2.16% globally. A strategy to protect a host from SARS-CoV-2 is by suppressing intracellular viral replication or preventing viral entry. We focused on the spike glycoprotein that is responsible for the entry of SARS-CoV-2 into the host cell. Recently, the US Food and Drug Administration/EU Medicines Agency authorized a vaccine and antibody to treat COVID-19 patients by emergency use approval in the absence of long-term clinical trials. Both commercial and academic efforts to develop preventive and therapeutic agents continue all over the world. In this review, we present a perspective on current reports about the spike glycoprotein of SARS-CoV-2 as a therapeutic target.

Stability of SARS-CoV-2 phylogenies

The SARS-CoV-2 pandemic has led to unprecedented, nearly real-time genetic tracing due to the rapid community sequencing response. Researchers immediately leveraged these data to infer the evolutionary relationships among viral samples and to study key biological questions, including whether host viral genome editing and recombination are features of SARS-CoV-2 evolution. This global sequencing effort is inherently decentralized and must rely on data collected by many labs using a wide variety of molecular and bioinformatic techniques. There is thus a strong possibility that systematic errors associated with lab-or protocol-specific practices affect some sequences in the repositories. We find that some recurrent mutations in reported SARS-CoV-2 genome sequences have been observed predominantly or exclusively by single labs, co-localize with commonly used primer binding sites and are more likely to affect the protein-coding sequences than other similarly recurrent mutations. We show that their inclusion can affect phylogenetic inference on scales relevant to local lineage tracing, and make it appear as though there has been an excess of recurrent mutation or recombination among viral lineages. We suggest how samples can be screened and problematic variants removed, and we plan to regularly inform the scientific community with our updated results as more SARS-CoV-2 genome sequences are shared (https://virological.org/t/issues-with-sars-cov-2-sequencing-data/473 and https://virological.org/t/masking-strategies-for-sars-cov-2-alignments/480). We also develop tools for comparing and visualizing differences among very large phylogenies and we show that consistent clade- and tree-based comparisons can be made between phylogenies produced by different groups. These will facilitate evolutionary inferences and comparisons among phylogenies produced for a wide array of purposes. Building on the SARS-CoV-2 Genome Browser at UCSC, we present a toolkit to compare, analyze and combine SARS-CoV-2 phylogenies, find and remove potential sequencing errors and establish a widely shared, stable clade structure for a more accurate scientific inference and discourse.

Mutations Observed in the SARS-CoV-2 Spike Glycoprotein and Their Effects in the Interaction of Virus with ACE-2 Receptor

Coronaviruses (CoVs) classified in the Coronaviridae family infect a very large spectrum of vertebrate group. Seven CoVs that cause human disease consist of Alpha-CoVs, which are HCoV-229E, and NL63 and beta-CoVs, which are MERS-CoV, SARS-CoV, HCoV-OC43, HCoV-HKU1, and SARS-CoV-2. SARS-CoV-2 is an enveloped, positive-polarity, single-stranded RNA virus responsible for a new Coronavirus disease 2019 (COVID-19). The mutagenic ability of the SARS-CoV-2 directs its evolution and genome variability, thus allowing viruses to escape from host immunity and develop drug resistance. Tracing viral mutations is also important for the development of new vaccines, antiviral drugs, and diagnostic systems. During replication in the host cell, genomic mutations occur in the virus and these mutations are transferred to new generations. For this reason, systematic monitoring of mutations in the SARS-CoV-2 genome allows observation of the national and international molecular epidemiology of the virus. SARS-CoV-2 spike (S) glycoprotein is vital in the binding of the virus to the host cell receptor that is angiotensin converting-enzyme 2 (ACE2), membrane fusion, vaccine studies and immune response to the virus. Therefore, mutations in the gene encoding the S glycoprotein and especially the possible variations in the receptor binding domain (RBD) in S gene are important issues to be emphasized. In this article, information about the mutations observed in the SARS-CoV-2 S glycoprotein and their possible effects are presented.

Characterization of local SARS-CoV-2 isolates and pathogenicity in IFNAR-/- mice

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) recently a global pandemic with unprecedented public health, economic and social impact. The development of effective mitigation strategies, therapeutics and vaccines relies on detailed genomic and biological characterization of the regional viruses. This study was carried out to isolate SARS-CoV-2 viruses circulating in Anatolia, and to investigate virus propagation in frequently-used cells and experimental animals. We obtained two SARS-CoV-2 viruses from nasopharngeal swabs of confirmed cases in Vero E6 cells, visualized the virions using atomic force and scanning electron microscopy and determined size distribution of the particles. Viral cytopathic effects on Vero E6 cells were initially observed at 72 h post-inoculation and reached 90% of the cells on the 5th day. The isolates displayed with similar infectivity titers, time course and infectious progeny yields. Genome sequencing revealed the viruses to be well-conserved, with less than 1% diversity compared to the prototype virus. The analysis of the viral genomes, along with the available 62 complete genomes from Anatolia, showed limited diversity (up to 0.2% on deduced amino acids) and no evidence of recombination. The most prominent sequence variation was observed on the spike protein, resulting in the substitution D614G, with a prevalence of 56.2%. The isolates produced non-fatal infection in the transgenic type I interferon knockout (IFNAR^-/-) mice, with varying neutralizing antibody titers. Hyperemia, regional consolidation and subpleural air accumulation was observed on necropsy, with similar histopathological and immunohistochemistry findings in the lungs, heart, stomach, intestines, liver, spleen and kidneys. Peak viral loads were detected in the lungs, with virus RNA present in the kidneys, jejunum, liver, spleen and heart. In conclusion, we characterized two local isolates, investigated in vitro growth dynamics in Vero E6 cells and identified IFNAR-/- mice as a potential animal model for SARS-CoV-2 experiments.

The New Coronavirus (SARS-CoV-2): A Comprehensive Review on Immunity and the Application of Bioinformatics and Molecular Modeling to the Discovery of Potential Anti-SARS-CoV-2 Agents

On March 11, 2020, the World Health Organization (WHO) officially declared the outbreak caused by the new coronavirus (SARS-CoV-2) a pandemic. The rapid spread of the disease surprised the scientific and medical community. Based on the latest reports, news, and scientific articles published, there is no doubt that the coronavirus has overloaded health systems globally. Practical actions against the recent emergence and rapid expansion of the SARS-CoV-2 require the development and use of tools for discovering new molecular anti-SARS-CoV-2 targets. Thus, this review presents bioinformatics and molecular modeling strategies that aim to assist in the discovery of potential anti-SARS-CoV-2 agents. Besides, we reviewed the relationship between SARS-CoV-2 and innate immunity, since understanding the structures involved in this infection can contribute to the development of new therapeutic targets. Bioinformatics is a technology that assists researchers in coping with diseases by investigating genetic sequencing and seeking structural models of potential molecular targets present in SARS-CoV2. The details provided in this review provide future points of consideration in the field of virology and medical sciences that will contribute to clarifying potential therapeutic targets for anti-SARS-CoV-2 and for understanding the molecular mechanisms responsible for the pathogenesis and virulence of SARS-CoV-2.

Investigation of compatibility of severe acute respiratory syndrome coronavirus 2 reverse transcriptase-PCR kits containing different gene targets during coronavirus disease 2019 pandemic

Aim:

In the diagnosis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), reverse transcriptase-PCR (RT-PCR) technique is often used. We evaluated the compatibility of SARS-CoV-2 RT-PCR kits containing different gene targets during the pandemic.

Materials & methods:

Samples were tested by Bio-Speddy^® (RdRp gene) and Diagnovital^® (RdRp + E genes). The correlation between two assays were determined by Deming regression analysis and chi-square analyses.

Results:

Diagnovital PCR kit showed amplification in a narrow Ct range and conveniently sharper exponential amplification curves than Bio-Speedy PCR kit. While the correlation between the findings of the two kits was apparent even with single gene target, this correlation increased when a secondary biomarker was added to the correlation calculations.

Conclusion:

We have observed high correlation between different PCR kits, however, using different PCR kits during the pandemic may provide a more accurate diagnosis of SARS-CoV-2, since despite correlation there are a number of patients showing contradicting diagnosis.