Bioinformatics Approaches for Determining the Functional Impact of Repetitive Elements on Non-coding RNAs.

With a large number of annotated non-coding RNAs (ncRNAs), repetitive sequences are found to constitute functional components (termed as repetitive elements) in ncRNAs that perform specific biological functions. Bioinformatics analysis is a powerful tool for improving our understanding of the role of repetitive elements in ncRNAs. This chapter summarizes recent findings that reveal the role of repetitive elements in ncRNAs. Furthermore, relevant bioinformatics approaches are systematically reviewed, which promises to provide valuable resources for studying the functional impact of repetitive elements on ncRNAs.

Glycerol kinase stimulates uncoupling protein 1 expression by regulating fatty acid metabolism in beige adipocytes.

Browning of adipose tissue is induced by specific stimuli such as cold exposure and consists of up-regulation of thermogenesis in white adipose tissue. Recently, it has emerged as an attractive target for managing obesity in humans. Here, we performed a comprehensive analysis to identify genes associated with browning in murine adipose tissue. We focused on glycerol kinase (GYK) because its mRNA expression pattern is highly correlated with that of uncoupling protein 1 (UCP1), which regulates the thermogenic capacity of adipocytes. Cold exposure-induced Ucp1 up-regulation in inguinal white adipose tissue (iWAT) was partially abolished by Gyk knockdown (KD) in vivo Consistently, the Gyk KD inhibited Ucp1 expression induced by treatment with the ß-adrenergic receptors (ßAR) agonist isoproterenol (Iso) in vitro and resulted in impaired uncoupled respiration. Gyk KD also suppressed Iso- and adenylate cyclase activator-induced transcriptional activation and phosphorylation of the cAMP response element-binding protein (CREB). However, we did not observe these effects with a cAMP analog. Therefore Gyk KD related to Iso-induced cAMP products. In Iso-treated Gyk KD adipocytes, stearoyl-CoA desaturase 1 (SCD1) was up-regulated, and monounsaturated fatty acids such as palmitoleic acid (POA) accumulated. Moreover, a SCD1 inhibitor treatment recovered the Gyk KD-induced Ucp1 down-regulation and POA treatment down-regulated Iso-activated Ucp1 Our findings suggest that Gyk stimulates Ucp1 expression via a mechanism that partially depends on the ßAR-cAMP-CREB pathway and Gyk-mediated regulation of fatty acid metabolism.

Long non-coding RNA 2310069B03Rik functions as a suppressor of Ucp1 expression under prolonged cold exposure in murine beige adipocytes.

Specific conditions, such as exposure to cold, can induce the production of brown-like adipocytes in white adipose tissue. These adipocytes express high levels of uncoupling protein 1 (UCP1) and energy expended by generating heat. Thus, these are a potential target for the prevention or treatment of obesity. The present study involved a comprehensive analysis of the adipose tissue to understand the relationship between long non-coding RNA (lncRNA) 2310069B03Rik and UCP1. Cold exposure increased both lncRNA 2310069B03Rik and Ucp1 expression in inguinal white adipose tissue (iWAT). However, overexpression of lncRNA 2310069B03Rik suppressed the Ucp1 mRNA expression and the promoter activity of UCP1 in the iWAT primary adipocytes. In addition, compared to the early induction of Ucp1 expression by cold stimulation, the induction of lncRNA 2310069B03Rik expression was later. These results suggest that lncRNA 2310069B03Rik functions as a suppression factor of Ucp1 expression.

Common genetic variants associated with Parkinson's disease display widespread signature of epigenetic plasticity.

Parkinson disease (PD) is characterized by a pivotal progressive loss of substantia nigra dopaminergic neurons and aggregation of α-synuclein protein encoded by the SNCA gene. Genome-wide association studies identified almost 100 sequence variants linked to PD in SNCA. However, the consequences of this genetic variability are rather unclear. Herein, our analysis on selective single nucleotide polymorphisms (SNPs) which are highly associated with the PD susceptibility revealed that several SNP sites attribute to the nucleosomes and overlay with bivalent regions poised to adopt either active or repressed chromatin states. We also identified large number of transcription factor (TF) binding sites associated with these variants. In addition, we located two docking sites in the intron-1 methylation prone region of SNCA which are required for the putative interactions with DNMT1. Taken together, our analysis reflects an additional layer of epigenomic contribution for the regulation of the SNCA gene in PD.

