Search

Although long non-coding RNAs (lncRNAs) are non-protein-coding transcripts by definition, recent studies have shown that a fraction of putative small open reading frames within lncRNAs are translated. However, the biological significance of these hidden polypeptides is still unclear. Here we identify and functionally characterize a novel polypeptide encoded by the lncRNA LINC00961. This polypeptide is conserved between human and mouse, is localized to the late endosome/lysosome and interacts with the lysosomal v-ATPase to negatively regulate mTORC1 activation. This regulation of mTORC1 is specific to activation of mTORC1 by amino acid stimulation, rather than by growth factors. Hence, we termed this polypeptide â small regulatory polypeptide of amino acid response' (SPAR). We show that the SPAR-encoding lncRNA is highly expressed in a subset of tissues and use CRISPR/Cas9 engineering to develop a SPAR-polypeptide-specific knockout mouse while maintaining expression of the host lncRNA. We find that the SPAR-encoding lncRNA is downregulated in skeletal muscle upon acute injury, and using this in vivo model we establish that SPAR downregulation enables efficient activation of mTORC1 and promotes muscle regeneration. Our data provide a mechanism by which mTORC1 activation may be finely regulated in a tissue-specific manner in response to injury, and a paradigm by which lncRNAs encoding small polypeptides can modulate general biological pathways and processes to facilitate tissue-specific requirements, consistent with their restricted and highly regulated expression profile.

Medullary thymic epithelial cells (mTECs) are essential for thymic negative selection to prevent autoimmunity. Previous studies show that mTEC development is dependent on the signal transducers TRAF6 and NIK. However, the downstream target genes of signals controlled by these molecules remain unknown. We performed a microarray analysis on mRNAs down-regulated by deficiencies in TRAF6 or functional NIK in an in vitro organ culture of fetal thymic stromata (2DG-FTOC). An in silico analysis of transcription factor binding sites in plausible promoter regions of differentially expressed genes suggests that STAT1 is involved in TRAF6- and NIK-dependent gene expression. Indeed, the signal of RANK, a TNF receptor family member that activates TRAF6 and NIK, induces the activation of STAT1 in 2DG-FTOC. Moreover, RANK signaling induces the up-regulation of interferon (IFN)-stimulated gene (ISG) expression, suggesting that the RANKL-dependent activation of STAT1 up-regulates ISG expression. The RANKL-dependent expression levels of ISGs were reduced but not completely abolished in interferon α receptor 1-deficient (Ifnar1(-/-)) 2DG-FTOC. Our data suggest that RANK signaling induces ISG expression in both type I interferon-independent and interferon-dependent mechanisms.

DataBase of Transcription Start Sites (DBTSS) is a database which contains precise positional information for transcription start sites (TSSs) of eukaryotic mRNAs. In this update, we included 330 million new tags generated by massively sequencing the 5′-end of oligo-cap selected cDNAs in humans and mice. The tags were collected from normal fetal or adult human tissues, including brain, thymus, liver, kidney and heart, from 6 human cell lines in 21 diverse growth conditions as well as from mouse NIH3T3 cell line: altogether 31 different cell types or culture conditions are represented. This unprecedented increase in depth of data now allows DBTSS to faithfully represent the dynamically changing landscape of TSSs in different cell types and conditions, during development and in the course of evolution. Differential usage of alternative 5′-ends across cell types and conditions can be viewed in a series of new interfaces. Promoter sequence information is now displayed in a comparative genomics viewer where evolutionary turnover of the TSSs can be evaluated. DBTSS can be accessed at http://dbtss.hgc.jp/.

Although the knowledge accumulated on the transcriptional regulations of eukaryotes is significant, the knowledge on their translational regulations remains limited. Thus, we performed a comprehensive detection of terminal oligo-pyrimidine (TOP), which is one of the well-characterized cis-regulatory motifs for translational controls located immediately downstream of the transcriptional start sites of mRNAs. Utilizing our precise 5'-end information of the full-length cDNAs, we could screen 1645 candidate TOP genes by position specific matrix search. Among them, not only 75 out of 78 ribosomal protein genes but also eight previously identified non-ribosomal-protein TOP genes were included. We further experimentally validated the translational activities of 83 TOP candidate genes. Clear translational regulations exerted on the stimulation of 12-O-tetradecanoyl-1-phorbol-13-acetate for at least 41 of them was observed, indicating that there should be a few hundreds of human genes which are subjected to regulation at translation levels via TOPs. Our result suggests that TOP genes code not only formerly characterized ribosomal proteins and translation-related proteins but also a wider variety of proteins, such as lysosome-related proteins and metabolism-related proteins, playing pivotal roles in gene expression controls in the majority of cellular mRNAs.

