Analysis of Circular RNAs Using the Web Tool CircInteractome.

Circular RNAs (circRNAs) are generated through nonlinear back splicing, during which the 5' and 3' ends are covalently joined. Consequently, the lack of free ends makes them very stable compared to their counterpart linear RNAs. By selectively interacting with microRNAs and RNA-binding proteins (RBPs), circRNAs have been shown to influence gene expression programs. We designed a web tool, CircInteractome, in order to (1) explore potential interactions of circRNAs with RBPs, (2) design specific divergent primers to detect circRNAs, (3) study tissue- and cell-specific circRNAs, (4) identify gene-specific circRNAs, (5) explore potential miRNAs interacting with circRNAs, and (6) design specific siRNAs to silence circRNAs. Here, we review the CircInteractome tool and explain recent updates to the site. The database is freely accessible at http://circinteractome.nia.nih.gov .

WIG1 is crucial for AGO2-mediated ACOT7 mRNA silencing via miRNA-dependent and -independent mechanisms.

RNA-binding proteins (RBPs) are involved in mRNA splicing, maturation, transport, translation, storage and turnover. Here, we identified ACOT7 mRNA as a novel target of human WIG1. ACOT7 mRNA decay was triggered by the microRNA miR-9 in a WIG1-dependent manner via classic recruitment of Argonaute 2 (AGO2). Interestingly, AGO2 was also recruited to ACOT7 mRNA in a WIG1-dependent manner in the absence of miR-9, which indicates an alternative model whereby WIG1 controls AGO2-mediated gene silencing. The WIG1-AGO2 complex attenuated translation initiation via an interaction with translation initiation factor 5B (eIF5B). These results were confirmed using a WIG1 tethering system based on the MS2 bacteriophage coat protein and a reporter construct containing an MS2-binding site, and by immunoprecipitation of WIG1 and detection of WIG1-associated proteins using liquid chromatography-tandem mass spectrometry. We also identified WIG1-binding motifs using photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation analyses. Altogether, our data indicate that WIG1 governs the miRNA-dependent and the miRNA-independent recruitment of AGO2 to lower the stability of and suppress the translation of ACOT7 mRNA.

High-purity circular RNA isolation method (RPAD) reveals vast collection of intronic circRNAs.

High-throughput RNA sequencing methods coupled with specialized bioinformatic analyses have recently uncovered tens of thousands of unique circular (circ)RNAs, but their complete sequences, genes of origin and functions are largely unknown. Given that circRNAs lack free ends and are thus relatively stable, their association with microRNAs (miRNAs) and RNA-binding proteins (RBPs) can influence gene expression programs. While exoribonuclease treatment is widely used to degrade linear RNAs and enrich circRNAs in RNA samples, it does not efficiently eliminate all linear RNAs. Here, we describe a novel method for the isolation of highly pure circRNA populations involving RNase R treatment followed by Polyadenylation and poly(A)+ RNA Depletion (RPAD), which removes linear RNA to near completion. High-throughput sequencing of RNA prepared using RPAD from human cervical carcinoma HeLa cells and mouse C2C12 myoblasts led to two surprising discoveries: (i) many exonic circRNA (EcircRNA) isoforms share an identical backsplice sequence but have different body sizes and sequences, and (ii) thousands of novel intronic circular RNAs (IcircRNAs) are expressed in cells. In sum, isolating high-purity circRNAs using the RPAD method can enable quantitative and qualitative analyses of circRNA types and sequence composition, paving the way for the elucidation of circRNA functions.

Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1.

HuR influences gene expression programs and hence cellular phenotypes by binding to hundreds of coding and noncoding linear RNAs. However, whether HuR binds to circular RNAs (circRNAs) and impacts on their function is unknown. Here, we have identified en masse circRNAs binding HuR in human cervical carcinoma HeLa cells. One of the most prominent HuR target circRNAs was hsa_circ_0031288, renamed CircPABPN1 as it arises from the PABPN1 pre-mRNA. Further analysis revealed that HuR did not influence CircPABPN1 abundance; interestingly, however, high levels of CircPABPN1 suppressed HuR binding to PABPN1 mRNA. Evaluation of PABPN1 mRNA polysomes indicated that PABPN1 translation was modulated positively by HuR and hence negatively by CircPABPN1. We propose that the extensive binding of CircPABPN1 to HuR prevents HuR binding to PABPN1 mRNA and lowers PABPN1 translation, providing the first example of competition between a circRNA and its cognate mRNA for an RBP that affects translation.