Comprehensive epigenome characterization reveals diverse transcriptional regulation across human vascular endothelial cells.

BACKGROUND: Endothelial cells (ECs) make up the innermost layer throughout the entire vasculature. Their phenotypes and physiological functions are initially regulated by developmental signals and extracellular stimuli. The underlying molecular mechanisms responsible for the diverse phenotypes of ECs from different organs are not well understood. RESULTS: To characterize the transcriptomic and epigenomic landscape in the vascular system, we cataloged gene expression and active histone marks in nine types of human ECs (generating 148 genome-wide datasets) and carried out a comprehensive analysis with chromatin interaction data. We developed a robust procedure for comparative epigenome analysis that circumvents variations at the level of the individual and technical noise derived from sample preparation under various conditions. Through this approach, we identified 3765 EC-specific enhancers, some of which were associated with disease-associated genetic variations. We also identified various candidate marker genes for each EC type. We found that the nine EC types can be divided into two subgroups, corresponding to those with upper-body origins and lower-body origins, based on their epigenomic landscape. Epigenomic variations were highly correlated with gene expression patterns, but also provided unique information. Most of the deferentially expressed genes and enhancers were cooperatively enriched in more than one EC type, suggesting that the distinct combinations of multiple genes play key roles in the diverse phenotypes across EC types. Notably, many homeobox genes were differentially expressed across EC types, and their expression was correlated with the relative position of each organ in the body. This reflects the developmental origins of ECs and their roles in angiogenesis, vasculogenesis and wound healing. CONCLUSIONS: This comprehensive analysis of epigenome characterization of EC types reveals diverse transcriptional regulation across human vascular systems. These datasets provide a valuable resource for understanding the vascular system and associated diseases.

Improvement of detection performance of fusion genes from RNA-seq data by clustering short reads.

Fusion genes are involved in cancer, and their detection using RNA-Seq is insufficient given the relatively short reading length. Therefore, we proposed a shifted short-read clustering (SSC) method, which focuses on overlapping reads from the same loci and extends them as a representative sequence. To verify their usefulness, we applied the SSC method to RNA-Seq data from four types of cell lines (BT-474, MCF-7, SKBR-3, and T-47D). As the slide width of the SSC method increased to one, two, five, or ten bases, the read length was extended from 201 bases to 217 (108%), 234 (116%), 282 (140%), or 317 (158%) bases, respectively. Furthermore, fusion genes were investigated using STAR-Fusion, a fusion gene detection tool, with and without the SSC method. When one base was shifted by the SSC method, the reads mapped to multiple loci decreased from 9.7% to 4.6%, and the sensitivity of the fusion gene was improved from 47% to 54% on average (BT-474: from 48% to 57%, MCF-7: 49% to 53%, SKBR-3: 50% to 57%, and T-47D: 43% to 50%) compared with original data. When the reads are shifted more, the positive predictive value was also improved. The SSC method could be an effective method for fusion gene detection.

The dipeptidyl peptidase-4 (DPP-4) inhibitor teneligliptin enhances brown adipose tissue function, thereby preventing obesity in mice.

To clarify the effects of a dipeptidyl peptidase-4 (DPP-4) inhibitor on whole-body energy metabolism, we treated mice fed a high-fat diet (HFD) with teneligliptin, a clinically available DPP-4 inhibitor. Teneligliptin significantly prevented HFD-induced obesity and obesity-associated metabolic disorders. It also increased oxygen consumption rate and upregulated uncoupling protein 1 (UCP1) expression in both brown adipose tissue (BAT) and inguinal white adipose tissue (iWAT), suggesting that it enhances BAT function. Soluble DPP-4 inhibited ß-adrenoreceptor-stimulated UCP1 expression in primary adipocytes, and this inhibition was prevented in the presence of teneligliptin, or an extracellular signal-related kinase inhibitor. These results indicate that soluble DPP-4 inhibits ß-adrenoreceptor-stimulated UCP1 induction and that chronic DPP-4 inhibitor treatment may prevent obesity through the activation of BAT function.

