Mapping the human miRNA interactome by CLASH reveals frequent noncanonical binding.

MicroRNAs (miRNAs) play key roles in gene regulation, but reliable bioinformatic or experimental identification of their targets remains difficult. To provide an unbiased view of human miRNA targets, we developed a technique for ligation and sequencing of miRNA-target RNA duplexes associated with human AGO1. Here, we report data sets of more than 18,000 high-confidence miRNA-mRNA interactions. The binding of most miRNAs includes the 5' seed region, but around 60% of seed interactions are noncanonical, containing bulged or mismatched nucleotides. Moreover, seed interactions are generally accompanied by specific, nonseed base pairing. 18% of miRNA-mRNA interactions involve the miRNA 3' end, with little evidence for 5' contacts, and some of these were functionally validated. Analyses of miRNA:mRNA base pairing showed that miRNA species systematically differ in their target RNA interactions, and strongly overrepresented motifs were found in the interaction sites of several miRNAs. We speculate that these affect the response of RISC to miRNA-target binding.

Nascent Transcript Folding Plays a Major Role in Determining RNA Polymerase Elongation Rates.

Transcription elongation rates influence RNA processing, but sequence-specific regulation is poorly understood. We addressed this in vivo, analyzing RNAPI in S. cerevisiae. Mapping RNAPI by Miller chromatin spreads or UV crosslinking revealed 5' enrichment and strikingly uneven local polymerase occupancy along the rDNA, indicating substantial variation in transcription speed. Two features of the nascent transcript correlated with RNAPI distribution: folding energy and GC content in the transcription bubble. In vitro experiments confirmed that strong RNA structures close to the polymerase promote forward translocation and limit backtracking, whereas high GC in the transcription bubble slows elongation. A mathematical model for RNAPI elongation confirmed the importance of nascent RNA folding in transcription. RNAPI from S. pombe was similarly sensitive to transcript folding, as were S. cerevisiae RNAPII and RNAPIII. For RNAPII, unstructured RNA, which favors slowed elongation, was associated with faster cotranscriptional splicing and proximal splice site use, indicating regulatory significance for transcript folding.

RNA Conformation Capture by Proximity Ligation.

RNA proximity ligation is a set of molecular biology techniques used to analyze the conformations and spatial proximity of RNA molecules within cells. A typical experiment starts with cross-linking of a biological sample using UV light or psoralen, followed by partial fragmentation of RNA, RNA-RNA ligation, library preparation, and high-throughput sequencing. In the past decade, proximity ligation has been used to study structures of individual RNAs, networks of interactions between small RNAs and their targets, and whole RNA-RNA interactomes, in models ranging from bacteria to animal tissues and whole animals. Here, we provide an overview of the field, highlight the main findings, review the recent experimental and computational developments, and provide troubleshooting advice for new users. In the final section, we draw parallels between DNA and RNA proximity ligation and speculate on possible future research directions.

Novel functional small RNAs are selectively loaded onto mammalian Ago1.

Argonaute (Ago) proteins function in RNA silencing as components of the RNA-induced silencing complex (RISC). In lower organisms, the small interfering RNA and miRNA pathways diverge due in part to sorting mechanisms that direct distinct small RNA (sRNA) duplexes onto specific Ago-RISCs. However, such sorting mechanisms appear to be lost in mammals. miRNAs appear not to distinguish among Ago1-4. To determine the effect of viral infection on the sorting system, we compared the content of deep-sequenced RNA extracted from immunoprecipitation experiments with the Ago1 and Ago2 proteins using Epstein-Barr virus (EBV)-infected cells. Consistent with previous observations, sequence tags derived from miRNA loci in EBV and humans globally associate in approximately equivalent amounts with Ago1 and Ago2. Interestingly, additional sRNAs, which have not been registered as miRNAs, were associated with Ago1. Among them, some unique sequence tags derived from tandem loci in the human genome associate exclusively with Ago1 but not, or rarely, with Ago2. This is supported by the observation that the expression of the unique sRNAs in the cells is highly dependent on Ago1 proteins. When we knocked down Ago1, the expression of the Ago1-specific sRNAs decreased dramatically. Most importantly, the Ago1-specific sRNAs bound to mRNAs and regulated target genes and were dramatically upregulated, depending on the EBV life cycle. Therefore, even in mammals, the sorting mechanism in the Ago1-4 family is functional. Moreover, the existence of Ago1-specific sRNAs implies vital roles in some aspects of mammalian biology.

