GISAID Database Research Articles

SARS-CoV-2 lineage B.6 was the major contributor to early pandemic transmission in Malaysia.

Malaysia had 10,219 confirmed cases of COVID-19 as of September 20, 2020. About 33% were associated with a Tablighi Jamaat religious mass gathering held in Kuala Lumpur between February 27 and March 3, 2020, which drove community transmission during Malaysia's second wave. We analysed genome sequences of SARS-CoV-2 from Malaysia to better understand the molecular epidemiology and spread. We obtained 58 SARS-CoV-2 whole genome sequences from patients in Kuala Lumpur and performed phylogenetic analyses on these and a further 57 Malaysian sequences available in the GISAID database. Nine different SARS-CoV-2 lineages (A, B, B.1, B.1.1, B.1.1.1, B.1.36, B.2, B.3 and B.6) were detected in Malaysia. The B.6 lineage was first reported a week after the Tablighi mass gathering and became predominant (65.2%) despite being relatively rare (1.4%) globally. Direct epidemiological links between lineage B.6 viruses and the mass gathering were identified. Increases in reported total cases, Tablighi-associated cases, and community-acquired B.6 lineage strains were temporally linked. Non-B.6 lineages were mainly travel-associated and showed limited onward transmission. There were also temporally correlated increases in B.6 sequences in other Southeast Asian countries, India and Australia, linked to participants returning from this event. Over 95% of global B.6 sequences originated from Asia Pacific. We also report a nsp3-C6310A substitution found in 47.3% of global B.6 sequences which was associated with reduced sensitivity using a commercial diagnostic real-time PCR assay. Lineage B.6 became the predominant cause of community transmission in Malaysia after likely introduction during a religious mass gathering. This event also contributed to spikes of lineage B.6 in other countries in the Asia-Pacific. Mass gatherings can be significant causes of local and global spread of COVID-19. Shared genomic surveillance can be used to identify SARS-CoV-2 transmission chains to aid prevention and control, and to monitor diagnostic molecular assays. Clinical Trial Registration: COVID-19 paper.

SARS-CoV-2 Among Military and Civilian Patients, Metro Manila, Philippines.

ABSTRACTIntroductionSevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) belonging to the family Coronaviridae and genus Betacoronavirus is the causative agent of COVID-19 disease and was first identified in Wuhan, China. SARS-CoV-2 spread globally with >28 million cases and 911,000 deaths recorded worldwide as of September 12, 2020. The Philippines reported the first case of community transmission on March 5, 2020, and despite the government imposing one of the longest and strictest lockdowns in Southeast Asia, the number of confirmed COVID-19 cases still surged with >250,000 cases and 4,000 deaths reported as of September 12, 2020. It is important to estimate the burden and impact of SARS-CoV-2 on the military population since this can affect the military readiness.Materials and MethodsNasopharyngeal and oropharyngeal swabs were collected and SARS-CoV-2 real-time RT-PCR testing was performed on the samples. Phylogenetic analysis was performed using sequences from 23 SARS-CoV-2-positive specimens from this study and sequences retrieved from GenBank and GISAID databases.ResultsFrom April 14 to August 15, 2020, a total of 12,432 specimens were tested with 763 (6%) unique individuals testing positive for SARS-CoV-2 by rRT-PCR. In the military population, majority of the patients who were tested (80%) and those who tested positive for SARS-CoV-2 (86%) were male. Military and civilian status was available for 7,672 patients with 515/5,042 (10%) positive among military patients and 248/2,630 (9%) positive among civilian patients. Both military and civilian populations had the highest case counts of SARS-CoV-2-positive cases in the 21- to 30- and 31- to 40-year-old age groups, while the 71- to 80-year-old age group had the highest proportion (18%) of SARS-CoV-2-positive cases. Sequencing analysis showed 19 different variants in the 23 genomes. Twenty of the 23 genomes were classified under clade GR/B1.1, 2 genomes were classified under clade GR/B1.1.28, and 1 genome was classified under Clade O/B.6. Twenty-two of the 23 sequences collected after June 25, 2020, contained the D614G mutation.ConclusionWe describe here the results of SARS-CoV-2 testing for military and civilian patients and personnel. The 21- to 30- and 31- to 40-year-old age groups had the highest case counts of SARS-CoV-2-positive cases. Sequencing results showed the presence of the D614G mutation in the spike protein in a majority of specimens collected from the end of June to July 2020.

