Emergence and spread of SARS-CoV-2 variants from farmed mink to humans and back during the epidemic in Denmark, June-November 2020.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) not only caused the COVID-19 pandemic but also had a major impact on farmed mink production in several European countries. In Denmark, the entire population of farmed mink (over 15 million animals) was culled in late 2020. During the period of June to November 2020, mink on 290 farms (out of about 1100 in the country) were shown to be infected with SARS-CoV-2. Genome sequencing identified changes in the virus within the mink and it is estimated that about 4000 people in Denmark became infected with these mink virus variants. However, the routes of transmission of the virus to, and from, the mink have been unclear. Phylogenetic analysis revealed the generation of multiple clusters of the virus within the mink. Detailed analysis of changes in the virus during replication in mink and, in parallel, in the human population in Denmark, during the same time period, has been performed here. The majority of cases in mink involved variants with the Y453F substitution and the H69/V70 deletion within the Spike (S) protein; these changes emerged early in the outbreak. However, further introductions of the virus, by variants lacking these changes, from the human population into mink also occurred. Based on phylogenetic analysis of viral genome data, we estimate, using a conservative approach, that about 17 separate examples of mink to human transmission occurred in Denmark but up to 59 such events (90% credible interval: (39-77)) were identified using parsimony to count cross-species jumps on transmission trees inferred using Bayesian methods. Using the latter approach, 136 jumps (90% credible interval: (117-164)) from humans to mink were found, which may underlie the farm-to-farm spread. Thus, transmission of SARS-CoV-2 from humans to mink, mink to mink, from mink to humans and between humans were all observed.

Evolutionary and Phylogenetic Dynamics of SARS-CoV-2 Variants: A Genetic Comparative Study of Taiyuan and Wuhan Cities of China.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a positive-sense, single-stranded RNA genome-containing virus which has infected millions of people all over the world. The virus has been mutating rapidly enough, resulting in the emergence of new variants and sub-variants which have reportedly been spread from Wuhan city in China, the epicenter of the virus, to the rest of China and all over the world. The occurrence of mutations in the viral genome, especially in the viral spike protein region, has resulted in the evolution of multiple variants and sub-variants which gives the virus the benefit of host immune evasion and thus renders modern-day vaccines and therapeutics ineffective. Therefore, there is a continuous need to study the genetic characteristics and evolutionary dynamics of the SARS-CoV-2 variants. Hence, in this study, a total of 832 complete genomes of SARS-CoV-2 variants from the cities of Taiyuan and Wuhan in China was genetically characterized and their phylogenetic and evolutionary dynamics studied using phylogenetics, genetic similarity, and phylogenetic network analyses. This study shows that the four most prevalent lineages in Taiyuan and Wuhan are as follows: the Omicron lineages EG.5.1.1, followed by HK.3, FY.3, and XBB.1.16 (Pangolin classification), and clades 23F (EG.5.1), followed by 23H (HK.3), 22F (XBB), and 23D (XBB.1.9) (Nextclade classification), and lineage B followed by the Omicron FY.3, lineage A, and Omicron FL.2.3 (Pangolin classification), and the clades 19A, followed by 22F (XBB), 23F (EG.5.1), and 23H (HK.3) (Nextclade classification), respectively. Furthermore, our genetic similarity analysis show that the SARS-CoV-2 clade 19A-B.4 from Wuhan (name starting with 412981) has the least genetic similarity of about 95.5% in the spike region of the genome as compared to the query sequence of Omicron XBB.2.3.2 from Taiyuan (name starting with 18495234), followed by the Omicron FR.1.4 from Taiyuan (name starting with 18495199) with ~97.2% similarity and Omicron DY.3 (name starting with 17485740) with ~97.9% similarity. The rest of the variants showed ≥98% similarity with the query sequence of Omicron XBB.2.3.2 from Taiyuan (name starting with 18495234). In addition, our recombination analysis results show that the SARS-CoV-2 variants have three statistically significant recombinant events which could have possibly resulted in the emergence of Omicron XBB.1.16 (recombination event 3), FY.3 (recombination event 5), and FL.2.4 (recombination event 7), suggesting some very important information regarding viral evolution. Also, our phylogenetic tree and network analyses show that there are a total of 14 clusters and more than 10,000 mutations which may have probably resulted in the emergence of cluster-I, followed by 47 mutations resulting in the emergence of cluster-II and so on. The clustering of the viral variants of both cities reveals significant information regarding the phylodynamics of the virus among them. The results of our temporal phylogenetic analysis suggest that the variants of Taiyuan have likely emerged as independent variants separate from the variants of Wuhan. This study, to the best of our knowledge, is the first ever genetic comparative study between Taiyuan and Wuhan cities in China. This study will help us better understand the virus and cope with the emergence and spread of new variants at a local as well as an international level, and keep the public health authorities informed for them to make better decisions in designing new viral vaccines and therapeutics. It will also help the outbreak investigators to better examine any future outbreak.

