Consequences of depleting TNRC6, AGO, and DROSHA proteins on expression of microRNAs.

The potential for microRNAs (miRNAs) to regulate gene expression remains incompletely understood. DROSHA initiates the biogenesis of miRNAs while variants of Argonaute (AGO) and trinucleotide repeat containing six (TNRC6) family proteins form complexes with miRNAs to facilitate RNA recognition and gene regulation. Here we investigate the fate of miRNAs in the absence of these critical RNAi protein factors. Knockout of DROSHA expression reduces levels of some miRNAs annotated in miRBase but not others. The identity of miRNAs with reduced expression matches the identity of miRNAs previously identified by experimental approaches. The MirGeneDB resource offers the closest alignment with experimental results. In contrast, the loss of TNRC6 proteins had much smaller effects on miRNA levels. Knocking out AGO proteins, which directly contact the mature miRNA, decreased expression of the miRNAs most strongly associated with AGO2 as determined from enhanced crosslinking immunoprecipitation (AGO2-eCLIP). Evaluation of miRNA binding to endogenously expressed AGO proteins revealed that miRNA:AGO association was similar for AGO1, AGO2, AGO3, and AGO4. Our data emphasize the need to evaluate annotated miRNAs based on approximate cellular abundance, DROSHA-dependence, and physical association with AGO when forming hypotheses related to their function.

Argonaute binding within human nuclear RNA and its impact on alternative splicing.

Mammalian RNA interference (RNAi) is often linked to the regulation of gene expression in the cytoplasm. Synthetic RNAs, however, can also act through the RNAi pathway to regulate transcription and splicing. While nuclear regulation by synthetic RNAs can be robust, a critical unanswered question is whether endogenous functions for nuclear RNAi exist in mammalian cells. Using enhanced crosslinking immunoprecipitation (eCLIP) in combination with RNA sequencing (RNA-seq) and multiple AGO knockout cell lines, we mapped AGO2 protein binding sites within nuclear RNA. The strongest AGO2 binding sites were mapped to micro RNAs (miRNAs). The most abundant miRNAs were distributed similarly between the cytoplasm and nucleus, providing no evidence for mechanisms that facilitate localization of miRNAs in one compartment versus the other. Beyond miRNAs, most statistically significant AGO2 binding was within introns. Splicing changes were confirmed by RT-PCR and recapitulated by synthetic miRNA mimics complementary to the sites of AGO2 binding. These data support the hypothesis that miRNAs can control gene splicing. While nuclear RNAi proteins have the potential to be natural regulatory mechanisms, careful study will be necessary to identify critical RNA drivers of normal physiology and disease.

Impact of scaffolding protein TNRC6 paralogs on gene expression and splicing.

TNRC6 is a scaffolding protein that bridges interactions between small RNAs, argonaute (AGO) protein, and effector proteins to control gene expression. There are three paralogs in mammalian cells, TNRC6A, TNRC6B, and TNRC6C These paralogs have ∼40% amino acid sequence identity and the extent of their unique or redundant functions is unclear. Here, we use knockout cell lines, enhanced crosslinking immunoprecipitation (eCLIP), and high-throughput RNA sequencing (RNA-seq) to explore the roles of TNRC6 paralogs in RNA-mediated control of gene expression. We find that the paralogs are largely functionally redundant and changes in levels of gene expression are well-correlated with those observed in AGO knockout cell lines. Splicing changes observed in AGO knockout cell lines are also observed in TNRC6 knockout cells. These data further define the roles of the TNRC6 isoforms as part of the RNA interference (RNAi) machinery.

Argonaute binding within 3'-untranslated regions poorly predicts gene repression.

