Rapid Diversification of FoxP2 in Teleosts through Gene Duplication in the Teleost-Specific Whole Genome Duplication Event

As one of the most conserved genes in vertebrates, FoxP2 is widely involved in a number of important physiological and developmental processes. We systematically studied the evolutionary history and functional adaptations of FoxP2 in teleosts. The duplicated FoxP2 genes (FoxP2a and FoxP2b), which were identified in teleosts using synteny and paralogon analysis on genome databases of eight organisms, were probably generated in the teleost-specific whole genome duplication event. A credible classification with FoxP2, FoxP2a and FoxP2b in phylogenetic reconstructions confirmed the teleost-specific FoxP2 duplication. The unavailability of FoxP2b in Danio rerio suggests that the gene was deleted through nonfunctionalization of the redundant copy after the Otocephala-Euteleostei split. Heterogeneity in evolutionary rates among clusters consisting of FoxP2 in Sarcopterygii (Cluster 1), FoxP2a in Teleostei (Cluster 2) and FoxP2b in Teleostei (Cluster 3), particularly between Clusters 2 and 3, reveals asymmetric functional divergence after the gene duplication. Hierarchical cluster analyses of hydrophobicity profiles demonstrated significant structural divergence among the three clusters with verification of subsequent stepwise discriminant analysis, in which FoxP2 of Leucoraja erinacea and Lepisosteus oculatus were classified into Cluster 1, whereas FoxP2b of Salmo salar was grouped into Cluster 2 rather than Cluster 3. The simulated thermodynamic stability variations of the forkhead box domain (monomer and homodimer) showed remarkable divergence in FoxP2, FoxP2a and FoxP2b clusters. Relaxed purifying selection and positive Darwinian selection probably were complementary driving forces for the accelerated evolution of FoxP2 in ray-finned fishes, especially for the adaptive evolution of FoxP2a and FoxP2b in teleosts subsequent to the teleost-specific gene duplication.     

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.     

Characterisation and differential expression during development of a duplicate Disabled-1 (Dab1) gene from zebrafish

We identified a new duplicated Dab1 gene (drDab1b) spanning around 25 kb of genomic DNA in zebrafish. Located in zebrafish chromosome 2, it is composed of 11 encoding exons and shows high sequence similarity to other Dab1 genes, including drDab1a, a zebrafish Dab1 gene previously characterised. drDab1b encodes by alternative splicing at least five different isoforms. Both drDab1a and drDab1b show differential gene expression levels in distinct adult tissues and during development. drDab1b is expressed in peripheral tissues (gills, heart, intestine, muscle), the immune system (blood, liver) and the central nervous system (CNS), whereas drDab1a is only expressed in gills, muscle and the CNS, suggesting a division of functions for two Dab1 genes in zebrafish adult tissues. RT-PCR analysis also reveals that both drDab1 genes show distinct developmental-specific expression patterns throughout development. drDab1b expression was higher than that of drDab1a, suggesting a major role of drDab1b in comparison with drDab1a during development and in different adult tissues. In addition, new putative Dab1 (a and/or b) from different teleost species were identified in silico and predicted protein products are compared with the previously characterised Dab1, demonstrating that the Dab1b group is more ancestral than their paralogue, the Dab1a group.     

Functional differentiation of bmp2a and bmp2b genes in zebrafish

Bone morphogenetic protein 2 plays an important role in the regulation of osteoblast proliferation and differentiation. Phylogenetic analysis showed that the bmp2 ortholog evolved from the same ancestral gene family in vertebrates and was duplicated in teleost, which were named bmp2a and bmp2b. The results of whole-mount in situ hybridization showed that the expression locations of bmp2a and bmp2b in zebrafish were different in different periods (24 hpf, 48 hpf, 72 hpf), which revealed potential functional differentiation between bmp2a and bmp2b. Phenotypic analysis showed that bmp2a mutations caused partial rib and vertebral deformities in zebrafish, while bmp2b^−/− embryos died massively after 12 hpf due to abnormal somite formation. We further explored the expression pattern changes of genes (bmp2a, bmp2b, smad1, fgf4, runx2b, alp) related to skeletal development at different developmental stages (20 dpf, 60 dpf, 90 dpf) in wild-type and bmp2a^−/− zebrafish. The results showed that the expression of runx2b in bmp2a^−/− was significantly downregulated at three stages and the expression of other genes were significantly downregulated at 90 dpf compared with wild-type zebrafish. The study revealed functional differentiation of bmp2a and bmp2b in zebrafish embryonic and skeletal development.     

