SARS-CoV-2-related bat viruses evade human intrinsic immunity but lack efficient transmission capacity.

Circulating bat coronaviruses represent a pandemic threat. However, our understanding of bat coronavirus pathogenesis and transmission potential is limited by the lack of phenotypically characterized strains. We created molecular clones for the two closest known relatives of SARS-CoV-2, BANAL-52 and BANAL-236. We demonstrated that BANAL-CoVs and SARS-CoV-2 have similar replication kinetics in human bronchial epithelial cells. However, BANAL-CoVs have impaired replication in human nasal epithelial cells and in the upper airway of mice. We also observed reduced pathogenesis in mice and diminished transmission in hamsters. Further, we observed that diverse bat coronaviruses evade interferon and downregulate major histocompatibility complex class I. Collectively, our study demonstrates that despite high genetic similarity across bat coronaviruses, prediction of pandemic potential of a virus necessitates functional characterization. Finally, the restriction of bat coronavirus replication in the upper airway highlights that transmission potential and innate immune restriction can be uncoupled in this high-risk family of emerging viruses.

The West Nile virus genome harbors essential riboregulatory elements with conserved and host-specific functional roles.

West Nile virus (WNV) is an arthropod-borne, positive-sense RNA virus that poses an increasing global threat due to warming climates and lack of effective therapeutics. Like other enzootic viruses, little is known about how host context affects the structure of the full-length RNA genome. Here, we report a complete secondary structure of the entire WNV genome within infected mammalian and arthropod cell lines. Our analysis affords structural insights into multiple, conserved aspects of flaviviral biology. We show that the WNV genome folds with minimal host dependence, and we prioritize well-folded regions for functional validation using structural homology between hosts as a guide. Using structure-disrupting, antisense locked nucleic acids, we then demonstrate that the WNV genome contains riboregulatory structures with conserved and host-specific functional roles. These results reveal promising RNA drug targets within flaviviral genomes, and they highlight the therapeutic potential of ASO-LNAs as both WNV-specific and pan-flaviviral therapeutic agents.

The SMN complex drives structural changes in human snRNAs to enable snRNP assembly.

Spliceosomal snRNPs are multicomponent particles that undergo a complex maturation pathway. Human Sm-class snRNAs are generated as 3'-end extended precursors, which are exported to the cytoplasm and assembled together with Sm proteins into core RNPs by the SMN complex. Here, we provide evidence that these pre-snRNA substrates contain compact, evolutionarily conserved secondary structures that overlap with the Sm binding site. These structural motifs in pre-snRNAs are predicted to interfere with Sm core assembly. We model structural rearrangements that lead to an open pre-snRNA conformation compatible with Sm protein interaction. The predicted rearrangement pathway is conserved in Metazoa and requires an external factor that initiates snRNA remodeling. We show that the essential helicase Gemin3, which is a component of the SMN complex, is crucial for snRNA structural rearrangements during snRNP maturation. The SMN complex thus facilitates ATP-driven structural changes in snRNAs that expose the Sm site and enable Sm protein binding.

Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium identifies target cells, alterations in gene expression, and cell state changes.

There are currently limited Food and Drug Administration (FDA)-approved drugs and vaccines for the treatment or prevention of Coronavirus Disease 2019 (COVID-19). Enhanced understanding of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection and pathogenesis is critical for the development of therapeutics. To provide insight into viral replication, cell tropism, and host-viral interactions of SARS-CoV-2, we performed single-cell (sc) RNA sequencing (RNA-seq) of experimentally infected human bronchial epithelial cells (HBECs) in air-liquid interface (ALI) cultures over a time course. This revealed novel polyadenylated viral transcripts and highlighted ciliated cells as a major target at the onset of infection, which we confirmed by electron and immunofluorescence microscopy. Over the course of infection, the cell tropism of SARS-CoV-2 expands to other epithelial cell types including basal and club cells. Infection induces cell-intrinsic expression of type I and type III interferons (IFNs) and interleukin (IL)-6 but not IL-1. This results in expression of interferon-stimulated genes (ISGs) in both infected and bystander cells. This provides a detailed characterization of genes, cell types, and cell state changes associated with SARS-CoV-2 infection in the human airway.

Comprehensive in vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms.

Severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2) is the positive-sense RNA virus that causes coronavirus disease 2019 (COVID-19). The genome of SARS-CoV-2 is unique among viral RNAs in its vast potential to form RNA structures, yet as much as 97% of its 30 kilobases have not been structurally explored. Here, we apply a novel long amplicon strategy to determine the secondary structure of the SARS-CoV-2 RNA genome at single-nucleotide resolution in infected cells. Our in-depth structural analysis reveals networks of well-folded RNA structures throughout Orf1ab and reveals aspects of SARS-CoV-2 genome architecture that distinguish it from other RNA viruses. Evolutionary analysis shows that several features of the SARS-CoV-2 genomic structure are conserved across ß-coronaviruses, and we pinpoint regions of well-folded RNA structure that merit downstream functional analysis. The native, secondary structure of SARS-CoV-2 presented here is a roadmap that will facilitate focused studies on the viral life cycle, facilitate primer design, and guide the identification of RNA drug targets against COVID-19.

The global and local distribution of RNA structure throughout the SARS-CoV-2 genome.

SARS-CoV-2 is the causative viral agent of COVID-19, the disease at the center of the current global pandemic. While knowledge of highly structured regions is integral for mechanistic insights into the viral infection cycle, very little is known about the location and folding stability of functional elements within the massive, ∼30kb SARS-CoV-2 RNA genome. In this study, we analyze the folding stability of this RNA genome relative to the structural landscape of other well-known viral RNAs. We present an in-silico pipeline to predict regions of high base pair content across long genomes and to pinpoint hotspots of well-defined RNA structures, a method that allows for direct comparisons of RNA structural complexity within the several domains in SARS-CoV-2 genome. We report that the SARS-CoV-2 genomic propensity for stable RNA folding is exceptional among RNA viruses, superseding even that of HCV, one of the most structured viral RNAs in nature. Furthermore, our analysis suggests varying levels of RNA structure across genomic functional regions, with accessory and structural ORFs containing the highest structural density in the viral genome. Finally, we take a step further to examine how individual RNA structures formed by these ORFs are affected by the differences in genomic and subgenomic contexts, which given the technical difficulty of experimentally separating cellular mixtures of sgRNA from gRNA, is a unique advantage of our in-silico pipeline. The resulting findings provide a useful roadmap for planning focused empirical studies of SARS-CoV-2 RNA biology, and a preliminary guide for exploring potential SARS-CoV-2 RNA drug targets.Importance The RNA genome of SARS-CoV-2 is among the largest and most complex viral genomes, and yet its RNA structural features remain relatively unexplored. Since RNA elements guide function in most RNA viruses, and they represent potential drug targets, it is essential to chart the architectural features of SARS-CoV-2 and pinpoint regions that merit focused study. Here we show that RNA folding stability of SARS-CoV-2 genome is exceptional among viral genomes and we develop a method to directly compare levels of predicted secondary structure across SARS-CoV-2 domains. Remarkably, we find that coding regions display the highest structural propensity in the genome, forming motifs that differ between the genomic and subgenomic contexts. Our approach provides an attractive strategy to rapidly screen for candidate structured regions based on base pairing potential and provides a readily interpretable roadmap to guide functional studies of RNA viruses and other pharmacologically relevant RNA transcripts.

Comprehensive in-vivo secondary structure of the SARS-CoV-2 genome reveals novel regulatory motifs and mechanisms.

SARS-CoV-2 is the positive-sense RNA virus that causes COVID-19, a disease that has triggered a major human health and economic crisis. The genome of SARS-CoV-2 is unique among viral RNAs in its vast potential to form stable RNA structures and yet, as much as 97% of its 30 kilobases have not been structurally explored in the context of a viral infection. Our limited knowledge of SARS-CoV-2 genomic architecture is a fundamental limitation to both our mechanistic understanding of coronavirus life cycle and the development of COVID-19 RNA-based therapeutics. Here, we apply a novel long amplicon strategy to determine for the first time the secondary structure of the SARS-CoV-2 RNA genome probed in infected cells. In addition to the conserved structural motifs at the viral termini, we report new structural features like a conformationally flexible programmed ribosomal frameshifting pseudoknot, and a host of novel RNA structures, each of which highlights the importance of studying viral structures in their native genomic context. Our in-depth structural analysis reveals extensive networks of well-folded RNA structures throughout Orf1ab and reveals new aspects of SARS-CoV-2 genome architecture that distinguish it from other single-stranded, positive-sense RNA viruses. Evolutionary analysis of RNA structures in SARS-CoV-2 shows that several features of its genomic structure are conserved across beta coronaviruses and we pinpoint individual regions of well-folded RNA structure that merit downstream functional analysis. The native, complete secondary structure of SAR-CoV-2 presented here is a roadmap that will facilitate focused studies on mechanisms of replication, translation and packaging, and guide the identification of new RNA drug targets against COVID-19.