SARS-CoV-2 spread across the Colombian-Venezuelan border.

IntroductionVenezuela and Colombia both adopted measures of containment early in response to the COVID-19 pandemic. However, Venezuela's ongoing humanitarian crisis has decimated its health care system, and forced millions of Venezuelans to flee through its porous border with Colombia. The extensive shared border, and illegal cross-border transit through improvised trails between the two countries are major challenges for public health authorities. We report the first SARS-CoV-2 genomes from Venezuela, and present a snapshot of the SARS-CoV-2 epidemiologic landscape in the Colombian-Venezuelan border region.MethodsWe sequenced and assembled viral genomes from total RNA extracted from nasopharyngeal (NP) clinical specimens using a custom reference-based analysis pipeline. Three assemblies obtained were subjected to typing using the Phylogenetic Assignment of Named Global Outbreak LINeages ‘Pangolin’ tool. A total of 376 publicly available SARS-CoV-2 genomes from South America were obtained from the GISAID database to perform comparative genomic analyses. Additionally, the Wuhan-1 strain was used as reference.ResultsWe found that two of the SARS-CoV-2 genomes from Venezuela belonged to the B1 lineage, and the third to the B.1.13 lineage. We observed a point mutation in the Spike protein gene (D614G substitution), previously reported to be associated with increased infectivity, in all three Venezuelan genomes. Additionally, three mutations (R203K/G204R substitution) were present in the nucleocapsid (N) gene of one Venezuelan genome.ConclusionsGenomic sequencing demonstrates similarity between SARS-CoV-2 lineages from Venezuela and viruses collected from patients in bordering areas in Colombia and from Brazil, consistent with cross-border transit despite administrative measures including lockdowns. The presence of mutations associated with increased infectivity in the 3 Venezuelan genomes we report and Colombian SARS-CoV-2 genomes from neighboring borders areas may pose additional challenges for control of SARS-CoV-2 spread in the complex epidemiological landscape in Latin American countries. Public health authorities should carefully follow the progress of the pandemic and its impact on displaced populations within the region.

Analysis of SARS‐CoV‐2 mutations in Mexico, Belize, and isolated regions of Guatemala and its implication in the diagnosis

The genomic sequences of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) worldwide are publicly available and are derived from studies due to the increase in the number of cases. The importance of study of mutations is related to the possible virulence and diagnosis of SARS‐CoV‐2. To identify circulating mutations present in SARS‐CoV‐2 genomic sequences in Mexico, Belize, and Guatemala to find out if the same strain spread to the south, and analyze the specificity of the primers used for diagnosis in these samples. Twenty three complete SARS‐CoV‐2 genomic sequences, available in the GISAID database from May 8 to September 11, 2020 were analyzed and aligned versus the genomic sequence reported in Wuhan, China (NC_045512.2), using Clustal Omega. Open reading frames were translated using the ExPASy Translate Tool and UCSF Chimera (v.1.12) for amino acid substitutions analysis. Finally, the sequences were aligned versus primers used in the diagnosis of COVID‐19. One hundred and eighty seven distinct variants were identified, of which 102 are missense, 66 synonymous and 19 noncoding. P4715L and P5828L substitutions in replicase polyprotein were found, as well as D614G in spike protein and L84S in ORF8 in Mexico, Belize, and Guatemala. The primers design by CDC of United States showed a positive E value. The genomic sequences of SARS‐CoV‐2 in Mexico, Belize, and Guatemala present similar mutations related to a virulent strain of greater infectivity, which could mean a greater capacity for inclusion in the host genome and be related to an increased spread of the virus in these countries, furthermore, its diagnosis would be affected.

