Spike mutation D614G alters SARS-CoV-2 fitness and neutralization susceptibility

Sequence Analysis of Hot Spot Regions of Spike and RNA-dependent-RNA polymerase (RdRp) Genes of SARS-CoV-2 in Kerman, Iran

Background

Mutations in the SARS-CoV-2 genome might influence pathogenicity, transmission rate, and evasion of the host immune system. Therefore, the purpose of the present study was to investigate the genetic alteration as well as assess their effects on the receptor binding domain (RBD) of the spike and the putative RNA binding site of the RdRp genes of SARS-CoV-2 using bioinformatics tools.

Materials and Method

In this cross-sectional study, 45 confirmed COVID-19 patients using qRT-PCR were included and divided into mild, severe, and critical groups based on the severity of the disease. RNA was extracted from nasopharyngeal swab samples using a commercial kit. RT-PCR was performed to amplify the target sequences of the spike and RdRp genes and sequence them by the Sanger method. Clustal OMEGA, MEGA 11 software, I-mutant tools, SWISS-MODEL, and HDOCK web servers were used for bioinformatics analyses.

Results

The mean age of the patients was 50.68±2.73. The results showed that four of six mutations (L452R, T478K, N501Y, and D614G) in RBD and three of eight in the putative RNA binding site (P314L, E1084D, V1883T) were missense. In the putative RNA binding site, another deletion was discovered. Among missense mutations, N501Y and V1883T were responsible for increasing structural stability, while others were responsible for decreasing it. The various homology models designed showed that these homologies were like the Wuhan model. The molecular docking analysis revealed that the T478K mutation in RBD had the highest binding affinity. In addition, 35 RBD samples (89.7%) and 33 putative RNA binding site samples (84.6%) were similar to the Delta variant.

Conclusion

Our results indicated that double mutations (T478K and N501Y) in the S protein might increase the binding affinity of SARS-CoV-2 to human ACE2 compared to the wild-type (WT) strain. Moreover, variations in the spike and RdRp genes might influence the stability of encoded proteins.

Potential of conserved antigenic sites in development of universal SARS-like coronavirus vaccines

Given pandemic risks of zoonotic SARS-CoV-2 variants and other SARS-like coronaviruses in the future, it is valuable to perform studies on conserved antigenic sites to design universal SARS-like coronavirus vaccines. By using antibodies obtained from convalescent COVID-19 patients, we succeeded in functional comparison of conserved antigenic sites at multiple aspects with each other, and even with SARS-CoV-2 unique antigenic sites, which promotes the cognition of process of humoral immune response to the conserved antigenic sites. The conserved antigenic sites between SARS-CoV-2 and SARS-CoV can effectively induce affinity maturation of cross-binding antibodies, finally resulting in broadly neutralizing antibodies against multiple variants of concern, which provides an important basis for universal vaccine design, however they are subdominant, putatively due to their lower accessibility relative to SARS-CoV-2 unique antigenic sites. Furthermore, we preliminarily design RBDs to improve the immunogenicity of these conserved antigenic sites. Our study focusing on conserved antigenic sites provides insights for promoting the development of universal SARS-like coronavirus vaccines, thereby enhancing our pandemic preparedness.

Bivalent SARS-CoV-2 mRNA vaccines increase breadth of neutralization and protect against the BA.5 Omicron variant

The emergence of SARS-CoV-2 variants in the Omicron lineage with large numbers of substitutions in the spike protein that can evade antibody neutralization has resulted in diminished vaccine efficacy and persistent transmission. One strategy to broaden vaccine-induced immunity is to administer bivalent vaccines that encode for spike proteins from both historical and newly-emerged variant strains. Here, we evaluated the immunogenicity and protective efficacy of two bivalent vaccines that recently were authorized for use in Europe and the United States and contain two mRNAs encoding Wuhan-1 and either BA.1 (mRNA-1273.214) or BA.4/5 (mRNA-1273.222) spike proteins. As a primary immunization series in BALB/c mice, both bivalent vaccines induced broader neutralizing antibody responses than the constituent monovalent vaccines (mRNA-1273 [Wuhan-1], mRNA-1273.529 [BA.1], and mRNA-1273-045 [BA.4/5]). When administered to K18-hACE2 transgenic mice as a booster at 7 months after the primary vaccination series with mRNA-1273, the bivalent vaccines induced greater breadth and magnitude of neutralizing antibodies compared to an mRNA-1273 booster. Moreover, the response in bivalent vaccine-boosted mice was associated with increased protection against BA.5 infection and inflammation in the lung. Thus, boosting with bivalent Omicron-based mRNA-1273.214 or mRNA-1273.222 vaccines enhances immunogenicity and protection against currently circulating SARS-CoV-2 strains.

A comprehensive evaluation of the immune system response and type-I Interferon signaling pathway in hospitalized COVID-19 patients

BACKGROUND

