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.

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

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.

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.

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.

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.

Identification of Guide-Intrinsic Determinants of Cas9 Specificity

Considerable effort has been devoted to developing a comprehensive understanding of CRISPR nuclease specificity. In silico predictions and multiple genome-wide cellular and biochemical approaches have revealed a basic understanding of the Cas9 specificity profile. However, none of these approaches has delivered a model that allows accurate prediction of a CRISPR nuclease's ability to cleave a site based entirely on the sequence of the guide RNA (gRNA) and the target. We describe a library-based biochemical assay that directly reports the cleavage efficiency of a particular Cas9-guide complex by measuring both uncleaved and cleaved target molecules over a wide range of mismatched library members. We applied our assay using libraries of targets to evaluate the specificity of Staphylococcus aureus Cas9 under a variety of experimental conditions. Surprisingly, our data show an unexpectedly high variation in the random gRNA:target DNA mismatch tolerance when cleaving with different gRNAs, indicating guide-intrinsic mismatch permissiveness and challenging the assumption of universal specificity models. We use data generated by our assay to create the first off-target, guide-specific cleavage models. The barcoded libraries of targets approach is rapid, highly modular, and capable of generating protein- and guide-specific models, as well as illuminating the biophysics of Cas9 binding versus cutting. These models may be useful in identifying potential off-targets, and the gRNA-intrinsic nature of mismatch tolerance argues for coupling these specificity models with orthogonal methods for a more complete assessment of gRNA specificity.

Sensitive detection of structural features and rearrangements in long, structured RNA molecules

Technical innovations in structural probing have drastically advanced the field of RNA structure analysis. These advances have led to parallel approaches developed in separate labs for analyzing RNA structure and dynamics. With the wealth of methodologies available, it can be difficult to determine which is best suited for a given application. Here, using a long, highly structured viral RNA as an example (the positive strand genome of Hepatitis C Virus), we present a semi-comprehensive analysis and describe the major approaches for analyzing the architecture of RNA that is modified with structure-sensitive probes. Additionally, we present an updated method for generating in vitro transcribed and folded RNA that maintains native secondary structures in long RNA molecules. We anticipate that the methods described here will streamline the use of current approaches and help investigators who are unfamiliar with structure probing, obviating the need for time-consuming and expensive optimization.

Pairwise library screen systematically interrogates Staphylococcus aureus Cas9 specificity in human cells

Therapeutic genome editing with Staphylococcus aureus Cas9 (SaCas9) requires a rigorous understanding of its potential off-target activity in the human genome. Here we report a high-throughput screening approach to measure SaCas9 genome editing variation in human cells across a large repertoire of 88,692 single guide RNAs (sgRNAs) paired with matched or mismatched target sites in a synthetic cassette. We incorporate randomized barcodes that enable whitelisting of correctly synthesized molecules for further downstream analysis, in order to circumvent the limitation of oligonucleotide synthesis errors. We find SaCas9 sgRNAs with 21-mer or 22-mer spacer sequences are generally more active, although high efficiency 20-mer spacers are markedly less tolerant of mismatches. Using this dataset, we developed an SaCas9 specificity model that performs robustly in ranking off-target sites. The barcoded pairwise library screen enabled high-fidelity recovery of guide-target relationships, providing a scalable framework for the investigation of CRISPR enzyme properties and general nucleic acid interactions.

Reductions in the mitochondrial ABC transporter Abcb10 affect the transcriptional profile of heme biosynthesis genes

ATP-binding cassette subfamily B member 10 (Abcb10) is a mitochondrial ATP-binding cassette (ABC) transporter that complexes with mitoferrin1 and ferrochelatase to enhance heme biosynthesis in developing red blood cells. Reductions in Abcb10 levels have been shown to reduce mitoferrin1 protein levels and iron import into mitochondria, resulting in reduced heme biosynthesis. As an ABC transporter, Abcb10 binds and hydrolyzes ATP, but its transported substrate is unknown. Here, we determined that decreases in Abcb10 did not result in protoporphyrin IX accumulation in morphant-treated zebrafish embryos or in differentiated Abcb10-specific shRNA murine Friend erythroleukemia (MEL) cells in which Abcb10 was specifically silenced with shRNA. We also found that the ATPase activity of Abcb10 is necessary for hemoglobinization in MEL cells, suggesting that the substrate transported by Abcb10 is important in mediating increased heme biosynthesis during erythroid development. Inhibition of 5-aminolevulinic acid dehydratase (EC 4.2.1.24) with succinylacetone resulted in both 5-aminolevulinic acid (ALA) accumulation in control and Abcb10-specific shRNA MEL cells, demonstrating that reductions in Abcb10 do not affect ALA export from mitochondria and indicating that Abcb10 does not transport ALA. Abcb10 silencing resulted in an alteration in the heme biosynthesis transcriptional profile due to repression by the transcriptional regulator Bach1, which could be partially rescued by overexpression of Alas2 or Gata1, providing a mechanistic explanation for why Abcb10 shRNA MEL cells exhibit reduced hemoglobinization. In conclusion, our findings rule out that Abcb10 transports ALA and indicate that Abcb10's ATP-hydrolysis activity is critical for hemoglobinization and that the substrate transported by Abcb10 provides a signal that optimizes hemoglobinization.

