Medakafish as a model system for vertebrate developmental genetics

Several teleosts, such as the zebrafish and the medakafish or medaka (Oryzias latipes), are used as vertebrate model systems in various fields of biology. The medaka is suitable for use in genomic studies because of its small genome size. Moreover, our recent results of small-scale mutagenesis in the medaka indicate that it is possible to identify mutations, the phenotypes of which could not be found in zebrafish mutants obtained by large-scale mutagenesis. An example is Oot (One-sided optic tectum), a maternal-effect mutation. In the Oot phenotype, bilateral symmetry is broken in the optic tectum in the early developmental stages, and either the left or right morphology is duplicated on both sides. Medaka inbred strains can be produced and used to study quantitative traits in vertebrate development. Data presented support the use of medaka as another important fish model for the study of vertebrate developmental genetics.     

The rice field eel as a model system for vertebrate sexual development

Complex developmental mechanisms of vertebrates are unraveled using comparative genomic approaches. Several teleosts, such as zebrafish, medaka and pufferfish, are used as genetic model systems because they are amenable to studies of gene function. The rice field eel, a freshwater fish, is emerging as a specific model system for studies of vertebrate sexual development because of its small genome size and naturally occurring sex reversal. Data presented here support the use of the rice field eel as another important fish model for comparative genome studies, especially in vertebrate sexual development. This model system is complementary rather than redundant.     

Teleost fish with specific genome duplication as unique models of vertebrate evolution

Whole-genome duplication (WGD) is believed to be one of the major evolutionary events that shaped the genome organization of vertebrates. Here, we review recent research on vertebrate genome evolution, specifically on WGD and its consequences for gene and genome evolution in teleost fishes. Recent genome analyses confirmed that all vertebrates experienced two rounds of WGD early in their evolution, and that teleosts experienced a subsequent additional third-round (3R)-WGD. The 3R-WGD was estimated to have occurred 320–400 million years ago in a teleost ancestor, but after its divergence from a common ancestor with living non-teleost actinopterygians (Bichir, Sturgeon, Bowfin, and Gar) based on the analyses of teleost-specific duplicate genes. This 3R-WGD was confirmed by synteny analysis and ancestral karyotype inference using the genome sequences of Tetraodon and medaka. Most of the tetrapods, on the other hand, have not experienced an additional WGD; however, they have experienced repeated chromosomal rearrangements throughout the whole genome. Therefore, different types of chromosomal events have characterized the genomes of teleosts and tetrapods, respectively. The 3R-WGD is useful to investigate the consequences of WGD because it is an evolutionarily recent WGD and thus teleost genomes retain many more WGD-derived duplicates and “traces” of their evolution. In addition, the remarkable morphological, physiological, and ecological diversity of teleosts may facilitate understanding of macrophenotypic evolution on the basis of genetic/genomic information. We highlight the teleosts with 3R-WGD as unique models for future studies on ecology and evolution taking advantage of emerging genomics technologies and systems biology environments.     

Temporal pattern of loss/persistence of duplicate genes involved in signal transduction and metabolic pathways after teleost-specific genome duplication

Background  

Recent genomic studies have revealed a teleost-specific third-round whole genome duplication (3R-WGD) event occurred in a common ancestor of teleost fishes. However, it is unclear how the genes duplicated in this event were lost or persisted during the diversification of teleosts, and therefore, how many of the duplicated genes contribute to the genetic differences among teleosts. This subject is also important for understanding the process of vertebrate evolution through WGD events. We applied a comparative evolutionary approach to this question by focusing on the genes involved in long-term potentiation, taste and olfactory transduction, and the tricarboxylic acid cycle, based on the whole genome sequences of four teleosts; zebrafish, medaka, stickleback, and green spotted puffer fish.     

Functional lysophosphatidic acid receptors expressed in Oryzias latipes