Characteristics of functional enrichment and gene expression level of human putative transcriptional target genes.

BACKGROUND: Transcriptional target genes show functional enrichment of genes. However, how many and how significantly transcriptional target genes include functional enrichments are still unclear. To address these issues, I predicted human transcriptional target genes using open chromatin regions, ChIP-seq data and DNA binding sequences of transcription factors in databases, and examined functional enrichment and gene expression level of putative transcriptional target genes. RESULTS: Gene Ontology annotations showed four times larger numbers of functional enrichments in putative transcriptional target genes than gene expression information alone, independent of transcriptional target genes. To compare the number of functional enrichments of putative transcriptional target genes between cells or search conditions, I normalized the number of functional enrichment by calculating its ratios in the total number of transcriptional target genes. With this analysis, native putative transcriptional target genes showed the largest normalized number of functional enrichments, compared with target genes including 5-60% of randomly selected genes. The normalized number of functional enrichments was changed according to the criteria of enhancer-promoter interactions such as distance from transcriptional start sites and orientation of CTCF-binding sites. Forward-reverse orientation of CTCF-binding sites showed significantly higher normalized number of functional enrichments than the other orientations. Journal papers showed that the top five frequent functional enrichments were related to the cellular functions in the three cell types. The median expression level of transcriptional target genes changed according to the criteria of enhancer-promoter assignments (i.e. interactions) and was correlated with the changes of the normalized number of functional enrichments of transcriptional target genes. CONCLUSIONS: Human putative transcriptional target genes showed significant functional enrichments. Functional enrichments were related to the cellular functions. The normalized number of functional enrichments of human putative transcriptional target genes changed according to the criteria of enhancer-promoter assignments and correlated with the median expression level of the target genes. These analyses and characters of human putative transcriptional target genes would be useful to examine the criteria of enhancer-promoter assignments and to predict the novel mechanisms and factors such as DNA binding proteins and DNA sequences of enhancer-promoter interactions.

Transcription factor IRF8 plays a critical role in the development of murine basophils and mast cells.

Basophils and mast cells play critical roles in host defense against pathogens and allergic disorders. However, the molecular mechanism by which these cells are generated is not completely understood. Here we demonstrate that interferon regulatory factor-8 (IRF8), a transcription factor essential for the development of several myeloid lineages, also regulates basophil and mast cell development. Irf8(-/-) mice displayed a severe reduction in basophil counts, which was accounted for by the absence of pre-basophil and mast cell progenitors (pre-BMPs). Although Irf8(-/-) mice retained peripheral tissue mast cells, remaining progenitors from Irf8(-/-) mice including granulocyte progenitors (GPs) were unable to efficiently generate either basophils or mast cells, indicating that IRF8 also contributes to the development of mast cells. IRF8 appeared to function at the GP stage, because IRF8 was expressed in GPs, but not in basophils, mast cells, and basophil/mast cell-restricted progenitor cells. Furthermore, we demonstrate that GATA2, a transcription factor known to promote basophil and mast cell differentiation, acts downstream of IRF8. These results shed light on the pathways and mechanism underlying the development of basophils and mast cells.

Essential role of the IRF8-KLF4 transcription factor cascade in murine monocyte differentiation.

Monocytes regulate host defenses, inflammation, and tissue homeostasis. The transcription factor interferon regulatory factor-8 (IRF8) stimulates monocyte/macrophage differentiation, yet genome-wide understanding of the differentiation program initiated by IRF8 is lacking. By combining chromatin immunoprecipitation sequencing with gene expression profiling, we show that during IRF8-dependent monocyte differentiation, IRF8 binding occurs at both promoter-proximal and promotor-distal regions together with the transcription factor PU.1 and is associated with gene induction. Many of the promoter-distal IRF8 binding sites show an increase in histone H3 lysine 4 monomethylation, a signature for enhancers. However, about half the IRF8-induced genes were not bound by IRF8, suggesting the involvement of downstream transcription factors. Analysis of DNA motifs in cis-regulatory elements of these indirect IRF8 target genes predicted that Krüppel-like factor-4 (KLF4)-essential for Ly6C(+) monocyte development-is one such factor. Indeed, monocyte development in Irf8(-/-) mice is as defective as that in Klf4(-/-) chimeric mice. Moreover, Irf8(-/-) monocyte-dendritic cell progenitors do not express Klf4 messenger RNA. Introduction of KLF4 into an Irf8(-/-) myeloid progenitor cell line induced a subset of IRF8 target genes and caused partial monocyte differentiation. Taken together, our present results uncover genome-wide behavior of IRF8 and identify an IRF8-KLF4 axis that operates during monocyte differentiation.