HuR and GRSF1 modulate the nuclear export and mitochondrial localization of the lncRNA RMRP.

Some mitochondrial long noncoding RNAs (lncRNAs) are encoded by nuclear DNA, but the mechanisms that mediate their transport to mitochondria are poorly characterized. Using affinity RNA pull-down followed by mass spectrometry analysis, we found two RNA-binding proteins (RBPs), HuR (human antigen R) and GRSF1 (G-rich RNA sequence-binding factor 1), that associated with the nuclear DNA-encoded lncRNA RMRP and mobilized it to mitochondria. In cultured human cells, HuR bound RMRP in the nucleus and mediated its CRM1 (chromosome region maintenance 1)-dependent export to the cytosol. After RMRP was imported into mitochondria, GRSF1 bound RMRP and increased its abundance in the matrix. Loss of GRSF1 lowered the mitochondrial levels of RMRP, in turn suppressing oxygen consumption rates and modestly reducing mitochondrial DNA replication priming. Our findings indicate that RBPs HuR and GRSF1 govern the cytoplasmic and mitochondrial localization of the lncRNA RMRP, which is encoded by nuclear DNA but has key functions in mitochondria.

Generation and gene expression profiling of 48 transcription-factor-inducible mouse embryonic stem cell lines.

Mouse embryonic stem cells (ESCs) can differentiate into a wide range - and possibly all cell types in vitro, and thus provide an ideal platform to study systematically the action of transcription factors (TFs) in cell differentiation. Previously, we have generated and analyzed 137 TF-inducible mouse ESC lines. As an extension of this "NIA Mouse ESC Bank," we generated and characterized 48 additional mouse ESC lines, in which single TFs in each line could be induced in a doxycycline-controllable manner. Together, with the previous ESC lines, the bank now comprises 185 TF-manipulable ESC lines (>10% of all mouse TFs). Global gene expression (transcriptome) profiling revealed that the induction of individual TFs in mouse ESCs for 48 hours shifts their transcriptomes toward specific differentiation fates (e.g., neural lineages by Myt1 Isl1, and St18; mesodermal lineages by Pitx1, Pitx2, Barhl2, and Lmx1a; white blood cells by Myb, Etv2, and Tbx6, and ovary by Pitx1, Pitx2, and Dmrtc2). These data also provide and lists of inferred target genes of each TF and possible functions of these TFs. The results demonstrate the utility of mouse ESC lines and their transcriptome data for understanding the mechanism of cell differentiation and the function of TFs.

Novel RNA-binding activity of MYF5 enhances Ccnd1/Cyclin D1 mRNA translation during myogenesis.

Skeletal muscle contains long multinucleated and contractile structures known as muscle fibers, which arise from the fusion of myoblasts into multinucleated myotubes during myogenesis. The myogenic regulatory factor (MRF) MYF5 is the earliest to be expressed during myogenesis and functions as a transcription factor in muscle progenitor cells (satellite cells) and myocytes. In mouse C2C12 myocytes, MYF5 is implicated in the initial steps of myoblast differentiation into myotubes. Here, using ribonucleoprotein immunoprecipitation (RIP) analysis, we discovered a novel function for MYF5 as an RNA-binding protein which associated with a subset of myoblast mRNAs. One prominent MYF5 target was Ccnd1 mRNA, which encodes the key cell cycle regulator CCND1 (Cyclin D1). Biotin-RNA pulldown, UV-crosslinking and gel shift experiments indicated that MYF5 was capable of binding the 3' untranslated region (UTR) and the coding region (CR) of Ccnd1 mRNA. Silencing MYF5 expression in proliferating myoblasts revealed that MYF5 promoted CCND1 translation and modestly increased transcription of Ccnd1 mRNA. Accordingly, overexpressing MYF5 in C2C12 cells upregulated CCND1 expression while silencing MYF5 reduced myoblast proliferation as well as differentiation of myoblasts into myotubes. Moreover, MYF5 silencing reduced myogenesis, while ectopically restoring CCND1 abundance partially rescued the decrease in myogenesis seen after MYF5 silencing. We propose that MYF5 enhances early myogenesis in part by coordinately elevating Ccnd1 transcription and Ccnd1 mRNA translation.

CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs.