The COVID-19 pandemic has become the world's main life-threatening challenge in the third decade of the twenty-first century. Numerous studies have been conducted on SARS-CoV2 virus structure and pathogenesis to find reliable treatments and vaccines. The present study aimed to evaluate the immune-phenotype and IFN-I signaling pathways of COVID-19 patients with mild and severe conditions.

MATERIAL AND METHODS

A total of 100 COVID-19 patients (50 with mild and 50 with severe conditions) were enrolled in this study. The frequency of CD4 + T, CD8 + T, Th17, Treg, and B lymphocytes beside NK cells was evaluated using flow cytometry. IFN-I downstream signaling molecules, including JAK-1, TYK-2, STAT-1, and STAT-2, and Interferon regulatory factors (IRF) 3 and 7 expressions at RNA and protein status were investigated using real-time PCR and western blotting techniques, respectively. Immune levels of cytokines (e.g., IL-1β, IL-6, IL-17, TNF-α, IL-2R, IL-10, IFN-α, and IFN-β) and the existence of anti-IFN-α autoantibodies were evaluated via enzyme-linked immunosorbent assay (ELISA).

RESULTS

Immune-phenotyping results showed a significant decrease in the absolute count of NK cells, CD4 + T, CD8 + T, and B lymphocytes in COVID-19 patients. The frequency of Th17 and Treg cells showed a remarkable increase and decrease, respectively. All signaling molecules of the IFN-I downstream pathway and IRFs (i.e., JAK-1, TYK-2, STAT-1, STAT-2, IRF-3, and IRF-7) showed very reduced expression levels in COVID-19 patients with the severe condition compared to healthy individuals at both RNA and protein levels. Of 50 patients with severe conditions, 14 had anti-IFN-α autoantibodies in sera. Meanwhile, this result was 2 and 0 for patients with mild symptoms and healthy controls, respectively.

CONCLUSION

Our results indicate a positive association of the existence of anti-IFN-α autoantibodies and immune cells dysregulation with the severity of illness in COVID-19 patients. However, comprehensive studies are necessary to find out more about this context. Video abstract.

A So-Far Overlooked Secondary Conformation State in the Binding Mode of SARS-CoV-2 Spike Protein to Human ACE2 and Its Conversion Rate Are Crucial for Estimating Infectivity Efficacy of the Underlying Virus Variant

Since its outbreak in 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread with high transmission efficiency across the world, putting health care as well as economic systems under pressure. During the course of the pandemic, the originally identified SARS-CoV-2 variant has been multiple times replaced by various mutant versions, which showed enhanced fitness due to increased infection and transmission rates. In order to find an explanation for why SARS-CoV-2 and its emerging mutated versions showed enhanced transmission efficiency compared with SARS-CoV (2002), an enhanced binding affinity of the spike protein to human angiotensin converting enzyme 2 (hACE2) has been proposed by crystal structure analysis and was identified in cell culture models. Kinetic analysis of the interaction of various spike protein constructs with hACE2 was considered to be best described by a Langmuir-based 1:1 stoichiometric interaction. However, we demonstrate in this report that the SARS-CoV-2 spike protein interaction with hACE2 is best described by a two-step interaction, which is defined by an initial binding event followed by a slower secondary rate transition that enhances the stability of the complex by a factor of ~190 (primary versus secondary state) with an overall equilibrium dissociation constant (K_D) of 0.20 nM. In addition, we show that the secondary rate transition is not only present in SARS-CoV-2 wild type ("wt"; Wuhan strain) but also found in the B.1.1.7 variant, where its transition rate is 5-fold increased. IMPORTANCE The current SARS-CoV-2 pandemic is characterized by the high infectivity of SARS-CoV-2 and its derived variants of concern (VOCs). It has been widely assumed that the reason for its increased cell entry compared with SARS-CoV (2002) is due to alterations in the viral spike protein, where single amino acid residue substitutions can increase affinity for hACE2. So far, the interaction of a single unit of the CoV-2 spike protein has been described using the 1:1 Langmuir interaction kinetic. However, we demonstrate here that there is a secondary state binding step that may be essential for novel VOCs in order to further increase their infectivity. These findings are important for quantitatively understanding the infection process of SARS-CoV-2 and characterization of emerging SARS-CoV-2 variants of spike proteins. Thus, they provide a tool for predicting the potential infectivity of the respective viral variants based on secondary rate transition and secondary complex stability.

Characterization of SARS-CoV-2 Spike mutations important for infection of mice and escape from human immune sera

Due to differences in human and murine angiotensin converting enzyme 2 (ACE-2) receptor, initially available SARS-CoV-2 isolates could not infect mice. Here we show that serial passaging of USA-WA1/2020 strain in mouse lungs results in "mouse-adapted" SARS-CoV-2 (MA-SARS-CoV-2) with mutations in S, M, and N genes, and a twelve-nucleotide insertion in the S gene. MA-SARS-CoV-2 infection causes mild disease, with more pronounced morbidity depending on genetic background and in aged and obese mice. Two mutations in the S gene associated with mouse adaptation (N501Y, H655Y) are present in SARS-CoV-2 variants of concern (VoCs). N501Y in the receptor binding domain of viruses of the B.1.1.7, B.1.351, P.1 and B.1.1.529 lineages (Alpha, Beta, Gamma and Omicron variants) is associated with high transmissibility and allows VoCs to infect wild type mice. We further show that S protein mutations of MA-SARS-CoV-2 do not affect neutralization efficiency by human convalescent and post vaccination sera.