[A Bioinformatics Method for Identification of Human Proteases Active against Viral Envelope Glycoproteins: A Case Study on the SARS-CoV-2 Spike Protein].

Many viruses, including SARS-CoV-2, the coronavirus responsible for the COVID-19 pandemic, enter host cells through a process of cell-viral membrane fusion that is activated by proteolytic enzymes. Typically, these enzymes are host cell proteases. Identifying the proteases that activate the virus is not a simple task but is important for the development of new antiviral drugs. In this study, we developed a bioinformatics method for identifying proteases that can cleave viral envelope glycoproteins. The proposed approach involves the use of predictive models for the substrate specificity of human proteases and the application of a structural analysis method for predicting the vulnerability of protein regions to proteolysis based on their 3D structures. Specificity models were constructed for 169 human proteases using information on their known substrates. A previously developed method for structural analysis of potential proteolysis sites was applied in parallel with specificity models. Validation of the proposed approach was performed on the SARS-CoV-2 spike protein, whose proteolysis sites have been well studied.

IFITM3 rs12252 polymorphism association with COVID-19 severity and mortality in a Brazilian sample: an update and a meta-analysis.

This study investigated the association between the IFITM3 rs12252 polymorphism and the severity and mortality of COVID-19 in hospitalized Brazilian patients. A total of 102 COVID-19 patients were included, and the outcomes of interest were defined as death and the need for mechanical ventilation. Genotypes were assessed using Taqman probes. No significant associations were found between the rs12252 polymorphism and COVID-19 outcomes in the original sample, both for death and the need for mechanical ventilation. A meta-analysis, incorporating previous studies that used death as a severity indicator, revealed no association in the allelic and C-recessive models. However, due to the rarity of the T allele and its absence in the sample, further replication studies in larger and more diverse populations are needed to clarify the role of rs12252 in COVID-19 prognosis.

Untangling the Evolution of the Receptor-Binding Motif of SARS-CoV-2.

The spike protein determines the host-range specificity of coronaviruses. In particular, the Receptor-Binding Motif in the spike protein from SARS-CoV-2 contains the amino acids involved in molecular recognition of the host Angiotensin Converting Enzyme 2. Therefore, to understand how SARS-CoV-2 acquired its capacity to infect humans it is necessary to reconstruct the evolution of this important motif. Early during the pandemic, it was proposed that the SARS-CoV-2 Receptor-Binding Domain was acquired via recombination with a pangolin infecting coronavirus. This proposal was challenged by an alternative explanation that suggested that the Receptor-Binding Domain from SARS-CoV-2 did not originated via recombination with a coronavirus from a pangolin. Instead, this alternative hypothesis proposed that the Receptor-Binding Motif from the bat coronavirus RaTG13, was acquired via recombination with an unidentified coronavirus. And as a consequence of this event, the Receptor-Binding Domain from the pangolin coronavirus appeared as phylogenetically closer to SARS-CoV-2. Recently, the genomes from coronaviruses from Cambodia (bat_RShST182/200) and Laos (BANAL-20-52/103/247) which are closely related to SARS-CoV-2 were reported. However, no detailed analysis of the evolution of the Receptor-Binding Motif from these coronaviruses was reported. Here we revisit the evolution of the Receptor-Binding Domain and Motif in the light of the novel coronavirus genome sequences. Specifically, we wanted to test whether the above coronaviruses from Cambodia and Laos were the source of the Receptor-Binding Domain from RaTG13. We found that the Receptor-Binding Motif from these coronaviruses is phylogenetically closer to SARS-CoV-2 than to RaTG13. Therefore, the source of the Receptor-Binding Domain from RaTG13 is still unidentified. In accordance with previous studies, our results are consistent with the hypothesis that the Receptor-Binding Motif from SARS-CoV-2 evolved by vertical inheritance from a bat-infecting population of coronaviruses.

