SARS-CoV-2 RNA stabilizes host mRNAs to elicit immunopathogenesis.

SARS-CoV-2 RNA interacts with host factors to suppress interferon responses and simultaneously induces cytokine release to drive the development of severe coronavirus disease 2019 (COVID-19). However, how SARS-CoV-2 hijacks host RNAs to elicit such imbalanced immune responses remains elusive. Here, we analyzed SARS-CoV-2 RNA in situ structures and interactions in infected cells and patient lung samples using RIC-seq. We discovered that SARS-CoV-2 RNA forms 2,095 potential duplexes with the 3' UTRs of 205 host mRNAs to increase their stability by recruiting RNA-binding protein YBX3 in A549 cells. Disrupting the SARS-CoV-2-to-host RNA duplex or knocking down YBX3 decreased host mRNA stability and reduced viral replication. Among SARS-CoV-2-stabilized host targets, NFKBIZ was crucial for promoting cytokine production and reducing interferon responses, probably contributing to cytokine storm induction. Our study uncovers the crucial roles of RNA-RNA interactions in the immunopathogenesis of RNA viruses such as SARS-CoV-2 and provides valuable host targets for drug development.

Capture RIC-seq reveals positional rules of PTBP1-associated RNA loops in splicing regulation.

RNA-binding proteins (RBPs) bind at different positions of the pre-mRNA molecules to promote or reduce the usage of a particular exon. Seeking to understand the working principle of these positional effects, we develop a capture RIC-seq (CRIC-seq) method to enrich specific RBP-associated in situ proximal RNA-RNA fragments for deep sequencing. We determine hnRNPA1-, SRSF1-, and PTBP1-associated proximal RNA-RNA contacts and regulatory mechanisms in HeLa cells. Unexpectedly, the 3D RNA map analysis shows that PTBP1-associated loops in individual introns preferentially promote cassette exon splicing by accelerating asymmetric intron removal, whereas the loops spanning across cassette exon primarily repress splicing. These "positional rules" can faithfully predict PTBP1-regulated splicing outcomes. We further demonstrate that cancer-related splicing quantitative trait loci can disrupt RNA loops by reducing PTBP1 binding on pre-mRNAs to cause aberrant splicing in tumors. Our study presents a powerful method for exploring the functions of RBP-associated RNA-RNA proximal contacts in gene regulation and disease.

LncRNA CCTT-mediated RNA-DNA and RNA-protein interactions facilitate the recruitment of CENP-C to centromeric DNA during kinetochore assembly.

Kinetochore assembly on centromeres is central for chromosome segregation, and defects in this process cause mitotic errors and aneuploidy. Besides the well-established protein network, emerging evidence suggests the involvement of regulatory RNA in kinetochore assembly; however, it has remained elusive about the identity of such RNA, let alone its mechanism of action in this critical process. Here, we report CCTT, a previously uncharacterized long non-coding RNA (lncRNA) transcribed from the arm of human chromosome 17, which plays a vital role in kinetochore assembly. We show that CCTT highly localizes to all centromeres via the formation of RNA-DNA triplex and specifically interacts with CENP-C to help engage this blueprint protein in centromeres, and consequently, CCTT loss triggers extensive mitotic errors and aneuploidy. These findings uncover a non-centromere-derived lncRNA that recruits CENP-C to centromeres and shed critical lights on the function of centromeric DNA sequences as anchor points for kinetochore assembly.

MicroRNA directly enhances mitochondrial translation during muscle differentiation.

MicroRNAs are well known to mediate translational repression and mRNA degradation in the cytoplasm. Various microRNAs have also been detected in membrane-compartmentalized organelles, but the functional significance has remained elusive. Here, we report that miR-1, a microRNA specifically induced during myogenesis, efficiently enters the mitochondria where it unexpectedly stimulates, rather than represses, the translation of specific mitochondrial genome-encoded transcripts. We show that this positive effect requires specific miR:mRNA base-pairing and Ago2, but not its functional partner GW182, which is excluded from the mitochondria. We provide evidence for the direct action of Ago2 in mitochondrial translation by crosslinking immunoprecipitation coupled with deep sequencing (CLIP-seq), functional rescue with mitochondria-targeted Ago2, and selective inhibition of the microRNA machinery in the cytoplasm. These findings unveil a positive function of microRNA in mitochondrial translation and suggest a highly coordinated myogenic program via miR-1-mediated translational stimulation in the mitochondria and repression in the cytoplasm.

CLP1 founder mutation links tRNA splicing and maturation to cerebellar development and neurodegeneration.