Circular RNAs (circRNAs) are widely expressed in animal cells, but their biogenesis and functions are poorly understood. CircRNAs have been shown to act as sponges for miRNAs and may also potentially sponge RNA-binding proteins (RBPs) and are thus predicted to function as robust posttranscriptional regulators of gene expression. The joint analysis of large-scale transcriptome data coupled with computational analyses represents a powerful approach to elucidate possible biological roles of ribonucleoprotein (RNP) complexes. Here, we present a new web tool, CircInteractome (circRNA interactome), for mapping RBP- and miRNA-binding sites on human circRNAs. CircInteractome searches public circRNA, miRNA, and RBP databases to provide bioinformatic analyses of binding sites on circRNAs and additionally analyzes miRNA and RBP sites on junction and junction-flanking sequences. CircInteractome also allows the user the ability to (1) identify potential circRNAs which can act as RBP sponges, (2) design junction-spanning primers for specific detection of circRNAs of interest, (3) design siRNAs for circRNA silencing, and (4) identify potential internal ribosomal entry sites (IRES). In sum, the web tool CircInteractome, freely accessible at http://circinteractome.nia.nih.gov, facilitates the analysis of circRNAs and circRNP biology.

Transient bursts of Zscan4 expression are accompanied by the rapid derepression of heterochromatin in mouse embryonic stem cells.

Mouse embryonic stem cells (mESCs) have a remarkable capacity to maintain normal genome stability and karyotype in culture. We previously showed that infrequent bursts of Zscan4 expression (Z4 events) are important for the maintenance of telomere length and genome stability in mESCs. However, the molecular details of Z4 events remain unclear. Here we show that Z4 events involve unexpected transcriptional derepression in heterochromatin regions that usually remain silent. During a Z4 event, we see rapid derepression and rerepression of heterochromatin leading to a burst of transcription that coincides with transient histone hyperacetylation and DNA demethylation, clustering of pericentromeric heterochromatin around the nucleolus, and accumulation of activating and repressive chromatin remodelling complexes. This heterochromatin-based transcriptional activity suggests that mESCs may maintain their extraordinary genome stability at least in part by transiently resetting their heterochromatin.

7SL RNA represses p53 translation by competing with HuR.

Noncoding RNAs (ncRNAs) and RNA-binding proteins are potent post-transcriptional regulators of gene expression. The ncRNA 7SL is upregulated in cancer cells, but its impact upon the phenotype of cancer cells is unknown. Here, we present evidence that 7SL forms a partial hybrid with the 3'-untranslated region (UTR) of TP53 mRNA, which encodes the tumor suppressor p53. The interaction of 7SL with TP53 mRNA reduced p53 translation, as determined by analyzing p53 expression levels, nascent p53 translation and TP53 mRNA association with polysomes. Silencing 7SL led to increased binding of HuR to TP53 mRNA, an interaction that led to the promotion of p53 translation and increased p53 abundance. We propose that the competition between 7SL and HuR for binding to TP53 3'UTR contributes to determining the magnitude of p53 translation, in turn affecting p53 levels and the growth-suppressive function of p53. Our findings suggest that targeting 7SL may be effective in the treatment of cancers with reduced p53 levels.

Chromatin properties of regulatory DNA probed by manipulation of transcription factors.

Transcription factors (TFs) bind to DNA and regulate the transcription of nearby genes. However, only a small fraction of TF binding sites have such regulatory effects. Here we search for the predictors of functional binding sites by carrying out a systematic computational screening of a variety of contextual factors (histone modifications, nuclear lamin-bindings, and cofactor bindings). We used regression analysis to test if contextual factors are associated with upregulation or downregulation of neighboring genes following the induction or knockdown of the 9 TFs in mouse embryonic stem (ES) cells. Functional TF binding sites appeared to be either active (i.e., bound by P300, CHD7, mediator, cohesin, and SWI/SNF) or repressed (i.e., with H3K27me3 histone marks and bound by Polycomb factors). Active binding sites mediated the downregulation of nearby genes upon knocking down the activating TFs or inducing repressors. Repressed TF binding sites mediated the upregulation of nearby genes (e.g., poised developmental regulators) upon inducing TFs. In addition, repressed binding sites mediated repressive effects of TFs, identified by the downregulation of target genes after the induction of TFs or by the upregulation of target genes after the knockdown of TFs. The contextual factors associated with functions of DNA-bound TFs were used to improve the identification of candidate target genes regulated by TFs.

SHAVE: shrinkage estimator measured for multiple visits increases power in GWAS of quantitative traits.

