Budrul Ahsan

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.

Mutations in COQ2 in familial and sporadic multiple-system atrophy the multiple-system atrophy research collaboration

BACKGROUND Multiple-system atrophy is an intractable neurodegenerative disease characterized by autonomic failure in addition to various combinations of parkinsonism, cerebellar ataxia, and pyramidal dysfunction. Although multiple-system atrophy is widely considered to be a nongenetic disorder, we previously identified multiplex families with this disease, which indicates the involvement of genetic components. METHODS In combination with linkage analysis, we performed whole-genome sequencing of a sample obtained from a member of a multiplex family in whom multiple-system atrophy had been diagnosed on autopsy. We also performed mutational analysis of samples from members of five other multiplex families and from a Japanese series (363 patients and two sets of controls, one of 520 persons and one of 2383 persons), a European series (223 patients and 315 controls), and a North American series (172 patients and 294 controls). On the basis of these analyses, we used a yeast complementation assay and measured enzyme activity of parahydroxybenzoate-polyprenyl transferase. This enzyme is encoded by the gene COQ2 and is essential for the biosynthesis of coenzyme Q10. Levels of coenzyme Q10 in lymphoblastoid cells and brain tissue were measured on high-performance liquid chromatography. RESULTS We identified a homozygous mutation (M78V-V343A/M78V-V343A) and compound heterozygous mutations (R337X/V343A) in COQ2 in two multiplex families. Furthermore, we found that a common variant (V343A) and multiple rare variants in COQ2, all of which are functionally impaired, are associated with sporadic multiple-system atrophy. The V343A variant was exclusively observed in the Japanese population. CONCLUSIONS Functionally impaired variants of COQ2 were associated with an increased risk of multiple-system atrophy in multiplex families and patients with sporadic disease, providing evidence of a role of impaired COQ2 activities in the pathogenesis of this disease. (Funded by the Japan Society for the Promotion of Science and others.).

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.

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.

High-Resolution Analysis of the 5′-End Transcriptome Using a Next Generation DNA Sequencer

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.

Variants associated with Gaucher disease in multiple system atrophy

Glucocerebrosidase gene (GBA) variants that cause Gaucher disease are associated with Parkinson disease (PD) and dementia with Lewy bodies (DLB). To investigate the role of GBA variants in multiple system atrophy (MSA), we analyzed GBA variants in a large case–control series.

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/

CSF1R mutations identified in three families with autosomal dominantly inherited leukoencephalopathy

Genetic and phenotypic heterogeneities are considerably high in adult‐onset leukoencephalopathy, in which comprehensive mutational analyses of the candidate genes by conventional methods are too laborious. We applied exome sequencing to conduct a comprehensive mutational analysis of genes for autosomal dominant leukoencephalopathies. Genomic DNA samples from four patients of three families with autosomal dominantly inherited adult‐onset leukodystrophy were subjected to exome sequencing. On the basis of the results, 21 patients with adult‐onset sporadic leukodystrophy and one patient with pathologically proven HDLS were additionally screened for CSF1R mutations. Exome sequencing identified heterozygous CSF1R mutations (p.I794T and p.R777W) in two families. I794T has recently been reported as a causative mutation for hereditary diffuse leukoencephalopathy with spheroids (HDLS), and R777W is a novel mutation. Although mutational analysis of CSF1R in 21 sporadic cases revealed no mutations, another novel CSF1R mutation, p.C653Y, was identified in one patient with autopsy‐proven HDSL. These variants were located in the PTK domain where the causative mutations cluster. Functional prediction of the mutant CSF1R as well as cross‐species conservation of the affected amino acids supports the notion that these variants are pathogenic for HDLS. Exome sequencing is useful for a comprehensive mutational analysis of causative genes for hereditary leukoencephalopathies, and CSF1R should be considered a candidate gene for patients with autosomal dominant leukoencephalopathies.

MachiBase: a Drosophila melanogaster 5′-end mRNA transcription database

MachiBase (http://machibase.gi.k.u-tokyo.ac.jp/) provides a comprehensive and freely accessible resource regarding Drosophila melanogaster 5′-end mRNA transcription at different developmental states, supporting studies on the variabilities of promoter transcriptional activities and gene-expression profiles in the fruitfly. The data were generated in conjunction with the recently developed high-throughput genome sequencer Illumina/Solexa using a newly developed 5′-end mRNA collection method.

Posterior column ataxia with retinitis pigmentosa in a Japanese family with a novel mutation in FLVCR1

Posterior column ataxia with retinitis pigmentosa (PCARP) is an autosomal recessive neurodegenerative disorder characterized by retinitis pigmentosa and sensory ataxia. Previous studies of PCARP in two families showed a linkage to 1q31–q32. However, detailed investigations on the clinical presentations as well as molecular genetics of PCARP have been limited. Here, we describe a Japanese consanguineous family with PCARP. Two affected siblings suffered from childhood-onset retinitis pigmentosa and slowly progressive sensory ataxia. They also showed mild mental retardation, which has not been described in patients with PCARP. Parametric linkage analysis using high-density single nucleotide polymorphism arrays supported a linkage to the same locus. Target capture and high-throughput sequencing technologies revealed a novel homozygous c.1477G>C (G493R) mutation in FLVCR1, which cosegregated with the disease. A recent study has identified three independent mutations in FLVCR1 in the original and other families. Our results further confirmed that PCARP is caused by mutations in FLVCR1.