Lysophosphatidic acid (LPA) signaling is known to play biological and pathophysiological roles in many types of animals. Medaka (Oryzias latipes) is an experimental fish that can be easily maintained, propagated, and analyzed, and whose genome has been completely sequenced. However, there is limited information available regarding medaka LPA receptors. Here, using information from the medaka genome database, we examine the genomic structures, expression, and functions of six LPA receptor genes, Lpar1–Lpar6. Our analyses reveal that the genomic structures of Lpar1 and Lpar4 are different from those deduced from the database. Functional analyses using a heterologous expression system demonstrate that all medaka LPA receptors except for LPA_5b respond to LPA treatment with cytoskeletal changes. These findings provide useful information on the structure and function of medaka LPA receptor genes, and identify medaka as a useful experimental model for exploration of the biological significance of LPA signaling.     

Liver receptor homologue-1 (LRH-1) activates the promoter of brain aromatase (cyp19a2) in a teleost fish, the medaka, Oryzias latipes

The medaka, Oryzias latipes, like other fish, have two distinct aromatase genes, the ovarian (cyp19a1) and brain (cyp19a2) forms. We previously reported that Ad4BP/SF-1, a member of the NR5A subfamily, plays an important role in the regulation of cyp19a1 expression in medaka ovarian follicles during vitellogenesis. In the present study, we investigated whether liver receptor homologue-1 (LRH-1), another NR5A subfamily member, is involved in the regulation of cyp19a2 expression in the medaka brain. In situ hybridization analysis revealed that LRH-1 was expressed in the hypothalamus, where it colocalized with aromatase (cyp19a2). We then showed by transient transfection assays that LRH-1 was able to increase expression of a cyp19a2 reporter gene in various mammalian cell lines, and that mutation of a putative LRH-1 binding site within the cyp19a2 promoter abolished this effect. Taken together, these findings suggest that LRH-1 plays a role in regulating cyp19a2 expression in the medaka brain. This is the first to demonstrate in vitro the activation of brain aromatase by LRH-1 in the vertebrate brain.     

Consequences of Lineage-Specific Gene Loss on Functional Evolution of Surviving Paralogs: ALDH1A and Retinoic Acid Signaling in Vertebrate Genomes

Genome duplications increase genetic diversity and may facilitate the evolution of gene subfunctions. Little attention, however, has focused on the evolutionary impact of lineage-specific gene loss. Here, we show that identifying lineage-specific gene loss after genome duplication is important for understanding the evolution of gene subfunctions in surviving paralogs and for improving functional connectivity among human and model organism genomes. We examine the general principles of gene loss following duplication, coupled with expression analysis of the retinaldehyde dehydrogenase Aldh1a gene family during retinoic acid signaling in eye development as a case study. Humans have three ALDH1A genes, but teleosts have just one or two. We used comparative genomics and conserved syntenies to identify loss of ohnologs (paralogs derived from genome duplication) and to clarify uncertain phylogenies. Analysis showed that Aldh1a1 and Aldh1a2 form a clade that is sister to Aldh1a3-related genes. Genome comparisons showed secondarily loss of aldh1a1 in teleosts, revealing that Aldh1a1 is not a tetrapod innovation and that aldh1a3 was recently lost in medaka, making it the first known vertebrate with a single aldh1a gene. Interestingly, results revealed asymmetric distribution of surviving ohnologs between co-orthologous teleost chromosome segments, suggesting that local genome architecture can influence ohnolog survival. We propose a model that reconstructs the chromosomal history of the Aldh1a family in the ancestral vertebrate genome, coupled with the evolution of gene functions in surviving Aldh1a ohnologs after R1, R2, and R3 genome duplications. Results provide evidence for early subfunctionalization and late subfunction-partitioning and suggest a mechanistic model based on altered regulation leading to heterochronic gene expression to explain the acquisition or modification of subfunctions by surviving ohnologs that preserve unaltered ancestral developmental programs in the face of gene loss.     

Immunoglobulin light chains in medaka (Oryzias latipes)