Chromosome-biased binding and gene regulation by the Caenorhabditis elegans DRM complex.

DRM is a conserved transcription factor complex that includes E2F/DP and pRB family proteins and plays important roles in development and cancer. Here we describe new aspects of DRM binding and function revealed through genome-wide analyses of the Caenorhabditis elegans DRM subunit LIN-54. We show that LIN-54 DNA-binding activity recruits DRM to promoters enriched for adjacent putative E2F/DP and LIN-54 binding sites, suggesting that these two DNA-binding moieties together direct DRM to its target genes. Chromatin immunoprecipitation and gene expression profiling reveals conserved roles for DRM in regulating genes involved in cell division, development, and reproduction. We find that LIN-54 promotes expression of reproduction genes in the germline, but prevents ectopic activation of germline-specific genes in embryonic soma. Strikingly, C. elegans DRM does not act uniformly throughout the genome: the DRM recruitment motif, DRM binding, and DRM-regulated embryonic genes are all under-represented on the X chromosome. However, germline genes down-regulated in lin-54 mutants are over-represented on the X chromosome. We discuss models for how loss of autosome-bound DRM may enhance germline X chromosome silencing. We propose that autosome-enriched binding of DRM arose in C. elegans as a consequence of germline X chromosome silencing and the evolutionary redistribution of germline-expressed and essential target genes to autosomes. Sex chromosome gene regulation may thus have profound evolutionary effects on genome organization and transcriptional regulatory networks.

Insertional mutagenesis by the Tol2 transposon-mediated enhancer trap approach generated mutations in two developmental genes: tcf7 and synembryn-like.

Gene trap and enhancer trap methods using transposon or retrovirus have been recently described in zebrafish. However, insertional mutants using these methods have not been reported. We report here development of an enhancer trap method by using the Tol2 transposable element and identification and characterization of insertional mutants. We created 73 fish lines that carried single copy insertions of an enhancer trap construct, which contained the zebrafish hsp70 promoter and the GFP gene, in their genome and expressed GFP in specific cells, tissues and organs, indicating that the hsp70 promoter is highly capable of responding to chromosomal enhancers. First, we analyzed genomic DNA surrounding these insertions. Fifty-one of them were mapped onto the current version of the genomic sequence and 43% (22/51) were located within transcribed regions, either exons or introns. Then, we crossed heterozygous fish carrying the same insertions and identified two insertions that caused recessive mutant phenotypes. One disrupted the tcf7 gene, which encodes a transcription factor of the Tcf/Lef family mediating Wnt signaling, and caused shorter and wavy median fin folds and pectoral fins. We knocked down Lef1, another member of the Tcf/Lef family also expressed in the fin bud, in the tcf7 mutant, and revealed functional redundancy of these factors and their essential role in establishment of the apical ectodermal ridge (AER). The other disrupted the synembryn-like gene (synbl), a homolog of the C. elegans synembryn gene, and caused embryonic lethality and small pigment spots. The pigment phenotype was rescued by application of forskolin, an activator of adenylyl cyclase, suggesting that the synbl gene activates the Galpha(S) pathway leading to activation of adenylyl cyclase. We thus demonstrated that the transposon-mediated enhancer trap approach can indeed create insertional mutations in developmental genes. Our present study provides a basis for the development of efficient transposon-mediated insertional mutagenesis in a vertebrate.

The evolutionary emergence of cell type-specific genes inferred from the gene expression analysis of Hydra.