Cis-regulatory characterization of sequence conservation surrounding the Hox4 genes

Hox genes are key regulators of anterior-posterior axis patterning and have a major role in hindbrain development. The zebrafish Hox4 paralogs have strong overlapping activities in hindbrain rhombomeres 7 and 8, in the spinal cord and in the pharyngeal arches. With the aim to predict enhancers that act on the hoxa4a, hoxb4a, hoxc4a and hoxd4a genes, we used sequence conservation around the Hox4 genes to analyze all fish:human conserved non-coding sequences by reporter assays in stable zebrafish transgenesis. Thirty-four elements were functionally tested in GFP reporter gene constructs and more than 100 F1 lines were analyzed to establish a correlation between sequence conservation and cis-regulatory function, constituting a catalog of Hox4 CNEs. Sixteen tissue-specific enhancers could be identified. Multiple alignments of the CNEs revealed paralogous cis-regulatory sequences, however, the CNE sequence similarities were found not to correlate with tissue specificity. To identify ancestral enhancers that direct Hox4 gene activity, genome sequence alignments of mammals, teleosts, horn shark and the cephalochordate amphioxus, which is the most basal extant chordate possessing a single prototypical Hox cluster, were performed. Three elements were identified and two of them exhibited regulatory activity in transgenic zebrafish, however revealing no specificity. Our data show that the approach to identify cis-regulatory sequences by genome sequence alignments and subsequent testing in zebrafish transgenesis can be used to define enhancers within the Hox clusters and that these have significantly diverged in their function during evolution.     

The Ndst Gene Family in Zebrafish: Role of Ndst1b in Pharyngeal Arch Formation

Heparan sulfate (HS) proteoglycans are ubiquitous components of the extracellular matrix and plasma membrane of metazoans. The sulfation pattern of the HS glycosaminoglycan chain is characteristic for each tissue and changes during development. The glucosaminyl N-deacetylase/N-sulfotransferase (NDST) enzymes catalyze N-deacetylation and N-sulfation during HS biosynthesis and have a key role in designing the sulfation pattern. We here report on the presence of five NDST genes in zebrafish. Zebrafish ndst1a, ndst1b, ndst2a and ndst2b represent duplicated mammalian orthologues of NDST1 and NDST2 that arose through teleost specific genome duplication. Interestingly, the single zebrafish orthologue ndst3, is equally similar to tetrapod Ndst3 and Ndst4. It is likely that a local duplication in the common ancestor of lobe-finned fish and tetrapods gave rise to these two genes. All zebrafish Ndst genes showed distinct but partially overlapping expression patterns during embryonic development. Morpholino knockdown of ndst1b resulted in delayed development, craniofacial cartilage abnormalities, shortened body and pectoral fin length, resembling some of the features of the Ndst1 mouse knockout.     

The ascl1a and dlx genes have a regulatory role in the development of GABAergic interneurons in the zebrafish diencephalon

During development of the mouse forebrain interneurons, the Dlx genes play a key role in a gene regulatory network (GRN) that leads to the GABAergic phenotype. Here, we have examined the regulatory relationships between the ascl1a, dlx, and gad1b genes in the zebrafish forebrain. Expression of ascl1a overlaps with dlx1a in the telencephalon and diencephalon during early forebrain development. The loss of Ascl1a function results in a loss of dlx expression, and subsequent losses of dlx5a and gad1b expression in the diencephalic prethalamus and hypothalamus. Loss of Dlx1a and Dlx2a function, and, to a lesser extent, of Dlx5a and Dlx6a, impairs gad1b expression in the prethalamus and hypothalamus. We conclude that dlx1a/2a act downstream of ascl1a but upstream of dlx5a/dlx6a and gad1b to activate GABAergic specification. This pathway is conserved in the diencephalon, but has diverged between mammals and teleosts in the telencephalon.     

Additional hox clusters in the zebrafish: divergent expression patterns belie equivalent activities of duplicate hoxB5 genes