The gene segments encoding antibodies have been studied in many capacities and represent some of the best-characterized gene families in traditional animal disease models (mice and humans). To date, multiple immunoglobulin light chain (IgL) isotypes have been found in vertebrates and it is unclear as to which isotypes might be more primordial in nature. Sequence data emerging from an array of fish genome projects is a valuable resource for discerning complex multigene assemblages in this critical branch point of vertebrate phylogeny. Herein, we have analyzed the genomic organization of medaka (Oryzias latipes) IgL gene segments based on recently released genome data. The medaka IgL locus located on chromosome 11 contains at least three clusters of IgL gene segments comprised of multiple gene assemblages of the kappa light chain isotype. These data suggest that medaka IgL gene segments may undergo both intra- and inter-cluster rearrangements as a means to generate additional diversity. Alignments of expressed sequence tags to concordant gene segments which revealed each of the three IgL clusters are expressed. Collectively, these data provide a genomic framework for IgL genes in medaka and indicate that Ig diversity in this species is achieved from at least three distinct chromosomal regions.     

Expression of Aromatase in Radial Glial Cells in the Brain of the Japanese Eel Provides Insight into the Evolution of the cyp191a Gene in Actinopterygians

The cyp19a1 gene that encodes aromatase, the only enzyme permitting conversion of C19 aromatizable androgens into estrogens, is present as a single copy in the genome of most vertebrate species, except in teleosts in which it has been duplicated. This study aimed at investigating the brain expression of a cyp19a1 gene expressed in both gonad and brain of Japanese eel, a basal teleost. By means of immunohistochemistry and in situ hybridization, we show that cyp19a1 is expressed only in radial glial cells of the brain and in pituitary cells. Treatments with salmon pituitary homogenates (female) or human chorionic gonadotrophin (male), known to turn on steroid production in immature eels, strongly stimulated cyp19a1 messenger and protein expression in radial glial cells and pituitary cells. Using double staining studies, we also showed that aromatase-expressing radial glial cells exhibit proliferative activity in both the brain and the pituitary. Altogether, these data indicate that brain and pituitary expression of Japanese eel cyp19a1 exhibits characteristics similar to those reported for the brain specific cyp19a1b gene in teleosts having duplicated cyp19a1 genes. This supports the hypothesis that, despite the fact that eels also underwent the teleost specific genome duplication, they have a single cyp19a1 expressed in both brain and gonad. Such data also suggest that the intriguing features of brain aromatase expression in teleost fishes were not gained after the whole genome duplication and may reflect properties of the cyp19a1 gene of ancestral Actinopterygians.     

Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis.

Myosin heavy chain genes (MYHs) are the most important functional domains of myosins, which are highly conserved throughout evolution. The human genome contains 15 MYHs, whereas the corresponding number in teleost appears to be much higher. Although teleosts comprise more than one-half of all vertebrate species, our knowledge of MYHs in teleosts is rather limited. A comprehensive analysis of the torafugu (Takifugu rubripes) genome database enabled us to detect at least 28 MYHs, almost twice as many as in humans. RT-PCR revealed that at least 16 torafugu MYH representatives (5 fast skeletal, 3 cardiac, 2 slow skeletal, 1 superfast, 2 smooth, and 3 nonmuscle types) are actually transcribed. Among these, MYH(M743-2) and MYH(M5) of fast and slow skeletal types, respectively, are expressed during development of torafugu embryos. Syntenic analysis reveals that torafugu fast skeletal MYHs are distributed across five genomic regions, three of which form clusters. Interestingly, while human fast skeletal MYHs form one cluster, its syntenic region in torafugu is duplicated, although each locus contains just a single MYH in torafugu. The results of the syntenic analysis were further confirmed by corresponding analysis of MYHs based on databases from Tetraodon, zebrafish, and medaka genomes. Phylogenetic analysis suggests that fast skeletal MYHs evolved independently in teleosts and tetrapods after fast skeletal MYHs had diverged from four ancestral MYHs.     

The evolution of teleost pigmentation and the fish‐specific genome duplication

