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.

Long-range RNA structures in the human transcriptome beyond evolutionarily conserved regions.

RNA structure has been increasingly recognized as a critical player in the biogenesis and turnover of many transcripts classes. In eukaryotes, the prediction of RNA structure by thermodynamic modeling meets fundamental limitations due to the large sizes and complex, discontinuous organization of eukaryotic genes. Signatures of functional RNA structures can be found by detecting compensatory substitutions in homologous sequences, but a comparative approach is applicable only within conserved sequence blocks. Here, we developed a computational pipeline called PHRIC, which is not limited to conserved regions and relies on RNA contacts derived from RNA in situ conformation sequencing (RIC-seq) experiments. It extracts pairs of short RNA fragments surrounded by nested clusters of RNA contacts and predicts long, nearly perfect complementary base pairings formed between these fragments. In application to a panel of RIC-seq experiments in seven human cell lines, PHRIC predicted ~12,000 stable long-range RNA structures with equilibrium free energy below -15 kcal/mol, the vast majority of which fall outside of regions annotated as conserved among vertebrates. These structures, nevertheless, show some level of sequence conservation and remarkable compensatory substitution patterns in other clades. Furthermore, we found that introns have a higher propensity to form stable long-range RNA structures between each other, and moreover that RNA structures tend to concentrate within the same intron rather than connect adjacent introns. These results for the first time extend the application of proximity ligation assays to RNA structure prediction beyond conserved regions.

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.

RNA in situ conformation sequencing reveals novel long-range RNA structures with impact on splicing.

Over recent years, long-range RNA structure has emerged as a factor that is fundamental to alternative splicing regulation. An increasing number of human disorders are now being associated with splicing defects; hence it is essential to develop methods that assess long-range RNA structure experimentally. RNA in situ conformation sequencing (RIC-seq) is a method that recapitulates RNA structure within physiological RNA-protein complexes. In this work, we juxtapose pairs of conserved complementary regions (PCCRs) that were predicted in silico with the results of RIC-seq experiments conducted in seven human cell lines. We show statistically that RIC-seq support of PCCRs correlates with their properties, such as equilibrium free energy, presence of compensatory substitutions, and occurrence of A-to-I RNA editing sites and forked eCLIP peaks. Exons enclosed in PCCRs that are supported by RIC-seq tend to have weaker splice sites and lower inclusion rates, which is indicative of post-transcriptional splicing regulation mediated by RNA structure. Based on these findings, we prioritize PCCRs according to their RIC-seq support and show, using antisense nucleotides and minigene mutagenesis, that PCCRs in two disease-associated human genes, PHF20L1 and CASK, and also PCCRs in their murine orthologs, impact alternative splicing. In sum, we demonstrate how RIC-seq experiments can be used to discover functional long-range RNA structures, and particularly those that regulate alternative splicing.

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.

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.

APAF1-Binding Long Noncoding RNA Promotes Tumor Growth and Multidrug Resistance in Gastric Cancer by Blocking Apoptosome Assembly.

Chemotherapeutics remain the first choice for advanced gastric cancers (GCs). However, drug resistance and unavoidable severe toxicity lead to chemotherapy failure and poor prognosis. Long noncoding RNAs (lncRNAs) play critical roles in tumor progression in many cancers, including GC. Here, through RNA screening, an apoptotic protease-activating factor 1 (APAF1)-binding lncRNA (ABL) that is significantly elevated in cancerous GC tissues and an independent prognostic factor for GC patients is identified. Moreover, ABL overexpression inhibits GC cell apoptosis and promotes GC cell survival and multidrug resistance in GC xenograft and organoid models. Mechanistically, ABL directly binds to the RNA-binding protein IGF2BP1 via its KH1/2 domain, and then IGF2BP1 further recognizes the METTL3-mediated m6A modification on ABL, which maintains ABL stability. In addition, ABL can bind to the WD1/WD2 domain of APAF1, which competitively prevent cytochrome c from interacting with APAF1, blocking apoptosome assembly and caspase-9/3 activation; these events lead to resistance to cell death in GC cells. Intriguingly, targeting ABL using encapsulated liposomal siRNA can significantly enhance the sensitivity of GC cells to chemotherapy. Collectively, the results suggest that ABL can be a potential prognostic biomarker and therapeutic target in GC.