Molecular basis of hippopotamus ACE2 binding to SARS-CoV-2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a wide range of hosts, including hippopotami, which are semi-aquatic mammals and phylogenetically closely related to Cetacea. In this study, we characterized the binding properties of hippopotamus angiotensin-converting enzyme 2 (hiACE2) to the spike (S) protein receptor binding domains (RBDs) of the SARS-CoV-2 prototype (PT) and variants of concern (VOCs). Furthermore, the cryo-electron microscopy (cryo-EM) structure of the SARS-CoV-2 PT S protein complexed with hiACE2 was resolved. Structural and mutational analyses revealed that L30 and F83, which are specific to hiACE2, played a crucial role in the hiACE2/SARS-CoV-2 RBD interaction. In addition, comparative and structural analysis of ACE2 orthologs suggested that the cetaceans may have the potential to be infected by SARS-CoV-2. These results provide crucial molecular insights into the susceptibility of hippopotami to SARS-CoV-2 and suggest the potential risk of SARS-CoV-2 VOCs spillover and the necessity for surveillance. IMPORTANCE: The hippopotami are the first semi-aquatic artiodactyl mammals wherein SARS-CoV-2 infection has been reported. Exploration of the invasion mechanism of SARS-CoV-2 will provide important information for the surveillance of SARS-CoV-2 in hippopotami, as well as other semi-aquatic mammals and cetaceans. Here, we found that hippopotamus ACE2 (hiACE2) could efficiently bind to the RBDs of the SARS-CoV-2 prototype (PT) and variants of concern (VOCs) and facilitate the transduction of SARS-CoV-2 PT and VOCs pseudoviruses into hiACE2-expressing cells. The cryo-EM structure of the SARS-CoV-2 PT S protein complexed with hiACE2 elucidated a few critical residues in the RBD/hiACE2 interface, especially L30 and F83 of hiACE2 which are unique to hiACE2 and contributed to the decreased binding affinity to PT RBD compared to human ACE2. Our work provides insight into cross-species transmission and highlights the necessity for monitoring host jumps and spillover events on SARS-CoV-2 in semi-aquatic/aquatic mammals.

A positive selection at binding site 501 in the B.1 lineage might have triggered the highly infectious sub-lineages of SARS-CoV-2.

The descendants of the B lineage are the most predominant variants among the SARS-CoV-2 virus due to the incorporation of new mutations augmenting the infectivity of the virus. There is a substantial increase in the transition transversion bias, nucleotide diversity and purifying selection on the spike protein in the descendants of the B lineage of the SARS-CoV-2 virus on a temporal scale. A strong bias for C-to-U substitutions is found in the genes encoding spike protein in this lineage. The positive selection has operated on the spike gene of B lineages and its sub-lineages. The B.1 lineage has undergone positive selection on site 501 of the receptor binding domain ultimately reflected in a key substitution N501Y in its three descendant lineages namely B.1.1.7, B.1.351 and P.1. The intensity of purifying selection on the multiple sites of the spike gene has increased substantially in the sub-lineages of B.1 in a timescale. The binding site 501 on the spike protein in B lineage is found to coevolve with other amino acid sites. This study sheds light on the evolutionary trajectory of the B lineage into highly infectious descendants in the recent past under the influence of positive and purifying selection exerted by natural immunity and vaccination of the host.

Bioinformatics and system biology approaches to determine the connection of SARS-CoV-2 infection and intrahepatic cholangiocarcinoma.

INTRODUCTION: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causal agent of coronavirus disease 2019 (COVID-19), has infected millions of individuals worldwide, which poses a severe threat to human health. COVID-19 is a systemic ailment affecting various tissues and organs, including the lungs and liver. Intrahepatic cholangiocarcinoma (ICC) is one of the most common liver cancer, and cancer patients are particularly at high risk of SARS-CoV-2 infection. Nonetheless, few studies have investigated the impact of COVID-19 on ICC patients. METHODS: With the methods of systems biology and bioinformatics, this study explored the link between COVID-19 and ICC, and searched for potential therapeutic drugs. RESULTS: This study identified a total of 70 common differentially expressed genes (DEGs) shared by both diseases, shedding light on their shared functionalities. Enrichment analysis pinpointed metabolism and immunity as the primary areas influenced by these common genes. Subsequently, through protein-protein interaction (PPI) network analysis, we identified SCD, ACSL5, ACAT2, HSD17B4, ALDOA, ACSS1, ACADSB, CYP51A1, PSAT1, and HKDC1 as hub genes. Additionally, 44 transcription factors (TFs) and 112 microRNAs (miRNAs) were forecasted to regulate the hub genes. Most importantly, several drug candidates (Periodate-oxidized adenosine, Desipramine, Quercetin, Perfluoroheptanoic acid, Tetrandrine, Pentadecafluorooctanoic acid, Benzo[a]pyrene, SARIN, Dorzolamide, 8-Bromo-cAMP) may prove effective in treating ICC and COVID-19. CONCLUSION: This study is expected to provide valuable references and potential drugs for future research and treatment of COVID-19 and ICC.

