Combining our full-length cDNA method and the massively parallel sequencing technology, we developed a simple method to collect precise positional information of transcriptional start sites (TSSs) together with digital information of the gene-expression levels in a high throughput manner. We applied this method to observe gene-expression changes in a colon cancer cell line cultured in normoxic and hypoxic conditions. We generated more than 100 million 36-base TSS-tag sequences and revealed comprehensive features of hypoxia responsive alterations in the transcriptional landscape of the human genome. The features include presence of inducible 'hot regions' in 54 genomic regions, 220 novel hypoxia inducible promoters that may drive non-protein-coding transcripts, 191 hypoxia responsive alternative promoters and detailed views of 120 novel as well as known hypoxia responsive genes. We further analyzed hypoxic response of different cells using additional 60 million TSS-tags and found that the degree of the gene-expression changes were different among cell lines, possibly reflecting cellular robustness against hypoxia. The novel dynamic figure of the human gene transcriptome will deepen our understanding of the transcriptional program of the human genome as well as bringing new insights into the biology of cancer cells in hypoxia.

Using full-length cDNA sequences, we compared alternative splicing (AS) in humans and mice. The alignment of the human and mouse genomes showed that 86% of 199 426 total exons in human AS variants were conserved in the mouse genome. Of the 20 392 total human AS variants, however, 59% consisted of all conserved exons. Comparing AS patterns between human and mouse transcripts revealed that only 431 transcripts from 189 loci were perfectly conserved AS variants. To exclude the possibility that the full-length human cDNAs used in the present study, especially those with retained introns, were cloning artefacts or prematurely spliced transcripts, we experimentally validated 34 such cases. Our results indicate that even retained-intron type transcripts are typically expressed in a highly controlled manner and interact with translating ribosomes. We found non-conserved AS exons to be predominantly outside the coding sequences (CDSs). This suggests that non-conserved exons in the CDSs of transcripts cause functional constraint. These findings should enhance our understanding of the relationship between AS and species specificity of human genes.