Despite two decades of study, the full scope of RNAi in mammalian cells has remained obscure. Here we combine: (i) Knockout of argonaute (AGO) variants; (ii) RNA sequencing analysis of gene expression changes and (iii) Enhanced Crosslinking Immunoprecipitation Sequencing (eCLIP-seq) using anti-AGO2 antibody to identify potential microRNA (miRNA) binding sites. We find that knocking out AGO1, AGO2 and AGO3 together are necessary to achieve full impact on steady state levels of mRNA. eCLIP-seq located AGO2 protein associations within 3'-untranslated regions. The standard mechanism of miRNA action would suggest that these associations should repress gene expression. Contrary to this expectation, associations between AGO and RNA are poorly correlated with gene repression in wild-type versus knockout cells. Many clusters are associated with increased steady state levels of mRNA in wild-type versus knock out cells, including the strongest cluster within the MYC 3'-UTR. Our results suggest that assumptions about miRNA action should be re-examined.

Analyzing pre-symptomatic tissue to gain insights into the molecular and mechanistic origins of late-onset degenerative trinucleotide repeat disease.

How genetic defects trigger the molecular changes that cause late-onset disease is important for understanding disease progression and therapeutic development. Fuchs' endothelial corneal dystrophy (FECD) is an RNA-mediated disease caused by a trinucleotide CTG expansion in an intron within the TCF4 gene. The mutant intronic CUG RNA is present at one-two copies per cell, posing a challenge to understand how a rare RNA can cause disease. Late-onset FECD is a uniquely advantageous model for studying how RNA triggers disease because: (i) Affected tissue is routinely removed during surgery; (ii) The expanded CTG mutation is one of the most prevalent disease-causing mutations, making it possible to obtain pre-symptomatic tissue from eye bank donors to probe how gene expression changes precede disease; and (iii) The affected tissue is a homogeneous single cell monolayer, facilitating accurate transcriptome analysis. Here, we use RNA sequencing (RNAseq) to compare tissue from individuals who are pre-symptomatic (Pre_S) to tissue from patients with late stage FECD (FECD_REP). The abundance of mutant repeat intronic RNA in Pre_S and FECD_REP tissue is elevated due to increased half-life in a corneal cells. In Pre_S tissue, changes in splicing and extracellular matrix gene expression foreshadow the changes observed in advanced disease and predict the activation of the fibrosis pathway and immune system seen in late-stage patients. The absolute magnitude of splicing changes is similar in pre-symptomatic and late stage tissue. Our data identify gene candidates for early drivers of disease and biomarkers that may represent diagnostic and therapeutic targets for FECD. We conclude that changes in alternative splicing and gene expression are observable decades prior to the diagnosis of late-onset trinucleotide repeat disease.

Human GW182 Paralogs Are the Central Organizers for RNA-Mediated Control of Transcription.

In the cytoplasm, small RNAs can control mammalian translation by regulating the stability of mRNA. In the nucleus, small RNAs can also control transcription and splicing. The mechanisms for RNA-mediated nuclear regulation are not understood and remain controversial, hindering the effective application of nuclear RNAi and investigation of its natural regulatory roles. Here, we reveal that the human GW182 paralogs TNRC6A/B/C are central organizing factors critical to RNA-mediated transcriptional activation. Mass spectrometry of purified nuclear lysates followed by experimental validation demonstrates that TNRC6A interacts with proteins involved in protein degradation, RNAi, the CCR4-NOT complex, the mediator complex, and histone-modifying complexes. Functional analysis implicates TNRC6A, NAT10, MED14, and WDR5 in RNA-mediated transcriptional activation. These findings describe protein complexes capable of bridging RNA-mediated sequence-specific recognition of noncoding RNA transcripts with the regulation of gene transcription.

Design and bioinformatics analysis of genome-wide CLIP experiments.