ISMI-VAE: A deep learning model for classifying disease cells using gene expression and SNV data.

Various studies have linked several diseases, including cancer and COVID-19, to single nucleotide variations (SNV). Although single-cell RNA sequencing (scRNA-seq) technology can provide SNV and gene expression data, few studies have integrated and analyzed these multimodal data. To address this issue, we introduce Interpretable Single-cell Multimodal Data Integration Based on Variational Autoencoder (ISMI-VAE). ISMI-VAE leverages latent variable models that utilize the characteristics of SNV and gene expression data to overcome high noise levels and uses deep learning techniques to integrate multimodal information, map them to a low-dimensional space, and classify disease cells. Moreover, ISMI-VAE introduces an attention mechanism to reflect feature importance and analyze genetic features that could potentially cause disease. Experimental results on three cancer data sets and one COVID-19 data set demonstrate that ISMI-VAE surpasses the baseline method in terms of both effectiveness and interpretability and can effectively identify disease-causing gene features.

A measles-vectored vaccine candidate expressing prefusion-stabilized SARS-CoV-2 spike protein brought to phase I/II clinical trials: protection of African green monkeys from COVID-19 disease.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged at the end of 2019 and is responsible for the largest human pandemic in 100 years. Thirty-four vaccines are currently approved for use worldwide, and approximately 67% of the world population has received a complete primary series of one, yet countries are dealing with new waves of infections, variant viruses continue to emerge, and breakthrough infections are frequent secondary to waning immunity. Here, we evaluate a measles virus (MV)-vectored vaccine expressing a stabilized prefusion SARS-CoV-2 spike (S) protein (MV-ATU3-S2PΔF2A; V591) with demonstrated immunogenicity in mouse models (see companion article [J. Brunet, Z. Choucha, M. Gransagne, H. Tabbal, M.-W. Ku et al., J Virol 98:e01693-23, 2024, https://doi.org/10.1128/jvi.01693-23]) in an established African green monkey model of disease. Animals were vaccinated with V591 or the control vaccine (an equivalent MV-vectored vaccine with an irrelevant antigen) intramuscularly using a prime/boost schedule, followed by challenge with an early pandemic isolate of SARS-CoV-2 at 56 days post-vaccination. Pre-challenge, only V591-vaccinated animals developed S-specific antibodies that had virus-neutralizing activity as well as S-specific T cells. Following the challenge, V591-vaccinated animals had lower infectious virus and viral (v) RNA loads in mucosal secretions and stopped shedding virus in these secretions earlier. vRNA loads were lower in these animals in respiratory and gastrointestinal tract tissues at necropsy. This correlated with a lower disease burden in the lungs as quantified by PET/CT at early and late time points post-challenge and by pathological analysis at necropsy.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the largest human pandemic in 100 years. Even though vaccines are currently available, countries are dealing with new waves of infections, variant viruses continue to emerge, breakthrough infections are frequent, and vaccine hesitancy persists. This study uses a safe and effective measles vaccine as a platform for vaccination against SARS-CoV-2. The candidate vaccine was used to vaccinate African green monkeys (AGMs). All vaccinated AGMs developed robust antigen-specific immune responses. After challenge, these AGMs produced less virus in mucosal secretions, for a shorter period, and had a reduced disease burden in the lungs compared to control animals. At necropsy, lower levels of viral RNA were detected in tissue samples from vaccinated animals, and the lungs of these animals lacked the histologic hallmarks of SARS-CoV-2 disease observed exclusively in the control AGMs.

Identification of a membrane-associated element (MAE) in the C-terminal region of SARS-CoV-2 nsp6 that is essential for viral replication.

The coronavirus disease 2019 (COVID-19) pandemic, caused by the novel coronavirus severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), has rapidly spread worldwide since its emergence in late 2019. Its ongoing evolution poses challenges for antiviral drug development. Coronavirus nsp6, a multiple-spanning transmembrane protein, participates in the biogenesis of the viral replication complex, which accommodates the viral replication-transcription complex. The roles of its structural domains in viral replication are not well studied. Herein, we predicted the structure of the SARS-CoV-2 nsp6 protein using AlphaFold2 and identified a highly folded C-terminal region (nsp6C) downstream of the transmembrane helices. The enhanced green fluorescent protein (EGFP)-fused nsp6C was found to cluster in the cytoplasm and associate with membranes. Functional mapping identified a minimal membrane-associated element (MAE) as the region from amino acids 237 to 276 (LGV-KLL), which is mainly composed of the α-helix H1 and the α-helix H2; the latter exhibits characteristics of an amphipathic helix (AH). Mutagenesis studies and membrane flotation experiments demonstrate that AH-like H2 is required for MAE-mediated membrane association. This MAE was functionally conserved across MERS-CoV, HCoV-OC43, HCoV-229E, HCoV-HKU1, and HCoV-NL63, all capable of mediating membrane association. In a SARS-CoV-2 replicon system, mutagenesis studies of H2 and replacements of H1 and H2 with their homologous counterparts demonstrated requirements of residues on both sides of the H2 and properly paired H1-H2 for MAE-mediated membrane association and viral replication. Notably, mutations I266A and K274A significantly attenuated viral replication without dramatically affecting membrane association, suggesting a dual role of the MAE in viral replication: mediating membrane association as well as participating in protein-protein interactions.IMPORTANCESevere acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) assembles a double-membrane vesicle (DMV) by the viral non-structural proteins for viral replication. Understanding the mechanisms of the DMV assembly is of paramount importance for antiviral development. Nsp6, a multiple-spanning transmembrane protein, plays an important role in the DMV biogenesis. Herein, we predicted the nsp6 structure of SARS-CoV-2 and other human coronaviruses using AlphaFold2 and identified a putative membrane-associated element (MAE) in the highly conserved C-terminal regions of nsp6. Experimentally, we verified a functionally conserved minimal MAE composed of two α-helices, the H1, and the amphipathic helix-like H2. Mutagenesis studies confirmed the requirement of H2 for MAE-mediated membrane association and viral replication and demonstrated a dual role of the MAE in viral replication, by mediating membrane association and participating in residue-specific interactions. This functionally conserved MAE may serve as a novel anti-viral target.

