Search

Despite the key role of the human ribosome in protein biosynthesis, little is known about the extent of sequence variation in ribosomal DNA (rDNA) or its pre-rRNA and rRNA products. We recovered ribosomal DNA segments from a single human chromosome 21 using transformation-associated recombination (TAR) cloning in yeast. Accurate long-read sequencing of 13 isolates covering ∼0.82 Mb of the chromosome 21 rDNA complement revealed substantial variation among tandem repeat rDNA copies, several palindromic structures and potential errors in the previous reference sequence. These clones revealed 101 variant positions in the 45S transcription unit and 235 in the intergenic spacer sequence. Approximately 60% of the 45S variants were confirmed in independent whole-genome or RNA-seq data, with 47 of these further observed in mature 18S/28S rRNA sequences. TAR cloning and long-read sequencing enabled the accurate reconstruction of multiple rDNA units and a new, high-quality 44 838 bp rDNA reference sequence, which we have annotated with variants detected from chromosome 21 of a single individual. The large number of variants observed reveal heterogeneity in human rDNA, opening up the possibility of corresponding variations in ribosome dynamics.

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.

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.

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.

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 .

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.

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-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.

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.

Retinoic acid (RA) is one of the most potent inducers of differentiation of mouse embryonic stem cells (ESCs). However, previous studies show that RA treatment of cells cultured in the presence of a leukemia inhibitory factor (LIF) also result in the upregulation of a gene called Zscan4, whose transient expression is a marker for undifferentiated ESCs. We explored the balance between these two seemingly antagonistic effects of RA. ESCs indeed differentiated in the presence of LIF after RA treatment, but colonies of undifferentiated ESCs eventually emerged from these differentiated cells - even in the presence of RA. These colonies, named secondary colonies, consist of three cell types: typical undifferentiated ESCs expressing pluripotency genes such as Pou5f1, Sox2, and Nanog; cells expressing Zscan4; and endodermal-like cells located at the periphery of the colony. The capacity to form secondary colonies was confirmed for all eight tested ESC lines. Cells from the secondary colonies - after transfer to the standard ESC medium - retained pluripotency, judged by their strong alkaline phosphatase (ALP) staining, typical colony morphology, gene expression profile, stable karyotype, capacity to differentiate into all three germ layers in embryoid body formation assays, and successful contribution to chimeras after injection into blastocysts. Based on flow cytometry analysis (FACS), the proportion of Zscan4-positive cells in secondary colonies was higher than in standard ESC colonies, which may explain the capacity of ESCs to resist the differentiating effects of RA and instead form secondary colonies of undifferentiated ESCs. This hypothesis is supported by cell-lineage tracing analysis, which showed that most cells in the secondary colonies were descendents of cells transiently expressing Zscan4.

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.

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.

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 catalog of mouse genes expressed in early embryos, embryonic and adult stem cells was assembled, including 250000 ESTs, representing approximately 39000 unique transcripts. The cDNA libraries, enriched in full-length clones, were condensed into the NIA 15 and 7.4K clone sets, freely distributed to the research community, providing a standard platform for expression studies using microarrays. They are essential tools for studying mammalian development and stem cell biology, and to provide hints about the differential nature of embryonic and adult stem cells.

A set of 7407 cDNA clones (NIA mouse 7.4K) was assembled from >20 cDNA libraries constructed mainly from early mouse embryos, including several stem cell libraries. The clone set was assembled from embryonic and newborn organ libraries consisting of ∼120,000 cDNA clones, which were initially re-arrayed into a set of ∼11,000 unique cDNA clones. A set of tubes was constructed from the racks in this set to prevent contamination and potential mishandling errors in all further re-arrays. Sequences from this set (11K) were analyzed further for quality and clone identity, and high-quality clones with verified identity were re-arrayed into the final set (7.4K). The set is freely available, and a corresponding database was built to provide comprehensive annotation for those clones with known identity or homology, and has been made available through an extensive Web site that includes many link-outs to external databases and analysis servers.[The sequence data from this study have been submitted to GenBank under accession nos. BQ550036–BQ563104.]

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.

Networks of transcription factors (TFs) are thought to determine and maintain the identity of cells. Here we systematically repressed each of 100 TFs with shRNA and carried out global gene expression profiling in mouse embryonic stem (ES) cells. Unexpectedly, only the repression of a handful of TFs significantly affected transcriptomes, which changed in two directions/trajectories: one trajectory by the repression of either Pou5f1 or Sox2; the other trajectory by the repression of either Esrrb, Sall4, Nanog, or Tcfap4. The data suggest that the trajectories of gene expression change are already preconfigured by the gene regulatory network and roughly correspond to extraembryonic and embryonic fates of cell differentiation, respectively. These data also indicate the robustness of the pluripotency gene network, as the transient repression of most TFs did not alter the transcriptomes.

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.

To examine transcription factor (TF) network(s), we created mouse ESC lines, in each of which 1 of 50 TFs tagged with a FLAG moiety is inserted into a ubiquitously controllable tetracycline-repressible locus. Of the 50 TFs, Cdx2 provoked the most extensive transcriptome perturbation in ESCs, followed by Esx1, Sox9, Tcf3, Klf4, and Gata3. ChIP-Seq revealed that CDX2 binds to promoters of upregulated target genes. By contrast, genes downregulated by CDX2 did not show CDX2 binding but were enriched with binding sites for POU5F1, SOX2, and NANOG. Genes with binding sites for these core TFs were also downregulated by the induction of at least 15 other TFs, suggesting a common initial step for ESC differentiation mediated by interference with the binding of core TFs to their target genes. These ESC lines provide a fundamental resource to study biological networks in ESCs and mice. © 2009 Elsevier Inc. All rights reserved.

An in-situ-synthesized 60-mer oligonucleotide microarray designed to detect transcripts from all mouse genes is presented. Exogenous RNA controls derived from yeast allow quantitative estimation of absolute endogenous transcript abundance, 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.

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.

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. © Oxford University Press 2004; all rights reserved.

Applications of microarray technologies to mouse embryology/genetics have been limited, due to the nonavailability of microarrays containing large numbers of embryonic genes and the gap between microgram quantities of RNA required by typical microarray methods and the miniscule amounts of tissue available to researchers. To overcome these problems, we have developed a microarray platform containing in situ-synthesized 60-mer oligonucleotide probes representing approximately 22,000 unique mouse transcripts, assembled primarily from sequences of stem cell and embryo cDNA libraries. We have optimized RNA labeling protocols and experimental designs to use as little as 2 ng total RNA reliably and reproducibly. At least 98% of the probes contained in the microarray correspond to clones in our publicly available collections, making cDNAs readily available for further experimentation on genes of interest. These characteristics, combined with the ability to profile very small samples, make this system a resource for stem cell and embryogenomics research.[Supplemental material is available online atwww.genome.org and at the NIA Mouse cDNA Project Web site,http://lgsun.grc.nia.nih.gov/cDNA/cDNA.html.]

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.

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/.