Measurement error and biological variability generate distortions in quantitative phenotypic data. In longitudinal studies with repeated measurements, the multiple measurements provide a route to reduce noise and correspondingly increase the strength of signals in genome-wide association studies (GWAS).To optimize noise correction, we have developed Shrunken Average (SHAVE), an approach using a Bayesian Shrinkage estimator. This estimator uses regression toward the mean for every individual as a function of (1) their average across visits; (2) their number of visits; and (3) the correlation between visits. Computer simulations support an increase in power, with results very similar to those expected by the assumptions of the model. The method was applied to a real data set for 14 anthropomorphic traits in ∼6000 individuals enrolled in the SardiNIA project, with up to three visits (measurements) for each participant. Results show that additional measurements have a large impact on the strength of GWAS signals, especially when participants have different number of visits, with SHAVE showing a clear increase in power relative to single visits. In addition, we have derived a relation to assess the improvement in power as a function of number of visits and correlation between visits. It can also be applied in the optimization of experimental designs or usage of measuring devices. SHAVE is fast and easy to run, written in R and freely available online.

Generation of mouse ES cell lines engineered for the forced induction of transcription factors.

Here we report the generation and characterization of 84 mouse ES cell lines with doxycycline-controllable transcription factors (TFs) which, together with the previous 53 lines, cover 7-10% of all TFs encoded in the mouse genome. Global gene expression profiles of all 137 lines after the induction of TFs for 48 hrs can associate each TF with the direction of ES cell differentiation, regulatory pathways, and mouse phenotypes. These cell lines and microarray data provide building blocks for a variety of future biomedical research applications as a community resource.

Prediction of evolutionarily conserved interologs in Mus musculus.

BACKGROUND: Identification of protein-protein interactions is an important first step to understand living systems. High-throughput experimental approaches have accumulated large amount of information on protein-protein interactions in human and other model organisms. Such interaction information has been successfully transferred to other species, in which the experimental data are limited. However, the annotation transfer method could yield false positive interologs due to the lack of conservation of interactions when applied to phylogenetically distant organisms. RESULTS: To address this issue, we used phylogenetic profile method to filter false positives in interologs based on the notion that evolutionary conserved interactions show similar patterns of occurrence along the genomes. The approach was applied to Mus musculus, in which the experimentally identified interactions are limited. We first inferred the protein-protein interactions in Mus musculus by using two approaches: i) identifying mouse orthologs of interacting proteins (interologs) based on the experimental protein-protein interaction data from other organisms; and ii) analyzing frequency of mouse ortholog co-occurrence in predicted operons of bacteria. We then filtered possible false-positives in the predicted interactions using the phylogenetic profiles. We found that this filtering method significantly increased the frequency of interacting protein-pairs coexpressed in the same cells/tissues in gene expression omnibus (GEO) database as well as the frequency of interacting protein-pairs shared the similar Gene Ontology (GO) terms for biological processes and cellular localizations. The data supports the notion that phylogenetic profile helps to reduce the number of false positives in interologs. CONCLUSION: We have developed protein-protein interaction database in mouse, which contains 41109 interologs. We have also developed a web interface to facilitate the use of database http://lgsun.grc.nia.nih.gov/mppi/.

CisView: a browser and database of cis-regulatory modules predicted in the mouse genome.