Conserved structures and dynamics in 5'-proximal regions of Betacoronavirus RNA genomes.

Betacoronaviruses are a genus within the Coronaviridae family of RNA viruses. They are capable of infecting vertebrates and causing epidemics as well as global pandemics in humans. Mitigating the threat posed by Betacoronaviruses requires an understanding of their molecular diversity. The development of novel antivirals hinges on understanding the key regulatory elements within the viral RNA genomes, in particular the 5'-proximal region, which is pivotal for viral protein synthesis. Using a combination of cryo-electron microscopy, atomic force microscopy, chemical probing, and computational modeling, we determined the structures of 5'-proximal regions in RNA genomes of Betacoronaviruses from four subgenera: OC43-CoV, SARS-CoV-2, MERS-CoV, and Rousettus bat-CoV. We obtained cryo-electron microscopy maps and determined atomic-resolution models for the stem-loop-5 (SL5) region at the translation start site and found that despite low sequence similarity and variable length of the helical elements it exhibits a remarkable structural conservation. Atomic force microscopy imaging revealed a common domain organization and a dynamic arrangement of structural elements connected with flexible linkers across all four Betacoronavirus subgenera. Together, these results reveal common features of a critical regulatory region shared between different Betacoronavirus RNA genomes, which may allow targeting of these RNAs by broad-spectrum antiviral therapeutics.

DinoKnot: Duplex Interaction of Nucleic Acids With PseudoKnots.