Research advances on interferon (IFN) response during BVDV infection.

Bovine viral diarrhea virus (BVDV) is an important pathogen responsible for significant economic loss to cattle. BVDV infection in pregnant cattle leads to fetal infection and reproductive losses, including early embryonic death, abortion, and stillbirth. Importantly, vaccinated heifers could not provide fetal protection against BVDV. It can be divided into two genotypes (BVDV-1 and BVDV-2) and two biotypes (cytopathic (CP) and non-cytopathic (NCP)). Infection with NCP-BVDV during gestation, the fetus becomes persistently infected (PI) and sheds BVDV throughout life, serving as the main source of infection for other cattle. BVDV potentially induces immunosuppression and aggravates bovine respiratory disease (BRD). Accordingly, BVDV infection results in a heterogeneous range of clinical signs and immune responses. Interferon (IFN) plays a vital role by mediating the innate immune response against antiviral infection through the Janus Kinase (JAK)-signal transducer and activator of transcription (STAT) pathway. BVDV infection can reportedly exert variable degrees of influence on IFN response. Interestingly, reports have suggested that IFN can exert a significant inhibitory effect on various viruses. Human IFN-α was used to restrain BVDV in vitro. In this article, we summarized the latest researches on IFN response during BVDV infection.

Recent advances in RNA structurome.

RNA structures are essential to support RNA functions and regulation in various biological processes. Recently, a range of novel technologies have been developed to decode genome-wide RNA structures and novel modes of functionality across a wide range of species. In this review, we summarize key strategies for probing the RNA structurome and discuss the pros and cons of representative technologies. In particular, these new technologies have been applied to dissect the structural landscape of the SARS-CoV-2 RNA genome. We also summarize the functionalities of RNA structures discovered in different regulatory layers-including RNA processing, transport, localization, and mRNA translation-across viruses, bacteria, animals, and plants. We review many versatile RNA structural elements in the context of different physiological and pathological processes (e.g., cell differentiation, stress response, and viral replication). Finally, we discuss future prospects for RNA structural studies to map the RNA structurome at higher resolution and at the single-molecule and single-cell level, and to decipher novel modes of RNA structures and functions for innovative applications.

Research Progress of H2S Donors Conjugate Drugs Based on ADTOH.

H2S is an endogenous gas signaling molecule and its multiple biological effects have been demonstrated. The abnormal level of H2S is closely related to the occurrence and development of many diseases, and H2S donors has important pharmacological implications. In recent years, H2S donors represented by ADTOH (5-(4-hydroxyphenyl)-3H-1,2-dithiole-3-thione) are often used to synthesize new 'conjugate' compounds that can release H2S and parent drugs. These hybrids retain the pharmacological activity of the parent drugs and H2S and have a synergistic effect. ADTOH and parent drug hybrids have become one of the important strategies for the development of H2S donor conjugate drugs. This review summarizes molecular hybrids between ADTOH and clinical drugs to provide new ideas for the study of H2S donor drug design.

The architecture of the SARS-CoV-2 RNA genome inside virion.

SARS-CoV-2 carries the largest single-stranded RNA genome and is the causal pathogen of the ongoing COVID-19 pandemic. How the SARS-CoV-2 RNA genome is folded in the virion remains unknown. To fill the knowledge gap and facilitate structure-based drug development, we develop a virion RNA in situ conformation sequencing technology, named vRIC-seq, for probing viral RNA genome structure unbiasedly. Using vRIC-seq data, we reconstruct the tertiary structure of the SARS-CoV-2 genome and reveal a surprisingly "unentangled globule" conformation. We uncover many long-range duplexes and higher-order junctions, both of which are under purifying selections and contribute to the sequential package of the SARS-CoV-2 genome. Unexpectedly, the D614G and the other two accompanying mutations may remodel duplexes into more stable forms. Lastly, the structure-guided design of potent small interfering RNAs can obliterate the SARS-CoV-2 in Vero cells. Overall, our work provides a framework for studying the genome structure, function, and dynamics of emerging deadly RNA viruses.