Hyb: a bioinformatics pipeline for the analysis of CLASH (crosslinking, ligation and sequencing of hybrids) data.

Associations between proteins and RNA-RNA duplexes are important in post-transcriptional regulation of gene expression. The CLASH (Cross-linking, Ligation and Sequencing of Hybrids) technique captures RNA-RNA interactions by physically joining two RNA molecules associated with a protein complex into a single chimeric RNA molecule. These events are relatively rare and considerable effort is needed to detect a small number of chimeric sequences amongst millions of non-chimeric cDNA reads resulting from a CLASH experiment. We present the "hyb" bioinformatics pipeline, which we developed to analyse high-throughput cDNA sequencing data from CLASH experiments. Although primarily designed for use with AGO CLASH data, hyb can also be used for the detection and annotation of chimeric reads in other high-throughput sequencing datasets. We examined the sensitivity and specificity of chimera detection in a test dataset using the BLAST, BLAST+, BLAT, pBLAT and Bowtie2 read alignment programs. We obtained the most reliable results in the shortest time using a combination of preprocessing with Flexbar and subsequent read-mapping using Bowtie2. The "hyb" software is distributed under the GNU GPL (General Public License) and can be downloaded from https://github.com/gkudla/hyb.

Murine cytomegalovirus encodes a miR-27 inhibitor disguised as a target.

Individual microRNAs (miRNAs) are rapidly down-regulated during conditions of cellular activation and infection, but factors mediating miRNA turnover are poorly understood. Infection of mouse cells with murine cytomegalovirus (MCMV) induces the rapid down-regulation of an antiviral cellular miRNA, miR-27. Here, we identify a transcript produced by MCMV that binds to miR-27 and mediates its degradation. UV-crosslinking and high-throughput sequencing [CRAC (UV-crosslinking and analysis of cDNA)] identified MCMV RNA segments associated with the miRNA-binding protein Argonaute 2 (Ago2). A cluster of hits mapped to a predicted miR-27-binding site in the 3'UTR of the previously uncharacterized ORF, m169. The expression kinetics of the m169 transcript correlated with degradation of miR-27 during infection, and m169 expression inhibited miR-27 functional activity in a reporter assay. siRNA knockdown of m169 demonstrated its requirement for miR-27 degradation following infection and did not affect other host miRNAs. Substitution of the miR-27-binding site in m169 to create complementarity to a different cellular miRNA, miR-24, resulted in down-regulation of only miR-24 following infection. The m169 transcript is cytoplasmic, capped, polyadenylated, and interacts with miRNA-27 through seed pairing: characteristic features of the normal messenger RNA (mRNA) targets of miRNAs. This virus-host interaction reveals a mode of miRNA regulation in which a mRNA directs the degradation of a miRNA. We speculate that RNA-mediated miRNA degradation could be a more general viral strategy for manipulating host cells.

A disease-linked lncRNA mutation in RNase MRP inhibits ribosome synthesis.

RMRP encodes a non-coding RNA forming the core of the RNase MRP ribonucleoprotein complex. Mutations cause Cartilage Hair Hypoplasia (CHH), characterized by skeletal abnormalities and impaired T cell activation. Yeast RNase MRP cleaves a specific site in the pre-ribosomal RNA (pre-rRNA) during ribosome synthesis. CRISPR-mediated disruption of RMRP in human cells lines caused growth arrest, with pre-rRNA accumulation. Here, we analyzed disease-relevant primary cells, showing that mutations in RMRP impair mouse T cell activation and delay pre-rRNA processing. Patient-derived human fibroblasts with CHH-linked mutations showed similar pre-rRNA processing delay. Human cells engineered with the most common CHH mutation (70AG in RMRP) show specifically impaired pre-rRNA processing, resulting in reduced mature rRNA and a reduced ratio of cytosolic to mitochondrial ribosomes. Moreover, the 70AG mutation caused a reduction in intact RNase MRP complexes. Together, these results indicate that CHH is a ribosomopathy.