Using genomics to improve preparedness and response of future epidemics or pandemics in Africa

Emerging infectious diseases, such as COVID-191 and Ebola virus disease,2 pose a major threat to public health. Despite the presence of Ebola virus disease in central and west Africa for several years,2 many African countries are yet to fully utilise genomics in public health. This situation is highlighted by the current COVID-19 pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2).1 By Oct 5, 2020, 9 months into the COVID-19 pandemic, 2980 SARS-CoV-2 genomes sampled from patients with COVID-19, from 21 of 54 African countries, had been deposited in the GISAID database.

Genomic Diversity and Hotspot Mutations in 30,983 SARS-CoV-2 Genomes: Moving Toward a Universal Vaccine for the "Confined Virus"?

The COVID-19 pandemic has been ongoing since its onset in late November 2019 in Wuhan, China. Understanding and monitoring the genetic evolution of the virus, its geographical characteristics, and its stability are particularly important for controlling the spread of the disease and especially for the development of a universal vaccine covering all circulating strains. From this perspective, we analyzed 30,983 complete SARS-CoV-2 genomes from 79 countries located in the six continents and collected from 24 December 2019, to 13 May 2020, according to the GISAID database. Our analysis revealed the presence of 3206 variant sites, with a uniform distribution of mutation types in different geographic areas. Remarkably, a low frequency of recurrent mutations has been observed; only 169 mutations (5.27%) had a prevalence greater than 1% of genomes. Nevertheless, fourteen non-synonymous hotspot mutations (>10%) have been identified at different locations along the viral genome; eight in ORF1ab polyprotein (in nsp2, nsp3, transmembrane domain, RdRp, helicase, exonuclease, and endoribonuclease), three in nucleocapsid protein, and one in each of three proteins: Spike, ORF3a, and ORF8. Moreover, 36 non-synonymous mutations were identified in the receptor-binding domain (RBD) of the spike protein with a low prevalence (<1%) across all genomes, of which only four could potentially enhance the binding of the SARS-CoV-2 spike protein to the human ACE2 receptor. These results along with intra-genomic divergence of SARS-CoV-2 could indicate that unlike the influenza virus or HIV viruses, SARS-CoV-2 has a low mutation rate which makes the development of an effective global vaccine very likely.

Geographical reconstruction of the SARS‐CoV‐2 outbreak in Lombardy (Italy) during the early phase

The first identification of autochthonous transmission of SARS‐CoV‐2 in Italy was documented by the Laboratory of Clinical Microbiology, Virology and Bioemergencies of L. Sacco Hospital (Milano, Italy) on 20th February 2020 in a 38 years old male patient, who was found positive for pneumonia at the Codogno Hospital. Thereafter Lombardy has reported the highest prevalence of COVID‐19 cases in the country, especially in Milano, Brescia and Bergamo provinces. The aim of this study was to assess the potential presence of different viral clusters belonging to the six main provinces involved in Lombardy COVID‐19 cases in order to highlight peculiar province‐dependent viral characteristics. A phylogenetic analysis was conducted on 20 full length genomes obtained from patients addressing to several Lombard hospitals from February 20th to April 4th, 2020, aligned with 41 Italian viral genome assemblies available on GISAID database as of 30th March, 2020: two main monophyletic clades, containing 8 and 53 isolates, respectively, were identified. Noteworthy, Bergamo isolates mapped inside the small clade harbouring M gene D3G mutation. The molecular clock analysis estimated a cluster divergence approximately one month before the first patient identification, supporting the hypothesis that different SARS‐CoV‐2 strains had spread worldwide at different times, but their presence became evident only in late February along with Italian epidemic emergence. Therefore, this epidemiological reconstruction suggests that virus initial circulation in Lombardy was ascribable to multiple introduction. The phylogenetic reconstruction robustness, however, will be improved when more genomic sequences are available, in order to guarantee a complete epidemiological surveillance.