Global profiling of RNA-binding protein target sites by LACE-seq.

RNA-binding proteins (RBPs) have essential functions during germline and early embryo development. However, current methods are unable to identify the in vivo targets of a RBP in these low-abundance cells. Here, by coupling RBP-mediated reverse transcription termination with linear amplification of complementary DNA ends and sequencing, we present the LACE-seq method for identifying RBP-regulated RNA networks at or near the single-oocyte level. We determined the binding sites and regulatory mechanisms for several RBPs, including Argonaute 2 (Ago2), Mili, Ddx4 and Ptbp1, in mature mouse oocytes. Unexpectedly, transcriptomics and proteomics analysis of Ago2-/- oocytes revealed that Ago2 interacts with endogenous small interfering RNAs (endo-siRNAs) to repress mRNA translation globally. Furthermore, the Ago2 and endo-siRNA complexes fine-tune the transcriptome by slicing long terminal repeat retrotransposon-derived chimeric transcripts. The precise mapping of RBP-binding sites in low-input cells opens the door to studying the roles of RBPs in embryonic development and reproductive diseases.

Global in situ profiling of RNA-RNA spatial interactions with RIC-seq.

Emerging evidence has demonstrated that RNA-RNA interactions are vital in controlling diverse biological processes, including transcription, RNA splicing and protein translation. RNA in situ conformation sequencing (RIC-seq) is a technique for capturing protein-mediated RNA-RNA proximal interactions globally in living cells at single-base resolution. Cells are first treated with formaldehyde to fix all the protein-mediated RNA-RNA interactions in situ. After cell permeabilization and micrococcal nuclease digestion, the proximally interacting RNAs are 3' end-labeled with pCp-biotin and subsequently ligated using T4 RNA ligase. The chimeric RNAs are then enriched and converted into libraries for paired-end sequencing. After deep sequencing, computational analysis yields interaction strength scores for every base on proximally interacting RNAs in the starting populations. The whole experimental procedure is designed to be completed within 6 d, followed by an additional 8 d for computational analysis. RIC-seq technology can unbiasedly detect intra- and intermolecular RNA-RNA interactions, thereby rendering it useful for reconstructing RNA higher-order structures and identifying direct noncoding RNA targets.

Enhancer RNA: biogenesis, function, and regulation.

Enhancers are noncoding DNA elements that are present upstream or downstream of a gene to control its spatial and temporal expression. Specific histone modifications, such as monomethylation on histone H3 lysine 4 (H3K4me1) and H3K27ac, have been widely used to assign enhancer regions in mammalian genomes. In recent years, emerging evidence suggests that active enhancers are bidirectionally transcribed to produce enhancer RNAs (eRNAs). This finding not only adds a new reliable feature to define enhancers but also raises a fundamental question of how eRNAs function to activate transcription. Although some believe that eRNAs are merely transcriptional byproducts, many studies have demonstrated that eRNAs execute crucial tasks in regulating chromatin conformation and transcription activation. In this review, we summarize the current understanding of eRNAs from their biogenesis, functions, and regulation to their pathological significance. Additionally, we discuss the challenges and possible mechanisms of eRNAs in regulated transcription.

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.

The RNA-binding protein ROD1/PTBP3 cotranscriptionally defines AID-loading sites to mediate antibody class switch in mammalian genomes.

Activation-induced cytidine deaminase (AID) mediates class switching by binding to a small fraction of single-stranded DNA (ssDNA) to diversify the antibody repertoire. The precise mechanism for highly selective AID targeting in the genome has remained elusive. Here, we report an RNA-binding protein, ROD1 (also known as PTBP3), that is both required and sufficient to define AID-binding sites genome-wide in activated B cells. ROD1 interacts with AID via an ultraconserved loop, which proves to be critical for the recruitment of AID to ssDNA using bi-directionally transcribed nascent RNAs as stepping stones. Strikingly, AID-specific mutations identified in human patients with hyper-IgM syndrome type 2 (HIGM2) completely disrupt the AID interacting surface with ROD1, thereby abolishing the recruitment of AID to immunoglobulin (Ig) loci. Together, our results suggest that bi-directionally transcribed RNA traps the RNA-binding protein ROD1, which serves as a guiding system for AID to load onto specific genomic loci to induce DNA rearrangement during immune responses.