N-glycosylation profiles of the SARS-CoV-2 spike D614G mutant and its ancestral protein characterized by advanced mass spectrometry

N-glycosylation plays an important role in the structure and function of membrane and secreted proteins. The spike protein on the surface of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, is heavily glycosylated and the major target for developing vaccines, therapeutic drugs and diagnostic tests. The first major SARS-CoV-2 variant carries a D614G substitution in the spike (S-D614G) that has been associated with altered conformation, enhanced ACE2 binding, and increased infectivity and transmission. In this report, we used mass spectrometry techniques to characterize and compare the N-glycosylation of the wild type (S-614D) or variant (S-614G) SARS-CoV-2 spike glycoproteins prepared under identical conditions. The data showed that half of the N-glycosylation sequons changed their distribution of glycans in the S-614G variant. The S-614G variant showed a decrease in the relative abundance of complex-type glycans (up to 45%) and an increase in oligomannose glycans (up to 33%) on all altered sequons. These changes led to a reduction in the overall complexity of the total N-glycosylation profile. All the glycosylation sites with altered patterns were in the spike head while the glycosylation of three sites in the stalk remained unchanged between S-614G and S-614D proteins.

Efficacy of Interferon-β in Moderate-to-Severe Hospitalised Cases of COVID-19: A Systematic Review and Meta-analysis

Background and Objective

Interferon-β, as with several other anti-viral agents, has been investigated as a treatment option for COVID-19 as a repurposed drug. The present study is a systematic review and meta-analysis of interferon-β to determine its efficacy among moderate-to-severe COVID-19 patients.

Methods

A systematic literature search was done using relevant terms for ‘COVID-19’ and ‘interferon-β’. Randomised controlled trials (RCT) evaluating the efficacy of interferon-β in COVID-19 were included. Data were extracted for outcome measures, namely mortality, time to clinical improvement and length of hospital stay. Random effects meta-analysis was performed using RevMan V.5.4.1 to calculate overall effect estimate as odds ratio/hazard ratio for categorical variables and mean difference for continuous variable.

Result

Eight RCTs were eligible for qualitative synthesis and seven for meta-analysis. The overall effect estimate (odds ratio [OR] 0.59; 95 % CI 0.91, 1.12) and (mean difference [MD] − 1.41; 95 % CI − 2.84, 0.02) indicated no statistically significant difference between effect of IFN-β and that of control on mortality and length of hospital stay, respectively. However, the overall effect estimate (hazard ratio [HR] 1.95; 95 % CI 1.36, 2.79) denoted a favourable effect of INF-β on reducing the time to clinical improvement in moderate-to-severe COVID-19 patients.

Conclusion

Addition of interferon-β to standard of care resulted in significant reduction in time to clinical improvement but no significant benefit in terms of reduction in mortality and length of hospital stay in moderate-to-severe cases of COVID-19.

The origin of human mutation in light of genomic data

Despite years of active research into the role of DNA repair and replication in mutagenesis, surprisingly little is known about the origin of spontaneous human mutation in the germ line. With the advent of high-throughput sequencing, genome-scale data have revealed statistical properties of mutagenesis in humans. These properties include variation of the mutation rate and spectrum along the genome at different scales in relation to epigenomic features and dependency on parental age. Moreover, mutations originated in mothers are less frequent than mutations originated in fathers and have a distinct genomic distribution. Statistical analyses that interpret these patterns in the context of known biochemistry can provide mechanistic models of mutagenesis in humans.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called ‘Tiled-ClickSeq’, which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5’UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Characterization and structural basis of a lethal mouse-adapted SARS-CoV-2

There is an urgent need for animal models to study SARS-CoV-2 pathogenicity. Here, we generate and characterize a novel mouse-adapted SARS-CoV-2 strain, MASCp36, that causes severe respiratory symptoms, and mortality. Our model exhibits age- and gender-related mortality akin to severe COVID-19. Deep sequencing identified three amino acid substitutions, N501Y, Q493H, and K417N, at the receptor binding domain (RBD) of MASCp36, during in vivo passaging. All three RBD mutations significantly enhance binding affinity to its endogenous receptor, ACE2. Cryo-electron microscopy analysis of human ACE2 (hACE2), or mouse ACE2 (mACE2), in complex with the RBD of MASCp36, at 3.1 to 3.7 Å resolution, reveals the molecular basis for the receptor-binding switch. N501Y and Q493H enhance the binding affinity to hACE2, whereas triple mutations at N501Y/Q493H/K417N decrease affinity and reduce infectivity of MASCp36. Our study provides a platform for studying SARS-CoV-2 pathogenesis, and unveils the molecular mechanism for its rapid adaptation and evolution.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

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.