Genetic grouping of SARS-CoV-2 coronavirus sequences using informative subtype markers for pandemic spread visualization.

We propose an efficient framework for genetic subtyping of SARS-CoV-2, the novel coronavirus that causes the COVID-19 pandemic. Efficient viral subtyping enables visualization and modeling of the geographic distribution and temporal dynamics of disease spread. Subtyping thereby advances the development of effective containment strategies and, potentially, therapeutic and vaccine strategies. However, identifying viral subtypes in real-time is challenging: SARS-CoV-2 is a novel virus, and the pandemic is rapidly expanding. Viral subtypes may be difficult to detect due to rapid evolution; founder effects are more significant than selection pressure; and the clustering threshold for subtyping is not standardized. We propose to identify mutational signatures of available SARS-CoV-2 sequences using a population-based approach: an entropy measure followed by frequency analysis. These signatures, Informative Subtype Markers (ISMs), define a compact set of nucleotide sites that characterize the most variable (and thus most informative) positions in the viral genomes sequenced from different individuals. Through ISM compression, we find that certain distant nucleotide variants covary, including non-coding and ORF1ab sites covarying with the D614G spike protein mutation which has become increasingly prevalent as the pandemic has spread. ISMs are also useful for downstream analyses, such as spatiotemporal visualization of viral dynamics. By analyzing sequence data available in the GISAID database, we validate the utility of ISM-based subtyping by comparing spatiotemporal analyses using ISMs to epidemiological studies of viral transmission in Asia, Europe, and the United States. In addition, we show the relationship of ISMs to phylogenetic reconstructions of SARS-CoV-2 evolution, and therefore, ISMs can play an important complementary role to phylogenetic tree-based analysis, such as is done in the Nextstrain project. The developed pipeline dynamically generates ISMs for newly added SARS-CoV-2 sequences and updates the visualization of pandemic spatiotemporal dynamics, and is available on Github at https://github.com/EESI/ISM (Jupyter notebook), https://github.com/EESI/ncov_ism (command line tool) and via an interactive website at https://covid19-ism.coe.drexel.edu/.

Variations in SARS-CoV-2 Spike Protein Cell Epitopes and Glycosylation Profiles During Global Transmission Course of COVID-19.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to an outbreak of a pandemic worldwide. For better understanding the viral spike (S) protein variations and its potential effects on the interaction with the host immune system and also in vaccine development, the cell epitopes, glycosylation profile and their changes during the global transmission course were characterized and compared with SARS-CoV for their glycosylation profile. We analyzed totally 7,813 sequences screened from 8,897 whole genome sequences on GISAID database up to April 26, and 18 S protein amino acid variations with relatively high frequency (≥10−3) were identified. A total of 228 sequences of variants had multiple variations, of note, most of them harboring the D614G mutation. Among the predicted 69 linear B cell epitopes, 175 discontinuous B cell epitopes and 41 cytotoxic T lymphocyte epitopes in the viral S protein, we found that the protein structure and its potential function of some sites changed, such as the linear epitope length shortened and discontinuous epitope disappeared of G476S. In addition, we detected 9 predicted N-glycosylation sites and 3 O-glycosylation sites unique to SARS-CoV-2, but no evidently observed variation of the glycan sites so far. Our findings provided an important snapshot of temporal and geographical distributions on SARS-CoV-2 S protein cell epitopes and glycosylation sites, which would be an essential basis for the selection of vaccine candidates.

Improved binding of SARS-CoV-2 Envelope protein to tight junction-associated PALS1 could play a key role in COVID-19 pathogenesis