Neurodegenerative diseases can occur so early as to affect neurodevelopment. From a cohort of more than 2,000 consanguineous families with childhood neurological disease, we identified a founder mutation in four independent pedigrees in cleavage and polyadenylation factor I subunit 1 (CLP1). CLP1 is a multifunctional kinase implicated in tRNA, mRNA, and siRNA maturation. Kinase activity of the CLP1 mutant protein was defective, and the tRNA endonuclease complex (TSEN) was destabilized, resulting in impaired pre-tRNA cleavage. Germline clp1 null zebrafish showed cerebellar neurodegeneration that was rescued by wild-type, but not mutant, human CLP1 expression. Patient-derived induced neurons displayed both depletion of mature tRNAs and accumulation of unspliced pre-tRNAs. Transfection of partially processed tRNA fragments into patient cells exacerbated an oxidative stress-induced reduction in cell survival. Our data link tRNA maturation to neuronal development and neurodegeneration through defective CLP1 function in humans.

Complementary Alu sequences mediate enhancer-promoter selectivity.

Enhancers determine spatiotemporal gene expression programs by engaging with long-range promoters1-4. However, it remains unknown how enhancers find their cognate promoters. We recently developed a RNA in situ conformation sequencing technology to identify enhancer-promoter connectivity using pairwise interacting enhancer RNAs and promoter-derived noncoding RNAs5,6. Here we apply this technology to generate high-confidence enhancer-promoter RNA interaction maps in six additional cell lines. Using these maps, we discover that 37.9% of the enhancer-promoter RNA interaction sites are overlapped with Alu sequences. These pairwise interacting Alu and non-Alu RNA sequences tend to be complementary and potentially form duplexes. Knockout of Alu elements compromises enhancer-promoter looping, whereas Alu insertion or CRISPR-dCasRx-mediated Alu tethering to unregulated promoter RNAs can create new loops to homologous enhancers. Mapping 535,404 noncoding risk variants back to the enhancer-promoter RNA interaction maps enabled us to construct variant-to-function maps for interpreting their molecular functions, including 15,318 deletions or insertions in 11,677 Alu elements that affect 6,497 protein-coding genes. We further demonstrate that polymorphic Alu insertion at the PTK2 enhancer can promote tumorigenesis. Our study uncovers a principle for determining enhancer-promoter pairing specificity and provides a framework to link noncoding risk variants to their molecular functions.

Direct conversion of fibroblasts to neurons by reprogramming PTB-regulated microRNA circuits.

The induction of pluripotency or trans-differentiation of one cell type to another can be accomplished with cell-lineage-specific transcription factors. Here, we report that repression of a single RNA binding polypyrimidine-tract-binding (PTB) protein, which occurs during normal brain development via the action of miR-124, is sufficient to induce trans-differentiation of fibroblasts into functional neurons. Besides its traditional role in regulated splicing, we show that PTB has a previously undocumented function in the regulation of microRNA functions, suppressing or enhancing microRNA targeting by competitive binding on target mRNA or altering local RNA secondary structure. A key event during neuronal induction is the relief of PTB-mediated blockage of microRNA action on multiple components of the REST complex, thereby derepressing a large array of neuronal genes, including miR-124 and multiple neuronal-specific transcription factors, in nonneuronal cells. This converts a negative feedback loop to a positive one to elicit cellular reprogramming to the neuronal lineage.

Climate-driven flyway changes and memory-based long-distance migration.

Millions of migratory birds occupy seasonally favourable breeding grounds in the Arctic1, but we know little about the formation, maintenance and future of the migration routes of Arctic birds and the genetic determinants of migratory distance. Here we established a continental-scale migration system that used satellite tracking to follow 56 peregrine falcons (Falco peregrinus) from 6 populations that breed in the Eurasian Arctic, and resequenced 35 genomes from 4 of these populations. The breeding populations used five migration routes across Eurasia, which were probably formed by longitudinal and latitudinal shifts in their breeding grounds during the transition from the Last Glacial Maximum to the Holocene epoch. Contemporary environmental divergence between the routes appears to maintain their distinctiveness. We found that the gene ADCY8 is associated with population-level differences in migratory distance. We investigated the regulatory mechanism of this gene, and found that long-term memory was the most likely selective agent for divergence in ADCY8 among the peregrine populations. Global warming is predicted to influence migration strategies and diminish the breeding ranges of peregrine populations of the Eurasian Arctic. Harnessing ecological interactions and evolutionary processes to study climate-driven changes in migration can facilitate the conservation of migratory birds.

Author Correction: Reversing a model of Parkinson's disease with in situ converted nigral neurons.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

RIC-seq for global in situ profiling of RNA-RNA spatial interactions.