Coincident rapid expansion of two SARS-CoV-2 lineages with enhanced infectivity in Nigeria

The emergence of new SARS-CoV-2 variants with enhanced transmissibility or decreased susceptibility to immune responses is a major threat to global efforts to end the coronavirus disease 2019 (COVID-19) pandemic. Disparities in viral genomic surveillance capabilities and efforts have resulted in gaps in our understanding of the viral population dynamics across the globe. Nigeria, despite having the largest population of any nation in Africa, has had relatively little SARS-CoV-2 sequence data made publicly available. Here we report the whole-genome sequences of 74 SARS-CoV-2 isolates collected from individuals in Oyo State, Nigeria in January 2021. Most isolates belonged to either the B.1.1.7 Alpha "variant of concern" or the B.1.525 Eta lineage, which is currently considered a "variant of interest" containing multiple spike protein mutations previously associated with enhanced transmissibility and possible immune escape. Nigeria has the highest reported frequency of the B.1.525 lineage globally with phylogenetic characteristics consistent with a recent monophyletic origin and rapid expansion. Spike protein from the B.1.525 lineage displayed both increased infectivity and decreased neutralization by convalescent sera compared to Spike proteins from other clades. These results, along with indications that the virus is outpacing the B.1.1.7 lineage in Nigeria, suggest that the B.1.525 lineage represents another "variant of concern" and further underline the importance of genomic surveillance in undersampled regions across the globe.

Identifying Primate ACE2 Variants That Confer Resistance to SARS-CoV-2

SARS-CoV-2 infects humans through the binding of viral S-protein (spike protein) to human angiotensin I converting enzyme 2 (ACE2). The structure of the ACE2-S-protein complex has been deciphered and we focused on the 27 ACE2 residues that bind to S-protein. From human sequence databases, we identified nine ACE2 variants at ACE2-S-protein binding sites. We used both experimental assays and protein structure analysis to evaluate the effect of each variant on the binding affinity of ACE2 to S-protein. We found one variant causing complete binding disruption, two and three variants, respectively, strongly and mildly reducing the binding affinity, and two variants strongly enhancing the binding affinity. We then collected the ACE2 gene sequences from 57 nonhuman primates. Among the 6 apes and 20 Old World monkeys (OWMs) studied, we found no new variants. In contrast, all 11 New World monkeys (NWMs) studied share four variants each causing a strong reduction in binding affinity, the Philippine tarsier also possesses three such variants, and 18 of the 19 prosimian species studied share one variant causing a strong reduction in binding affinity. Moreover, one OWM and three prosimian variants increased binding affinity by >50%. Based on these findings, we proposed that the common ancestor of primates was strongly resistant to and that of NWMs was completely resistant to SARS-CoV-2 and so is the Philippine tarsier, whereas apes and OWMs, like most humans, are susceptible. This study increases our understanding of the differences in susceptibility to SARS-CoV-2 infection among primates.

Scrutinizing Coronaviruses Using Publicly Available Bioinformatic Tools: The Viral Structural Proteins as a Case Study

Since early 2020, the world suffers from a new beta-coronavirus, called SARS-CoV-2, that has devastating effects globally due to its associated disease, Covid-19. Until today, Covid-19, which not only causes life-threatening lung infections but also impairs various other organs and tissues, has killed hundreds of thousands of people and caused irreparable damage to many others. Since the very onset of the pandemic, huge efforts were made worldwide to fully understand this virus and numerous studies were, and still are, published. Many of these deal with structural analyses of the viral spike glycoprotein and with vaccine development, antibodies and antiviral molecules or immunomodulators that are assumed to become essential tools in the struggle against the virus. This paper summarizes knowledge on the properties of the four structural proteins (spike protein S, membrane protein M, envelope protein E and nucleocapsid protein N) of the SARS-CoV-2 virus and its relatives, SARS-CoV and MERS-CoV, that emerged few years earlier. Moreover, attention is paid to ways to analyze such proteins using freely available bioinformatic tools and, more importantly, to bring these proteins alive by looking at them on a computer/laptop screen with the easy-to-use but highly performant and interactive molecular graphics program DeepView. It is hoped that this paper will stimulate non-bioinformaticians and non-specialists in structural biology to scrutinize these and other macromolecules and as such will contribute to establishing procedures to fight these and maybe other forthcoming viruses.