Cell lineages of cnidarians including Hydra represent the fundamental cell types of metazoans and provides us a unique opportunity to study the evolutionary diversification of cell type in the animal kingdom. Hydra contains epithelial cells as well as a multipotent interstitial cell (I-cell) that gives rise to nematocytes, nerve cells, gland cells, and germ-line cells. We used cDNA microarrays to identify cell type-specific genes by comparing gene expression in normal Hydra with animals lacking the I-cell lineage, so-called epithelial Hydra. We then performed in situ hybridization to localize expression to specific cell types. Eighty-six genes were shown to be expressed in specific cell types of the I-cell lineage. An additional 29 genes were expressed in epithelial cells and were down-regulated in epithelial animals lacking I-cells. Based on the above information, we constructed a database (http://hydra.lab.nig.ac.jp/hydra/), which describes the expression patterns of cell type-specific genes in Hydra. Most genes expressed specifically in either I-cells or epithelial cells have homologues in higher metazoans. By comparison, most nematocyte-specific genes and approximately half of the gland cell- and nerve cell-specific genes are unique to the cnidarian lineage. Because nematocytes, gland cells, and nerve cells appeared along with the emergence of cnidarians, this suggests that lineage-specific genes arose in cnidarians in conjunction with the evolution of new cell types required by the cnidarians.

Transcriptional interferences in cis natural antisense transcripts of humans and mice.

For a significant fraction of mRNAs, their expression is regulated by other RNAs, including cis natural antisense transcripts (cis-NATs) that are complementary mRNAs transcribed from opposite strands of DNA at the same genomic locus. The regulatory mechanism of mRNA expression by cis-NATs is unknown, although a few possible explanations have been proposed. To understand this regulatory mechanism, we conducted a large-scale analysis of the currently available data and examined how the overlapping arrangements of cis-NATs affect their expression level. Here, we show that for both human and mouse the expression level of cis-NATs decreases as the length of the overlapping region increases. In particular, the proportions of the highly expressed cis-NATs in all cis-NATs examined were approximately 36 and 47% for human and mouse, respectively, when the overlapping region was <200 bp. However, both proportions decreased to virtually zero when the overlapping regions were >2000 bp in length. Moreover, the distribution of the expression level of cis-NATs changes according to different types of the overlapping pattern of cis-NATs in the genome. These results are consistent with the transcriptional collision model for the regulatory mechanism of gene expression by cis-NATs.

A genome-wide and nonredundant mouse transcription factor database.

Here we describe the development of a genome-wide and nonredundant mouse transcription factor database and its viewer (http://genome.gsc.riken.gp/TFdb/). We systematically selected transcription factors with DNA-binding properties and their regulators on the basis of their LocusLink and Gene Ontology annotations. We also incorporated into our database information regarding the corresponding available cDNA clones and their structural properties. Because of these features, our database is unique and may provide useful information for systematic genome-wide studies of transcriptional regulation.

Identification of region-specific transcription factor genes in the adult mouse brain by medium-scale real-time RT-PCR.

We established a medium-scale real-time RT-PCR system focusing on transcription factors and applied it to their expression profiles in the adult mouse 11 brain regions (http://genome.gsc.riken.jp/qRT-PCR/). Almost 90% of the examined genes showed significant expression in at least one region. We successfully extracted 179 region-specific genes by clustering analysis. Interestingly, the transcription factors involved in the development of the pituitary were still expressed in the adult pituitary, suggesting that they also play important roles in maintenance of the pituitary. These results provide unique molecular markers that may account for the molecular basis of the unique functions of specific brain regions.

Targeting a complex transcriptome: the construction of the mouse full-length cDNA encyclopedia.

We report the construction of the mouse full-length cDNA encyclopedia,the most extensive view of a complex transcriptome,on the basis of preparing and sequencing 246 libraries. Before cloning,cDNAs were enriched in full-length by Cap-Trapper,and in most cases,aggressively subtracted/normalized. We have produced 1,442,236 successful 3'-end sequences clustered into 171,144 groups, from which 60,770 clones were fully sequenced cDNAs annotated in the FANTOM-2 annotation. We have also produced 547,149 5' end reads,which clustered into 124,258 groups. Altogether, these cDNAs were further grouped in 70,000 transcriptional units (TU),which represent the best coverage of a transcriptome so far. By monitoring the extent of normalization/subtraction, we define the tentative equivalent coverage (TEC),which was estimated to be equivalent to >12,000,000 ESTs derived from standard libraries. High coverage explains discrepancies between the very large numbers of clusters (and TUs) of this project,which also include non-protein-coding RNAs,and the lower gene number estimation of genome annotations. Altogether,5'-end clusters identify regions that are potential promoters for 8637 known genes and 5'-end clusters suggest the presence of almost 63,000 transcriptional starting points. An estimate of the frequency of polyadenylation signals suggests that at least half of the singletons in the EST set represent real mRNAs. Clones accounting for about half of the predicted TUs await further sequencing. The continued high-discovery rate suggests that the task of transcriptome discovery is not yet complete.

Antisense transcripts with FANTOM2 clone set and their implications for gene regulation.

We have used the FANTOM2 mouse cDNA set (60,770 clones), public mRNA data, and mouse genome sequence data to identify 2481 pairs of sense-antisense transcripts and 899 further pairs of nonantisense bidirectional transcription based upon genomic mapping. The analysis greatly expands the number of known examples of sense-antisense transcript and nonantisense bidirectional transcription pairs in mammals. The FANTOM2 cDNA set appears to contain substantially large numbers of noncoding transcripts suitable for antisense transcript analysis. The average proportion of loci encoding sense-antisense transcript and nonantisense bidirectional transcription pairs on autosomes was 15.1 and 5.4%, respectively. Those on the X chromosome were 6.3 and 4.2%, respectively. Sense-antisense transcript pairs, rather than nonantisense bidirectional transcription pairs, may be less prevalent on the X chromosome, possibly due to X chromosome inactivation. Sense and antisense transcripts tended to be isolated from the same libraries, where nonantisense bidirectional transcription pairs were not apparently coregulated. The existence of large numbers of natural antisense transcripts implies that the regulation of gene expression by antisense transcripts is more common that previously recognized. The viewer showing mapping patterns of sense-antisense transcript pairs and nonantisense bidirectional transcription pairs on the genome and other related statistical data is available on our Web site.

Antisense transcripts with rice full-length cDNAs.

BACKGROUND: Natural antisense transcripts control gene expression through post-transcriptional gene silencing by annealing to the complementary sequence of the sense transcript. Because many genome and mRNA sequences have become available recently, genome-wide searches for sense-antisense transcripts have been reported, but few plant sense-antisense transcript pairs have been studied. The Rice Full-Length cDNA Sequencing Project has enabled computational searching of a large number of plant sense-antisense transcript pairs. RESULTS: We identified sense-antisense transcript pairs from 32,127 full-length rice cDNA sequences produced by this project and public rice mRNA sequences by aligning the cDNA sequences with rice genome sequences. We discovered 687 bidirectional transcript pairs in rice, including sense-antisense transcript pairs. Both sense and antisense strands of 342 pairs (50%) showed homology to at least one expressed sequence tag other than that of the pair. Microarray analysis showed 82 pairs (32%) out of 258 pairs on the microarray were more highly expressed than the median expression intensity of 21,938 rice transcriptional units. Both sense and antisense strands of 594 pairs (86%) had coding potential. CONCLUSIONS: The large number of plant sense-antisense transcript pairs suggests that gene regulation by antisense transcripts occurs in plants and not only in animals. On the basis of our results, experiments should be carried out to analyze the function of plant antisense transcripts.

Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana.

The NAC domain was originally characterized from consensus sequences from petunia NAM and from Arabidopsis ATAF1, ATAF2, and CUC2. Genes containing the NAC domain (NAC family genes) are plant-specific transcriptional regulators and are expressed in various developmental stages and tissues. We performed a comprehensive analysis of NAC family genes in Oryza sativa (a monocot) and Arabidopsis thaliana (a dicot). We found 75 predicted NAC proteins in full-length cDNA data sets of O. sativa (28,469 clones) and 105 in putative genes (28,581 sequences) from the A. thaliana genome. NAC domains from both predicted and known NAC family proteins were classified into two groups and 18 subgroups by sequence similarity. There were a few differences in amino acid sequences in the NAC domains between O. sativa and A. thaliana. In addition, we found 13 common sequence motifs from transcriptional activation regions in the C-terminal regions of predicted NAC proteins. These motifs probably diverged having correlations with NAC domain structures. We discuss the relationship between the structure and function of the NAC family proteins in light of our results and the published data. Our results will aid further functional analysis of NAC family genes.