Asao Fujiyama

The amphioxus genome and the evolution of the chordate karyotype.

Lancelets (‘amphioxus’) are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic ∼520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.

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.

A map of rice genome variation reveals the origin of cultivated rice

Crop domestications are long-term selection experiments that have greatly advanced human civilization. The domestication of cultivated rice (Oryza sativa L.) ranks as one of the most important developments in history. However, its origins and domestication processes are controversial and have long been debated. Here we generate genome sequences from 446 geographically diverse accessions of the wild rice species Oryza rufipogon, the immediate ancestral progenitor of cultivated rice, and from 1,083 cultivated indica and japonica varieties to construct a comprehensive map of rice genome variation. In the search for signatures of selection, we identify 55 selective sweeps that have occurred during domestication. In-depth analyses of the domestication sweeps and genome-wide patterns reveal that Oryza sativa japonica rice was first domesticated from a specific population of O. rufipogon around the middle area of the Pearl River in southern China, and that Oryza sativa indica rice was subsequently developed from crosses between japonica rice and local wild rice as the initial cultivars spread into South East and South Asia. The domestication-associated traits are analysed through high-resolution genetic mapping. This study provides an important resource for rice breeding and an effective genomics approach for crop domestication research.

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.

Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads

Although many de novo genome assembly projects have recently been conducted using high-throughput sequencers, assembling highly heterozygous diploid genomes is a substantial challenge due to the increased complexity of the de Bruijn graph structure predominantly used. To address the increasing demand for sequencing of nonmodel and/or wild-type samples, in most cases inbred lines or fosmid-based hierarchical sequencing methods are used to overcome such problems. However, these methods are costly and time consuming, forfeiting the advantages of massive parallel sequencing. Here, we describe a novel de novo assembler, Platanus, that can effectively manage high-throughput data from heterozygous samples. Platanus assembles DNA fragments (reads) into contigs by constructing de Bruijn graphs with automatically optimized k-mer sizes followed by the scaffolding of contigs based on paired-end information. The complicated graph structures that result from the heterozygosity are simplified during not only the contig assembly step but also the scaffolding step. We evaluated the assembly results on eukaryotic samples with various levels of heterozygosity. Compared with other assemblers, Platanus yields assembly results that have a larger scaffold NG50 length without any accompanying loss of accuracy in both simulated and real data. In addition, Platanus recorded the largest scaffold NG50 values for two of the three low-heterozygosity species used in the de novo assembly contest, Assemblathon 2. Platanus therefore provides a novel and efficient approach for the assembly of gigabase-sized highly heterozygous genomes and is an attractive alternative to the existing assemblers designed for genomes of lower heterozygosity.

Role for piRNAs and Noncoding RNA in de Novo DNA Methylation of the Imprinted Mouse Rasgrf1 Locus

Small, noncoding PIWI-interacting RNAs regulate the imprinting of a mouse gene. Genomic imprinting causes parental origin–specific monoallelic gene expression through differential DNA methylation established in the parental germ line. However, the mechanisms underlying how specific sequences are selectively methylated are not fully understood. We have found that the components of the PIWI-interacting RNA (piRNA) pathway are required for de novo methylation of the differentially methylated region (DMR) of the imprinted mouse Rasgrf1 locus, but not other paternally imprinted loci. A retrotransposon sequence within a noncoding RNA spanning the DMR was targeted by piRNAs generated from a different locus. A direct repeat in the DMR, which is required for the methylation and imprinting of Rasgrf1, served as a promoter for this RNA. We propose a model in which piRNAs and a target RNA direct the sequence-specific methylation of Rasgrf1.

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.

Cloning and structural analyses of hepatitis B virus DNAs, subtype adr

Entire genomes of hepatitis B virus (subtype adr) have been cloned. The nucleotide sequence data were compared with other sequences of HBV genome including: adw [Valenzuela et al. (1981) in Animal Virus Genetics. Fields et al. eds. Academic Press, Inc., NY. pp. 57-70], ayw [Galibert et al. (1979) Nature, 281, 646-650], and adyw [Pasek et al. (1979) Nature 282, 575-579]. Four open coding frames for polypeptides larger than 6,000 dalton were found to be conserved and were highly compressed by overlapping with each other in one strand (L-strand). Sites of initiation of the S gene and termination of the P gene were not conserved. No conserved coding frame was found on the opposite strand (S strand). Amino acid sequences of six surface antigen (HBsAg) peptides, including subtypes adr, adw, and ayw, are deduced from the DNA sequences, and the substitution of amino acid residues which are consistent with the change of subtypes are demonstrated.

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

DNA sequence and comparative analysis of chimpanzee chromosome 22

Human–chimpanzee comparative genome research is essential for narrowing down genetic changes involved in the acquisition of unique human features, such as highly developed cognitive functions, bipedalism or the use of complex language. Here, we report the high-quality DNA sequence of 33.3 megabases of chimpanzee chromosome 22. By comparing the whole sequence with the human counterpart, chromosome 21, we found that 1.44% of the chromosome consists of single-base substitutions in addition to nearly 68,000 insertions or deletions. These differences are sufficient to generate changes in most of the proteins. Indeed, 83% of the 231 coding sequences, including functionally important genes, show differences at the amino acid sequence level. Furthermore, we demonstrate different expansion of particular subfamilies of retrotransposons between the lineages, suggesting different impacts of retrotranspositions on human and chimpanzee evolution. The genomic changes after speciation and their biological consequences seem more complex than originally hypothesized.Human–chimpanzee comparative genome research is essential for narrowing down genetic changes involved in the acquisition of unique human features, such as highly developed cognitive functions, bipedalism or the use of complex language. Here, we report the high-quality DNA sequence of 33.3 megabases of chimpanzee chromosome 22. By comparing the whole sequence with the human counterpart, chromosome 21, we found that 1.44% of the chromosome consists of single-base substitutions in addition to nearly 68,000 insertions or deletions. These differences are sufficient to generate changes in most of the proteins. Indeed, 83% of the 231 coding sequences, including functionally important genes, show differences at the amino acid sequence level. Furthermore, we demonstrate different expansion of particular subfamilies of retrotransposons between the lineages, suggesting different impacts of retrotranspositions on human and chimpanzee evolution. The genomic changes after speciation and their biological consequences seem more complex than originally hypothesized.