Teleosts comprise more than half of all vertebrate species and have adapted to a variety of marine and freshwater habitats 1 . Their genome evolution and diversification are important subjects for the understanding of vertebrate evolution. Although draft genome sequences of two pufferfishes have been published 2,3 , analysis of more fish genomes is desirable. Here we report a high-quality draft genome sequence of a small egg-laying freshwater teleost, medaka (Oryzias latipes). Medaka is native to East Asia and an excellent model system for a wide range of biology, including ecotoxicology, carcinogenesis, sex determination 4-6 and developmental genetics 7 . In the assembled medaka genome (700 megabases), which is less than half of the zebrafish genome, we predicted 20,141 genes, including 2,900 new genes, using 59-end serial analysis of gene expression tag information. We found single nucleotide polymorphisms (SNPs) at an average rate of 3.42% between the two inbred strains derived from two regional populations; this is the highest SNP rate seen in any vertebrate species. Analyses based on the dense SNP information show a strict genetic separation of 4 million years (Myr) between the two populations, and suggest that differential selective pressures acted on specific gene categories. Four-way comparisons with the human, pufferfish (Tetraodon), zebrafish and medaka genomes revealed that eight major interchromosomal rearrangements took place in a remarkably short period of 50 Myr after the whole-genome duplication event in the teleost ancestor and afterwards, intriguingly, the medaka genome preserved its ancestral karyotype for more than 300 Myr.We applied the whole-genome shotgun approach to an inbred strain, , derived from the southern Japanese population, as the main target. A total of 13.8 million reads amounting to approximately 10.6-fold genome coverage were obtained from the shotgun plasmid, fosmid and bacterial artificial chromosome (BAC) libraries. A newly developed RAMEN assembler was used to process the shotgun reads to generate contigs and scaffolds. The N50 values (50% of nucleotides in an assembly are in scaffolds-or contigs-longer than or equal to the N50 value) are ,1.41 megabases (Mb) for scaffolds and ,9.8 kilobases (Kb) for contigs. The total length of the contigs reached 700.4 Mb, which, from now on, we refer to as the medaka genome size.To construct ultracontigs, the scaffolds were integrated with the medaka genetic map by using SNP markers. For this purpose, we further obtained about 2.8-fold coverage of shotgun reads from another inbred strain HNI (refs 9, 10), which is derived from the northern Japanese population. The reads were assembled by RAMEN to scaffolds covering 648 Mb. Aligning the HNI contigs with the HdrR genome using BLASTZ 11 , we identified 16.4 million SNPs as well as 1.40 million insertions and 1.45 million deletions in non-repetitive regions (Supplementary Table 2). We selected 2,401 SNPs and genetically mapped them onto medaka chromosomes using a backcross panel between the...

We observed here that the expression of B lymphocyte chemokine (BLC/CXCL13) was markedly enhanced in the thymus and kidney in aged (NZB × NZW)F1 (BWF1) mice developing lupus nephritis, but not in similarly aged NZB and NZW mice. BLC-positive cells were present in the cellular infiltrates in the target organs with a reticular pattern of staining. CD11b+CD11c+ dendritic cells were increased in the thymus and spleen in aged BWF1 mice and identified as the major cell source for BLC. CD4+ T cells as well as B cells were dramatically increased in the thymus in aged BWF1 mice, whereas no increase was observed in aged NZB and NZW mice. B1/B2 ratio in the thymus was significantly higher than those in the spleen and peripheral blood in aged BWF1 mice. Interestingly, BLC showed preferential chemotactic activity for B1 cells derived from several mouse strains, including nonautoimmune mice. Cell surface CXCR5 expression on B1 cells was significantly higher than that on B2 cells. Thus, aberrant high expression of BLC by myeloid dendritic cells in the target organs in aged BWF1 mice may play a pivotal role in breaking immune tolerance in the thymus and in recruiting autoantibody-producing B cells in the development of murine lupus.

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.

We performed a genome-wide analysis of transcriptional start sites (TSSs) in human genes by multifaceted use of a massively parallel sequencer. By analyzing 800 million sequences that were obtained from various types of transcriptome analyses, we characterized 140 million TSS tags in 12 human cell types. Despite the large number of TSS clusters (TSCs), the number of TSCs was observed to decrease sharply with increasing expression levels. Highly expressed TSCs exhibited several characteristic features: Nucleosome-seq analysis revealed highly ordered nucleosome structures, ChIP-seq analysis detected clear RNA polymerase II binding signals in their surrounding regions, evaluations of previously sequenced and newly shotgunsequenced complete cDNA sequences showed that they encode preferable transcripts for protein translation, and RNA-seq analysis of polysome-incorporated RNAs yielded direct evidence that those transcripts are actually translated into proteins. We also demonstrate that integrative interpretation of transcriptome data is essential for the selection of putative alternative promoter TSCs, two of which also have protein consequences. Furthermore, discriminative chromatin features that separate TSCs at different expression levels were found for both genic TSCs and intergenic TSCs. The collected integrative information should provide a useful basis for future biological characterization of TSCs.

Pulmonary fibrosis is characterized by accumulation of activated fibroblasts that produce excessive amounts of extracellular matrix components such as collagen type I. However, the dynamics and activation signatures of fibroblasts during fibrogenesis remain poorly understood, especially in vivo. We examined changes in lung tissue cell populations and in the phenotype of activated fibroblasts after acute injury in a model of bleomycin-induced pulmonary fibrosis. Despite clustering of collagen type I-producing fibroblasts in fibrotic regions, flow cytometry-based quantitative analysis of whole lungs revealed that the number of fibroblasts in the lungs remained constant. At the peak of inflammation, fibroblast proliferation and apoptosis were both increased, suggesting that the clustering was not merely a result of proliferation, but also of fibroblast migration from nearby alveolar walls. Parabiosis experiments demonstrated that fibroblasts were not supplied from the circulation. Comprehensive gene expression analysis of freshly isolated fibroblasts revealed a detailed activation signature associated with fibrogenesis, including changes in genes responsible for migration and extracellular matrix construction. The Spp1 gene, which encodes osteopontin, was highly up-regulated and was an identifying characteristic of activated fibroblasts present at the sites of remodeling. Osteopontin may serve as a useful marker of profibrotic fibroblasts. These results provide insights into the cellular and molecular mechanisms underlying pulmonary fibrosis and provide a foundation for development of specific antifibrotic therapies.