Implications of the Novel Mutations in the SARS-CoV-2 Genome for Transmission, Disease Severity, and the Vaccine Development

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative virus of the coronavirus disease 2019 (COVID-19), has been identified in China in late December 2019. SARS-CoV-2 is an enveloped, positive-sense, single-stranded RNA betacoronavirus of the Coronaviridae family. Coronaviruses have genetic proofreading mechanism that corrects copying mistakes and thus SARS-CoV-2 genetic diversity is extremely low. Despite lower mutation rate of the virus, researchers have detected a total of 12,706 mutations in the SARS-CoV-2 genome, the majority of which were single nucleotide polymorphisms. Sequencing data revealed that the SARS-CoV-2 accumulates two-single nucleotide mutations per month in its genome. Recently, an amino acid aspartate (D) to glycine (G) (D614G) mutation due to an adenine to guanine nucleotide change at position 23,403 at the 614th amino-acid position of the spike protein in the original reference genotype has been identified. The SARS-CoV-2 viruses that carry the spike protein D614G mutation have become dominant variant around the world. The D614G mutation has been found to be associated with 3 other mutations in the spike protein. Clinical and pseudovirus experimental studies have demonstrated that the spike protein D614G mutation alters the virus phenotype. However, the impact of the mutation on the rate of transmission between people, disease severity and the vaccine and therapeutic development remains unclear. Three variants of SARS-CoV-2 have recently been identified. They are B.1.1.7 (UK) variant, B.1.351 (N501Y.V2, South African) variant and B.1.1.28 (Brazilian) variant. Epidemiological data suggest that they have a higher transmissibility than the original variant. There are reports that some vaccines are less efficacious against the B.1.351 variant. This review article discusses the effects of novel mutations in the SARS-CoV-2 genome on transmission, clinical outcomes and vaccine development.

European context of the diversity and phylogenetic position of SARS-CoV-2 sequences from Polish COVID-19 patients

To provide a comprehensive analysis of the SARS-CoV-2 sequence diversity in Poland in the European context. All publicly available (n = 115; GISAID database) whole-genome SARS-Cov-2 sequences from Polish samples, including those obtained during coronavirus testing performed in our COVID-19 Lab, were examined. Multiple sequence alignment of Polish isolates, phylogenetic analysis (ML tree), and multidimensional scaling (based on the pairwise DNA distances) were complemented by the comparison of the coronavirus clades frequency and diversity in the subset of over 5000 European GISAID sequences. Approximately seventy-seven percent of isolates in the European dataset carried frequent and ubiquitously found haplotypes; the remaining haplotype diversity was population-specific and resulted from population-specific mutations, homoplasies, and recombinations. Coronavirus strains circulating in Poland represented the variability found in other European countries. The prevalence of clades circulating in Poland was shifted in favor of GR, both in terms of the diversity (number of distinct haplotypes) and the frequency (number of isolates) of the clade. Polish-specific haplotypes were rare and could be explained by changes affecting common European strains. The analysis of the whole viral genomes allowed detection of several tight clusters of isolates, presumably reflecting local outbreaks. New mutations, homoplasies, and, to a smaller extent, recombinations increase SARS-CoV-2 haplotype diversity, but the majority of these variants do not increase in frequency and remains rare and population-specific. The spectrum of SARS-CoV-2 haplotypes in the Polish dataset reflects many independent transfers from a variety of sources, followed by many local outbreaks. The prevalence of the sequences belonging to the GR clade among Polish isolates is consistent with the European trend of the GR clade frequency increase.

Evolution of SARS-CoV-2: Review of Mutations, Role of the Host Immune System

Since the reporting of the first cases of coronavirus in China and the publication of the first sequence of SARS-CoV-2 in December 2019, the virus has undergone numerous mutations. In Europe, the spring outbreak (March-April) was followed by a drop in the number of cases and deaths. The disease may have evolved into a milder form. The increase in PCR-positive cases in late summer 2020 did not lead to the expected increase in hospitalizations, ICU admissions, and deaths, based on the severity of the disease in the spring. This difference in disease severity could be due to factors independent of the virus or to the evolution of the virus. This review attempts to identify the mutations that have appeared since the beginning of the pandemic and their role in the temporal evolution of the pandemic. There are a cell and humoral type cross-reactivity in a large part of the population to common cold coronaviruses (HCoVs) and SARS-CoV-2. Evolutionarily important mutations and deletions have emerged in the SARS-CoV-2 genes encoding proteins that interact with the host immune system. In addition, one of the major mutations (in viral polymerase) is logically associated with a higher frequency of mutations throughout the genome. This frequency fluctuates over time and shows a peak at the time when the epidemic was most active. The rate of mutations in proteins involved in the relationship to the immune system continues to increase after the first outbreak. The cross-reactivity on the 1 hand and the viral mutations observed on the other hand could explain the evolution of the pandemic until the summer of 2020, partly due to the evolution of the virus in relation to the host immune system. The immunization campaign began in December 2020: concerns are emerging about a possible escape of the circulating variants vaccines in early 2021. These variants could also escape immunity acquired through infection with the 2020 strains.

Unraveling the stability landscape of mutations in the SARS-CoV-2 receptor-binding domain