Although recent studies have revealed that the majority of human genes are subject to regulation of alternative promoters, the biological relevance of this phenomenon remains unclear. We have also demonstrated that roughly half of the human RefSeq genes examined contain putative alternative promoters (PAPs). Here we report large-scale comparative studies of PAPs between human and mouse counterpart genes. Detailed sequence comparison of the 17,245 putative promoter regions (PPRs) in 5463 PAP-containing human genes revealed that PPRs in only a minor fraction of genes (807 genes) showed clear evolutionary conservation as one or more pairs. Also, we found that there were substantial qualitative differences between conserved and non-conserved PPRs, with the latter class being AT-rich PPRs of relative minor usage, enriched in repetitive elements and sometimes producing transcripts that encode small or no proteins. Systematic luciferase assays of these PPRs revealed that both classes of PPRs did have promoter activity, but that their strength ranges were significantly different. Furthermore, we demonstrate that these characteristic features of the non-conserved PPRs are shared with the PPRs of previously discovered putative non-protein coding transcripts. Taken together, our data suggest that there are two distinct classes of promoters in humans, with the latter class of promoters emerging frequently during evolution.

DBTSS was first constructed in 2002 based on precise, experimentally determined 5' end clones. Several major updates and additions have been made since the last report. First, the number of human clones has drastically increased, going from 190,964 to 1,359,000. Second, information about potential alternative promoters is presented because the number of 5' end clones is now sufficient to determine several promoters for one gene. Namely, we defined putative promoter groups by clustering transcription start sites (TSSs) separated by500 bases. A total of 8308 human genes and 4276 mouse genes were found to have putative multiple promoters. Third, DBTSS provides detailed sequence comparisons of user-specified TSSs. Finally, we have added TSS information for zebrafish, malaria and schyzon (a red algae model organism). DBTSS is accessible at http://dbtss.hgc.jp.

Ciona intestinalis is a useful model organism to analyze chordate development and genetics. However, unlike vertebrates, it shares a unique mechanism called trans-splicing with lower eukaryotes. In the computational analysis of trans-splicing in C. intestinalis we report here, we discovered that although the amount of non-trans-spliced and trans-spliced genes is usually equivalent, the expression ratio between the two groups varies significantly with tissues and developmental stages. Among the seven tissues studied, the observed ratios ranged from 2.53 in "gonad" to 19.53 in "endostyle", and during development they increased from 1.68 at the "egg" stage to 7.55 at the "juvenile" stage. We further hypothesize that this enrichment in trans-spliced mRNAs in early developmental stages might be related to the abundance of trans-spliced mRNAs in "gonad". Our analysis indicates that in C. intestinalis, although there may not exist strong fundamental requirements for genes to be trans-spliced, the populations of non-trans-spliced and trans-spliced genes are likely to be spatially and temporally regulated differently.

DBTSS is a database of transcriptional start sites, based on our unique collection of precise, experimentally determined 5′-end sequences of full-length cDNAs. Since its first release in 2002, several major updates have been made. In this update, we expanded the human transcriptional start site dataset by 19 million uniquely mapped, and RefSeq-associated, 5′-end sequences, which were generated by a newly introduced Solexa sequencer. Moreover, in order to provide means for interpreting those massive TSS data, we implemented two new analytical tools: one for connecting expression information with predicted transcription factor binding sites; the other for examining evolutionary conservation or species-specificity of promoters and transcripts, which can be browsed by our own comparative genome viewer. With the expanded dataset and the enhanced functionalities, DBTSS provides a unique platform that enables in-depth transcriptome analyses. DBTSS is accessible at http://dbtss.hgc.jp/.

Epigenetic regulation including DNA methylation plays an important role in several differentiation processes. We profiled global DNA methylation in the neural differentiation of P19 embryonic carcinoma cells using a microarray-based method called MIAMI. We found a genome-wide demethylation of genes. This suggests demethylation rather than methylation is important in neural differentiation.

