Masahiro Kasahara

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

Chromatin-associated periodicity in genetic variation downstream of transcriptional start sites

Might DNA sequence variation reflect germline genetic activity and underlying chromatin structure? We investigated this question using medaka (Japanese killifish, Oryzias latipes), by comparing the genomic sequences of two strains (Hd-rR and HNI) and by mapping ∼37.3 million nucleosome cores from Hd-rR blastulae and 11,654 representative transcription start sites from six embryonic stages. We observed a distinctive ∼200–base pair (bp) periodic pattern of genetic variation downstream of transcription start sites; the rate of insertions and deletions longer than 1 bp peaked at positions of approximately +200, +400, and +600 bp, whereas the point mutation rate showed corresponding valleys. This ∼200-bp periodicity was correlated with the chromatin structure, with nucleosome occupancy minimized at positions 0, +200, +400, and +600 bp. These data exemplify the potential for genetic activity (transcription) and chromatin structure to contribute to molding the DNA sequence on an evolutionary time scale.

Large Scale Full-Length cDNA Sequencing Reveals a Unique Genomic Landscape in a Lepidopteran Model Insect, Bombyx mori

The establishment of a complete genomic sequence of silkworm, the model species of Lepidoptera, laid a foundation for its functional genomics. A more complete annotation of the genome will benefit functional and comparative studies and accelerate extensive industrial applications for this insect. To realize these goals, we embarked upon a large-scale full-length cDNA collection from 21 full-length cDNA libraries derived from 14 tissues of the domesticated silkworm and performed full sequencing by primer walking for 11,104 full-length cDNAs. The large average intron size was 1904 bp, resulting from a high accumulation of transposons. Using gene models predicted by GLEAN and published mRNAs, we identified 16,823 gene loci on the silkworm genome assembly. Orthology analysis of 153 species, including 11 insects, revealed that among three Lepidoptera including Monarch and Heliconius butterflies, the 403 largest silkworm-specific genes were composed mainly of protective immunity, hormone-related, and characteristic structural proteins. Analysis of testis-/ovary-specific genes revealed distinctive features of sexual dimorphism, including depletion of ovary-specific genes on the Z chromosome in contrast to an enrichment of testis-specific genes. More than 40% of genes expressed in specific tissues mapped in tissue-specific chromosomal clusters. The newly obtained FL-cDNA sequences enabled us to annotate the genome of this lepidopteran model insect more accurately, enhancing genomic and functional studies of Lepidoptera and comparative analyses with other insect orders, and yielding new insights into the evolution and organization of lepidopteran-specific genes.

Activation of a cyanobacterial adenylate cyclase, CyaC, by autophosphorylation and a subsequent phosphotransfer reaction.

The CyaC protein, a cyanobacterial adenylate cyclase, has a unique primary structure composed of the catalytic domain of adenylate cyclase and the conserved domains of bacterial two-component regulatory systems, one transmitter domain and two receiver domains. In the present work, CyaC was produced inEscherichia coli as a histidine-tagged recombinant protein and purified to homogeneity. CyaC showed ability to autophosphorylatein vitro with the γ-phosphate of [γ-32P]ATP. CyaC derivatives were constructed by site-directed mutagenesis in which the highly conserved phosphorylation sites in the transmitter domain (His572) and receiver domains (Asp60 or Asp895) were replaced by glutamine and alanine residues, respectively. After autophosphorylation of the CyaC derivatives, the chemical stabilities of the phosphoryl groups bound to the derivatives were determined. It was found that His572 is the initial phosphorylation site and that the phosphoryl group once bound to His572 is transferred to Asp895. The enzyme activities of the CyaC derivatives defective in His572 or Asp895 were considerably reduced. Asp895 is phosphorylated by acetyl [32P]phosphate, a small phosphoryl molecule, but Asp60 is not. Acetyl phosphate stimulates adenylate cyclase activity only when Asp895 is intact. These results suggest that the phosphorylation of Asp895 is essential for the activation of adenylate cyclase and that Asp60 functions differently from Asp895 in regulating the enzyme activity.

UTGB/medaka: genomic resource database for medaka biology

Medaka (Oryzias latipes) is a small egg-laying freshwater teleost native to East Asia that has become an excellent model system for developmental genetics and evolutionary biology. The draft medaka genome sequence (700 Mb) was reported in June 2007, and its substantial genomic resources have been opened to the public through the University of Tokyo Genome Browser Medaka (UTGB/medaka) database. This database provides basic genomic information, such as predicted genes, expressed sequence tags (ESTs), guanine/cytosine (GC) content, repeats and comparative genomics, as well as unique data resources including (i) 2473 genetic markers and experimentally confirmed PCR primers that amplify these markers, (ii) 142 414 bacterial artificial chromosome (BAC) and 217 344 fosmid end sequences that amount to 15.0- and 11.1-fold clone coverage of the entire genome, respectively, and were used for draft genome assembly, (iii) 16 519 460 single nucleotide polymorphisms (SNPs), and 2 859 905 insertions/deletions detected between two medaka inbred strain genomes and (iv) 841 235 5′-end serial analyses of gene-expression (SAGE) tags that identified 344 266 transcription start sites on the genome. UTGB/medaka is available at: http://medaka.utgenome.org/