The Envelope (E) protein of SARS-CoV-2 is the most enigmatic protein among the four structural ones. Most of its current knowledge is based on the direct comparison to the SARS E protein, initially mistakenly undervalued and subsequently proved to be a key factor in the ER-Golgi localization and in tight junction disruption.We compared the genomic sequences of E protein of SARS-CoV-2, SARS-CoV and the closely related genomes of bats and pangolins obtained from the GISAID and GenBank databases. When compared to the known SARS E protein, we observed a significant difference in amino acid sequence in the C-terminal end of SARS-CoV-2 E protein.Subsequently, in silico modelling analyses of E proteins conformation and docking provide evidences of a strengthened binding of SARS-CoV-2 E protein with the tight junction-associated PALS1 protein. Based on our computational evidences and on data related to SARS-CoV, we believe that SARS-CoV-2 E protein interferes more stably with PALS1 leading to an enhanced epithelial barrier disruption, amplifying the inflammatory processes, and promoting tissue remodelling. These findings raise a warning on the underestimated role of the E protein in the pathogenic mechanism and open the route to detailed experimental investigations.

Impact of lockdown on Covid-19 case fatality rate and viral mutations spread in 7 countries in Europe and North America

BackgroundSevere acute respiratory syndrome CoV-2 (SARS-CoV-2) caused the first coronavirus disease 2019 (COVID-19) outbreak in China and has become a public health emergency of international concern. SARS-CoV-2 outbreak has been declared a pandemic by WHO on March 11th, 2020 and the same month several Countries put in place different lockdown restrictions and testing strategies in order to contain the spread of the virus.MethodsThe calculation of the Case Fatality Rate of SARS-CoV-2 in the Countries selected was made by using the data available at https://github.com/owid/covi-19-data/tree/master/public/data. Case fatality rate was calculated as the ratio between the death cases due to COVID-19, over the total number of SARS-CoV-2 reported cases 14 days before. Standard Case Fatality Rate values were normalized by the Country-specific ρ factor, i.e. the number of PCR tests/1 million inhabitants over the number of reported cases/1 million inhabitants. Case-fatality rates between Countries were compared using proportion test. Post-hoc analysis in the case of more than two groups was performed using pairwise comparison of proportions and p value was adjusted using Holm method. We also analyzed 487 genomic sequences from the GISAID database derived from patients infected by SARS-CoV-2 from January 2020 to April 2020 in Italy, Spain, Germany, France, Sweden, UK and USA. SARS-CoV-2 reference genome was obtained from the GenBank database (NC_045512.2). Genomes alignment was performed using Muscle and Jalview software. We, then, calculated the Case Fatality Rate of SARS-CoV-2 in the Countries selected.ResultsIn this study we analyse how different lockdown strategies and PCR testing capability adopted by Italy, France, Germany, Spain, Sweden, UK and USA have influenced the Case Fatality Rate and the viral mutations spread. We calculated case fatality rates by dividing the death number of a specific day by the number of patients with confirmed COVID-19 infection observed 14 days before and normalized by a ρ factor which takes into account the diagnostic PCR testing capability of each Country and the number of positive cases detected. We notice the stabilization of a clear pattern of mutations at sites nt241, nt3037, nt14408 and nt23403. A novel nonsynonymous SARS-CoV-2 mutation in the spike protein (nt24368) has been found in genomes sequenced in Sweden, which enacted a soft lockdown strategy.ConclusionsStrict lockdown strategies together with a wide diagnostic PCR testing of the population were correlated with a relevant decline of the case fatality rate in different Countries. The emergence of specific patterns of mutations concomitant with the decline in case fatality rate needs further confirmation and their biological significance remains unclear.

Emerging of a SARS-CoV-2 viral strain with a deletion in nsp1