Interaction of nucleic acid molecules is essential for their functional roles in the cell and their applications in biotechnology. While simple duplex interactions have been studied before, the problem of efficiently predicting the minimum free energy structure of more complex interactions with possibly pseudoknotted structures remains a challenge. In this work, we introduce a novel and efficient algorithm for prediction of Duplex Interaction of Nucleic acids with pseudoKnots, DinoKnot follows the hierarchical folding hypothesis to predict the secondary structure of two interacting nucleic acid strands (both homo- and hetero-dimers). DinoKnot utilizes the structure of molecules before interaction as a guide to find their duplex structure allowing for possible base pair competitions. To showcase DinoKnots's capabilities we evaluated its predicted structures against (1) experimental results for SARS-CoV-2 genome and nine primer-probe sets, (2) a clinically verified example of a mutation affecting detection, and (3) a known nucleic acid interaction involving a pseudoknot. In addition, we compared our results against our closest competition, RNAcofold, further highlighting DinoKnot's strengths. We believe DinoKnot can be utilized for various applications including screening new variants for potential detection issues and supporting existing applications involving DNA/RNA interactions, adding structural considerations to the interaction to elicit functional information.

Global patterns of phylogenetic diversity and transmission of bat coronavirus.

Bats are reservoirs for multiple coronaviruses (CoVs). However, the phylogenetic diversity and transmission of global bat-borne CoVs remain poorly understood. Here, we performed a Bayesian phylogeographic analysis based on 3,594 bat CoV RdRp gene sequences to study the phylogenetic diversity and transmission of bat-borne CoVs and the underlying driving factors. We found that host-switching events occurred more frequently for α-CoVs than for ß-CoVs, and the latter was highly constrained by bat phylogeny. Bat species in the families Molossidae, Rhinolophidae, Miniopteridae, and Vespertilionidae had larger contributions to the cross-species transmission of bat CoVs. Regions of eastern and southern Africa, southern South America, Western Europe, and Southeast Asia were more frequently involved in cross-region transmission events of bat CoVs than other regions. Phylogenetic and geographic distances were the most important factors limiting CoV transmission. Bat taxa and global geographic hotspots associated with bat CoV phylogenetic diversity were identified, and bat species richness, mean annual temperature, global agricultural cropland, and human population density were strongly correlated with the phylogenetic diversity of bat CoVs. These findings provide insight into bat CoV evolution and ecological transmission among bat taxa. The identified hotspots of bat CoV evolution and transmission will guide early warnings of bat-borne CoV zoonotic diseases.

Human coronaviruses: The emergence of SARS-CoV-2 and management of COVID-19.

To date, a total of seven human coronaviruses (HCoVs) have been identified, all of which are important respiratory pathogens. Recently, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19), has led to a global pandemic causing millions of infections and deaths. Here, we summarize the discovery and fundamental virology of HCoVs, discuss their zoonotic transmission and highlight the weak species barrier of SARS-CoV-2. We also discuss the possible origins of SARS-CoV-2 variants of concern identified to date and discuss the experimental challenges in characterizing mutations of interest and propose methods to circumvent them. As the COVID-19 treatment and prevention landscape rapidly evolves, we summarize current therapeutics and vaccines, and their implications on SARS-CoV-2 variants. Finally, we explore how interspecies transmission of SARS-CoV-2 may drive the emergence of novel strains, how disease severity may evolve and how COVID-19 will likely continue to burden healthcare systems globally.

Genomic Comparisons of Alphacoronaviruses and Betacoronaviruses from Korean Bats.

Coronaviruses are well known as a diverse family of viruses that affect a wide range of hosts. Since the outbreak of severe acute respiratory syndrome, a variety of bat-associated coronaviruses have been identified in many countries. However, they do not represent all the specific geographic locations of their hosts. In this study, full-length genomes representing newly identified bat coronaviruses in South Korea were obtained using an RNA sequencing approach. The analysis, based on genome structure, conserved replicase domains, spike gene, and nucleocapsid genes revealed that bat Alphacoronaviruses are from three different viral species. Among them, the newly identified B20-97 strain may represent a new putative species, closely related to PEDV. In addition, the newly-identified MERS-related coronavirus exhibited shared genomic nucleotide identities of less than 76.4% with other Merbecoviruses. Recombination analysis and multiple alignments of spike and RBD amino acid sequences suggested that this strain underwent recombination events and could possibly use hDPP4 molecules as its receptor. The bat SARS-related CoV B20-50 is unlikely to be able to use hACE2 as its receptor and lack of an open reading frame in ORF8 gene region. Our results illustrate the diversity of coronaviruses in Korean bats and their evolutionary relationships. The evolution of the bat coronaviruses related ORF8 accessory gene is also discussed.

Analysis of SARS-CoV-2 synonymous codon usage evolution throughout the COVID-19 pandemic.

