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.

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.

Fusion genes make for one of the mechanisms of tumorigenesis. The identification of fusion genes by RNA-Seq has attracted attention. Various methods for detecting fusion genes have been proposed, but their accuracy is not sufficient. One of the causes of this problem is the relatively short reading length in RNA-Seq data. Therefore, before mapping RNA-Seq data, we proposed a method, which is based on shifted short-read clustering (SSC), to identify shifted reads of the same origin and extend them as representative sequences. As a result, we assumed that the percentage of uniquely mapped reads would be increased, and the detection rates of the fusion genes could be improved. To verify these hypotheses, we applied the SSC method to RNA-Seq data from three cell-lines (BT-474, MCF-7, and SKBR-3). When only one base was shifted, the average read lengths of BT-474, MCF-7, and SKBR-3 were extended from 201 to 223 bases (111%), 201 to 214 bases (106%), and 201 to 213 bases (106%), respectively. Furthermore, the effectiveness of the SSC method is demonstrated by comparing the performances of a fusion gene detection tool's results, STAR-Fusion, with and without the SSC method of the reads. The percentage of uniquely mapped reads of BT-474, MCF-7, and SKBR-3 were improved from 88% to 93%, 88% to 94%, and 92% to 95%, respectively. Finally, the fusion gene detection rates of BT-474, MCF-7, and SKBR-3 were increased from 48% to 57%, 49% to 53%, and 50% to 53% respectively. The SSC method is considered to be an effective method not only for improving the percentage of uniquely mapped reads but also for fusion gene detection.

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

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

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.

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 G(alpha 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.

Here we report the new features and improvements in our latest release of the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/), a comprehensive annotation resource for human genes and transcripts. H-InvDB, originally developed as an integrated database of the human transcriptome based on extensive annotation of large sets of full-length cDNA (FLcDNA) clones, now provides annotation for 120 558 human mRNAs extracted from the International Nucleotide Sequence Databases (INSD), in addition to 54 978 human FLcDNAs, in the latest release H-InvDB_4.6. We mapped those human transcripts onto the human genome sequences (NCBI build 36.1) and determined 34 699 human gene clusters, which could define 34 057 (98.1%) protein-coding and 642 (1.9%) non-protein-coding loci; 858 (2.5%) transcribed loci overlapped with predicted pseudogenes. For all these transcripts and genes, we provide comprehensive annotation including gene structures, gene functions, alternative splicing variants, functional non-protein-coding RNAs, functional domains, predicted sub cellular localizations, metabolic pathways, predictions of protein 3D structure, mapping of SNPs and microsatellite repeat motifs, co-localization with orphan diseases, gene expression profiles, orthologous genes, proteinprotein interactions (PPI) and annotation for gene families. The current H-InvDB annotation resources consist of two main views: Transcript view and Locus view and eight sub-databases: the DiseaseInfo Viewer, H-ANGEL, the Clustering Viewer, G-integra, the TOPO Viewer, Evola, the PPI view and the Gene family/group.

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 (1-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 1-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 1-cell lineage. An additional 29 genes were expressed in epithelial cells and were down-regulated in epithelial animals lacking 1-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-specif ic genes in Hydra. Most genes expressed specifically in either 1-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-specif ic 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 cniclarians.

For a significant fraction of mRNA-s, 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 similar to 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.

We present here the annotation of the complete genome of rice Oryza sativa L. ssp. japonica cultivar Nipponbare. All functional annotations for proteins and non-protein-coding RNA (npRNA) candidates were manually curated. Functions were identified or inferred in 19,969 (70%) of the proteins, and 131 possible npRNAs (including 58 antisense transcripts) were found. Almost 5000 annotated protein-coding genes were found to be disrupted in insertional mutant lines, which will accelerate future experimental validation of the annotations. The rice loci were determined by using cDNA sequences obtained from rice and other representative cereals. Our conservative estimate based on these loci and an extrapolation suggested that the gene number of rice is ∼32,000, which is smaller than previous estimates. We conducted comparative analyses between rice and Arabidopsis thaliana and found that both genomes possessed several lineage-specific genes, which might account for the observed differences between these species, while they had similar sets of predicted functional domains among the protein sequences. A system to control translational efficiency seems to be conserved across large evolutionary distances. Moreover, the evolutionary process of protein-coding genes was examined. Our results suggest that natural selection may have played a role for duplicated genes in both species, so that duplication was suppressed or favored in a manner that depended on the function of a gene. ©2007 by Cold Spring Harbor Laboratory Press.

Here we describe the development of a genome-wide and nonredundant mouse transcription factor database and its viewer (http://genome.gsc.riken.jp/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. (C) 2004 Elsevier Inc. All rights reserved.

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. (C) 2004 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.

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.

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

We collected and completely sequenced 28,469 full-length complementary DNA clones from Oryza sativa L. ssp. japonica cv. Nipponbare. Through homology searches of publicly available sequence data, we assigned tentative protein functions to 21,596 clones (75.86%). Mapping of the cDNA clones to genomic DNA revealed that there are 19,000 to 20,500 transcription units in the rice genome. Protein informatics analysis against the InterPro database revealed the existence of proteins presented in rice but not in Arabidopsis. Sixty-four percent of our cDNAs are homologous to Arabidopsis proteins.

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, S'-end clusters identify regions that are potential promoters for 8637 known genes and S'-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.

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.

Only a small proportion of the mouse genome is transcribed into mature messenger RNA transcripts. There is an international collaborative effort to identify all full-length mRNA transcripts from the mouse, and to ensure that each is represented in a physical collection of clones. Here we report the manual annotation of 60,770 full-length mouse complementary DNA sequences. These are clustered into 33,409 'transcriptional units', contributing 90.1% of a newly established mouse transcriptome database. Of these transcriptional units, 4,258 are new protein-coding and 11,665 are new non-coding messages, indicating that non-coding RNA is a major component of the transcriptome. 41% of all transcriptional units showed evidence of alternative splicing. In protein-coding transcripts, 79% of splice variations altered the protein product. Whole-transcriptome analyses resulted in the identification of 2,431 sense-antisense pairs. The present work, completely supported by physical clones, provides the most comprehensive survey of a mammalian transcriptome so far, and is a valuable resource for functional genomics.

We describe a computer-based method that selects representative clones for full-length sequencing in a full-length cDNA project. Our method classifies end sequences using two kinds of criteria, grouping, and clustering. Grouping places together variant cDNAs, family genes, and cDNAs with sequencing errors. Clustering separates those cDNA clones into distinct Clusters. The full-length sequences of the clones selected by grouping are determined preferentially, and then the sequences selected by clustering are determined. Grouping reduced the number of rice cDNA clones for full-length sequencing to 21% and mouse cDNA clones to 25%. Rice full-length sequences selected by grouping showed a 1.07-fold redundancy. Mouse full-length sequences showed a 1.04-fold redundancy, which can be reduced by similar to30% from the selection using our previous method. To estimate the coverage of unique genes, we used FANTOM (Functional Annotation of RIKEN Mouse cDNA Clones) clusters (the RIKEN Genome Exploration Research Group 2001). Grouping covered almost all unique genes (93% of FANTOM clusters), and clustering covered all genes. Therefore, our method is useful for the selection of appropriate representative clones for full-length Sequencing, thereby greatly reducing the cost, labor, and time necessary for this process.

We have developed a method to overcome sequencing problems caused by the presence of homopolymer stretches, such as polyA/T in cDNA libraries. PolvA tails are shortened by cleaving before cDNA cloning with type IIS restriction enzymes, such as GsuI, placed next to the oligo-dT used to prime the polyA tails of mRNAs. We constructed four rice Cap-Trapper-selected. full-length normalized cDNA libraries, of which the average residual polyA tail was 4 bases or shorter in most of the clones analyzed. Because of the removal of homopolymeric stretches, libraries prepared with this method can be used for direct sequencing and transcriptional sequencing without the slippage observed for libraries prepared with currently available methods, thus improving sequencing accuracy, operations, and throughput.