BackgroundThe new Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), which was first detected in Wuhan (China) in December of 2019 is responsible for the current global pandemic. Phylogenetic analysis revealed that it is similar to other betacoronaviruses, such as SARS-CoV and Middle-Eastern Respiratory Syndrome, MERS-CoV. Its genome is ∼ 30 kb in length and contains two large overlapping polyproteins, ORF1a and ORF1ab that encode for several structural and non-structural proteins. The non-structural protein 1 (nsp1) is arguably the most important pathogenic determinant, and previous studies on SARS-CoV indicate that it is both involved in viral replication and hampering the innate immune system response. Detailed experiments of site-specific mutagenesis and in vitro reconstitution studies determined that the mechanisms of action are mediated by (a) the presence of specific amino acid residues of nsp1 and (b) the interaction between the protein and the host’s small ribosomal unit. In fact, substitution of certain amino acids resulted in reduction of its negative effects.MethodsA total of 17,928 genome sequences were obtained from the GISAID database (December 2019 to July 2020) from patients infected by SARS-CoV-2 from different areas around the world. Genomes alignment was performed using MAFFT (REFF) and the nsp1 genomic regions were identified using BioEdit and verified using BLAST. Nsp1 protein of SARS-CoV-2 with and without deletion have been subsequently modelled using I-TASSER.ResultsWe identified SARS-CoV-2 genome sequences, from several Countries, carrying a previously unknown deletion of 9 nucleotides in position 686-694, corresponding to the AA position 241-243 (KSF). This deletion was found in different geographical areas. Structural prediction modelling suggests an effect on the C-terminal tail structure.ConclusionsModelling analysis of a newly identified deletion of 3 amino acids (KSF) of SARS-CoV-2 nsp1 suggests that this deletion could affect the structure of the C-terminal region of the protein, important for regulation of viral replication and negative effect on host’s gene expression. In addition, substitution of the two amino acids (KS) from nsp1 of SARS-CoV was previously reported to revert loss of interferon-alpha expression. The deletion that we describe indicates that SARS-CoV-2 is undergoing profound genomic changes. It is important to: (i) confirm the spreading of this particular viral strain, and potentially of strains with other deletions in the nsp1 protein, both in the population of asymptomatic and pauci-symptomatic subjects, and (ii) correlate these changes in nsp1 with potential decreased viral pathogenicity.

Impact of Early Pandemic Stage Mutations on Molecular Dynamics of SARS-CoV-2 Mpro.

A new coronavirus (SARS-CoV-2) is a global threat to world health and economy. Its dimeric main protease (Mpro), which is required for the proteolytic cleavage of viral precursor proteins, is a good candidate for drug development owing to its conservation and the absence of a human homolog. Improving our understanding of Mpro behavior can accelerate the discovery of effective therapies to reduce mortality. All-atom molecular dynamics (MD) simulations (100 ns) of 50 mutant Mpro dimers obtained from filtered sequences from the GISAID database were analyzed using root-mean-square deviation, root-mean-square fluctuation, Rg, averaged betweenness centrality, and geometry calculations. The results showed that SARS-CoV-2 Mpro essentially behaves in a similar manner to its SAR-CoV homolog. However, we report the following new findings from the variants: (1) Residues GLY15, VAL157, and PRO184 have mutated more than once in SARS CoV-2; (2) the D48E variant has lead to a novel “TSEEMLN”” loop at the binding pocket; (3) inactive apo Mpro does not show signs of dissociation in 100 ns MD; (4) a non-canonical pose for PHE140 widens the substrate binding surface; (5) dual allosteric pockets coinciding with various stabilizing and functional components of the substrate binding pocket were found to display correlated compaction dynamics; (6) high betweenness centrality values for residues 17 and 128 in all Mpro samples suggest their high importance in dimer stability—one such consequence has been observed for the M17I mutation whereby one of the N-fingers was highly unstable. (7) Independent coarse-grained Monte Carlo simulations suggest a relationship between the rigidity/mutability and enzymatic function. Our entire approach combining database preparation, variant retrieval, homology modeling, dynamic residue network (DRN), relevant conformation retrieval from 1-D kernel density estimates from reaction coordinates to other existing approaches of structural analysis, and data visualization within the coronaviral Mpro is also novel and is applicable to other coronaviral proteins.

Dataset of mutational analysis, miRNAs targeting SARS-CoV-2 genes and host gene expression in SARS-CoV and SARS-CoV-2 infections.