Full-Parasites (http://fullmal.hgc.jp/) is a transcriptome database of apicomplexa parasites, which include Plasmodium and Toxoplasma species. The latest version of Full-Parasites contains a total of 105,786 EST sequences from 12 parasites, of which 5925 full-length cDNAs have been completely sequenced. Full-Parasites also contain more than 30 million transcription start sites (TSS) for Plasmodium falciparum (Pf) and Toxoplasma gondii (Tg), which were identified using our novel oligo-capping-based protocol. Various types of cDNA data resources were interconnected with our original database functionalities. Specifically, in this update, we have included two unique RNA-Seq data sets consisting of 730 million mapped RNA-Seq tags. One is a dataset of 16 time-lapse experiments of cultured bradyzoite differentiation for Tg. The other dataset includes 31 clinical samples of Pf. Parasite RNA was extracted together with host human RNA, and the extracted mixed RNA was used for RNA sequencing, with the expectation that gene expression information from the host and parasite would be simultaneously represented. By providing the largest unique full-length cDNA and dynamic transcriptome data, Full-Parasites is useful for understanding host-parasite interactions and will help to eventually elucidate how monophyletic organisms have evolved to become parasites by adopting complex life cycles.

Understanding the genome sequence-specific positioning of nucleosomes is essential to understand various cellular processes, such as transcriptional regulation and replication. As a typical example, the 10-bp periodicity of AA/TT and GC dinucleotides has been reported in several species, but it is still unclear whether this feature can be observed in the whole genomes of all eukaryotes.

Massively parallel, tag-based sequencing systems, such as the SOLiD system, hold the promise of revolutionizing the study of whole genome gene expression due to the number of data points that can be generated in a simple and cost-effective manner. We describe the development of a 5'-end transcriptome workflow for the SOLiD system and demonstrate the advantages in sensitivity and dynamic range offered by this tag-based application over traditional approaches for the study of whole genome gene expression. 5'-end transcriptome analysis was used to study whole genome gene expression within a colon cancer cell line, HT-29, treated with the DNA methyltransferase inhibitor, 5-aza-2'-deoxycytidine (5Aza). More than 20 million 25-base 5'-end tags were obtained from untreated and 5Aza-treated cells and matched to sequences within the human genome. Seventy three percent of the mapped unique tags were associated with RefSeq cDNA sequences, corresponding to approximately 14,000 different protein-coding genes in this single cell type. The level of expression of these genes ranged from 0.02 to 4,704 transcripts per cell. The sensitivity of a single sequence run of the SOLiD platform was 100-1,000 fold greater than that observed from 5'end SAGE data generated from the analysis of 70,000 tags obtained by Sanger sequencing. The high-resolution 5'end gene expression profiling presented in this study will not only provide novel insight into the transcriptional machinery but should also serve as a basis for a better understanding of cell biology.

High-temperature fermentation technology with thermotolerant microbes has been expected to reduce the cost of bioconversion of cellulosic biomass to fuels or chemicals. Thermotolerant Kluyveromyces marxianus possesses intrinsic abilities to ferment and assimilate a wide variety of substrates including xylose and to efficiently produce proteins. These capabilities have been found to exceed those of the traditional ethanol producer Saccharomyces cerevisiae or lignocellulose-bioconvertible ethanologenic Scheffersomyces stipitis.

Chromosome condensation is a hallmark of mitosis in eukaryotes and is a prerequisite for faithful segregation of genetic material to daughter cells. Here we show that condensin, which is essential for assembling condensed chromosomes, helps to preclude the detrimental effects of gene transcription on mitotic condensation. ChIP-seq profiling reveals that the fission yeast condensin preferentially binds to active protein-coding genes in a transcription-dependent manner during mitosis. Pharmacological and genetic attenuation of transcription largely rescue bulk chromosome segregation defects observed in condensin mutants. We also demonstrate that condensin is associated with and reduces unwound DNA segments generated by transcription, providing a direct link between an in vitro activity of condensin and its in vivo function. The human condensin isoform condensin I also binds to unwound DNA regions at the transcription start sites of active genes, implying that our findings uncover a fundamental feature of condensin complexes.

Ecdysteroidogenesis is essential for arthropod development and reproduction. Although the importance of ecdysteroids has been demonstrated, there is little information on the sites and enzymes for synthesis of ecdysteroids from Chelicerates. Ecdysteroid functions have been well studied in the soft tick Ornithodoros moubata, making this species an excellent candidate for elucidating ecdysteroidogenesis in Chelicerates. Results showed that O. moubata has at least two ecdysteroidogenic enzymes, Spook (OmSpo) and Shade (OmShd). RNAi showed both enzymes were required for ecdysteroidogenesis. Enzymatic assays demonstrated OmShd has the conserved functions of ecdysone 20-hydroxylase. OmSpo showed specific expression in the ovaries of final nymphal and adult stages, indicating O. moubata utilizes the ovary as an ecdysteroidogenic tissue instead of specific tissues as seen in other arthropods. On the other hand, OmShd expression was observed in various tissues including the midgut, indicating functional ecdysteroids can be produced in these tissues. In nymphal stages, expression of both OmSpo and OmShd peaked before molting corresponding with high ecdysteroid titers in the hemolymph. In fed adult females, OmSpo expression peaked at 8-10 days after engorgement, while OmShd expression peaked immediately after engorgement. Mated females showed more frequent surges of OmShd than virgin females. These results indicate that the regulation of synthesis of ecdysteroids differs in nymphs and adult females, and mating modifies adult female ecdysteroidogenesis. This is the first report to focus on synthesis of ecdysteroids in ticks and provides essential knowledge for understanding the evolution of ecdysteroidogenesis in arthropods.

Neospora caninum is a protozoan parasite that causes neurological disorders in dogs and cattle. It can cause nonsuppurative meningoencephalitis and a variety of neuronal symptoms are observed, particularly in dogs. However, the pathogenic mechanism, including the relationship between the parasite distribution and the clinical signs, is unclear. In this study, to understand the pathogenic mechanism of neosporosis, parasite distribution and lesions were assessed in the brain of mice infected with N. caninum (strain Nc-1). Host gene expression was also analyzed with RNA sequencing (RNA-Seq). The histopathological lesions in the frontal lobe and the medulla oblongata were significantly more severe in symptomatic mice than in asymptomatic mice, although no association between the severity of the lesions and parasite numbers was found. In infected mice, the expression of 772 mouse brain genes was upregulated. A GOstat analysis predicted that the upregulated genes were involved in the host immune response. Genes whose expression correlated positively and negatively with parasite numbers were involved in the host immune response, and neuronal morphogenesis and lipid metabolic processes, respectively. These results suggest that changes in the gene expression profile associated with neuronal functions as well as immune responses can contribute to the pathogenesis in N. caninum-infected animals.

Teleost intestine is crucial for seawater acclimation by sensing osmolality of imbibed seawater and regulating drinking and water/ion absorption. Regulatory genes for transforming intestinal function have not been identified. A transcriptomic approach was used to search for such genes in the intestine of euryhaline medaka.

DBTSS (http://dbtss.hgc.jp/) was originally constructed as a collection of uniquely determined transcriptional start sites (TSSs) in humans and some other species in 2002. Since then, it has been regularly updated and in recent updates epigenetic information has also been incorporated because such information is useful for characterizing the biological relevance of these TSSs/downstream genes. In the newest release, Release 9, we further integrated public and original single nucleotide variation (SNV) data into our database. For our original data, we generated SNV data from genomic analyses of various cancer types, including 97 lung adenocarcinomas and 57 lung small cell carcinomas from Japanese patients as well as 26 cell lines of lung cancer origin. In addition, we obtained publically available SNV data from other cancer types and germline variations in total of 11,322 individuals. With these updates, users can examine the association between sequence variation pattern in clinical lung cancers with its corresponding TSS-seq, RNA-seq, ChIP-seq and BS-seq data. Consequently, DBTSS is no longer a mere storage site for TSS information but has evolved into an integrative platform of a variety of genome activity data.

Mutually exclusive selection of one exon in a cluster of exons is a rare form of alternative pre-mRNA splicing, yet suggests strict regulation. However, the repertoires of regulation mechanisms for the mutually exclusive (ME) splicing in vivo are still unknown. Here, we experimentally explore putative ME exons in C. elegans to demonstrate that 29 ME exon clusters in 27 genes are actually selected in a mutually exclusive manner. Twenty-two of the clusters consist of homologous ME exons. Five clusters have too short intervening introns to be excised between the ME exons. Fidelity of ME splicing relies at least in part on nonsense-mediated mRNA decay for 14 clusters. These results thus characterize all the repertoires of ME splicing in this organism.

We detected and characterized the binding sites of the representative Rest complex components Rest, Sin3A, and Lsd1. We compared their binding patterns in mouse embryonic stem (ES) cells and epiblast stem (EpiS) cells. We found few Rest sites unique to the EpiS cells. The ES-unique site features were distinct from those of the common sites, namely, the signal intensities were weaker, and the characteristic gene function categories differed. Our analyses showed that the Rest binding sites do not always overlap with the Sin3A and Lsd1 binding sites. The Sin3A binding pattern differed remarkably between the ES and EpiS cells and was accompanied by significant changes in acetylated-histone patterns in the surrounding regions. A series of transcriptome analyses in the same cell types unexpectedly showed that the putative target gene transcript levels were not dramatically different despite dynamic changes in the Rest complex binding patterns and chromatin statuses, which suggests that Rest is not the sole determinant of repression at its targets. Nevertheless, we identified putative Rest targets with explicitly enhanced transcription upon Rest knock-down in 143 and 60 common and ES-unique Rest target genes, respectively. Among such sites, several genes are involved in ES cell proliferation. In addition, we also found that long, intergenic non-coding RNAs were apparent Rest targets and shared similar features with the protein-coding target genes. Interestingly, such non-coding target genes showed less conservation through evolution than protein-coding targets. As a result of differences in the components and targets of the Rest complex, its functional roles may differ in ES and EpiS cells.

The Amoebozoa represent a clade of unicellular amoeboid organisms that display a wide variety of lifestyles, including free-living and parasitic species. For example, the social amoeba Dictyostelium discoideum has the ability to aggregate into a multicellular fruiting body upon starvation, while the pathogenic amoeba Entamoeba histolytica is a parasite of humans. Globins are small heme proteins that are present in almost all extant organisms. Although several genomes of amoebozoan species have been sequenced, little is known about the phyletic distribution of globin genes within this phylum. Only two flavohemoglobins (FHbs) of D. discoideum have been reported and characterized previously while the genomes of Entamoeba species are apparently devoid of globin genes. We investigated eleven amoebozoan species for the presence of globin genes by genomic and phylogenetic in silico analyses. Additional FHb genes were identified in the genomes of four social amoebas and the true slime mold Physarum polycephalum. Moreover, a single-domain globin (SDFgb) of Hartmannella vermiformis, as well as two truncated hemoglobins (trHbs) of Acanthamoeba castellanii were identified. Phylogenetic evidence suggests that these globin genes were independently acquired via horizontal gene transfer from some ancestral bacteria. Furthermore, the phylogenetic tree of amoebozoan FHbs indicates that they do not share a common ancestry and that a transfer of FHbs from bacteria to amoeba occurred multiple times.

Genome-wide distribution of the majority of H2A and H3 variants (H2A, H2AX, H2AZ, macroH2A, H3.1, H3.2 and H3.3) was simultaneously investigated in mouse embryonic stem cells by chromatin immunoprecipitation sequencing. Around the transcription start site, histone variant distribution differed between genes possessing promoters of high and low CpG density, regardless of their expression levels. In the intergenic regions, regulatory elements were enriched in H2A.Z and H3.3, whereas repeat elements were abundant in H2A and macroH2A, and H3.1, respectively. Analysis of H2A and H3 variant combinations composing nucleosomes revealed that the H2A.Z and H3.3 combinations were present at a higher frequency throughout the genome than the other combinations, suggesting that H2A.Z and H3.3 associate preferentially with each other to comprise the nucleosomes independently of genome region. Finally, we found that chromatin was unstable only in regions where it was enriched in both H2A.Z and H3.3, but strongly quantified stable in regions in which only H3.3 was abundant. Therefore, histone variant composition is an important determinant of chromatin structure, which is associated with specific genomic functions.

Alveolar echinococcosis is a worldwide zoonosis of great public health concern. Analysis of genome data for Echinococcus multilocularis has identified antigen families that can be used in diagnostic assays and vaccine development. However, little gene expression data is available for antigens of the egg and early larval stages. To address this information gap, we used a Next-Generation Sequencing approach to investigate three different stages (non-activated and activated oncospheres, and early stage metacestodes) of E. multilocularis (Nemuro strain). Transcriptome data analysis revealed that some diagnostic antigen gp50 isoforms and the antigen Eg95 family dominated in activated oncospheres, and the antigen B family dominated in early stage metacestodes. Furthermore, heat shock proteins and antigen II/3 are constantly expressed in the three stages. The expression pattern of various known antigens in E. multilocularis may give fundamental information for choosing candidate genes used in diagnosis and vaccine development.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

Although long non-coding RNAs (lncRNAs) have been associated with a variety of cancers, the interplay between lncRNAs and androgen receptor signaling in prostate cancer is still unclear. We identified an androgen-dependent lncRNA, POTEF-AS1, whose expression was regulated by androgen receptor in two androgen-dependent cells by using directional RNA sequencing analysis. POTEF-AS1 promoted cell growth, repressed genes related to the Toll-like receptor signaling and apoptosis pathways, and inhibited apoptosis in docetaxel-treated LNCaP cells. These findings suggest that POTEF-AS1 would play a key role in the progression of prostate cancer by repressing Toll-like receptor signaling.

Functionally polarized CD4+ T helper (Th) cells such as Th1, Th2 and Th17 cells are central to the regulation of acquired immunity. However, the molecular mechanisms governing the maintenance of the polarized functions of Th cells remain unclear. GATA3, a master regulator of Th2 cell differentiation, initiates the expressions of Th2 cytokine genes and other Th2-specific genes. GATA3 also plays important roles in maintaining Th2 cell function and in continuous chromatin remodeling of Th2 cytokine gene loci. However, it is unclear whether continuous expression of GATA3 is required to maintain the expression of various other Th2-specific genes. In this report, genome-wide DNA gene expression profiling revealed that GATA3 expression is critical for the expression of a certain set of Th2-specific genes. We demonstrated that GATA3 dependency is reduced for some Th2-specific genes in fully developed Th2 cells compared to that observed in effector Th2 cells, whereas it is unchanged for other genes. Moreover, effects of a loss of GATA3 expression in Th2 cells on the expression of cytokine and cytokine receptor genes were examined in detail. A critical role of GATA3 in the regulation of Th2-specific gene expression is confirmed in in vivo generated antigen-specific memory Th2 cells. Therefore, GATA3 is required for the continuous expression of the majority of Th2-specific genes involved in maintaining the Th2 cell identity.

The innate immune response is primarily mediated by the Toll-like receptors functioning through the MyD88-dependent and TRIF-dependent pathways. Despite being widely studied, it is not yet completely understood and systems-level analyses have been lacking. In this study, we identified a high-probability network of genes activated during the innate immune response using a novel approach to analyze time-course gene expression profiles of activated immune cells in combination with a large gene regulatory and protein-protein interaction network. We classified the immune response into three consecutive time-dependent stages and identified the most probable paths between genes showing a significant change in expression at each stage. The resultant network contained several novel and known regulators of the innate immune response, many of which did not show any observable change in expression at the sampled time points. The response network shows the dominance of genes from specific functional classes during different stages of the immune response. It also suggests a role for the protein phosphatase 2a catalytic subunit α in the regulation of the immunoproteasome during the late phase of the response. In order to clarify the differences between the MyD88-dependent and TRIF-dependent pathways in the innate immune response, time-course gene expression profiles from MyD88-knockout and TRIF-knockout dendritic cells were analyzed. Their response networks suggest the dominance of the MyD88-dependent pathway in the innate immune response, and an association of the circadian regulators and immunoproteasomal degradation with the TRIF-dependent pathway. The response network presented here provides the most probable associations between genes expressed in the early and the late phases of the innate immune response, while taking into account the intermediate regulators. We propose that the method described here can also be used in the identification of time-dependent gene sub-networks in other biological systems.