We present the second version of Melina, a web-based tool for promoter analysis. Melina II shows potential DNA motifs in promoter regions with a combination of several available programs, Consensus, MEME, Gibbs sampler, MDscan and Weeder, as well as several parameter settings. It allows running a maximum of four programs simultaneously, and comparing their results with graphical representations. In addition, users can build a weight matrix from a predicted motif and apply it to upstream sequences of several typical genomes (human, mouse, S. cerevisiae, E. coli, B. subtilis or A. thaliana) or to public motif databases (JASPAR or DBTBS) in order to find similar motifs. Melina II is a client/server system developed by using Adobe (Macromedia) Flash and is accessible over the web at http://melina.hgc.jp.

It has been envisaged that CpG islands are often observed near the transcriptional start sites (TSS) of housekeeping genes. However, neither the precise positions of CpG islands relative to TSS of genes nor the correlation between the presence of the CpG islands and the expression specificity of these genes is well-understood. Using thousands of sequences with known TSS in human and mouse, we found that there is a clear peak in the distribution of CpG islands around TSS in the genes of these two species. Thus, we classified human (mouse) genes into 6600 (2948) CpG+ genes and 2619 (1830) CpG- ones, based on the presence of a CpG island within the -100: +100 region. We estimated the degree of each gene being a housekeeper by the number of cDNA libraries where its ESTs were detected. Then, the tendency that a gene lacking CpG islands around its TSS is expressed with a higher degree of tissue specificity turned out to be evolutionarily conserved. We also confirmed this tendency by analyzing the gene ontology annotation of classified genes. Since no such clear correlation was found in the control data (mRNAs, pre-mRNAs, and chromosome banding pattern), we concluded that the effect of a CpG island near the TSS should be more important than the global GC content of the region where the gene resides.

We report the generation and initial characterization of a large-scale collection of sequences of putative promoter regions (PPRs) of human and mouse genes. Based on our unique collection of 400,225 and 580,209 human and mouse full-length cDNAs, we determined exact transcriptional start sites (TSSs). Using positional information of the TSSs, we could retrieve adjacent sequences as PPRs for 8,793 and 6,875 human and mouse genes, respectively. The positions of the PPRs were 4 kb upstream to previously reported 5'-ends of cDNAs on average, demonstrating that full-length cDNA information is indispensable for this purpose. Among those PPRs supported by experimentally validated TSSs, 3,324 could be paired as mutually homologous genes between human and mouse and were used for the comprehensive comparative studies. The sequence identities in the proximal regions of the TSSs were 45% on average, and 22,794 putative transcription factor binding sites that are conserved between human and mouse were identified. The data resource created in the present work and the results of the sequences' initial characterization should lay the firm foundation for deciphering the transcriptional modulations of human genes. All the data were deposited and made available through a database for comparative studies, DBTSS.

Comparative sequence analysis was carried out for the regions adjacent to experimentally validated transcriptional start sites (TSSs), using 3324 pairs of human and mouse genes. We aligned the upstream putative promoter sequences over the 1-kb proximal regions and found that the sequence conservation could not be further extended at, on average, 510 bp upstream positions of the TSSs. This discontinuous manner of the sequence conservation revealed a “block” structure in about one-third of the putative promoter regions. Consistently, we also observed that G+C content and CpG frequency were significantly different inside and outside the blocks. Within the blocks, the sequence identity was uniformly 65% regardless of their length. About 90% of the previously characterized transcription factor binding sites were located within those blocks. In 46% of the blocks, the 5′ ends were bounded by interspersed repetitive elements, some of which may have nucleated the genomic rearrangements. The length of the blocks was shortest in the promoters of genes encoding transcription factors and of genes whose expression patterns are brain specific, which suggests that the evolutional diversifications in the transcriptional modulations should be the most marked in these populations of genes.