The RNA-Binding Landscape of HAX1 Protein Indicates Its Involvement in Translation and Ribosome Assembly.

HAX1 is a human protein with no known homologues or structural domains. Mutations in the HAX1 gene cause severe congenital neutropenia through mechanisms that are poorly understood. Previous studies reported the RNA-binding capacity of HAX1, but the role of this binding in physiology and pathology remains unexplained. Here, we report the transcriptome-wide characterization of HAX1 RNA targets using RIP-seq and CRAC, indicating that HAX1 binds transcripts involved in translation, ribosome biogenesis, and rRNA processing. Using CRISPR knockouts, we find that HAX1 RNA targets partially overlap with transcripts downregulated in HAX1 KO, implying a role in mRNA stabilization. Gene ontology analysis demonstrated that genes differentially expressed in HAX1 KO (including genes involved in ribosome biogenesis and translation) are also enriched in a subset of genes whose expression correlates with HAX1 expression in four analyzed neoplasms. The functional connection to ribosome biogenesis was also demonstrated by gradient sedimentation ribosome profiles, which revealed differences in the small subunit:monosome ratio in HAX1 WT/KO. We speculate that changes in HAX1 expression may be important for the etiology of HAX1-linked diseases through dysregulation of translation.

Integrative vectors for regulated expression of SARS-CoV-2 proteins implicated in RNA metabolism.

Infection with SARS-CoV-2 is expected to result in substantial reorganization of host cell RNA metabolism. We identified 14 proteins that were predicted to interact with host RNAs or RNA binding proteins, based on published data for SARS-CoV and SARS-CoV-2. Here, we describe a series of affinity-tagged and codon-optimized expression constructs for each of these 14 proteins. Each viral gene was separately tagged at the N-terminus with Flag-His 8, the C-terminus with His 8-Flag, or left untagged. The resulting constructs were stably integrated into the HEK293 Flp-In T-REx genome. Each viral gene was expressed under the control of an inducible Tet-On promoter, allowing expression levels to be tuned to match physiological conditions during infection. Expression time courses were successfully generated for most of the fusion proteins and quantified by western blot. A few fusion proteins were poorly expressed, whereas others, including Nsp1, Nsp12, and N protein, were toxic unless care was taken to minimize background expression. All plasmids can be obtained from Addgene and cell lines are available. We anticipate that availability of these resources will facilitate a more detailed understanding of coronavirus molecular biology.

High guanine and cytosine content increases mRNA levels in mammalian cells.

Mammalian genes are highly heterogeneous with respect to their nucleotide composition, but the functional consequences of this heterogeneity are not clear. In the previous studies, weak positive or negative correlations have been found between the silent-site guanine and cytosine (GC) content and expression of mammalian genes. However, previous studies disregarded differences in the genomic context of genes, which could potentially obscure any correlation between GC content and expression. In the present work, we directly compared the expression of GC-rich and GC-poor genes placed in the context of identical promoters and UTR sequences. We performed transient and stable transfections of mammalian cells with GC-rich and GC-poor versions of Hsp70, green fluorescent protein, and IL2 genes. The GC-rich genes were expressed several-fold to over a 100-fold more efficiently than their GC-poor counterparts. This effect was not due to different translation rates of GC-rich and GC-poor mRNA. On the contrary, the efficient expression of GC-rich genes resulted from their increased steady-state mRNA levels. mRNA degradation rates were not correlated with GC content, suggesting that efficient transcription or mRNA processing is responsible for the high expression of GC-rich genes. We conclude that silent-site GC content correlates with gene expression efficiency in mammalian cells.

Identification of miRNA-Target RNA Interactions Using CLASH.

We present a detailed protocol for the experimental identification of miRNA-target RNA interaction sites using cross-linking, ligation, and sequencing of hybrids (CLASH). The basis of the technique is the purification of UV-stabilized Argonaute (AGO)-RNA complexes assembled in living cells, with subsequent ligation of AGO-associated RNA-RNA duplexes to form chimeric RNAs. Following cDNA synthesis, DNA library preparation and high-throughput sequencing, interacting RNA molecules are unambiguously identified as chimeric reads in bioinformatic analysis of sequencing data. CLASH potentially recovers any RNA duplex that is bound by RNA-binding protein, so modified approaches would be suitable for the identification of many other inter- and intramolecular RNA-RNA interactions. Since CLASH analysis is independent of bioinformatic predictions it allows the identification and analysis of RNA targeting rules in an unbiased way.