Predicting and Interpreting the Structure of Type IV Pilus of Electricigens by Molecular Dynamics Simulations.

Nanowires that transfer electrons to extracellular acceptors are important in organic matter degradation and nutrient cycling in the environment. Geobacter pili of the group of Type IV pilus are regarded as nanowire-like biological structures. However, determination of the structure of pili remains challenging due to the insolubility of monomers, presence of surface appendages, heterogeneity of the assembly, and low-resolution of electron microscopy techniques. Our previous study provided a method to predict structures for Type IV pili. In this work, we improved on our previous method using molecular dynamics simulations to optimize structures of Neisseria gonorrhoeae (GC), Neisseria meningitidis and Geobacter uraniireducens pilus. Comparison between the predicted structures for GC and Neisseria meningitidis pilus and their native structures revealed that proposed method could predict Type IV pilus successfully. According to the predicted structures, the structural basis for conductivity in G.uraniireducens pili was attributed to the three N-terminal aromatic amino acids. The aromatics were interspersed within the regions of charged amino acids, which may influence the configuration of the aromatic contacts and the rate of electron transfer. These results will supplement experimental research into the mechanism of long-rang electron transport along pili of electricigens.

Efficient identification of SNPs in pooled DNA samples using a dual mononucleotide addition-based sequencing method.

Identifying single nucleotide polymorphism (SNPs) from pooled samples is critical for many studies and applications. SNPs determined by next-generation sequencing results may suffer from errors in both base calling and read mapping. Taking advantage of dual mononucleotide addition-based pyrosequencing, we present Epds, a method to efficiently identify SNPs from pooled DNA samples. On the basis of only five patterns of non-synchronistic extensions between the wild and mutant sequences using dual mononucleotide addition-based pyrosequencing, we employed an enumerative algorithm to infer the mutant locus and estimate the proportion of mutant sequence. According to the profiles resulting from three runs with distinct dual mononucleotide additions, Epds could recover the mutant bases. Results showed that our method had a false-positive rate of less than 3%. Series of simulations revealed that Epds outperformed the current method (PSM) in many situations. Finally, experiments based on profiles produced by real sequencing proved that our method could be successfully applied for the identification of mutants from pooled samples. The software for implementing this method and the experimental data are available at http://bioinfo.seu.edu.cn/Epds .

Ehapp2: Estimate haplotype frequencies from pooled sequencing data with prior database information.

To reduce the cost of large-scale re-sequencing, multiple individuals are pooled together and sequenced called pooled sequencing. Pooled sequencing could provide a cost-effective alternative to sequencing individuals separately. To facilitate the application of pooled sequencing in haplotype-based diseases association analysis, the critical procedure is to accurately estimate haplotype frequencies from pooled samples. Here we present Ehapp2 for estimating haplotype frequencies from pooled sequencing data by utilizing a database which provides prior information of known haplotypes. We first translate the problem of estimating frequency for each haplotype into finding a sparse solution for a system of linear equations, where the NNREG algorithm is employed to achieve the solution. Simulation experiments reveal that Ehapp2 is robust to sequencing errors and able to estimate the frequencies of haplotypes with less than 3% average relative difference for pooled sequencing of mixture of real Drosophila haplotypes with 50× total coverage even when the sequencing error rate is as high as 0.05. Owing to the strategy that proportions for local haplotypes spanning multiple SNPs are accurately calculated first, Ehapp2 retains excellent estimation for recombinant haplotypes resulting from chromosomal crossover. Comparisons with present methods reveal that Ehapp2 is state-of-the-art for many sequencing study designs and more suitable for current massive parallel sequencing.