The interaction between the receptor-binding domain (RBD) of the SARS-CoV-2 spike glycoprotein and the ACE2 enzyme is believed to be the entry point of the virus into various cells in the body, including the lungs, heart, liver, and kidneys. The current focus of several therapeutic design efforts explores attempts at affecting the binding potential between the two proteins to limit the activity of the virus and disease progression. In this work, we analyze the stability of the spike protein under all possible single-point mutations in the RBD and computationally explore mutations that can affect the binding with the ACE2 enzyme. We unravel the mutation landscape of the receptor region and assess the toxicity potential of single and multi-point mutations, generating insights for future vaccine efforts on mutations that might further stabilize the spike protein and increase its infectivity. We developed a tool, called SpikeMutator, to construct full atomic protein structures of the mutant spike proteins and shared a database of 3800 single-point mutant structures. We analyzed the recent 65,000 reported spike sequences across the globe and observed the emergence of stable multi-point mutant structures. Using the landscape, we searched through 7.5 million possible 2-point mutation combinations and report that the (R355D K424E) mutation produces one of the strongest spike proteins that therapeutic efforts should investigate for the sake of developing effective vaccines.

Developing RT-LAMP Assays for Detection of SARS-CoV-2 in Saliva

The coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 has killed millions of people worldwide. The current crisis has created an unprecedented demand for rapid test of SARS-CoV-2 infection. Reverse transcription loop-mediated isothermal amplification (RT-LAMP) is a fast and convenient method to amplify and identify the transcripts of a targeted pathogen. However, the sensitivity and specificity of RT-LAMP were generally regarded as inferior when compared with the gold standard RT-qPCR. To address this issue, we combined bioinformatic and experimental analyses to improve the assay performance for COVID-19 diagnosis. First, we developed an improved algorithm to design LAMP primers targeting the nucleocapsid (N), membrane (M), and spike (S) genes of SARS-CoV-2. Next, we rigorously validated these new assays for their efficacy and specificity. Further, we demonstrated that multiplexed RT-LAMP assays could directly detect as low as a few copies of SARS-CoV-2 RNA in saliva, without the need of RNA isolation. Importantly, further testing using saliva samples from COVID-19 patients indicated that the new RT-LAMP assays were in total agreement in sensitivity and specificity with standard RT-qPCR. In summary, our new LAMP primer design algorithm along with the validated assays provide a fast and reliable method for the diagnosis of COVID-19 cases.

Neurotropic Viruses, Astrocytes, and COVID-19

At the end of 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was discovered in China, causing a new coronavirus disease, termed COVID-19 by the WHO on February 11, 2020. At the time of this paper (January 31, 2021), more than 100 million cases have been recorded, which have claimed over 2 million lives worldwide. The most important clinical presentation of COVID-19 is severe pneumonia; however, many patients present various neurological symptoms, ranging from loss of olfaction, nausea, dizziness, and headache to encephalopathy and stroke, with a high prevalence of inflammatory central nervous system (CNS) syndromes. SARS-CoV-2 may also target the respiratory center in the brainstem and cause silent hypoxemia. However, the neurotropic mechanism(s) by which SARS-CoV-2 affects the CNS remain(s) unclear. In this paper, we first address the involvement of astrocytes in COVID-19 and then elucidate the present knowledge on SARS-CoV-2 as a neurotropic virus as well as several other neurotropic flaviviruses (with a particular emphasis on the West Nile virus, tick-borne encephalitis virus, and Zika virus) to highlight the neurotropic mechanisms that target astroglial cells in the CNS. These key homeostasis-providing cells in the CNS exhibit many functions that act as a favorable milieu for virus replication and possibly a favorable environment for SARS-CoV-2 as well. The role of astrocytes in COVID-19 pathology, related to aging and neurodegenerative disorders, and environmental factors, is discussed. Understanding these mechanisms is key to better understanding the pathophysiology of COVID-19 and for developing new strategies to mitigate the neurotropic manifestations of COVID-19.

Combination of Angiotensin (1-7) Agonists and Convalescent Plasma as a New Strategy to Overcome Angiotensin Converting Enzyme 2 (ACE2) Inhibition for the Treatment of COVID-19

Coronavirus disease-2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is currently the most concerning health problem worldwide. SARS-CoV-2 infects cells by binding to angiotensin-converting enzyme 2 (ACE2). It is believed that the differential response to SARS-CoV-2 is correlated with the differential expression of ACE2. Several reports proposed the use of ACE2 pharmacological inhibitors and ACE2 antibodies to block viral entry. However, ACE2 inhibition is associated with lung and cardiovascular pathology and would probably increase the pathogenesis of COVID-19. Therefore, utilizing ACE2 soluble analogs to block viral entry while rescuing ACE2 activity has been proposed. Despite their protective effects, such analogs can form a circulating reservoir of the virus, thus accelerating its spread in the body. Levels of ACE2 are reduced following viral infection, possibly due to increased viral entry and lysis of ACE2 positive cells. Downregulation of ACE2/Ang (1-7) axis is associated with Ang II upregulation. Of note, while Ang (1-7) exerts protective effects on the lung and cardiovasculature, Ang II elicits pro-inflammatory and pro-fibrotic detrimental effects by binding to the angiotensin type 1 receptor (AT1R). Indeed, AT1R blockers (ARBs) can alleviate the harmful effects associated with Ang II upregulation while increasing ACE2 expression and thus the risk of viral infection. Therefore, Ang (1-7) agonists seem to be a better treatment option. Another approach is the transfusion of convalescent plasma from recovered patients with deteriorated symptoms. Indeed, this appears to be promising due to the neutralizing capacity of anti-COVID-19 antibodies. In light of these considerations, we encourage the adoption of Ang (1-7) agonists and convalescent plasma conjugated therapy for the treatment of COVID-19 patients. This therapeutic regimen is expected to be a safer choice since it possesses the proven ability to neutralize the virus while ensuring lung and cardiovascular protection through modulation of the inflammatory response.