Teleost fishes have evolved a unique complexity and diversity of pigmentation and colour patterning that is unmatched among vertebrates. Teleost colouration is mediated by five different major types of neural‐crest derived pigment cells, while tetrapods have a smaller repertoire of such chromatophores. The genetic basis of teleost colouration has been mainly uncovered by the cloning of pigmentation genes in mutants of zebrafish Danio rerio and medaka Oryzias latipes. Many of these teleost pigmentation genes were already known as key players in mammalian pigmentation, suggesting partial conservation of the corresponding developmental programme among vertebrates. Strikingly, teleost fishes have additional copies of many pigmentation genes compared with tetrapods, mainly as a result of a whole‐genome duplication that occurred 320–350 million years ago at the base of the teleost lineage, the so‐called fish‐specific genome duplication. Furthermore, teleosts have retained several duplicated pigmentation genes from earlier rounds of genome duplication in the vertebrate lineage, which were lost in other vertebrate groups. It was hypothesized that divergent evolution of such duplicated genes may have played an important role in pigmentation diversity and complexity in teleost fishes, which therefore not only provide important insights into the evolution of the vertebrate pigmentary system but also allow us to study the significance of genome duplications for vertebrate biodiversity.     

The DNA sequence of medaka chromosome LG22

We report the genomic DNA sequence of a single chromosome (linkage group 22; LG22) of the small teleost fish medaka (Oryzias latipes) as a first whole chromosome sequence from a non-mammalian vertebrate. The order and orientation of 633 protein-coding genes were deduced from 18,803,338 bp of DNA sequence, providing the opportunity to analyze chromosome evolution of vertebrate genomes by direct comparison with the human genome. The average number of genes in the "conserved gene cluster" (CGC), a strict definition of "synteny" at the sequence basis, between medaka and human was 1.6. These and other data suggest that approximately 38.8% of pair-wise gene relationships would have been broken from their common ancestor in the human and medaka lineages and further imply that approx 20,000 (15,520-23,280) breaks would have occurred from the entire genome of the common ancestor. These breaks were generated mainly by intra-chromosomal shufflings at a specific era in the vertebrate lineage. These precise comparative genomics allowed us to identify the pieces of ancient chromosomes of the common vertebrate ancestor and estimate chromosomal evolution in the vertebrate lineage.     

Different evolutionary fates of recently integrated human and chimpanzee LINE-1 retrotransposons

The long interspersed element-1 (LINE-1 or L1) is a highly successful retrotransposon in mammals. L1 elements have continued to actively propagate subsequent to the human–chimpanzee divergence,  6 million years ago, resulting in species-specific inserts. Here, we report a detailed characterization of chimpanzee-specific L1 subfamily diversity and a comparison with their human-specific counterparts. Our results indicate that L1 elements have experienced different evolutionary fates in humans and chimpanzees within the past  6 million years. Although the species-specific L1 copy numbers are on the same order in both species (1200–2000 copies), the number of retrotransposition-competent elements appears to be much higher in the human genome than in the chimpanzee genome. Also, while human L1 subfamilies belong to the same lineage, we identified two lineages of recently integrated L1 subfamilies in the chimpanzee genome. The two lineages seem to have coexisted for several million years, but only one shows evidence of expansion within the past three million years. These differential evolutionary paths may be the result of random variation, or the product of competition between L1 subfamily lineages. Our results suggest that the coexistence of several L1 subfamily lineages within a species may be resolved in a very short evolutionary period of time, perhaps in just a few million years. Therefore, the chimpanzee genome constitutes an excellent model in which to analyze the evolutionary dynamics of L1 retrotransposons.     

Cytochrome P450 1 genes in early deuterostomes (tunicates and sea urchins) and vertebrates (chicken and frog): origin and diversification of the CYP1 gene family