To facilitate the analysis of gene regulatory regions of the mouse genome, we developed a CisView (http://lgsun.grc.nia.nih.gov/cisview), a browser and database of genome-wide potential transcription factor binding sites (TFBSs) that were identified using 134 position-weight matrices and 219 sequence patterns from various sources and were presented with the information about sequence conservation, neighboring genes and their structures, GO annotations, protein domains, DNA repeats and CpG islands. Analysis of the distribution of TFBSs revealed that many TFBSs (N = 145) were over-represented near transcription start sites. We also identified potential cis-regulatory modules (CRMs) defined as clusters of conserved TFBSs in the entire mouse genome. Out of 739 074 CRMs, 157 442 had a significantly higher regulatory potential score than semi-random sequences generated with a 3rd-order Markov process. The CisView browser provides a user-friendly computer environment for studying transcription regulation on a whole-genome scale and can also be used for interpreting microarray experiments and identifying putative targets of transcription factors.

Transcript copy number estimation using a mouse whole-genome oligonucleotide microarray.

The ability to quantitatively measure the expression of all genes in a given tissue or cell with a single assay is an exciting promise of gene-expression profiling technology. An in situ-synthesized 60-mer oligonucleotide microarray designed to detect transcripts from all mouse genes was validated, as well as a set of exogenous RNA controls derived from the yeast genome (made freely available without restriction), which allow quantitative estimation of absolute endogenous transcript abundance.

Genome-wide assembly and analysis of alternative transcripts in mouse.

To build a mouse gene index with the most comprehensive coverage of alternative transcription/splicing (ATS), we developed an algorithm and a fully automated computational pipeline for transcript assembly from expressed sequences aligned to the genome. We identified 191,946 genomic loci, which included 27,497 protein-coding genes and 11,906 additional gene candidates (e.g., nonprotein-coding, but multiexon). Comparison of the resulting gene index with TIGR, UniGene, DoTS, and ESTGenes databases revealed that it had a greater number of transcripts, a greater average number of exons and introns with proper splicing sites per gene, and longer ORFs. The 27,497 protein-coding genes had 77,138 transcripts, i.e., 2.8 transcripts per gene on average. Close examination of transcripts led to a combinatorial table of 23 types of ATS units, only nine of which were previously described, i.e., 14 types of alternative splicing, seven types of alternative starts, and two types of alternative termination. The 47%, 18%, and 14% of 20,323 multiexon protein-coding genes with proper splice sites had alternative splicings, alternative starts, and alternative terminations, respectively. The gene index with the comprehensive ATS will provide a useful platform for analyzing the nature and mechanism of ATS, as well as for designing the accurate exon-based DNA microarrays. The sequence data from this study have been submitted to GenBank under accession numbers: CK329321-CK334090; CF891695-CF906652; CF906741-CF916750; CK334091-CK347104; CK387035-CK393993; CN660032-CN690720; CN690721-CN725493.

A web-based tool for principal component and significance analysis of microarray data.

UNLABELLED: We have developed a program for microarray data analysis, which features the false discovery rate for testing statistical significance and the principal component analysis using the singular value decomposition method for detecting the global trends of gene-expression patterns. Additional features include analysis of variance with multiple methods for error variance adjustment, correction of cross-channel correlation for two-color microarrays, identification of genes specific to each cluster of tissue samples, biplot of tissues and corresponding tissue-specific genes, clustering of genes that are correlated with each principal component (PC), three-dimensional graphics based on virtual reality modeling language and sharing of PC between different experiments. The software also supports parameter adjustment, gene search and graphical output of results. The software is implemented as a web tool and thus the speed of analysis does not depend on the power of a client computer. AVAILABILITY: The tool can be used on-line or downloaded at http://lgsun.grc.nia.nih.gov/ANOVA/

The status, quality, and expansion of the NIH full-length cDNA project: the Mammalian Gene Collection (MGC).

The National Institutes of Health's Mammalian Gene Collection (MGC) project was designed to generate and sequence a publicly accessible cDNA resource containing a complete open reading frame (ORF) for every human and mouse gene. The project initially used a random strategy to select clones from a large number of cDNA libraries from diverse tissues. Candidate clones were chosen based on 5'-EST sequences, and then fully sequenced to high accuracy and analyzed by algorithms developed for this project. Currently, more than 11,000 human and 10,000 mouse genes are represented in MGC by at least one clone with a full ORF. The random selection approach is now reaching a saturation point, and a transition to protocols targeted at the missing transcripts is now required to complete the mouse and human collections. Comparison of the sequence of the MGC clones to reference genome sequences reveals that most cDNA clones are of very high sequence quality, although it is likely that some cDNAs may carry missense variants as a consequence of experimental artifact, such as PCR, cloning, or reverse transcriptase errors. Recently, a rat cDNA component was added to the project, and ongoing frog (Xenopus) and zebrafish (Danio) cDNA projects were expanded to take advantage of the high-throughput MGC pipeline.

Age-associated alteration of gene expression patterns in mouse oocytes.

Decreasing oocyte competence with maternal aging is a major factor in human infertility. To investigate the age-dependent molecular changes in a mouse model, we compared the expression profiles of metaphase II oocytes collected from 5- to 6-week-old mice with those collected from 42- to 45-week-old mice using the NIA 22K 60-mer oligo microarray. Among approximately 11,000 genes whose transcripts were detected in oocytes, about 5% (530) showed statistically significant expression changes, excluding the possibility of global decline in transcript abundance. Consistent with the generally accepted view of aging, the differentially expressed genes included ones involved in mitochondrial function and oxidative stress. However, the expression of other genes involved in chromatin structure, DNA methylation, genome stability and RNA helicases was also altered, suggesting the existence of additional mechanisms for aging. Among the transcripts decreased with aging, we identified and characterized a group of new oocyte-specific genes, members of the human NACHT, leucine-rich repeat and PYD-containing (NALP) gene family. These results have implications for aging research as well as for clinical ooplasmic donation to rejuvenate aging oocytes.