The identification of host-miRNAs targeting mutated virus genes is crucial to understand the miRNA mediated host-defense mechanism in virus infections. To understand the mechanism in COVID-19 infections, we collected genome sequences of SARS-CoV-2 with its metadata from the GISAID database (submitted till April 2020) and identified mutational changes in the sequences. The dataset consists of genes with mutation event count and entropy scores. We predicted host-miRNAs targeting the genes in the genomes and compared it with that in related viral species. We have identified 2284 miRNAs targeting MERS genomes, 2074 miRNAs targeting SARS genomes, and 1599 miRNAs targeting SARS-CoV-2 genomes, identified using the miRNA target prediction software miRanda. The host miRNAs targeting SARS-CoV-2 genes were further validated to be anti-viral miRNAs and their role in respiratory diseases through a literature survey, which helped in the identification of 42 conserved antiviral miRNAs. The data could be used to validate the anti-viral role of the predicted miRNAs and design miRNA-based therapeutics against SARS-CoV-2.

Building the National SARS-CoV-2 Laboratory Diagnostic Capacity in Taiwan.

Through early and proactive laboratory testing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes novel coronavirus 2019, Taiwan has demonstrated an efficient and rapid control response to contain the outbreak. Two days after the World Health Organization announced the complete viral genome sequence, the national laboratory of the Taiwan Centers for Disease Control developed a specific real-time reverse transcription polymerase chain reaction (real-time RT-PCR) test for SARS-CoV-2. The national laboratory network was further strengthened through the recruitment of medical centers and regional hospitals distributed throughout most geographical regions of the country. Ultimately, a network of 60 laboratories with a capacity of 7,342 real-time RT-PCR tests per day was established. Between January 14 and August 5, 2020, a total of 158,772 tests were conducted, corresponding to 120,487 cases. Test results were obtained within 24 hours, enabling an efficient and rapid control response.

Genome-wide analysis of SARS-CoV-2 virus strains circulating worldwide implicates\xa0heterogeneity

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), a novel evolutionary divergent RNA virus, is responsible for the present devastating COVID-19 pandemic. To explore the genomic signatures, we comprehensively analyzed 2,492 complete and/or near-complete genome sequences of SARS-CoV-2 strains reported from across the globe to the GISAID database up to 30 March 2020. Genome-wide annotations revealed 1,516 nucleotide-level variations at different positions throughout the entire genome of SARS-CoV-2. Moreover, nucleotide (nt) deletion analysis found twelve deletion sites throughout the genome other than previously reported deletions at coding sequence of the ORF8 (open reading frame), spike, and ORF7a proteins, specifically in polyprotein ORF1ab (n = 9), ORF10 (n = 1), and 3´-UTR (n = 2). Evidence from the systematic gene-level mutational and protein profile analyses revealed a large number of amino acid (aa) substitutions (n = 744), demonstrating the viral proteins heterogeneous. Notably, residues of receptor-binding domain (RBD) showing crucial interactions with angiotensin-converting enzyme 2 (ACE2) and cross-reacting neutralizing antibody were found to be conserved among the analyzed virus strains, except for replacement of lysine with arginine at 378th position of the cryptic epitope of a Shanghai isolate, hCoV-19/Shanghai/SH0007/2020 (EPI_ISL_416320). Furthermore, our results of the preliminary epidemiological data on SARS-CoV-2 infections revealed that frequency of aa mutations were relatively higher in the SARS-CoV-2 genome sequences of Europe (43.07%) followed by Asia (38.09%), and North America (29.64%) while case fatality rates remained higher in the European temperate countries, such as Italy, Spain, Netherlands, France, England and Belgium. Thus, the present method of genome annotation employed at this early pandemic stage could be a promising tool for monitoring and tracking the continuously evolving pandemic situation, the associated genetic variants, and their implications for the development of effective control and prophylaxis strategies.

Genetic variants and source of introduction of SARS-CoV-2 in South America.