Functional importance of the D614G mutation in the SARS-CoV-2 spike protein

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an enveloped virus which binds its cellular receptor angiotensin-converting enzyme 2 (ACE2) and enters hosts cells through the action of its spike (S) glycoprotein displayed on the surface of the virion. Compared to the reference strain of SARS-CoV-2, the majority of currently circulating isolates possess an S protein variant characterized by an aspartic acid-to-glycine substitution at amino acid position 614 (D614G). Residue 614 lies outside the receptor binding domain (RBD) and the mutation does not alter the affinity of monomeric S protein for ACE2. However, S(G614), compared to S(D614), mediates more efficient ACE2-mediated transduction of cells by S-pseudotyped vectors and more efficient infection of cells and animals by live SARS-CoV-2. This review summarizes and synthesizes the epidemiological and functional observations of the D614G spike mutation, with focus on the biochemical and cell-biological impact of this mutation and its consequences for S protein function. We further discuss the significance of these recent findings in the context of the current global pandemic.

The D614G mutations in the SARS-CoV-2 spike protein: Implications for viral infectivity, disease severity and vaccine design

The development of the SARS-CoV-2 pandemic has prompted an extensive worldwide sequencing effort to characterise the geographical spread and molecular evolution of the virus. A point mutation in the spike protein, D614G, emerged as the virus spread from Asia into Europe and the USA, and has rapidly become the dominant form worldwide. Here we review how the D614G variant was identified and discuss recent evidence about the effect of the mutation on the characteristics of the virus, clinical outcome of infection and host immune response.

Effect of mutations in the SARS-CoV-2 spike protein on protein stability, cleavage, and cell-cell fusion function

The SARS-CoV-2 spike protein (S) is the sole viral protein responsible for both viral binding to a host cell and the membrane fusion event needed for cell entry. In addition to facilitating fusion needed for viral entry, S can also drive cell-cell fusion, a pathogenic effect observed in the lungs of SARS-CoV-2 infected patients. While several studies have investigated S requirements involved in viral particle entry, examination of S stability and factors involved in S cell-cell fusion remain limited. We demonstrate that S must be processed at the S1/S2 border in order to mediate cell-cell fusion, and that mutations at potential cleavage sites within the S2 subunit alter S processing at the S1/S2 border, thus preventing cell-cell fusion. We also identify residues within the internal fusion peptide and the cytoplasmic tail that modulate S cell-cell fusion. Additionally, we examine S stability and protein cleavage kinetics in a variety of mammalian cell lines, including a bat cell line related to the likely reservoir species for SARS-CoV-2, and provide evidence that proteolytic processing alters the stability of the S trimer. This work therefore offers insight into S stability, proteolytic processing, and factors that mediate S cell-cell fusion, all of which help give a more comprehensive understanding of this highly sought-after therapeutic target.

Paired heavy and light chain signatures contribute to potent SARS-CoV-2 neutralization in public antibody responses

Understanding protective mechanisms of antibody recognition can inform vaccine and therapeutic strategies against SARS-CoV-2. We discovered a new antibody, 910-30, that targets the SARS-CoV-2 ACE2 receptor binding site as a member of a public antibody response encoded by IGHV3-53/IGHV3-66 genes. We performed sequence and structural analyses to explore how antibody features correlate with SARS-CoV-2 neutralization. Cryo-EM structures of 910-30 bound to the SARS-CoV-2 spike trimer revealed its binding interactions and ability to disassemble spike. Despite heavy chain sequence similarity, biophysical analyses of IGHV3-53/3-66 antibodies highlighted the importance of native heavy:light pairings for ACE2 binding competition and for SARS-CoV-2 neutralization. We defined paired heavy:light sequence signatures and determined antibody precursor prevalence to be ~1 in 44,000 human B cells, consistent with public antibody identification in several convalescent COVID-19 patients. These data reveal key structural and functional neutralization features in the IGHV3-53/3-66 public antibody class to accelerate antibody-based medical interventions against SARS-CoV-2.

Highlights