Cytochrome P450 family 1 (CYP1) proteins are important in a large number of toxicological processes. CYP1A and CYP1B genes are well known in mammals, but the evolutionary history of the CYP1 family as a whole is obscure; that history may provide insight into endogenous functions of CYP1 enzymes. Here, we identify CYP1-like genes in early deuterostomes (tunicates and echinoderms), and several new CYP1 genes in vertebrates (chicken, Gallus gallus and frog, Xenopus tropicalis). Profile hidden Markov models (HMMs) generated from vertebrate CYP1A and CYP1B protein sequences were used to identify 5 potential CYP1 homologs in the tunicate Ciona intestinalis genome. The C. intestinalis genes were cloned and sequenced, confirming the predicted sequences. Orthologs of 4 of these genes were found in the Ciona savignyi genome. Bayesian phylogenetic analyses group the tunicate genes in the CYP1 family, provisionally in 2 new subfamilies, CYP1E and CYP1F, which fall in the CYP1A and CYP1B/1C clades. Bayesian and maximum likelihood analyses predict functional divergence between the tunicate and vertebrate CYP1s, and regions within CYP substrate recognition sites were found to differ significantly in position-specific substitution rates between tunicates and vertebrates. Subsequently, 10 CYP1-like genes were found in the echinoderm Strongylocentrotus purpuratus (sea urchin) genome. Several of the tunicate and echinoderm CYP1-like genes are expressed during development. Canonical xenobiotic response elements are present in the upstream genomic sequences of most tunicate and sea urchin CYP1s, and both groups are predicted to possess an aryl hydrocarbon receptor (AHR), suggesting possible regulatory linkage of AHR and these CYPs. The CYP1 family has undergone multiple rounds of gene duplication followed by functional divergence, with at least one gene lost in mammals. This study provides new insight into the origin and evolution of CYP1 genes.     

53BP1‐mediated recruitment of RASSF1A to ribosomal DNA breaks promotes local ATM signaling

DNA lesions occur across the genome and constitute a threat to cell viability; however, damage at specific genomic loci has a relatively greater impact on overall genome stability. The ribosomal RNA gene repeats (rDNA) are emerging fragile sites. Recent progress in understanding how the rDNA damage response is organized has highlighted a key role of adaptor proteins. Here, we show that the scaffold tumor suppressor RASSF1A is recruited to rDNA breaks. RASSF1A recruitment to double‐strand breaks is mediated by 53BP1 and depends on RASSF1A phosphorylation at Serine 131 by ATM kinase. Employing targeted rDNA damage, we uncover that RASSF1A recruitment promotes local ATM signaling. RASSF1A silencing, a common epigenetic event during malignant transformation, results in persistent breaks, rDNA copy number alterations and decreased cell viability. Overall, we identify a novel role for RASSF1A at rDNA break sites, provide mechanistic insight into how the DNA damage response is organized in a chromatin context, and provide further evidence for how silencing of the RASSF1A tumor suppressor contributes to genome instability.     

Five Nkx5 genes show differential expression patterns in anlagen of sensory organs in medaka: insight into the evolution of the gene family

We report the identification and characterisation of five different Nkx5-related genes in medaka fish (Oryzias latipes). They constitute homologues of genes previously isolated in higher vertebrates, Nkx5--1, Nkx5--2, Hmx1/Nkx5--3 and SOHo-1, and were named accordingly: OlNkx5--1.1, OlNkx5--2, OlNkx5--3 and OlSOHo. For the Nkx5--1 gene a new, second homologue, OlNkx5--1.2, was isolated. In medaka, Nkx5 and SOHo genes are differentially expressed in three developing sensory organs: eye, ear and lateral line and later in defined brain regions. Phylogenetic analyses of the entire Nkx5 family revealed that four paralogous Nkx5 groups, Nkx5--1, Nkx5--2, Hmx1/Nkx5--3/GH6 and SOHo, are present in vertebrates. Only some of the Nkx5 family members have been identified in singular vertebrate species so far. Here we present, for the first time, the isolation of representatives of each Nkx5 subgroup in one species, the medaka fish. Based on similarities in sequence and expression patterns, and genomic organisation we propose a model of the evolutionary history of the Nkx5 family. The model predicts that the four vertebrate Nkx5 genes arose by a tandem duplication, followed by chromosomal duplication. The two Nkx5--1 genes identified so far exclusively in medaka most probably result from an additional genome duplication in the fish lineage.