SUMMARY The evolution of metazoan body plans has involved changes to the Hox genes, which are involved in patterning the body axis and display striking evolutionary conservation of structure and expression. Invertebrates contain a single Hox cluster whereas tetrapods possess four clusters. The zebrafish has seven unlinked hox clusters, a finding that is difficult to reconcile with the notion that genomic complexity, reflected by Hox cluster number, and morphological complexity are causally linked, as the body plan of the zebrafish is not obviously more complex than that of the mouse or human. Why have the additional hox genes in zebrafish been conserved? To address the role of these additional zebrafish hox genes, we have examined the duplicate hoxB5 genes, hoxB5a, and hoxB5b. Conservation of gene duplicates can occur when one gene acquires a new function (neofunctionalization), or when the ancestral function is divided between the two duplicates (subfunctionalization). hoxB5a and hoxB5b are expressed in distinct domains, and their combined expression domain is strikingly similar to that of single Hoxb5 genes in other species. The biochemical functions encoded by the two genes were studied by overexpression, which resulted in identical developmental defects in the anterior hindbrain and cranial neural crest, suggesting strongly that hoxB5a and hoxB5b have equivalent biochemical properties with respect to early development. From these studies, we conclude that conservation of hoxB5a and hoxB5b is likely the result of division of the ancestral Hoxb5 function between the two genes, without significant changes in biochemical activity. These results suggest a resolution to the conundrum of the extra hox genes and clusters in the zebrafish, since if any of the additional hox genes in the zebrafish are similarly subfunctionalized, they are unlikely to supply novel genetic functions. Thus, the morphological complexity potentially conferred by the majority of additional zebrafish hox clusters may not be substantially greater than that conferred by the four tetrapod clusters.     

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal–gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.     

Comparative genomic analysis of teleost fish <Emphasis Type="Italic">bmal</Emphasis> genes

Bmal1 (Brain and muscle ARNT like 1) gene is a key circadian clock gene. Tetrapods also have the second Bmal gene, Bmal2. Fruit fly has only one bmal1/cycle gene. Interrogation of the five teleost fish genome sequences coupled with phylogenetic and splice site analyses found that zebrafish have two bmal1 genes, bmal1a and bmal1b, and bmal2a; Japanese pufferfish (fugu), green spotted pufferfish (tetraodon) and Japanese medaka fish each have two bmal2 genes, bmal2a and bmal2b, and bmal1a; and three-spine stickleback have bmal1a and bmal2b. Syntenic analysis further indicated that zebrafish bmal1a/bmal1b, and fugu, tetraodon and medaka bmal2a/bmal2b are ancient duplicates. Although the dN/dS ratios of these four fish bmal duplicates are all <1, implicating they have been under purifying selection, the Tajima relative rate test showed that fugu, tetraodon and medaka bmal2a/bmal2b have asymmetric evolutionary rates, suggesting that one of these duplicates have been subject to positive selection or relaxed functional constraint. These results support the notion that teleost fish bmal genes were derived from the fish-specific genome duplication (FSGD), divergent resolution following the duplication led to retaining different ancient bmal duplicates in different fishes, which could have shaped the evolution of the complex teleost fish timekeeping mechanisms. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Zebrafish β-adrenergic receptor mRNA expression and control of pigmentation

Beta adrenergic receptors (β-ARs) are members of the G-protein-coupled receptor superfamily and mediate various physiological processes in many species. The expression patterns and functions of β-ARs in zebrafish are, however, largely unknown. We have identified zebrafish β-AR orthologs, which we have designated as adrb1, adrb2a, adrb2b, adrb3a and adrb3b. adrb1 was found to be expressed in the heart and brain. Expression of adrb2a predominated in the brain and skin, whereas adrb2b was found to be highly expressed in muscle, pancreas and liver. Both adrb3a and adrb3b were exclusively expressed in blood. Knock-down of these β-ARs by morpholino oligonucleotides revealed a functional importance of adrb2a in pigmentation. Expression of atp5a1 and atp5b, genes that encode subunits of F1F0-ATPase, which is known to be involved in pigmentation, was significantly increased by knock-down of adrb2a. Our data suggest that adrb2a may regulate pigmentation, partly by modulating F1F0-ATPase.     

prox1b Activity is essential in zebrafish lymphangiogenesis