SARS-CoV-2, the seventh coronavirus known to infect humans, can cause severe life-threatening respiratory pathologies. To better understand SARS-CoV-2 evolution, genome-wide analyses have been made, including the general characterization of its codons usage profile. Here we present a bioinformatic analysis of the evolution of SARS-CoV-2 codon usage over time using complete genomes collected since December 2019. Our results show that SARS-CoV-2 codon usage pattern is antagonistic to, and it is getting farther away from that of the human host. Further, a selection of deoptimized codons over time, which was accompanied by a decrease in both the codon adaptation index and the effective number of codons, was observed. All together, these findings suggest that SARS-CoV-2 could be evolving, at least from the perspective of the synonymous codon usage, to become less pathogenic.

First detection of Hedgehog coronavirus 1 in Poland.

Hedgehogs are common in the majority of European countries and are known to host various pathogens, including viruses. The recent discovery of MERS-related coronaviruses (CoVs) in hedgehogs from Germany, France, the UK, China, and Italy suggests that hedgehogs may represent a wild reservoir of betacoronaviruses. This study reports the first detection and characterization of novel betacoronovirus, subgenus Merbecovirus in wild hedgehogs in Poland. The CoV RNA was detected in 10 out of 40 hedgehogs' rectal swabs and in 1 out of 18 samples of the lung. No viral RNA was identified in the duodenum and kidney. There was no significant relationship between clinical status, gender, hedgehogs' age, and coronaviral RNA detection. Phylogenetic analysis showed that CoVs detected in our study grouped together with other representatives of Hedgehog coronavirus 1 species identified in Western Europe. Our findings provide further evidence that hedgehogs are a natural reservoir of Merbecovirus. Considering the high mutation rate of CoVs and their potential for crossing interspecies barriers, the proper management of hedgehogs admitted to wildlife rehabilitation centres is needed. It cannot be excluded that merbecovirus strains detected in hedgehogs may recombine with other CoVs leading to new viruses with potential for interspecies transmission.

The unresolved question on COVID-19 virus origin: The three cards game?

The ongoing discussion about the real origin of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) feeds acrimonious debates. Where did SARS-CoV-2 come from? Was SARS-CoV-2 transmitted in the wild from an animal to a person before exploding in Wuhan or was it an engineered virus that escaped from research or a laboratory in Wuhan? Right now, we still don't know enough whether SARS-CoV-2 is human-made or not, and lab-leak theories remain essentially speculative. Many recent studies have pointed out several plausible scenarios. Anyhow, currently, even if suspicions by some about the possibility of lab-leak hypothesis still remain, the consensus view is that the pandemic probably started from a natural source and, to determine the real origin of the SARS-CoV-2 virus, further research is needed.

LinearTurboFold: Linear-time global prediction of conserved structures for RNA homologs with applications to SARS-CoV-2.

The constant emergence of COVID-19 variants reduces the effectiveness of existing vaccines and test kits. Therefore, it is critical to identify conserved structures in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes as potential targets for variant-proof diagnostics and therapeutics. However, the algorithms to predict these conserved structures, which simultaneously fold and align multiple RNA homologs, scale at best cubically with sequence length and are thus infeasible for coronaviruses, which possess the longest genomes (∼30,000 nt) among RNA viruses. As a result, existing efforts on modeling SARS-CoV-2 structures resort to single-sequence folding as well as local folding methods with short window sizes, which inevitably neglect long-range interactions that are crucial in RNA functions. Here we present LinearTurboFold, an efficient algorithm for folding RNA homologs that scales linearly with sequence length, enabling unprecedented global structural analysis on SARS-CoV-2. Surprisingly, on a group of SARS-CoV-2 and SARS-related genomes, LinearTurboFold's purely in silico prediction not only is close to experimentally guided models for local structures, but also goes far beyond them by capturing the end-to-end pairs between 5' and 3' untranslated regions (UTRs) (∼29,800 nt apart) that match perfectly with a purely experimental work. Furthermore, LinearTurboFold identifies undiscovered conserved structures and conserved accessible regions as potential targets for designing efficient and mutation-insensitive small-molecule drugs, antisense oligonucleotides, small interfering RNAs (siRNAs), CRISPR-Cas13 guide RNAs, and RT-PCR primers. LinearTurboFold is a general technique that can also be applied to other RNA viruses and full-length genome studies and will be a useful tool in fighting the current and future pandemics.