Highly structured RNA molecules usually interact with each other, and associate with various RNA-binding proteins, to regulate critical biological processes. However, RNA structures and interactions in intact cells remain largely unknown. Here, by coupling proximity ligation mediated by RNA-binding proteins with deep sequencing, we report an RNA in situ conformation sequencing (RIC-seq) technology for the global profiling of intra- and intermolecular RNA-RNA interactions. This technique not only recapitulates known RNA secondary structures and tertiary interactions, but also facilitates the generation of three-dimensional (3D) interaction maps of RNA in human cells. Using these maps, we identify noncoding RNA targets globally, and discern RNA topological domains and trans-interacting hubs. We reveal that the functional connectivity of enhancers and promoters can be assigned using their pairwise-interacting RNAs. Furthermore, we show that CCAT1-5L-a super-enhancer hub RNA-interacts with the RNA-binding protein hnRNPK, as well as RNA derived from the MYC promoter and enhancer, to boost MYC transcription by modulating chromatin looping. Our study demonstrates the power and applicability of RIC-seq in discovering the 3D structures, interactions and regulatory roles of RNA.

Reversing a model of Parkinson's disease with in situ converted nigral neurons.

Parkinson's disease is characterized by loss of dopamine neurons in the substantia nigra1. Similar to other major neurodegenerative disorders, there are no disease-modifying treatments for Parkinson's disease. While most treatment strategies aim to prevent neuronal loss or protect vulnerable neuronal circuits, a potential alternative is to replace lost neurons to reconstruct disrupted circuits2. Here we report an efficient one-step conversion of isolated mouse and human astrocytes to functional neurons by depleting the RNA-binding protein PTB (also known as PTBP1). Applying this approach to the mouse brain, we demonstrate progressive conversion of astrocytes to new neurons that innervate into and repopulate endogenous neural circuits. Astrocytes from different brain regions are converted to different neuronal subtypes. Using a chemically induced model of Parkinson's disease in mouse, we show conversion of midbrain astrocytes to dopaminergic neurons, which provide axons to reconstruct the nigrostriatal circuit. Notably, re-innervation of striatum is accompanied by restoration of dopamine levels and rescue of motor deficits. A similar reversal of disease phenotype is also accomplished by converting astrocytes to neurons using antisense oligonucleotides to transiently suppress PTB. These findings identify a potentially powerful and clinically feasible approach to treating neurodegeneration by replacing lost neurons.

Targeting intratumor heterogeneity suppresses colorectal cancer chemoresistance and metastasis.

Intratumor heterogeneity (ITH) is a barrier to effective therapy. However, it is largely unknown how ITH is established at the onset of tumor progression, such as in colorectal cancer (CRC). Here, we integrate single-cell RNA-seq and functional validation to show that asymmetric division of CRC stem-like cells (CCSC) is critical for early ITH establishment. We find that CCSC-derived xenografts contain seven cell subtypes, including CCSCs, that dynamically change during CRC xenograft progression. Furthermore, three of the subtypes are generated by asymmetric division of CCSCs. They are functionally distinct and appear at the early stage of xenografts. In particular, we identify a chemoresistant and an invasive subtype, and investigate the regulators that control their generation. Finally, we show that targeting the regulators influences cell subtype composition and CRC progression. Our findings demonstrate that asymmetric division of CCSCs contributes to the early establishment of ITH. Targeting asymmetric division may alter ITH and benefit CRC therapy.

NOVA1 prevents overactivation of the unfolded protein response and facilitates chromatin access during human white adipogenesis.

The molecular mechanism underlying white adipogenesis in humans has not been fully elucidated beyond the transcriptional level. Here, we found that the RNA-binding protein NOVA1 is required for the adipogenic differentiation of human mesenchymal stem cells. By thoroughly exploring the interactions between NOVA1 and its binding RNA, we proved that NOVA1 deficiency resulted in the aberrant splicing of DNAJC10 with an in-frame premature stop codon, reduced DNAJC10 expression at the protein level and hyperactivation of the unfolded protein response (UPR). Moreover, NOVA1 knockdown abrogated the down-regulation of NCOR2 during adipogenesis and up-regulated the 47b+ splicing isoform, which led to decreased chromatin accessibility at the loci of lipid metabolism genes. Interestingly, these effects on human adipogenesis could not be recapitulated in mice. Further analysis of multispecies genomes and transcriptomes indicated that NOVA1-targeted RNA splicing is evolutionarily regulated. Our findings provide evidence for human-specific roles of NOVA1 in coordinating splicing and cell organelle functions during white adipogenesis.

Sex-lethal regulates back-splicing and generation of the sex-differentially expressed circular RNAs.