The past decades have witnessed a surge of discoveries revealing RNA regulation as a central player in cellular processes. RNAs are regulated by RNA-binding proteins (RBPs) at all post-transcriptional stages, including splicing, transportation, stabilization and translation. Defects in the functions of these RBPs underlie a broad spectrum of human pathologies. Systematic identification of RBP functional targets is among the key biomedical research questions and provides a new direction for drug discovery. The advent of cross-linking immunoprecipitation coupled with high-throughput sequencing (genome-wide CLIP) technology has recently enabled the investigation of genome-wide RBP-RNA binding at single base-pair resolution. This technology has evolved through the development of three distinct versions: HITS-CLIP, PAR-CLIP and iCLIP. Meanwhile, numerous bioinformatics pipelines for handling the genome-wide CLIP data have also been developed. In this review, we discuss the genome-wide CLIP technology and focus on bioinformatics analysis. Specifically, we compare the strengths and weaknesses, as well as the scopes, of various bioinformatics tools. To assist readers in choosing optimal procedures for their analysis, we also review experimental design and procedures that affect bioinformatics analyses.

Intramolecular circularization increases efficiency of RNA sequencing and enables CLIP-Seq of nuclear RNA from human cells.

RNA sequencing (RNA-Seq) is a powerful tool for analyzing the identity of cellular RNAs but is often limited by the amount of material available for analysis. In spite of extensive efforts employing existing protocols, we observed that it was not possible to obtain useful sequencing libraries from nuclear RNA derived from cultured human cells after crosslinking and immunoprecipitation (CLIP). Here, we report a method for obtaining strand-specific small RNA libraries for RNA sequencing that requires picograms of RNA. We employ an intramolecular circularization step that increases the efficiency of library preparation and avoids the need for intermolecular ligations of adaptor sequences. Other key features include random priming for full-length cDNA synthesis and gel-free library purification. Using our method, we generated CLIP-Seq libraries from nuclear RNA that had been UV-crosslinked and immunoprecipitated with anti-Argonaute 2 (Ago2) antibody. Computational protocols were developed to enable analysis of raw sequencing data and we observe substantial differences between recognition by Ago2 of RNA species in the nucleus relative to the cytoplasm. This RNA self-circularization approach to RNA sequencing (RC-Seq) allows data to be obtained using small amounts of input RNA that cannot be sequenced by standard methods.

RNAi factors are present and active in human cell nuclei.

RNAi is widely appreciated as a powerful regulator of mRNA translation in the cytoplasm of mammalian cells. However, the presence and activity of RNAi factors in the mammalian nucleus has been the subject of considerable debate. Here, we show that Argonaute-2 (Ago2) and RNAi factors Dicer, TRBP, and TRNC6A/GW182 are in the human nucleus and associate together in multiprotein complexes. Small RNAs can silence nuclear RNA and guide site-specific cleavage of the targeted RNA, demonstrating that RNAi can function in the human nucleus. Nuclear Dicer is active and miRNAs are bound to nuclear Ago2, consistent with the existence of nuclear miRNA pathways. Notably, we do not detect loading of duplex small RNAs in nuclear extracts and known loading factors are absent. These results extend RNAi into the mammalian nucleus and suggest that regulation of RNAi via small RNA loading of Ago2 differs between the cytoplasm and the nucleus.

Promoter RNA links transcriptional regulation of inflammatory pathway genes.

Although many long non-coding RNAs (lncRNAs) have been discovered, their function and their association with RNAi factors in the nucleus have remained obscure. Here, we identify RNA transcripts that overlap the cyclooxygenase-2 (COX-2) promoter and contain two adjacent binding sites for an endogenous miRNA, miR-589. We find that miR-589 binds the promoter RNA and activates COX-2 transcription. In addition to miR-589, fully complementary duplex RNAs that target the COX-2 promoter transcript activate COX-2 transcription. Activation by small RNA requires RNAi factors argonaute-2 (AGO2) and GW182, but does not require AGO2-mediated cleavage of the promoter RNA. Instead, the promoter RNA functions as a scaffold. Binding of AGO2 protein/small RNA complexes to the promoter RNA triggers gene activation. Gene looping allows interactions between the promoters of COX-2 and phospholipase A2 (PLA2G4A), an adjacent pro-inflammatory pathway gene that produces arachidonic acid, the substrate for COX-2 protein. miR-589 and fully complementary small RNAs regulate both COX-2 and PLA2G4A gene expression, revealing an unexpected connection between key steps of the eicosanoid signaling pathway. The work demonstrates the potential for RNA to coordinate locus-dependent assembly of related genes to form functional operons through cis-looping.