Background

The lymphatic vascular system, draining interstitial fluids from most tissues and organs, exerts crucial functions in several physiological and pathological processes. Lymphatic system development depends on Prox1, the first marker to be expressed in the endothelial cells of the cardinal vein from where lymph vessels originate. Prox1 ortholog in the optically clear, easily manipulated zebrafish model has been previously isolated and its contribution to lymphangiogenesis has been clarified. Because of a round of genome duplication occurred at the base of teleosts radiation, several zebrafish genes have been retained in duplicate through evolution. We investigated for the presence of additional prox1 genes and determined their role in zebrafish lymphangiogenesis.

Methodology/Principal Findings

We isolated a second ortholog, named prox1b, and analyzed its expression during development by whole mount in situ hybridization (WISH). We detected strong prox1b expression in the endothelium of the posterior cardinal vein (PCV) from where lymphatic precursors originate. To analyze prox1b involvement in lymphangiogenesis we utilized the fli1:GFP transgenics and followed the formation of the toracic duct (TD), the primary lymph vessel in fish, after prox1b knockdown. Our findings clearly demonstrated that the absence of prox1b activity severely hampers the formation of the TD.

Conclusions/Significance

This work provides substantial progress toward the understanding of zebrafish lymphangiogenesis. In light of the features shared by the lymphatic systems of zebrafish and higher vertebrates, the establishment of such lymphatic model will provide a powerful tool to study, for instance, disorders of body fluid homeostasis, inflammation and cancer metastasis, and may ultimately contribute to novel therapies.     

The Origin and Diversity of Cpt1 Genes in Vertebrate Species

The Carnitine palmitoyltransferase I (Cpt1) gene family plays a crucial role in energy homeostasis since it is required for the occurrence of fatty acid β-oxidation in the mitochondria. The exact gene repertoire in different vertebrate lineages is variable. Presently, four genes are documented: Cpt1a, also known as Cpt1a1, Cpt1a2; Cpt1b and Cpt1c. The later is considered a mammalian innovation resulting from a gene duplication event in the ancestor of mammals, after the divergence of sauropsids. In contrast, Cpt1a2 has been found exclusively in teleosts. Here, we reassess the overall evolutionary relationships of Cpt1 genes using a combination of approaches, including the survey of the gene repertoire in basal gnathostome lineages. Through molecular phylogenetics and synteny studies, we find that Cpt1c is most likely a rapidly evolving orthologue of Cpt1a2. Thus, Cpt1c is present in other lineages such as cartilaginous fish, reptiles, amphibians and the coelacanth. We show that genome duplications (2R) and variable rates of sequence evolution contribute to the history of Cpt1 genes in vertebrates. Finally, we propose that loss of Cpt1b is the likely cause for the unusual energy metabolism of elasmobranch.     

Phylogenetic Timing of the Fish-Specific Genome Duplication Correlates with the Diversification of Teleost Fish

For many genes, ray-finned fish (Actinopterygii) have two paralogous copies, where only one ortholog is present in tetrapods. The discovery of an additional, almost-complete set of Hox clusters in teleosts (zebrafish, pufferfish, medaka, and cichlid) but not in basal actinopterygian lineages (Polypterus) led to the formulation of the fish-specific genome duplication hypothesis. The phylogenetic timing of this genome duplication during the evolution of ray-finned fish is unknown, since only a few species of basal fish lineages have been investigated so far. In this study, three nuclear genes (fzd8, sox11, tyrosinase) were sequenced from sturgeons (Acipenseriformes), gars (Semionotiformes), bony tongues (Osteoglossomorpha), and a tenpounder (Elopomorpha). For these three genes, two copies have been described previously teleosts (e.g., zebrafish, pufferfish), but only one orthologous copy is found in tetrapods. Individual gene trees for these three genes and a concatenated dataset support the hypothesis that the fish-specific genome duplication event took place after the split of the Acipenseriformes and the Semionotiformes from the lineage leading to teleost fish but before the divergence of Osteoglossiformes. If these three genes were duplicated during the proposed fish-specific genome duplication event, then this event separates the species-poor early-branching lineages from the species-rich teleost lineage. The additional number of genes resulting from this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish.[Reviewing Editor: Martin Kreitman]