A molecular study of IGHV3-53/3-66 public antibody responses reveals critical heavy and light chain features for potent neutralizationCryo-EM analyses detail the structure of a novel public antibody class member, antibody 910-30, in complex with SARS-CoV-2 spike trimerCryo-EM data reveal that 910-30 can both bind assembled trimer and can disassemble the SARS-CoV-2 spikeSequence-structure-function signatures defined for IGHV3-53/3-66 class antibodies including both heavy and light chainsIGHV3-53/3-66 class precursors have a prevalence of 1:44,000 B cells in healthy human antibody repertoires.

Isolation and characterization of an early SARS-CoV-2 isolate from the 2020 epidemic in Medellín, Colombia

Introduction: SARS-CoV-2 has been identified as the new coronavirus causing an outbreak of acute respiratory disease in China in December, 2019. This disease, currently named COVID-19, has been declared as a pandemic by the World Health Organization (WHO). The first case of COVID-19 in Colombia was reported on March 6, 2020. Here we characterize an early SARS-CoV-2 isolate from the pandemic recovered in April, 2020. Objective: To describe the isolation and characterization of an early SARS-CoV-2 isolate from the epidemic in Colombia. Materials and methods: A nasopharyngeal specimen from a COVID-19 positive patient was inoculated on different cell lines. To confirm the presence of SARS-CoV-2 on cultures we used qRT-PCR, indirect immunofluorescence assay, transmission and scanning electron microscopy, and next-generation sequencing. Results: We determined the isolation of SARS-CoV-2 in Vero-E6 cells by the appearance of the cytopathic effect three days post-infection and confirmed it by the positive results in the qRT-PCR and the immunofluorescence with convalescent serum. Transmission and scanning electron microscopy images obtained from infected cells showed the presence of structures compatible with SARS-CoV-2. Finally, a complete genome sequence obtained by next-generation sequencing allowed classifying the isolate as B.1.5 lineage. Conclusion: The evidence presented in this article confirms the first isolation of SARSCoV-2 in Colombia. In addition, it shows that this strain behaves in cell culture in a similar way to that reported in the literature for other isolates and that its genetic composition is consistent with the predominant variant in the world. Finally, points out the importance of viral isolation for the detection of neutralizing antibodies, for the genotypic and phenotypic characterization of the strain and for testing compounds with antiviral potential.

Molecular Architecture of Early Dissemination and Massive Second Wave of the SARS-CoV-2 Virus in a Major Metropolitan Area

We sequenced the genomes of 5,085 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains causing two coronavirus disease 2019 (COVID-19) disease waves in metropolitan Houston, TX, an ethnically diverse region with 7 million residents. The genomes were from viruses recovered in the earliest recognized phase of the pandemic in Houston and from viruses recovered in an ongoing massive second wave of infections. The virus was originally introduced into Houston many times independently. Virtually all strains in the second wave have a Gly614 amino acid replacement in the spike protein, a polymorphism that has been linked to increased transmission and infectivity. Patients infected with the Gly614 variant strains had significantly higher virus loads in the nasopharynx on initial diagnosis. We found little evidence of a significant relationship between virus genotype and altered virulence, stressing the linkage between disease severity, underlying medical conditions, and host genetics. Some regions of the spike protein-the primary target of global vaccine efforts-are replete with amino acid replacements, perhaps indicating the action of selection. We exploited the genomic data to generate defined single amino acid replacements in the receptor binding domain of spike protein that, importantly, produced decreased recognition by the neutralizing monoclonal antibody CR3022. Our report represents the first analysis of the molecular architecture of SARS-CoV-2 in two infection waves in a major metropolitan region. The findings will help us to understand the origin, composition, and trajectory of future infection waves and the potential effect of the host immune response and therapeutic maneuvers on SARS-CoV-2 evolution.IMPORTANCE There is concern about second and subsequent waves of COVID-19 caused by the SARS-CoV-2 coronavirus occurring in communities globally that had an initial disease wave. Metropolitan Houston, TX, with a population of 7 million, is experiencing a massive second disease wave that began in late May 2020. To understand SARS-CoV-2 molecular population genomic architecture and evolution and the relationship between virus genotypes and patient features, we sequenced the genomes of 5,085 SARS-CoV-2 strains from these two waves. Our report provides the first molecular characterization of SARS-CoV-2 strains causing two distinct COVID-19 disease waves.

Structural impact on SARS-CoV-2 spike protein by D614G substitution

Substitution for aspartic acid by glycine at position 614 in the spike (S) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the ongoing pandemic, appears to facilitate rapid viral spread. The G614 variant has now replaced the D614-carrying virus as the dominant circulating strain. We report here cryo-EM structures of a full-length S trimer carrying G614, which adopts three distinct prefusion conformations differing primarily by the position of one receptor-binding domain (RBD). A loop disordered in the D614 S trimer wedges between domains within a protomer in the G614 spike. This added interaction appears to prevent premature dissociation of the G614 trimer, effectively increasing the number of functional spikes and enhancing infectivity. The loop transition may also modulate structural rearrangements of S protein required for membrane fusion. These findings extend our understanding of viral entry and suggest an improved immunogen for vaccine development.