The C 2 H 2 zinc finger motif found in many transcription factors is thought to be for nucleic acid binding and/or dimerization. Nearly 1% of eukaryote genes are estimated to encode this motif. We investigated the gene family encoding this motif in the Mus musculus mRNA by molecular indexing, a technique used to select a subpopulation of cDNA by ligation of adapters to cDNA fragments digested by a class IIS restriction enzyme(s). In place of oligo-dT primers in the original method, a polymerase chain reaction primer designed based on the conserved sequence of the C 2 H 2 zinc finger protein stranded cDNA was made from various mouse tissue mRNAs, digested with Fok I and Bsm AI, and joined with adapters. Amplification of the cDNA with an adapter primer and zinc finger-specific primer yielded products enriched in zinc finger protein genes. Fragments were separated by subsequent denaturing polyacrylamide gel electrophoresis, and characterized by DNA sequencing. Consequently, 259 C 2 H 2 zinc finger motif sequences were obtained, among which 166 were novel. Combined with the reported sequences, these mouse motif sequences were compared with those of other species such as Saccharomyces cerevisiae and Caenorhabditis elegans . Some of the amino acids in the motif sequence showed strong bias among species. Most of the novel sequences were supposed to be DNA-binding according to the surface potential of predicted tertiary structures.

Background CpG islands are observed in mammals and other vertebrates, generally escape DNA methylation, and tend to occur in the promoters of widely expressed genes. Another class of promoter has lower G+C and CpG contents, and is thought to be involved in the spatiotemporal regulation of gene expression. Non-vertebrate deuterostomes are reported to have a single class of promoter with high-frequency CpG dinucleotides, suggesting that this is the original type of promoter. However, the limited annotation of these genes has impeded the large-scale analysis of their promoters. Results To determine the origins of the two classes of vertebrate promoters, we chose Ciona intestinalis, an invertebrate that is evolutionarily close to the vertebrates, and identified its transcription start sites genome-wide using a next-generation sequencer. We indeed observed a high CpG content around the transcription start sites, but their levels in the promoters and background sequences differed much less than in mammals. The CpG-rich stretches were also fairly restricted, so they appeared more similar to mammalian CpG-poor promoters. Conclusions From these data, we infer that CpG islands are not sufficiently ancient to be found in invertebrates. They probably appeared early in vertebrate evolution via some active mechanism and have since been maintained as part of vertebrate promoters.

SummaryThe transcription factor Foxp3 is essential for the development of regulatory T (Treg) cells, yet its expression is insufficient for establishing the Treg cell lineage. Here we showed that Treg cell development was achieved by the combination of two independent processes, i.e., the expression of Foxp3 and the establishment of Treg cell-specific CpG hypomethylation pattern. Both events were induced by T cell receptor stimulation. The Treg cell-type CpG hypomethylation began in the thymus and continued to proceed in the periphery and could be fully established without Foxp3. The hypomethylation was required for Foxp3+ T cells to acquire Treg cell-type gene expression, lineage stability, and full suppressive activity. Thus, those T cells in which the two events have concurrently occurred are developmentally set into the Treg cell lineage. This model explains how Treg cell fate and plasticity is controlled and can be exploited to generate functionally stable Treg cells.

The high similarity of tunicates and vertebrates during their development coupled with the transparency of tunicate larvae, their well-studied cell lineages and the availability of simple and efficient transgenesis methods makes of this subphylum an ideal system for the investigation of vertebrate physiological and developmental processes. Recently, the sequencing of two different Ciona genomes has lead to the identification of numerous genes. In order to better understand the regulation of these genes, a database was created containing information on regulation of tunicate genes collected from literature. It includes for instance information regarding the minimal promoter length, the transcription factors involved and their binding sites, as well as the localization of the gene expression. Additionally, binding sites for characterized transcription factors were predicted based on published in vitro recognition sites. Comparison of the promoters of homologous genes in different species is also provided to allow identification of conserved cis elements. At the time of writing, information about 184 promoters, containing 73 identified binding sites and >2000 newly predicted binding sites is available. This database is accessible at http://dbtgr.hgc.jp.

Background Interspecies sequence comparison is a powerful tool to extract functional or evolutionary information from the genomes of organisms. A number of studies have compared protein sequences or promoter sequences between mammals, which provided many insights into genomics. However, the correlation between protein conservation and promoter conservation remains controversial. Results We examined promoter conservation as well as protein conservation for 6,901 human and mouse orthologous genes, and observed a very weak correlation between them. We further investigated their relationship by decomposing it based on functional categories, and identified categories with significant tendencies. Remarkably, the 'ribosome' category showed significantly low promoter conservation, despite its high protein conservation, and the 'extracellular matrix' category showed significantly high promoter conservation, in spite of its low protein conservation. Conclusion Our results show the relation of gene function to protein conservation and promoter conservation, and revealed that there seem to be nonparallel components between protein and promoter sequence evolution.