Mapping the miRNA interactome by cross-linking ligation and sequencing of hybrids (CLASH).

RNA-RNA interactions have critical roles in many cellular processes, but studying them is difficult and laborious. Here we describe an experimental procedure, termed cross-linking ligation and sequencing of hybrids (CLASH), which allows high-throughput identification of sites of RNA-RNA interaction. During CLASH, a tagged bait protein is UV-cross-linked in cell cultures to stabilize RNA interactions, and it is purified under denaturing conditions. RNAs associated with the bait protein are partially truncated, and the ends of RNA duplexes are ligated together. After linker addition, cDNA library preparation and high-throughput sequencing, the ligated duplexes give rise to chimeric cDNAs, which unambiguously identify RNA-RNA interaction sites independent of bioinformatic predictions. This protocol is optimized for studying miRNA targets bound by Argonaute (AGO) proteins, but it should be easily adapted for other RNA-binding proteins and classes of RNA. The protocol requires ∼5 d to complete, excluding the time required for high-throughput sequencing and bioinformatic analyses.

Gene conversion and GC-content evolution in mammalian Hsp70.

To investigate the mechanisms regulating the nucleotide usage in mammalian genes, we analyzed the sequences of three physically linked Hsp70 paralogs in human and mouse. We report that the sequences of HSPA1A and HSPA1B genes are almost identical, whereas the HSPA1L gene contains some regions very similar to HSPA1A/B and some regions with much higher divergence. Phylogenetic analysis reveals that gene conversion has homogenized the entire coding regions of HSPA1A/B and several fragments of HSPA1L. The regions undergoing conversion are all very GC rich, contrarily to the regions not subject to conversion. The pattern of nucleotide substitution in mammalian orthologs suggests that the mechanism increasing the GC content is still functioning. To test the possibility that the high GC content facilitates the expression of Hsp70 during heat-shock, we performed in vitro translation experiments. We failed to detect any effect of GC content on the translation efficiency at high temperatures. Taken together, our data strongly support the biased gene conversion hypothesis of GC-content evolution.

Hsp90 chaperones wild-type p53 tumor suppressor protein.

Immortalized human fibroblasts were used to investigate the putative interactions of the Hsp90 molecular chaperone with the wild-type p53 tumor suppressor protein. We show that geldanamycin or radicicol, specific inhibitors of Hsp90, diminish specific wild-type p53 binding to the p21 promoter sequence. Consequently, these inhibitors decrease p21 mRNA levels, which lead to a reduction in cellular p21/Waf1 protein, known to induce cell cycle arrest. In control experiments, we show that neither geldanamycin nor radicicol affect p53 mRNA levels. A minor decrease in p53 protein level following the treatment of human fibroblasts with the inhibitors suggests the potential involvement of Hsp90 in the stabilization of wild-type p53. To support our in vivo findings, we used a reconstituted system with highly purified recombinant proteins to examine the effects of Hsp90 on wild-type p53 binding to the p21 promoter sequence. The human recombinant Hsp90 alpha-isoform as well as bovine brain Hsp90 were purified to homogeneity. Both of these molecular chaperones displayed ATPase activity and the ability to refold heat-inactivated luciferase in a geldanamycin- and radicicol-sensitive manner, suggesting that post-translational modifications are not involved in the modulation of Hsp90alpha activity. We show that the incubation of recombinant p53 at 37 degrees C decreases the level of its wild-type conformation and strongly inhibits the in vitro binding of p53 to the p21 promoter sequence. Interestingly, Hsp90 in an ATP-dependent manner can positively modulate p53 DNA binding after incubation at physiological temperature of 37 degrees C. Other recombinant human chaperones from Hsp70 and Hsp40 families were not able to efficiently substitute Hsp90 in this reaction. Consistent with our in vivo results, geldanamycin can suppress Hsp90 ability to regulate in vitro p53 DNA binding to the promoter sequence. In summary, the results presented in this article state that chaperone activity of Hsp90 is important for the transcriptional activity of genotypically wild-type p53.