Mutations in the iron-sulfur cluster biogenesis protein HSCB cause congenital sideroblastic anemia.

The congenital sideroblastic anemias (CSAs) can be caused by primary defects in mitochondrial iron-sulfur (Fe-S) cluster biogenesis. HSCB (heat shock cognate B), which encodes a mitochondrial cochaperone, also known as HSC20 (heat shock cognate protein 20), is the partner of mitochondrial heat shock protein A9 (HSPA9). Together with glutaredoxin 5 (GLRX5), HSCB and HSPA9 facilitate the transfer of nascent 2-iron, 2-sulfur clusters to recipient mitochondrial proteins. Mutations in both HSPA9 and GLRX5 have previously been associated with CSA. Therefore, we hypothesized that mutations in HSCB could also cause CSA. We screened patients with genetically undefined CSA and identified a frameshift mutation and a rare promoter variant in HSCB in a female patient with non-syndromic CSA. We found that HSCB expression was decreased in patient-derived fibroblasts and K562 erythroleukemia cells engineered to have the patient-specific promoter variant. Furthermore, gene knockdown and deletion experiments performed in K562 cells, zebrafish, and mice demonstrate that loss of HSCB results in impaired Fe-S cluster biogenesis, a defect in RBC hemoglobinization, and the development of siderocytes and more broadly perturbs hematopoiesis in vivo. These results further affirm the involvement of Fe-S cluster biogenesis in erythropoiesis and hematopoiesis and define HSCB as a CSA gene.

Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium.

SARS-CoV-2, the causative agent of COVID-19, has tragically burdened individuals and institutions around the world. There are currently no approved drugs or vaccines for the treatment or prevention of COVID-19. Enhanced understanding of SARS-CoV-2 infection and pathogenesis is critical for the development of therapeutics. To reveal insight into viral replication, cell tropism, and host-viral interactions of SARS-CoV-2 we performed single-cell RNA sequencing of experimentally infected human bronchial epithelial cells (HBECs) in air-liquid interface cultures over a time-course. This revealed novel polyadenylated viral transcripts and highlighted ciliated cells as a major target of infection, which we confirmed by electron microscopy. Over the course of infection, cell tropism of SARS-CoV-2 expands to other epithelial cell types including basal and club cells. Infection induces cell-intrinsic expression of type I and type III IFNs and IL6 but not IL1. This results in expression of interferon-stimulated genes in both infected and bystander cells. We observe similar gene expression changes from a COVID-19 patient ex vivo. In addition, we developed a new computational method termed CONditional DENSity Embedding (CONDENSE) to characterize and compare temporal gene dynamics in response to infection, which revealed genes relating to endothelin, angio-genesis, interferon, and inflammation-causing signaling pathways. In this study, we conducted an in-depth analysis of SARS-CoV-2 infection in HBECs and a COVID-19 patient and revealed genes, cell types, and cell state changes associated with infection.

Sequencing and Structure Probing of Long RNAs Using MarathonRT: A Next-Generation Reverse Transcriptase.

Reverse transcriptase (RT) enzymes are indispensable tools for interrogating diverse aspects of RNA metabolism and transcriptome composition. Due to the growing interest in sequence and structural complexity of long RNA molecules, processive RT enzymes are now required for preserving linkage and information content in mixed populations of transcripts, and the low-processivity RT enzymes that are commercially available cannot meet this need. MarathonRT is encoded within a eubacterial group II intron, and it has been shown to efficiently copy highly structured long RNA molecules in a single pass. In this work, we systematically characterize MarathonRT as a tool enzyme and optimize its performance in a variety of applications that include single-cycle reverse transcription of long RNAs, dimethyl sulfate mutational profiling (DMS-MaP), selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP), using ultra-long amplicons and the detection of natural RNA base modifications. By diversifying MarathonRT reaction protocols, we provide an upgraded suite of tools for cutting-edge RNA research and clinical application.