Restored iron transport by a small molecule promotes absorption and hemoglobinization in animals

Multiple human diseases ensue from a hereditary or acquired deficiency of iron-transporting protein function that diminishes transmembrane iron flux in distinct sites and directions. Because other iron-transport proteins remain active, labile iron gradients build up across the corresponding protein-deficient membranes. Here we report that a small-molecule natural product, hinokitiol, can harness such gradients to restore iron transport into, within, and/or out of cells. The same compound promotes gut iron absorption in DMT1-deficient rats and ferroportin-deficient mice, as well as hemoglobinization in DMT1- and mitoferrin-deficient zebrafish. These findings illuminate a general mechanistic framework for small molecule-mediated site- and direction-selective restoration of iron transport. They also suggest that small molecules that partially mimic the function of missing protein transporters of iron, and possibly other ions, may have potential in treating human diseases.

Mitochondrial ClpX Activates a Key Enzyme for Heme Biosynthesis and Erythropoiesis

The mitochondrion maintains and regulates its proteome with chaperones primarily inherited from its bacterial endosymbiont ancestor. Among these chaperones is the AAA+ unfoldase ClpX, an important regulator of prokaryotic physiology with poorly defined function in the eukaryotic mitochondrion. We observed phenotypic similarity in S. cerevisiae genetic interaction data between mitochondrial ClpX (mtClpX) and genes contributing to heme biosynthesis, an essential mitochondrial function. Metabolomic analysis revealed that 5-aminolevulinic acid (ALA), the first heme precursor, is 5-fold reduced in yeast lacking mtClpX activity and that total heme is reduced by half. mtClpX directly stimulates ALA synthase in vitro by catalyzing incorporation of its cofactor, pyridoxal phosphate. This activity is conserved in mammalian homologs; additionally, mtClpX depletion impairs vertebrate erythropoiesis, which requires massive upregulation of heme biosynthesis to supply hemoglobin. mtClpX, therefore, is a widely conserved stimulator of an essential biosynthetic pathway and uses a previously unrecognized mechanism for AAA+ unfoldases.

The mTORC1/4E-BP pathway coordinates hemoglobin production with L-leucine availability

In multicellular organisms, the mechanisms by which diverse cell types acquire distinct amino acids and how cellular function adapts to their availability are fundamental questions in biology. We found that increased neutral essential amino acid (NEAA) uptake was a critical component of erythropoiesis. As red blood cells matured, expression of the amino acid transporter gene Lat3 increased, which increased NEAA import. Inadequate NEAA uptake by pharmacologic inhibition or RNAi-mediated knockdown of LAT3 triggered a specific reduction in hemoglobin production in zebrafish embryos and murine erythroid cells through the mTORC1 (mammalian target of rapamycin complex 1)/4E-BP (eukaryotic translation initiation factor 4E-binding protein) pathway. CRISPR-mediated deletion of members of the 4E-BP family in murine erythroid cells rendered them resistant to mTORC1 and LAT3 inhibition and restored hemoglobin production. These results identify a developmental role for LAT3 in red blood cells and demonstrate that mTORC1 serves as a homeostatic sensor that couples hemoglobin production at the translational level to sufficient uptake of NEAAs, particularly L-leucine.

Iron regulatory protein-1 protects against mitoferrin-1-deficient porphyria

Mitochondrial iron is essential for the biosynthesis of heme and iron-sulfur ([Fe-S]) clusters in mammalian cells. In developing erythrocytes, iron is imported into the mitochondria by MFRN1 (mitoferrin-1, SLC25A37). Although loss of MFRN1 in zebrafish and mice leads to profound anemia, mutant animals showed no overt signs of porphyria, suggesting that mitochondrial iron deficiency does not result in an accumulation of protoporphyrins. Here, we developed a gene trap model to provide in vitro and in vivo evidence that iron regulatory protein-1 (IRP1) inhibits protoporphyrin accumulation. Mfrn1(+/gt);Irp1(-/-) erythroid cells exhibit a significant increase in protoporphyrin levels. IRP1 attenuates protoporphyrin biosynthesis by binding to the 5'-iron response element (IRE) of alas2 mRNA, inhibiting its translation. Ectopic expression of alas2 harboring a mutant IRE, preventing IRP1 binding, in Mfrn1(gt/gt) cells mimics Irp1 deficiency. Together, our data support a model whereby impaired mitochondrial [Fe-S] cluster biogenesis in Mfrn1(gt/gt) cells results in elevated IRP1 RNA-binding that attenuates ALAS2 mRNA translation and protoporphyrin accumulation.

Macrocytic anemia and mitochondriopathy resulting from a defect in sideroflexin 4

We used exome sequencing to identify mutations in sideroflexin 4 (SFXN4) in two children with mitochondrial disease (the more severe case also presented with macrocytic anemia). SFXN4 is an uncharacterized mitochondrial protein that localizes to the mitochondrial inner membrane. sfxn4 knockdown in zebrafish recapitulated the mitochondrial respiratory defect observed in both individuals and the macrocytic anemia with megaloblastic features of the more severe case. In vitro and in vivo complementation studies with fibroblasts from the affected individuals and zebrafish demonstrated the requirement of SFXN4 for mitochondrial respiratory homeostasis and erythropoiesis. Our findings establish mutations in SFXN4 as a cause of mitochondriopathy and macrocytic anemia.