Genome of the pitcher plant Cephalotus reveals genetic changes associated with carnivory

Carnivorous plants exploit animals as a nutritional source and have inspired long-standing questions about the origin and evolution of carnivory-related traits. To investigate the molecular bases of carnivory, we sequenced the genome of the heterophyllous pitcher plant Cephalotus follicularis, in which we succeeded in regulating the developmental switch between carnivorous and non-carnivorous leaves. Transcriptome comparison of the two leaf types and gene repertoire analysis identified genetic changes associated with prey attraction, capture, digestion and nutrient absorption. Analysis of digestive fluid proteins from C. follicularis and three other carnivorous plants with independent carnivorous origins revealed repeated co-options of stress-responsive protein lineages coupled with convergent amino acid substitutions to acquire digestive physiology. These results imply constraints on the available routes to evolve plant carnivory.

FMN binding and photochemical properties of plant putative photoreceptors containing two LOV domains, LOV/LOV proteins.

LOV domains function as blue light-sensing modules in various photoreceptors in plants, fungi, algae, and bacteria. A LOV/LOV protein (LLP) has been found from Arabidopsis thaliana (AtLLP) as a two LOV domain-containing protein. However, its function remains unknown. We isolated cDNA clones coding for an LLP homolog from tomato (Solanum lycopersicum) and two homologs from the moss Physcomitrella patens. The tomato LLP (SlLLP) contains two LOV domains (LOV1 and LOV2 domains), as in AtLLP. Most of the amino acids required for association with chromophore are conserved in both LOV domains, except that the amino acid at the position equivalent to the cysteine essential for cysteinyl adduct formation is glycine in the LOV1 domain as in AtLLP. When expressed in Escherichia coli, SlLLP binds FMN and undergoes a self-contained photocycle upon irradiation of blue light. Analyses using mutant SlLLPs revealed that SlLLP binds FMN in both LOV domains, although the LOV1 domain does not show spectral changes on irradiation. However, when Gly66 in the LOV1 domain, which is located at the position equivalent to the essential cysteine of LOV domains, is replaced by cysteine, the mutated LOV1 domain shows light-induced spectral changes. In addition, all four LOV domains of P. patens LLPs (PpLLP1 and PpLLP2) show the typical features of LOV domains, including the reactive cysteine in each. This study shows that plants have a new LOV domain-containing protein family with the typical biochemical and photochemical properties of other LOV domain-containing proteins such as the phototropins.

Cost-Effective Sequencing of Full-Length cDNA Clones Powered by a De Novo -Reference Hybrid Assembly

Background Sequencing full-length cDNA clones is important to determine gene structures including alternative splice forms, and provides valuable resources for experimental analyses to reveal the biological functions of coded proteins. However, previous approaches for sequencing cDNA clones were expensive or time-consuming, and therefore, a fast and efficient sequencing approach was demanded. Methodology We developed a program, MuSICA 2, that assembles millions of short (36-nucleotide) reads collected from a single flow cell lane of Illumina Genome Analyzer to shotgun-sequence ∼800 human full-length cDNA clones. MuSICA 2 performs a hybrid assembly in which an external de novo assembler is run first and the result is then improved by reference alignment of shotgun reads. We compared the MuSICA 2 assembly with 200 pooled full-length cDNA clones finished independently by the conventional primer-walking using Sanger sequencers. The exon-intron structure of the coding sequence was correct for more than 95% of the clones with coding sequence annotation when we excluded cDNA clones insufficiently represented in the shotgun library due to PCR failure (42 out of 200 clones excluded), and the nucleotide-level accuracy of coding sequences of those correct clones was over 99.99%. We also applied MuSICA 2 to full-length cDNA clones from Toxoplasma gondii, to confirm that its ability was competent even for non-human species. Conclusions The entire sequencing and shotgun assembly takes less than 1 week and the consumables cost only ∼US

An adenylyl cyclase with a phosphodiesterase domain in basal plants with a motile sperm system

3 per clone, demonstrating a significant advantage over previous approaches.

Introducing difference recurrence relations for faster semi-global alignment of long sequences

Adenylyl cyclase (AC), which produces the signalling molecule cAMP, has numerous important cellular functions in diverse organisms from prokaryotes to eukaryotes. Here we report the identification and characterization of an AC gene from the liverwort Marchantia polymorpha. The encoded protein has both a C-terminal AC catalytic domain similar to those of class III ACs and an N-terminal cyclic nucleotide phosphodiesterase (PDE) domain that degrades cyclic nucleotides, thus we designated the gene MpCAPE (COMBINED AC with PDE). Biochemical analyses of recombinant proteins showed that MpCAPE has both AC and PDE activities. In MpCAPE-promoter-GUS lines, GUS activity was specifically detected in the male sexual organ, the antheridium, suggesting MpCAPE and thus cAMP signalling may be involved in the male reproductive process. CAPE orthologues are distributed only in basal land plants and charophytes that use motile sperm as the male gamete. CAPE is a subclass of class III AC and may be important in male organ and cell development in basal plants.