After more than 4 months of the COVID‐19 pandemics with genomic information of SARS‐CoV‐2 around the globe, there are more than 1000 complete genomes of this virus. We used 691 genomes from the GISAID database. Several studies have been reporting mutations and hotspots according to viral evolution. Our work intends to show and compare positions that have variants in 30 complete viral genomes from South American countries. We classified strains according to point alterations and portray the source where strains came into this region. Most viruses entered South America from Europe, followed by Oceania. Only Chilean isolates demonstrated a relationship with Asian isolates. Some changes in South American genomes are near to specific domains related to viral replication or the S protein. Our work contributes to the global understanding of which sort of strains are spreading throughout South America, and the differences among them according to the first isolates introduced to this region.

Mortality in COVID-19 disease patients: Correlating the association of major histocompatibility complex (MHC) with severe acute respiratory syndrome 2 (SARS-CoV-2) variants

Genetic factors such as the HLA type of patients may play a role in regard to disease severity and clinical outcome of patients with COVID-19. Taking the data deposited in the GISAID database, we made predictions using the IEDB analysis resource (TepiTool) to gauge how variants in the SARS-CoV-2 genome may change peptide binding to the most frequent MHC-class I and -II alleles in Africa, Asia and Europe. We caracterized how a single mutation in the wildtype sequence of of SARS-CoV-2 could influence the peptide binding of SARS-CoV-2 variants to MHC class II, but not to MHC class I alleles. Assuming the ORF8 (L84S) mutation is biologically significant, selective pressure from MHC class II alleles may select for viral varients and subsequently shape the quality and quantity of cellular immune responses aginast SARS-CoV-2. MHC 4-digit typing along with viral sequence analysis should be considered in studies examining clinical outcomes in patients with COVID-19.

Identification of the nucleotide substitutions in 62 SARS-CoV-2 sequences from Turkey.

A previously unknown coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been shown to cause coronavirus disease 2019 (COVID-19) pandemic. The first case of COVID-19 in Turkey has been declared in March 11th, 2020 and from there on, more than 150,000 people in the country have been diagnosed with the disease. In this study, 62 viral sequences from Turkey, which have been uploaded to GISAID database, were analyzed by means of their nucleotide substitutions in comparison to the reference SARS-CoV-2 genome from Wuhan. Our results indicate that the viral isolates from Turkey harbor some common mutations with the viral strains from Europe, Oceania, North America and Asia. When the mutations were evaluated, C3037T, C14408T and A23403G were found to be the most common nucleotide substitutions among the viral isolates in Turkey, which are mostly seen as linked mutations and are part of a haplotype observed high in Europe.

Multiple assays in a real-time RT-PCR SARS-CoV-2 panel can mitigate the risk of loss of sensitivity by new genomic variants during the COVID-19 outbreak

ObjectivesIn this study, five SARS-CoV-2 PCR assay panels were evaluated against the accumulated genetic variability of the virus to assess the effect on sensitivity of the individual assays. Design or methodsAs of week 21, 2020, the complete set of available SARS-CoV-2 genomes from GISAID and GenBank databases were used in this study. SARS-CoV-2 primer sequences from publicly available panels (WHO, CDC, NMDC, and HKU) and QIAstat-Dx were included in the alignment, and accumulated genetic variability affecting any oligonucleotide annealing was annotated. ResultsA total of 11,627 (34.38%) genomes included single mutations affecting annealing of any PCR assay. Variations in 8,773 (25.94%) genomes were considered as high risk, whereas additional 2,854 (8.43%) genomes presented low frequent single mutations and were predicted to yield no impact on sensitivity. In case of the QIAstat-Dx SARS-CoV-2 Panel, 99.11% of the genomes matched with a 100% coverage all oligonucleotides, and critical variations were tested in vitro corroborating no loss of sensitivity. ConclusionsThis analysis stresses the importance of targeting more than one region in the viral genome for SARS-CoV-2 detection to mitigate the risk of loss of sensitivity due to the unknown mutation rate during this SARS-CoV-2 outbreak.