Mechanism of allele-selective inhibition of huntingtin expression by duplex RNAs that target CAG repeats: function through the RNAi pathway.

Huntington's disease is an incurable neurodegenerative disorder caused by expansion of a CAG trinucleotide repeat within one allele of the huntingtin (HTT) gene. Agents that block expression of mutant HTT and preserve expression of wild-type HTT target the cause of the disease and are an alternative for therapy. We have previously demonstrated that mismatch-containing duplex RNAs complementary to the expanded trinucleotide repeat are potent and allele-selective inhibitors of mutant HTT expression, but the mechanism of allele selectivity was not explored. We now report that anti-CAG duplex RNA preferentially recruits argonaute 2 (AGO2) to mutant rather than wild-type HTT mRNA. Efficient inhibition of mutant HTT protein expression requires less AGO2 than needed for inhibiting wild-type expression. In contrast, inhibiting the expression of mutant HTT protein is highly sensitive to reduced expression of GW182 (TNRC6A) and its two paralogs, a protein family associated with miRNA action. Allele-selective inhibition may involve cooperative binding of multiple protein-RNA complexes to the expanded repeat. These data suggest that allele-selective inhibition proceeds through an RNA interference pathway similar to that used by miRNAs and that discrimination between mutant and wild-type alleles of HTT mRNA is highly sensitive to the pool of AGO2 and GW182 family proteins inside cells.

Transcriptional silencing by hairpin RNAs complementary to a gene promoter.

Double-stranded RNAs can target gene promoters and inhibit transcription. To date, most research has focused on synthetic RNA duplexes. Transcriptional silencing by hairpin RNAs would facilitate a better understanding of endogenous RNA-mediated regulation of transcription within cells. Here we examine transcriptional silencing of progesterone receptor (PR) expression by hairpin RNAs. We identify the guide strand as the strand complementary to an antisense transcript at the PR promoter and that hairpin RNAs are active transcriptional silencing agents. The sequence of the hairpin loop affects activity, with the highest activity achieved when the loop has the potential for full complementarity to the antisense transcript target. Introduction of centrally mismatched bases relative to the target transcript does not prevent transcriptional silencing unless the mismatches are present on both the guide and passenger strands. These data demonstrate that hairpin RNAs can cause transcriptional silencing and offer insights into the mechanism of gene modulation by RNAs that target gene promoters.

Involvement of argonaute proteins in gene silencing and activation by RNAs complementary to a non-coding transcript at the progesterone receptor promoter.

Double-stranded RNAs that are complementary to non-coding transcripts at gene promoters can activate or inhibit gene expression in mammalian cells. Understanding the mechanism for modulating gene expression by promoter-targeted antigene RNAs (agRNAs) will require identification of the proteins involved in recognition. Previous reports have implicated argonaute (AGO) proteins, but identifications have differed with involvement of AGO1, AGO2, or both AGO1 and AGO2 being reported by different studies. The roles of AGO3 and AGO4 have not been investigated. Here, we examine the role of AGO 1-4 in gene silencing and activation of the progesterone receptor (PR) gene. Expression of AGO2 is necessary for efficient gene silencing or activation and AGO2 is recruited to the non-coding transcript that overlaps the promoter during both gene silencing and activation. Expression of AGO1, AGO3 and AGO4 are not necessary for gene silencing or activation nor are AGO1, AGO3, or AGO4 recruited to the target non-coding transcript during gene activation. These data indicate that AGO2 is the primary AGO variant involved in modulating expression of PR by agRNAs.

Transcriptional regulation by small RNAs at sequences downstream from 3' gene termini.