Conversely to canonical splicing, back-splicing connects the upstream 3' splice site (SS) with a downstream 5'SS and generates exonic circular RNAs (circRNAs) that are widely identified and have regulatory functions in eukaryotic gene expression. However, sex-specific back-splicing in Drosophila has not been investigated and its regulation remains unclear. Here, we performed multiple RNA analyses of a variety sex-specific Drosophila samples and identified over ten thousand circular RNAs, in which hundreds are sex-differentially and -specifically back-spliced. Intriguingly, we found that expression of SXL, an RNA-binding protein encoded by Sex-lethal (Sxl), the master Drosophila sex-determination gene that is only spliced into functional proteins in females, promoted back-splicing of many female-differential circRNAs in the male S2 cells, whereas expression of a SXL mutant (SXLRRM) did not promote those events. Using a monoclonal antibody, we further obtained the transcriptome-wide RNA-binding sites of SXL through PAR-CLIP. After splicing assay of mini-genes with mutations in the SXL-binding sites, we revealed that SXL-binding on flanking exons and introns of pre-mRNAs facilitates back-splicing, whereas SXL-binding on the circRNA exons inhibits back-splicing. This study provides strong evidence that SXL has a regulatory role in back-splicing to generate sex-specific and -differential circRNAs, as well as in the initiation of sex-determination cascade through canonical forward-splicing.

RBFox2 Binds Nascent RNA to Globally Regulate Polycomb Complex 2 Targeting in Mammalian Genomes.

Increasing evidence suggests that diverse RNA binding proteins (RBPs) interact with regulatory RNAs to regulate transcription. RBFox2 is a well-characterized pre-mRNA splicing regulator, but we now encounter an unexpected paradigm where depletion of this RBP induces widespread increase in nascent RNA production in diverse cell types. Chromatin immunoprecipitation sequencing (ChIP-seq) reveals extensive interaction of RBFox2 with chromatin in a nascent RNA-dependent manner. Bayesian network analysis connects RBFox2 to Polycomb complex 2 (PRC2) and H3K27me3, and biochemical experiments demonstrate the ability of RBFox2 to directly interact with PRC2. Strikingly, RBFox2 inactivation eradicates PRC2 targeting on the majority of bivalent gene promoters and leads to transcriptional de-repression. Together, these findings uncover a mechanism underlying the enigmatic association of PRC2 with numerous active genes, highlight the importance of gene body sequences to gauge transcriptional output, and suggest nascent RNAs as critical signals for transcriptional feedback control to maintain homeostatic gene expression in mammalian genomes.

RBFox2-miR-34a-Jph2 axis contributes to cardiac decompensation during heart failure.

Heart performance relies on highly coordinated excitation-contraction (EC) coupling, and defects in this critical process may be exacerbated by additional genetic defects and/or environmental insults to cause eventual heart failure. Here we report a regulatory pathway consisting of the RNA binding protein RBFox2, a stress-induced microRNA miR-34a, and the essential EC coupler JPH2. In this pathway, initial cardiac defects diminish RBFox2 expression, which induces transcriptional repression of miR-34a, and elevated miR-34a targets Jph2 to impair EC coupling, which further manifests heart dysfunction, leading to progressive heart failure. The key contribution of miR-34a to this process is further established by administrating its mimic, which is sufficient to induce cardiac defects, and by using its antagomir to alleviate RBFox2 depletion-induced heart dysfunction. These findings elucidate a potential feed-forward mechanism to account for a critical transition to cardiac decompensation and suggest a potential therapeutic avenue against heart failure.

WNT7A and PAX6 define corneal epithelium homeostasis and pathogenesis.

The surface of the cornea consists of a unique type of non-keratinized epithelial cells arranged in an orderly fashion, and this is essential for vision by maintaining transparency for light transmission. Cornea epithelial cells (CECs) undergo continuous renewal from limbal stem or progenitor cells (LSCs), and deficiency in LSCs or corneal epithelium--which turns cornea into a non-transparent, keratinized skin-like epithelium--causes corneal surface disease that leads to blindness in millions of people worldwide. How LSCs are maintained and differentiated into corneal epithelium in healthy individuals and which key molecular events are defective in patients have been largely unknown. Here we report establishment of an in vitro feeder-cell-free LSC expansion and three-dimensional corneal differentiation protocol in which we found that the transcription factors p63 (tumour protein 63) and PAX6 (paired box protein PAX6) act together to specify LSCs, and WNT7A controls corneal epithelium differentiation through PAX6. Loss of WNT7A or PAX6 induces LSCs into skin-like epithelium, a critical defect tightly linked to common human corneal diseases. Notably, transduction of PAX6 in skin epithelial stem cells is sufficient to convert them to LSC-like cells, and upon transplantation onto eyes in a rabbit corneal injury model, these reprogrammed cells are able to replenish CECs and repair damaged corneal surface. These findings suggest a central role of the WNT7A-PAX6 axis in corneal epithelial cell fate determination, and point to a new strategy for treating corneal surface diseases.