OryzaGenome2.1: Database of Diverse Genotypes in Wild Oryza Species.

BACKGROUND: OryzaGenome ( http://viewer.shigen.info/oryzagenome21detail/index.xhtml ), a feature within Oryzabase ( https://shigen.nig.ac.jp/rice/oryzabase/ ), is a genomic database for wild Oryza species that provides comparative and evolutionary genomics approaches for the rice research community. RESULTS: Here we release OryzaGenome2.1, the first major update of OryzaGenome. The main feature in this version is the inclusion of newly sequenced genotypes and their meta-information, giving a total of 217 accessions of 19 wild Oryza species (O. rufipogon, O. barthii, O. longistaminata, O. meridionalis, O. glumaepatula, O. punctata, O. minuta, O. officinalis, O. rhizomatis, O. eichingeri, O. latifolia, O. alta, O. grandiglumis, O. australiensis, O. brachyantha, O. granulata, O. meyeriana, O. ridleyi, and O. longiglumis). These 19 wild species belong to 9 genome types (AA, BB, CC, BBCC, CCDD, EE, FF, GG, and HHJJ), representing wide genomic diversity in the genus. Using the genotype information, we analyzed the genome diversity of Oryza species. Other features of OryzaGenome facilitate the use of information on single nucleotide polymorphisms (SNPs) between O. sativa and its wild progenitor O. rufipogon in rice research, including breeding as well as basic science. For example, we provide Variant Call Format (VCF) files for genome-wide SNPs of 33 O. rufipogon accessions against the O. sativa reference genome, IRGSP1.0. In addition, we provide a new SNP Effect Table function, allowing users to identify SNPs or small insertion/deletion polymorphisms in the 33 O. rufipogon accessions and to search for the effect of these polymorphisms on protein function if they reside in the coding region (e.g., are missense or nonsense mutations). Furthermore, the SNP Viewer for 446 O. rufipogon accessions was updated by implementing new tracks for possible selective sweep regions and highly mutated regions that were potentially exposed to selective pressures during the process of domestication. CONCLUSION: OryzaGenome2.1 focuses on comparative genomic analysis of diverse wild Oryza accessions collected around the world and on the development of resources to speed up the identification of critical trait-related genes, especially from O. rufipogon. It aims to promote the use of genotype information from wild accessions in rice breeding and potential future crop improvements. Diverse genotypes will be a key resource for evolutionary studies in Oryza, including polyploid biology.

Evolution and Diversity of the Wild Rice Oryza officinalis Complex, across Continents, Genome Types, and Ploidy Levels.

The Oryza officinalis complex is the largest species group in Oryza, with more than nine species from four continents, and is a tertiary gene pool that can be exploited in breeding programs for the improvement of cultivated rice. Most diploid and tetraploid members of this group have a C genome. Using a new reference C genome for the diploid species O. officinalis, and draft genomes for two other C genome diploid species Oryza eichingeri and Oryza rhizomatis, we examine the influence of transposable elements on genome structure and provide a detailed phylogeny and evolutionary history of the Oryza C genomes. The O. officinalis genome is 1.6 times larger than the A genome of cultivated Oryza sativa, mostly due to proliferation of Gypsy type long-terminal repeat transposable elements, but overall syntenic relationships are maintained with other Oryza genomes (A, B, and F). Draft genome assemblies of the two other C genome diploid species, Oryza eichingeri and Oryza rhizomatis, and short-read resequencing of a series of other C genome species and accessions reveal that after the divergence of the C genome progenitor, there was still a substantial degree of variation within the C genome species through proliferation and loss of both DNA and long-terminal repeat transposable elements. We provide a detailed phylogeny and evolutionary history of the Oryza C genomes and a genomic resource for the exploitation of the Oryza tertiary gene pool.

Comparative analysis of seven types of superoxide dismutases for their ability to respond to oxidative stress in Bombyx mori.

Insects are well adapted to changing environmental conditions. They have unique systems for eliminating reactive oxygen species (ROS). Superoxide dismutase (SOD) is a key enzyme that plays a primary role in removing ROS. Bombyx mori is a lepidopteran insect, whose body size is larger than the model insect Drosophila melanogaster, which enabled us to more easily examine gene expression at the tissue level. We searched B. mori SOD (BmSOD) genes using genome database, and we analyzed their function under different type of oxidative stress. Consequently, we identified four new types of BmSODs in addition to the three types already known. Two of the seven types had a unique domain architecture that has not been discovered previously in the SOD family, and they were expressed in different tissues and developmental stages. Furthermore, these BmSODs responded differently to several kinds of stressors. Our results showed that the seven types of BmSODs are likely to play different roles in B. mori; therefore, B. mori could be used to distinguish the functions of each SOD for resistance to oxidative stress that changes with the environmental conditions.

Effects of androgen and oestrogen on IGF pathways controlling phallus growth.

The development of the mammalian phallus involves hormone-dependent mesenchymal-epithelial signalling mechanisms that contribute to urethral closure and regulation of phallus elongation and growth. In marsupials, most differentiation and growth of the phallus occurs post-natally, making them amenable to direct hormone treatment. Expression of IGFs, FGFs, EFNB2, MAFB, DLX5 and AP-1 mRNAs in the phallus at day 50 post-partum (pp) were altered after treatment of tammar wallaby young from day 20 to 40 pp with androgen, oestrogen or after castration at day 25 pp. However, the most interesting changes occurred in the IGF pathway genes. Androgen treatment upregulated IGF1 in female phalluses and oestrogen treatment upregulated IGF1 in male phalluses, but it was downregulated by castration. IGFBP3 was higher in female phalluses and downregulated by androgen. IGF1 expression was higher in all untreated male than in female phalluses from day 50 to 150 pp, but IGFBP3 had the reverse pattern. At day 90 pp, when urethral closure in males is progressing and male phallus growth is accelerating. IGF1 and PCNA protein were only detected in the male urorectal septum, suggesting for the first time that closure and elongation may involve IGF1 activation of cell proliferation specifically in male phalluses. These effects of sex steroids on gene expression and on the IGF1 signalling pathway in particular, suggest that the developing phallus may be especially susceptible to perturbation by exogenous hormones.

Hormone-responsive genes in the SHH and WNT/ß-catenin signaling pathways influence urethral closure and phallus growth.

Environmental endocrine disruptors (EEDs) that affect androgen or estrogen activity may disrupt gene regulation during phallus development to cause hypospadias or a masculinized clitoris. We treated developing male tammar wallabies with estrogen and females with androgen from day 20-40 postpartum (pp) during the androgen imprinting window of sensitivity. Estrogen inhibited phallus elongation but had no effect on urethral closure and did not significantly depress testicular androgen synthesis. Androgen treatment in females did not promote phallus elongation but initiated urethral closure. Phalluses were collected for transcriptome sequencing at day 50 pp when they first become sexually dimorphic to examine changes in two signaling pathways, sonic hedgehog (SHH) and wingless-type MMTV integration site family (WNT)/ß-catenin. SHH mRNA and ß-catenin were predominantly expressed in the urethral epithelium in the tammar phallus, as in eutherian mammals. Estrogen treatment and castration of males induced an upregulation of SHH, while androgen treatment downregulated SHH. These effects appear to be direct since we detected putative estrogen receptor α (ERα) and androgen receptor (AR) binding sites near SHH. WNT5A, like SHH, was downregulated by androgen, while WNT4 was upregulated in female phalluses after androgen treatment. After estrogen treatment, WIF1 and WNT7A were both downregulated in male phalluses. After castration, WNT9A was upregulated. These results suggest that SHH and WNT pathways are regulated by both estrogen and androgen to direct the proliferation and elongation of the phallus during differentiation. Their response to exogenous hormones makes these genes potential targets of EEDs in the etiology of abnormal phallus development including hypospadias.

Expansions of intronic TTTCA and TTTTA repeats in benign adult familial myoclonic epilepsy.

Epilepsy is a common neurological disorder, and mutations in genes encoding ion channels or neurotransmitter receptors are frequent causes of monogenic forms of epilepsy. Here we show that abnormal expansions of TTTCA and TTTTA repeats in intron 4 of SAMD12 cause benign adult familial myoclonic epilepsy (BAFME). Single-molecule, real-time sequencing of BAC clones and nanopore sequencing of genomic DNA identified two repeat configurations in SAMD12. Intriguingly, in two families with a clinical diagnosis of BAFME in which no repeat expansions in SAMD12 were observed, we identified similar expansions of TTTCA and TTTTA repeats in introns of TNRC6A and RAPGEF2, indicating that expansions of the same repeat motifs are involved in the pathogenesis of BAFME regardless of the genes in which the expanded repeats are located. This discovery that expansions of noncoding repeats lead to neuronal dysfunction responsible for myoclonic tremor and epilepsy extends the understanding of diseases with such repeat expansion.

Construction of a simple evaluation system for the intestinal absorption of an orally administered medicine using Bombyx mori larvae.

Human intestinal absorption is estimated using a human colon carcinoma cell line (Caco-2) cells from human colorectal adenocarcinoma, intestinal perfusion, or a mammalian model. These current evaluation systems are limited in their ability to estimate human intestinal absorption. In addition, in vivo evaluation systems using laboratory animals such as mice and rats entail animal ethics problems, and it is difficult to screen compounds on a large scale at the drug discovery stage. Thus, we propose the use of Bombyx mori larvae for evaluation of intestinal absorption of compounds as an alternative system in this study. First, to compare the characteristics among Caco-2 cells, human intestine, and B. mori larval midgut, we analyzed their RNA-seq data, and we found 26 drug transporters common to humans and B. mori. Next, we quantitatively developed an oral administration technique in B. mori and established a method using silkworm B. mori larvae that can easily estimate the intestinal permeability of compounds. Consequently, we could determine the dose and technique for oral administration in B. mori larvae. We also developed a B. mori model to evaluate the intestinal permeability of orally administered. Our constructed evaluation system will be useful for evaluating intestinal permeability in medical drug development.

Acetylation of histone H4 lysine 5 and 12 is required for CENP-A deposition into centromeres.

Centromeres are specified epigenetically through the deposition of the centromere-specific histone H3 variant CENP-A. However, how additional epigenetic features are involved in centromere specification is unknown. Here, we find that histone H4 Lys5 and Lys12 acetylation (H4K5ac and H4K12ac) primarily occur within the pre-nucleosomal CENP-A-H4-HJURP (CENP-A chaperone) complex, before centromere deposition. We show that H4K5ac and H4K12ac are mediated by the RbAp46/48-Hat1 complex and that RbAp48-deficient DT40 cells fail to recruit HJURP to centromeres and do not incorporate new CENP-A at centromeres. However, C-terminally-truncated HJURP, that does not bind CENP-A, does localize to centromeres in RbAp48-deficient cells. Acetylation-dead H4 mutations cause mis-localization of the CENP-A-H4 complex to non-centromeric chromatin. Crucially, CENP-A with acetylation-mimetic H4 was assembled specifically into centromeres even in RbAp48-deficient DT40 cells. We conclude that H4K5ac and H4K12ac, mediated by RbAp46/48, facilitates efficient CENP-A deposition into centromeres.

Single-neuron and genetic correlates of autistic behavior in macaque.

Atypical neurodevelopment in autism spectrum disorder is a mystery, defying explanation despite increasing attention. We report on a Japanese macaque that spontaneously exhibited autistic traits, namely, impaired social ability as well as restricted and repetitive behaviors, along with our single-neuron and genomic analyses. Its social ability was measured in a turn-taking task, where two monkeys monitor each other's actions for adaptive behavioral planning. In its brain, the medial frontal neurons responding to others' actions, abundant in the controls, were almost nonexistent. In its genes, whole-exome sequencing and copy number variation analyses identified rare coding variants linked to human neuropsychiatric disorders in 5-hydroxytryptamine (serotonin) receptor 2C (HTR2C) and adenosine triphosphate (ATP)-binding cassette subfamily A13 (ABCA13). This combination of systems neuroscience and cognitive genomics in macaques suggests a new, phenotype-to-genotype approach to studying mental disorders.

The Gonium pectorale genome demonstrates co-option of cell cycle regulation during the evolution of multicellularity.

The transition to multicellularity has occurred numerous times in all domains of life, yet its initial steps are poorly understood. The volvocine green algae are a tractable system for understanding the genetic basis of multicellularity including the initial formation of cooperative cell groups. Here we report the genome sequence of the undifferentiated colonial alga, Gonium pectorale, where group formation evolved by co-option of the retinoblastoma cell cycle regulatory pathway. Significantly, expression of the Gonium retinoblastoma cell cycle regulator in unicellular Chlamydomonas causes it to become colonial. The presence of these changes in undifferentiated Gonium indicates extensive group-level adaptation during the initial step in the evolution of multicellularity. These results emphasize an early and formative step in the evolution of multicellularity, the evolution of cell cycle regulation, one that may shed light on the evolutionary history of other multicellular innovations and evolutionary transitions.

Rapid turnover of antimicrobial-type cysteine-rich protein genes in closely related Oryza genomes.

Defensive and reproductive protein genes undergo rapid evolution. Small, cysteine-rich secreted peptides (CRPs) act as antimicrobial agents and function in plant intercellular signaling and are over-represented among reproductively expressed proteins. Because of their roles in defense, reproduction and development and their presence in multigene families, CRP variation can have major consequences for plant phenotypic and functional diversification. We surveyed the CRP genes of six closely related Oryza genomes comprising Oryza sativa ssp. japonica and ssp. indica, Oryza glaberrima and three accessions of Oryza rufipogon to observe patterns of evolution in these gene families and the effects of variation on their gene expression. These Oryza genomes, like other plant genomes, have accumulated large reservoirs of CRP sequences, comprising 26 groups totaling between 676 and 843 genes, in contrast to antimicrobial CRPs in animal genomes. Despite the close evolutionary relationships between the genomes, we observed rapid changes in number and structure among CRP gene families. Many CRP sequences are in gene clusters generated by local duplications, have undergone rapid turnover and are more likely to be silent or specifically expressed. By contrast, conserved CRP genes are more likely to be highly and broadly expressed. Variable CRP genes created by repeated duplication, gene modification and inactivation can gain new functions and expression patterns in newly evolved gene copies. For the CRP proteins, the process of gain/loss by deletion or duplication at gene clusters seems to be an important mechanism in evolution of the gene families, which also contributes to their expression evolution.

A biochemical landscape of A-to-I RNA editing in the human brain transcriptome.

Inosine is an abundant RNA modification in the human transcriptome and is essential for many biological processes in modulating gene expression at the post-transcriptional level. Adenosine deaminases acting on RNA (ADARs) catalyze the hydrolytic deamination of adenosines to inosines (A-to-I editing) in double-stranded regions. We previously established a biochemical method called "inosine chemical erasing" (ICE) to directly identify inosines on RNA strands with high reliability. Here, we have applied the ICE method combined with deep sequencing (ICE-seq) to conduct an unbiased genome-wide screening of A-to-I editing sites in the transcriptome of human adult brain. Taken together with the sites identified by the conventional ICE method, we mapped 19,791 novel sites and newly found 1258 edited mRNAs, including 66 novel sites in coding regions, 41 of which cause altered amino acid assignment. ICE-seq detected novel editing sites in various repeat elements as well as in short hairpins. Gene ontology analysis revealed that these edited mRNAs are associated with transcription, energy metabolism, and neurological disorders, providing new insights into various aspects of human brain functions.

ERBB4 mutations that disrupt the neuregulin-ErbB4 pathway cause amyotrophic lateral sclerosis type 19.

Amyotrophic lateral sclerosis (ALS) is a devastating neurological disorder characterized by the degeneration of motor neurons and typically results in death within 3-5 years from onset. Familial ALS (FALS) comprises 5%-10% of ALS cases, and the identification of genes associated with FALS is indispensable to elucidating the molecular pathogenesis. We identified a Japanese family affected by late-onset, autosomal-dominant ALS in which mutations in genes known to be associated with FALS were excluded. A whole- genome sequencing and parametric linkage analysis under the assumption of an autosomal-dominant mode of inheritance with incomplete penetrance revealed the mutation c.2780G>A (p. Arg927Gln) in ERBB4. An extensive mutational analysis revealed the same mutation in a Canadian individual with familial ALS and a de novo mutation, c.3823C>T (p. Arg1275Trp), in a Japanese simplex case. These amino acid substitutions involve amino acids highly conserved among species, are predicted as probably damaging, and are located within a tyrosine kinase domain (p. Arg927Gln) or a C-terminal domain (p. Arg1275Trp), both of which mediate essential functions of ErbB4 as a receptor tyrosine kinase. Functional analysis revealed that these mutations led to a reduced autophosphorylation of ErbB4 upon neuregulin-1 (NRG-1) stimulation. Clinical presentations of the individuals with mutations were characterized by the involvement of both upper and lower motor neurons, a lack of obvious cognitive dysfunction, and relatively slow progression. This study indicates that disruption of the neuregulin-ErbB4 pathway is involved in the pathogenesis of ALS and potentially paves the way for the development of innovative therapeutic strategies such using NRGs or their agonists to upregulate ErbB4 functions.

Analysis of DNA repair and protection in the Tardigrade Ramazzottius varieornatus and Hypsibius dujardini after exposure to UVC radiation.

Tardigrades inhabiting terrestrial environments exhibit extraordinary resistance to ionizing radiation and UV radiation although little is known about the mechanisms underlying the resistance. We found that the terrestrial tardigrade Ramazzottius varieornatus is able to tolerate massive doses of UVC irradiation by both being protected from forming UVC-induced thymine dimers in DNA in a desiccated, anhydrobiotic state as well as repairing the dimers that do form in the hydrated animals. In R. varieornatus accumulation of thymine dimers in DNA induced by irradiation with 2.5 kJ/m(2) of UVC radiation disappeared 18 h after the exposure when the animals were exposed to fluorescent light but not in the dark. Much higher UV radiation tolerance was observed in desiccated anhydrobiotic R. varieornatus compared to hydrated specimens of this species. On the other hand, the freshwater tardigrade species Hypsibius dujardini that was used as control, showed much weaker tolerance to UVC radiation than R. varieornatus, and it did not contain a putative phrA gene sequence. The anhydrobiotes of R. varieornatus accumulated much less UVC-induced thymine dimers in DNA than hydrated one. It suggests that anhydrobiosis efficiently avoids DNA damage accumulation in R. varieornatus and confers better UV radiation tolerance on this species. Thus we propose that UV radiation tolerance in tardigrades is due to the both high capacities of DNA damage repair and DNA protection, a two-pronged survival strategy.

Genome-wide analysis of DNA methylation and expression of microRNAs in breast cancer cells.

DNA methylation of promoters is linked to transcriptional silencing of protein-coding genes, and its alteration plays important roles in cancer formation. For example, hypermethylation of tumor suppressor genes has been seen in some cancers. Alteration of methylation in the promoters of microRNAs (miRNAs) has also been linked to transcriptional changes in cancers; however, no systematic studies of methylation and transcription of miRNAs have been reported. In the present study, to clarify the relation between DNA methylation and transcription of miRNAs, next-generation sequencing and microarrays were used to analyze the methylation and expression of miRNAs, protein-coding genes, other non-coding RNAs (ncRNAs), and pseudogenes in the human breast cancer cell lines MCF7 and the adriamycin (ADR) resistant cell line MCF7/ADR. DNA methylation in the proximal promoter of miRNAs is tightly linked to transcriptional silencing, as it is with protein-coding genes. In protein-coding genes, highly expressed genes have CpG-rich proximal promoters whereas weakly expressed genes do not. This is only rarely observed in other gene categories, including miRNAs. The present study highlights the epigenetic similarities and differences between miRNA and protein-coding genes.

The medaka zic1/zic4 mutant provides molecular insights into teleost caudal fin evolution.

Teleosts have an asymmetrical caudal fin skeleton formed by the upward bending of the caudal-most portion of the body axis, the ural region. This homocercal type of caudal fin ensures powerful and complex locomotion and is regarded as one of the most important innovations for teleosts during adaptive radiation in an aquatic environment. However, the mechanisms that create asymmetric caudal fin remain largely unknown. The spontaneous medaka (teleost fish) mutant, Double anal fin (Da), exhibits a unique symmetrical caudal skeleton that resembles the diphycercal type seen in Polypterus and Coelacanth. We performed a detailed analysis of the Da mutant to obtain molecular insight into caudal fin morphogenesis. We first demonstrate that a large transposon, inserted into the enhancer region of the zic1 and zic4 genes (zic1/zic4) in Da, is associated with the mesoderm-specific loss of their transcription. We then show that zic1/zic4 are strongly expressed in the dorsal part of the ural mesenchyme and thereby induce asymmetric caudal fin development in wild-type embryos, whereas their expression is lost in Da. Comparative analysis further indicates that the dorsal mesoderm expression of zic1/zic4 is conserved in teleosts, highlighting the crucial role of zic1/zic4 in caudal fin morphogenesis.

Using the Acropora digitifera genome to understand coral responses to environmental change.

Despite the enormous ecological and economic importance of coral reefs, the keystone organisms in their establishment, the scleractinian corals, increasingly face a range of anthropogenic challenges including ocean acidification and seawater temperature rise. To understand better the molecular mechanisms underlying coral biology, here we decoded the approximately 420-megabase genome of Acropora digitifera using next-generation sequencing technology. This genome contains approximately 23,700 gene models. Molecular phylogenetics indicate that the coral and the sea anemone Nematostella vectensis diverged approximately 500 million years ago, considerably earlier than the time over which modern corals are represented in the fossil record (∼240 million years ago). Despite the long evolutionary history of the endosymbiosis, no evidence was found for horizontal transfer of genes from symbiont to host. However, unlike several other corals, Acropora seems to lack an enzyme essential for cysteine biosynthesis, implying dependency of this coral on its symbionts for this amino acid. Corals inhabit environments where they are frequently exposed to high levels of solar radiation, and analysis of the Acropora genome data indicates that the coral host can independently carry out de novo synthesis of mycosporine-like amino acids, which are potent ultraviolet-protective compounds. In addition, the coral innate immunity repertoire is notably more complex than that of the sea anemone, indicating that some of these genes may have roles in symbiosis or coloniality. A number of genes with putative roles in calcification were identified, and several of these are restricted to corals. The coral genome provides a platform for understanding the molecular basis of symbiosis and responses to environmental changes.

Identification of RNF213 as a susceptibility gene for moyamoya disease and its possible role in vascular development.

BACKGROUND: Moyamoya disease is an idiopathic vascular disorder of intracranial arteries. Its susceptibility locus has been mapped to 17q25.3 in Japanese families, but the susceptibility gene is unknown. METHODOLOGY/PRINCIPAL FINDINGS: Genome-wide linkage analysis in eight three-generation families with moyamoya disease revealed linkage to 17q25.3 (P<10(-4)). Fine mapping demonstrated a 1.5-Mb disease locus bounded by D17S1806 and rs2280147. We conducted exome analysis of the eight index cases in these families, with results filtered through Ng criteria. There was a variant of p.N321S in PCMTD1 and p.R4810K in RNF213 in the 1.5-Mb locus of the eight index cases. The p.N321S variant in PCMTD1 could not be confirmed by the Sanger method. Sequencing RNF213 in 42 index cases confirmed p.R4810K and revealed it to be the only unregistered variant. Genotyping 39 SNPs around RNF213 revealed a founder haplotype transmitted in 42 families. Sequencing the 260-kb region covering the founder haplotype in one index case did not show any coding variants except p.R4810K. A case-control study demonstrated strong association of p.R4810K with moyamoya disease in East Asian populations (251 cases and 707 controls) with an odds ratio of 111.8 (P = 10(-119)). Sequencing of RNF213 in East Asian cases revealed additional novel variants: p.D4863N, p.E4950D, p.A5021V, p.D5160E, and p.E5176G. Among Caucasian cases, variants p.N3962D, p.D4013N, p.R4062Q and p.P4608S were identified. RNF213 encodes a 591-kDa cytosolic protein that possesses two functional domains: a Walker motif and a RING finger domain. These exhibit ATPase and ubiquitin ligase activities. Although the mutant alleles (p.R4810K or p.D4013N in the RING domain) did not affect transcription levels or ubiquitination activity, knockdown of RNF213 in zebrafish caused irregular wall formation in trunk arteries and abnormal sprouting vessels. CONCLUSIONS/SIGNIFICANCE: We provide evidence suggesting, for the first time, the involvement of RNF213 in genetic susceptibility to moyamoya disease.

Sequence analysis of the genome of an oil-bearing tree, Jatropha curcas L.

The whole genome of Jatropha curcas was sequenced, using a combination of the conventional Sanger method and new-generation multiplex sequencing methods. Total length of the non-redundant sequences thus obtained was 285 858 490 bp consisting of 120 586 contigs and 29 831 singlets. They accounted for ~95% of the gene-containing regions with the average G + C content was 34.3%. A total of 40 929 complete and partial structures of protein encoding genes have been deduced. Comparison with genes of other plant species indicated that 1529 (4%) of the putative protein-encoding genes are specific to the Euphorbiaceae family. A high degree of microsynteny was observed with the genome of castor bean and, to a lesser extent, with those of soybean and Arabidopsis thaliana. In parallel with genome sequencing, cDNAs derived from leaf and callus tissues were subjected to pyrosequencing, and a total of 21 225 unigene data have been generated. Polymorphism analysis using microsatellite markers developed from the genomic sequence data obtained was performed with 12 J. curcas lines collected from various parts of the world to estimate their genetic diversity. The genomic sequence and accompanying information presented here are expected to serve as valuable resources for the acceleration of fundamental and applied research with J. curcas, especially in the fields of environment-related research such as biofuel production. Further information on the genomic sequences and DNA markers is available at http://www.kazusa.or.jp/jatropha/.

The medaka draft genome and insights into vertebrate genome evolution.

Teleosts comprise more than half of all vertebrate species and have adapted to a variety of marine and freshwater habitats. Their genome evolution and diversification are important subjects for the understanding of vertebrate evolution. Although draft genome sequences of two pufferfishes have been published, 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 and developmental genetics. In the assembled medaka genome (700 megabases), which is less than half of the zebrafish genome, we predicted 20,141 genes, including approximately 2,900 new genes, using 5'-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 approximately 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.