Transcriptome studies reveal many noncoding transcripts overlapping 3' gene termini. The function of these transcripts is unknown. Here we have characterized transcription at the progesterone receptor (PR) locus and identified noncoding transcripts that overlap the 3' end of the gene. Small RNAs complementary to sequences beyond the 3' terminus of PR mRNA modulated expression of PR, recruited argonaute 2 to a 3' noncoding transcript, altered occupancy of RNA polymerase II, induced chromatin changes at the PR promoter and affected responses to physiological stimuli. We found that the promoter and 3' terminal regions of the PR locus are in close proximity, providing a potential mechanism for RNA-mediated control of transcription over long genomic distances. These results extend the potential for small RNAs to regulate transcription to target sequences beyond the 3' termini of mRNA.

A pseudocatenane structure formed between DNA and A cyclic bisintercalator.

Targeting double-stranded DNA with small molecules remains an active area of basic research. Herein is described a cyclic DNA bisintercalator that is based on two naphthalene diimide (NDI) intercalating units tethered by one linking element specific for binding in the minor groove and the other linking element specific for binding in the major groove. DNase I footprinting revealed a strong preference for binding the sequence 5'-GGTACC-3'. NMR structural studies of the complex with d(CGGTACCG)(2) verified a pseudocatenane structure in which the NDI units reside four base pairs apart, with one linker segment located in the minor groove and the other in the major groove consistent with the linker designs. To the best of our knowledge, this is the first structurally well-characterized pseudocatenane complex between a sequence specific cyclic bisintercalator and intact DNA.

Structural characterization of a rigidified threading bisintercalator.

NMR spectroscopy was used to explore the sequence-specific interaction of DNA with a new threading bisintercalator (C1) consisting of two intercalating 1,4,5,8-naphthalenetetracarboxylic diimide (NDI) units connected by a rigid, tricyclic spiro linker. A structural model of C1 complexed to d(CGGTACCG)(2) was calculated using distance constraints obtained from solution NMR data. The model was also supported by the results from residual dipolar coupling (RDC) measurements obtained using Pf1-phage as a cosolvent. According to the model, the central cyclohexane ring of the linker connecting the two NDI units lies flat in the minor groove of DNA. Linker length, hydrogen bonding, steric, and hydrophobic factors all appear to contribute to the observed sequence specificity of binding. These results serve to illustrate the versatility of threading polyintercalation given that, in a previous study, a ligand consisiting of two NDI units joined by a flexible peptide linker was shown to bind sequence specifically within the major groove of this same sequence of DNA.

Screening of threading bis-intercalators binding to duplex DNA by electrospray ionization tandem mass spectrometry.

The DNA binding of novel threading bis-intercalators V1, trans-D1, and cis-C1, which contain two naphthalene diimide (NDI) intercalation units connected by a scaffold, was evaluated using electrospray ionization mass spectrometry (ESI-MS) and DNAse footprinting techniques. ESI-MS experiments confirmed that V1, the ligand containing the -Gly3-Lys- peptide scaffold, binds to a DNA duplex containing the 5'-GGTACC-3' specific binding site identified in previous NMR-based studies. The ligand formed complexes with a ligand/DNA binding stoichiometry of 1:1, even when there was excess ligand in solution. Trans-D1 and cis-C1 are new ligands containing a rigid spiro-tricyclic scaffold in the trans- and cis- orientations, respectively. Preliminary DNAse footprinting experiments identified possible specific binding sites of 5'-CAGTGA-5' for trans-D1 and 5'-GGTACC-3' for cis-C1. ESI-MS experiments revealed that both ligands bound to DNA duplexes containing the respective specific binding sequences, with cis-C1 exhibiting the most extensive binding based on a higher fraction of bound DNA value. Cis-C1 formed complexes with a dominant 1:1 binding stoichiometry, whereas trans-D1 was able to form 2:1 complexes at ligand/DNA molar ratios >or=1 which is suggestive of nonspecific binding. Collisional activated dissociation (CAD) experiments indicate that DNA complexes containing V1, trans-D1, and cis-C1 have a unique fragmentation pathway, which was also observed for complexes containing the commercially available bis-intercalator echinomycin, as a result of similar binding interactions, marked by intercalation in addition to hydrogen bonding by the scaffold with the DNA major or minor groove.