Genome evolution and meiotic maps by massively parallel DNA sequencing: spotted gar, an outgroup for the teleost genome duplication

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

The gar is a fish... is a bird... is a mammal?

Teleosts have emerged as important model organisms, yet their ancestrally duplicated genomes sometimes complicate developmental genetic analyses and comparisons to humans. A new genome sequence of spotted gar, a fish related to teleosts but lacking a duplicated genome, now helps to bridge human and teleost biology.

Epilepsy, Behavioral Abnormalities, and Physiological Comorbidities in Syntaxin-Binding Protein 1 (STXBP1) Mutant Zebrafish

Mutations in the synaptic machinery gene syntaxin-binding protein 1, STXBP1 (also known as MUNC18-1), are linked to childhood epilepsies and other neurodevelopmental disorders. Zebrafish STXBP1 homologs (stxbp1a and stxbp1b) have highly conserved sequence and are prominently expressed in the larval zebrafish brain. To understand the functions of stxbp1a and stxbp1b, we generated loss-of-function mutations using CRISPR/Cas9 gene editing and studied brain electrical activity, behavior, development, heart physiology, metabolism, and survival in larval zebrafish. Homozygous stxbp1a mutants exhibited a profound lack of movement, low electrical brain activity, low heart rate, decreased glucose and mitochondrial metabolism, and early fatality compared to controls. On the other hand, homozygous stxbp1b mutants had spontaneous electrographic seizures, and reduced locomotor activity response to a movement-inducing “dark-flash” visual stimulus, despite showing normal metabolism, heart rate, survival, and baseline locomotor activity. Our findings in these newly generated mutant lines of zebrafish suggest that zebrafish recapitulate clinical phenotypes associated with human syntaxin-binding protein 1 mutations.

Expression analysis and characterization of an autosome-localized tesk1 gene in half-smooth tongue sole (Cynoglossus semilaevis)

Testis-specific protein kinase 1 (tesk1) represents a conserved gene family functioning in many cellular processes. In this study, we cloned and characterized an autosome-localized tesk1 gene (Altesk1) from Cynoglossus semilaevis. The open reading frame consists of 2088 nucleotides and encodes a 665 amino acid polypeptide. Phylogenetic analyses show that vertebrate Tesk1s are divided into two clusters based on protein length and AlTesk1 belongs to "long-type" group. Semi-quantitative PCR reveals that Altesk1 is predominantly expressed in ovary, despite of relatively low detection in some other tissues. Among different development stages, Altesk1 transcripts are only observed in ovary samples of 210-day and 1-year fish. In situ hybridization analyses have further confirmed its major localization in oocyte cells. Comparison of methylation patterns in different sexual genotypes reveals the low methylation level of Altesk1 promoter in female, which is consistent with Altesk1 high expression level in female. Taken together, this is the first time that tesk1 gene has been found to show female-biased expression and in view of this, we postulate that AlTesk1 might be involved in some cellular processes specific in ovary, e.g. oogenesis.

Rab32 and Rab38 genes in chordate pigmentation: an evolutionary perspective

Background

The regulation of cellular membrane trafficking in all eukaryotes is a very complex mechanism, mostly regulated by the Rab family proteins. Among all membrane-enclosed organelles, melanosomes are the cellular site for synthesis, storage and transport of melanin granules, making them an excellent model for studies on organelle biogenesis and motility. Specific Rab proteins, as Rab32 and Rab38, have been shown to play a key role in melanosome biogenesis. We analysed the Rab32 and Rab38 genes in the teleost zebrafish and in the cephalochordate amphioxus, gaining insight on their evolutionary history following gene and genome duplications.

Results

We studied the molecular evolution of Rab supergroup III in deuterostomes by phylogenetic reconstruction, intron and synteny conservation. We discovered a novel amino acid stretch, named FALK, shared by three related classes belonging to Rab supergroup III: Rab7L1, Rab32LO and Rab32/Rab38. Among these, we demonstrated that the Rab32LO class, already present in the last common eukaryotic ancestor, was lost in urochordates and vertebrates. Synteny shows that one zebrafish gene, Rab38a, which is expressed in pigmented cells, retained the linkage with tyrosinase, a protein essential for pigmentation. Moreover, the chromosomal linkage of Rab32 or Rab38 with a member of the glutamate receptor metabotropic (Grm) family has been retained in all analysed gnathostomes, suggesting a conserved microsynteny in the vertebrate ancestor. Expression patterns of Rab32 and Rab38 genes in zebrafish, and Rab32/38 in amphioxus, indicate their involvement in development of pigmented cells and notochord.

Conclusions

Phylogenetic, intron conservation and synteny analyses point towards an evolutionary scenario based on a duplication of a single invertebrate Rab32/38 gene giving rise to vertebrate Rab32 and Rab38. The expression patterns of Rab38 paralogues highlight sub-functionalization event. Finally, the discovery of a chromosomal linkage between the Rab32 or Rab38 gene with a Grm opens new perspectives on possible conserved bystander gene regulation across the vertebrate evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0596-1) contains supplementary material, which is available to authorized users.

Allatostatin-type A, kisspeptin and galanin GPCRs and putative ligands as candidate regulatory factors of mantle function

Allatostatin-type A (AST-A), kisspeptin (KISS) and galanin (GAL) G-protein coupled receptor (GPCR) systems share a common ancestral origin in arthropods and the vertebrates where they regulate metabolism and reproduction. The molluscs are the second most diverse phylum in the animal kingdom, they occupy an important phylogenetic position, and their genome is more similar to deuterostomes than the arthropods and nematodes and thus they are good models for studies of gene family evolution and function. This mini-review intends to extend the current knowledge about AST-A, KISS and GAL GPCR system evolution and their putative function in the mollusc mantle. Comparative evolutionary analysis of the target GPCR systems was established by identifying homologues in genomes and tissue transcriptome datasets available for molluscs and comparing them to those of other metazoan systems. Studies in arthropods have revealed the existence of the AST-A system but the loss of homologues of the KISS and GAL systems. Homologues of the insect AST-AR and vertebrate KISSR genes were found in molluscs but putative GALR genes were absent. Receptor gene number suggested that members of this family have suffered lineage specific evolution during the molluscan radiation. In molluscs, orthologues of the insect AST-A peptides were not identified but buccalin peptides that are structurally related were identified and are putative receptor agonists. The identification of AST-AR and KISSR genes in molluscs strengthens the hypotheses that in metazoans members of the AST-AR subfamily share evolutionary proximity with KISSRs. The variable number of receptors and large repertoire of buccalin peptides may be indicative of the functional diversity of the AST-AR/KISSR systems in molluscs. The identification of AST-A and KISS receptors and ligands in the mantle transcriptome indicates that in molluscs they may have acquired a novel function and may play a role in shell development or sensory detection in the mantle.

Exploring a Nonmodel Teleost Genome Through RAD Sequencing-Linkage Mapping in Common Pandora, Pagellus erythrinus and Comparative Genomic Analysis

Common pandora (Pagellus erythrinus) is a benthopelagic marine fish belonging to the teleost family Sparidae, and a newly recruited species in Mediterranean aquaculture. The paucity of genetic information relating to sparids, despite their growing economic value for aquaculture, provides the impetus for exploring the genomics of this fish group. Genomic tool development, such as genetic linkage maps provision, lays the groundwork for linking genotype to phenotype, allowing fine-mapping of loci responsible for beneficial traits. In this study, we applied ddRAD methodology to identify polymorphic markers in a full-sib family of common pandora. Employing the Illumina MiSeq platform, we sampled and sequenced a size-selected genomic fraction of 99 individuals, which led to the identification of 920 polymorphic loci. Downstream mapping analysis resulted in the construction of 24 robust linkage groups, corresponding to the karyotype of the species. The common pandora linkage map showed varying degrees of conserved synteny with four other teleost genomes, namely the European seabass (Dicentrarchus labrax), Nile tilapia (Oreochromis niloticus), stickleback (Gasterosteus aculeatus), and medaka (Oryzias latipes), suggesting a conserved genomic evolution in Sparidae. Our work exploits the possibilities of genotyping by sequencing to gain novel insights into genome structure and evolution. Such information will boost the study of cultured species and will set the foundation for a deeper understanding of the complex evolutionary history of teleosts.

Gene map of large yellow croaker (Larimichthys crocea) provides insights into teleost genome evolution and conserved regions associated with growth

The genetic map of a species is essential for its whole genome assembly and can be applied to the mapping of important traits. In this study, we performed RNA-seq for a family of large yellow croakers (Larimichthys crocea) and constructed a high-density genetic map. In this map, 24 linkage groups comprised 3,448 polymorphic SNP markers. Approximately 72.4% (2,495) of the markers were located in protein-coding regions. Comparison of the croaker genome with those of five model fish species revealed that the croaker genome structure was closer to that of the medaka than to the remaining four genomes. Because the medaka genome preserves the teleost ancestral karyotype, this result indicated that the croaker genome might also maintain the teleost ancestral genome structure. The analysis also revealed different genome rearrangements across teleosts. QTL mapping and association analysis consistently identified growth-related QTL regions and associated genes. Orthologs of the associated genes in other species were demonstrated to regulate development, indicating that these genes might regulate development and growth in croaker. This gene map will enable us to construct the croaker genome for comparative studies and to provide an important resource for selective breeding of croaker.

Recombination patterns reveal information about centromere location on linkage maps

Linkage mapping is often used to identify genes associated with phenotypic traits and for aiding genome assemblies. Still, many emerging maps do not locate centromeres - an essential component of the genomic landscape. Here, we demonstrate that for genomes with strong chiasma interference, approximate centromere placement is possible by phasing the same data used to generate linkage maps. Assuming one obligate crossover per chromosome arm, information about centromere location can be revealed by tracking the accumulated recombination frequency along linkage groups, similar to half-tetrad analyses. We validate the method on a linkage map for sockeye salmon (Oncorhynchus nerka) with known centromeric regions. Further tests suggest that the method will work well in other salmonids and other eukaryotes. However, the method performed weakly when applied to a male linkage map (rainbow trout; O. mykiss) characterized by low and unevenly distributed recombination - a general feature of male meiosis in many species. Further, a high frequency of double crossovers along chromosome arms in barley reduced resolution for locating centromeric regions on most linkage groups. Despite these limitations, our method should work well for high-density maps in species with strong recombination interference and will enrich many existing and future mapping resources.

Correlated expression of retrocopies and parental genes in zebrafish

Previous studies of the function and evolution of retrocopies in plants, Drosophila and non-mammalian chordates provided new insights into the origin of novel genes. However, little is known about retrocopies and their parental genes in teleosts, and it remains obscure whether there is any correlation between them. The present study aimed to characterize the spatial and temporal expression profiles of retrogenes and their parental genes based on RNA-Seq data from Danio rerio embryos and tissues from adult. Using a modified pipeline, 306 retrocopies were identified in the zebrafish genome, most of which exhibited ancient retroposition, and 76 of these showed a Ks < 2.0. Expression of a retrocopy is generally expected to present no correlation with its parental gene, as regulatory regions are not part of the retroposition event. Here, this assumption was tested based on RNA-Seq data from eight stages and thirteen tissue types of zebrafish. However, the result suggested that retrocopies displayed correlated expression with their parental genes. The level of correlation was found to decrease during embryogenesis, but to increase slightly within a tissue using Ks as the proxy for the divergence time. Tissue specificity was also observed: retrocopies were found to be expressed at a more specific level compared with their parental genes. Unlike Drosophila, which has sex chromosomes, zebrafish do not show testis-biased expression. Our study elaborated temporal and spatial patterns of expression of retrocopies in zebrafish, examined the correlation between retrocopies and parental genes and analyzed potential source of regulated elements of retrocopies, which lay a foundation for further functional study of retrocopies.

Divergent evolution of two corticotropin-releasing hormone (CRH) genes in teleost fishes

Genome duplication, thought to have happened twice early in vertebrate evolution and a third time in teleost fishes, gives rise to gene paralogs that can evolve subfunctions or neofunctions via sequence and regulatory changes. To explore the evolution and functions of corticotropin-releasing hormone (CRH), we searched sequenced teleost genomes for CRH paralogs. Our phylogenetic and synteny analyses indicate that two CRH genes, crha and crhb, evolved via duplication of crh1 early in the teleost lineage. We examined the expression of crha and crhb in two teleost species from different orders: an African cichlid, Burton's mouthbrooder, (Astatotilapia burtoni; Order Perciformes) and zebrafish (Danio rerio; Order Cypriniformes). Furthermore, we compared expression of the teleost crha and crhb genes with the crh1 gene of an outgroup to the teleost clade: the spotted gar (Lepisosteus oculatus). In situ hybridization for crha and crhb mRNA in brains and eyes revealed distinct expression patterns for crha in different teleost species. In the cichlid, crha mRNA was found in the retina but not in the brain. In zebrafish, however, crha mRNA was not found in the retina, but was detected in the brain, restricted to the ventral hypothalamus. Spotted gar crh1 was found in the retina as well as the brain, suggesting that the ancestor of teleost fishes likely had a crh1 gene expressed in both retina and brain. Thus, genome duplication may have freed crha from constraints, allowing it to evolve distinct sequences, expression patterns, and likely unique functions in different lineages.

High-Density Genetic Mapping with Interspecific Hybrids of Two Sea Urchins, Strongylocentrotus nudus and S. intermedius, by RAD Sequencing

Sea urchins have long been used as research model organisms for developmental biology and evolutionary studies. Some of them are also important aquaculture species in East Asia. In this work, we report the construction of RAD-tag based high-density genetic maps by genotyping F₁ interspecific hybrids derived from a crossing between a female sea urchin Strongylocentrotus nudus and a male Strongylocentrotus intermedius. With polymorphisms present in these two wild individuals, we constructed a female meiotic map containing 3,080 markers for S. nudus, and a male meiotic map for S. intermedius which contains 1,577 markers. Using the linkage maps, we were able to anchor a total of 1,591 scaffolds (495.9 Mb) accounting for 60.8% of the genome assembly of Strongylocentrotus purpuratus. A genome-wide scan resulted in the identification of one putative QTL for body size which spanned from 25.3 cM to 30.3 cM. This study showed the efficiency of RAD-Seq based high-density genetic map construction using F₁ progenies for species with no prior genomic information. The genetic maps are essential for QTL mapping and are useful as framework to order and orientate contiguous scaffolds from sea urchin genome assembly. The integration of the genetic map with genome assembly would provide an unprecedented opportunity to conduct QTL analysis, comparative genomics, and population genetics studies.

Molecular characterization, origin, and evolution of teleost p68 gene family: Insights from Japanese flounder, Paralichthys olivaceus

Two rounds of whole-genome duplication occurred in the common ancestor of vertebrates. Later, a third round genome duplication occurred in the teleost fishes. As a prototype member of DEAD-box RNA helicases, the function of p68 helicase in development has been well investigated in human, however, limited information is available regarding the regulatory function of this gene in the development of teleosts. In this study, being an important farmed fish in North China, Japanese flounder (Paralichthys olivaceus) was used as model fish to investigate the role of p68 gene in teleost development. Two p68 genes were first identified from Japanese flounder. Molecular characterization of them was performed by analyzing the exon-intron boundaries. Then, we confirmed that such two teleost p68 genes originated from teleost-specific genome duplication through phylogenetic and synteny analyses. Additionally, comparative analyses of amino acid sequences, variation in selective pressure, and expression profiles of p68 genes revealed probable sub-functionalization fate of teleost p68 genes after the duplication. Therefore, this study supplements the evolutionary properties of teleost p68 gene family and provides the groundwork for further studying the regulatory function of p68 genes in the development of teleosts.

Evolutionary Fate of the Androgen Receptor-Signaling Pathway in Ray-Finned Fishes with a Special Focus on Cichlids

The emergence of the steroid system is coupled to the evolution of multicellular animals. In vertebrates in particular, the steroid receptor repertoire has been shaped by genome duplications characteristic to this lineage. Here, we investigate for the first time the composition of the androgen receptor–signaling pathway in ray-finned fish genomes by focusing in particular on duplicates that emerged from the teleost-specific whole-genome duplication. We trace lineage- and species-specific duplications and gene losses for the genomic and nongenomic pathway of androgen signaling and subsequently investigate the sequence evolution of these genes. In one particular fish lineage, the cichlids, we find evidence for differing selection pressures acting on teleost-specific whole-genome duplication paralogs at a derived evolutionary stage. We then look into the expression of these duplicated genes in four cichlid species from Lake Tanganyika indicating, once more, rapid changes in expression patterns in closely related fish species. We focus on a particular case, the cichlid specific duplication of the rac1 GTPase, which shows possible signs of a neofunctionalization event.

Evolution of Vertebrate Adam Genes; Duplication of Testicular Adams from Ancient Adam9/9-like Loci

Members of the disintegrin metalloproteinase (ADAM) family have important functions in regulating cell-cell and cell-matrix interactions as well as cell signaling. There are two major types of ADAMs: the somatic ADAMs (sADAMs) that have a significant presence in somatic tissues, and the testicular ADAMs (tADAMs) that are expressed predominantly in the testis. Genes encoding tADAMs can be further divided into two groups: group I (intronless) and group II (intron-containing). To date, tAdams have only been reported in placental mammals, and their evolutionary origin and relationship to sAdams remain largely unknown. Using phylogenetic and syntenic tools, we analyzed the Adam genes in various vertebrates ranging from fishes to placental mammals. Our analyses reveal duplication and loss of some sAdams in certain vertebrate species. In particular, there exists an Adam9-like gene in non-mammalian vertebrates but not mammals. We also identified putative group I and group II tAdams in all amniote species that have been examined. These tAdam homologues are more closely related to Adams 9 and 9-like than to other sAdams. In all amniote species examined, group II tAdams lie in close vicinity to Adam9 and hence likely arose from tandem duplication, whereas group I tAdams likely originated through retroposition because of their lack of introns. Clusters of multiple group I tAdams are also common, suggesting tandem duplication after retroposition. Therefore, Adam9/9-like and some of the derived tAdam loci are likely preferred targets for tandem duplication and/or retroposition. Consistent with this hypothesis, we identified a young retroposed gene that duplicated recently from Adam9 in the opossum. As a result of gene duplication, some tAdams were pseudogenized in certain species, whereas others acquired new expression patterns and functions. The rapid duplication of Adam genes has a major contribution to the diversity of ADAMs in various vertebrate species.

Physiological response of alligator gar juveniles (Atractosteus spatula) exposed to sub-lethal doses of pollutants

Alligator gar populations have declined because of overfishing, habitat loss and pollution. Over time, the exposure to different pollutants have affected these fishes as a consequence of their high trophic level, bottom-dwelling habits and long life span. In order to evaluate the physiological effects of pollutants on alligator gar, juveniles (6, 12 and 24 months) were exposed to sub-lethal doses of diazinon, β-naphthoflavone (BNF) and 17 β-estradiol (E2) by intraperitoneal injection. After 2 days of exposure, liver samples were taken to determine the activities of acetylcholinesterase, butyrylcholinesterase and carboxylesterase; alkaline and acid phosphatases (ALP and ACP); ethoxyresorufin o-deethylase (EROD); glutathione s-transferase (GST); superoxide dismutase (SOD), and vitellogenin (VTG) concentration. Two additional bioassays consisting on the exposure of compounds through water or food were performed and after 4 and 28 days, respectively, biomarkers were determined. All esterases were inhibited in organisms exposed to diazinon as well as in 6-months gar exposed to E2 and BNF. In contrast, ALP activity increased in gar exposed to diazinon and E2, while ACP activity did not show any variations. No EROD activity was registered after exposure to the different pollutants, despite being one of the most sensitive and common detoxification biomarkers used for fishes. GST activity reduction was detected when gar were exposed to E2 and BNF, while SOD activity increased after exposure to diazinon and E2. Finally, VTG levels were higher in animals exposed to E2 compared to other treatments. Overall, these results suggest that alligator gar juveniles have a low biotransformation metabolism and show that they are especially sensitive to those pollutants affecting the nervous system.

Evolution and functional diversity of aquaporins

In this review, we provide a brief synopsis of the evolution and functional diversity of the aquaporin gene superfamily in prokaryotic and eukaryotic organisms. Based upon the latest data, we discuss the expanding list of molecules shown to permeate the central pore of aquaporins, and the unexpected diversity of water channel genes in Archaea and Bacteria. We further provide new insight into the origin by horizontal gene transfer of plant glycerol-transporting aquaporins (NIPs), and the functional co-option and gene replacement of insect glycerol transporters. Finally, we discuss the origins of four major grades of aquaporins in Eukaryota, together with the increasing repertoires of aquaporins in vertebrates.

The sea lamprey meiotic map improves resolution of ancient vertebrate genome duplications

It is generally accepted that many genes present in vertebrate genomes owe their origin to two whole-genome duplications that occurred deep in the ancestry of the vertebrate lineage. However, details regarding the timing and outcome of these duplications are not well resolved. We present high-density meiotic and comparative genomic maps for the sea lamprey (Petromyzon marinus), a representative of an ancient lineage that diverged from all other vertebrates ∼550 million years ago. Linkage analyses yielded a total of 95 linkage groups, similar to the estimated number of germline chromosomes (1n ∼ 99), spanning a total of 5570.25 cM. Comparative mapping data yield strong support for the hypothesis that a single whole-genome duplication occurred in the basal vertebrate lineage, but do not strongly support a hypothetical second event. Rather, these comparative maps reveal several evolutionarily independent segmental duplications occurring over the last 600+ million years of chordate evolution. This refined history of vertebrate genome duplication should permit more precise investigations of vertebrate evolution.

RAD mapping reveals an evolving, polymorphic and fuzzy boundary of a plant pseudoautosomal region

How loss of genetic exchanges (recombination) evolves between sex chromosomes is a long-standing question. Suppressed recombination may evolve when a sexually antagonistic (SA) polymorphism occurs in a partially sex-linked 'pseudoautosomal' region (or 'PAR'), maintaining allele frequency differences between the two sexes, and creating selection for closer linkage with the fully sex-linked region of the Y chromosome in XY systems, or the W in ZW sex chromosome systems. Most evidence consistent with the SA polymorphism hypothesis is currently indirect, and more studies of the genetics and population genetics of PAR genes are clearly needed. The sex chromosomes of the plant Silene latifolia are suitable for such studies, as they evolved recently and the loss of recombination could still be ongoing. Here, we used RAD sequencing to genetically map sequences in this plant, which has a large genome (c. 3 gigabases) and no available whole-genome sequence. We mapped 83 genes on the sex chromosomes, and comparative mapping in the related species S. vulgaris supports previous evidence for additions to an ancestral PAR and identified at least 12 PAR genes. We describe evidence that recombination rates have been reduced in meiosis of both sexes, and differences in recombination between S. latifolia families suggest ongoing recombination suppression. Large allele frequency differences between the sexes were found at several loci closely linked to the PAR boundary, and genes in different regions of the PAR showed striking sequence diversity patterns that help illuminate the evolution of the PAR.

C2H2 zinc finger proteins of the SP/KLF, Wilms tumor, EGR, Huckebein, and Klumpfuss families in metazoans and beyond

Specificity proteins (SPs) and Krüppel-Like Factors (KLFs) are C2H2-type zinc finger transcription factors that play essential roles in differentiation, development, proliferation and cell death. SP/KLF proteins, similarly to Wilms tumor protein 1 (WT1), Early Growth Response (EGR), Huckebein, and Klumpfuss, prefer to bind GC-rich sequences such as GC-box and CACCC-box (GT-box). We searched various genomes and transcriptomes of metazoans and single-cell holozoans for members of these families. Seven groups of KLFs (KLFA-G) and three groups of SPs (SPA-C) were identified in the three lineages of Bilateria (Deuterostomia, Ecdysozoa, and Lophotrochozoa). The last ancestor of jawed vertebrates was inferred to have at least 18 KLFs (group A: KLF1/2/4/17, group B: KLF3/8/12; group C: KLF5/5l; group D: KLF6/7; group E: KLF9/13/16; group F: KLF10/KLF11; group G: KLF15/15l) and 10 SPs (group A: SP1/2/3/4; group B: SP5/5l; group C: SP6/7/8/9), since they were found in both cartilaginous and boned fishes. Placental mammals have added KLF14 (group E) and KLF18 (group A), and lost KLF5l (KLF5-like) and KLF15l (KLF15-like). Multiple KLF members were found in basal metazoans (Ctenophora, Porifera, Placozoa, and Cnidaria). Ctenophora has the least number of KLFs and no SPs, which could be attributed to its proposed sister group relationship to other metazoans or gene loss. While SP, EGR and Klumpfuss were only detected in metazoans, KLF, WT1, and Huckebein are present in nonmetazoan holozoans. Of the seven metazoan KLF groups, only KLFG, represented by KLF15 in human, was found in nonmetazoans. In addition, two nonmetazoan groups of KLFs are present in Choanoflagellatea and Filasterea. WT1 could be evolutionarily the earliest among these GC/GT-box-binding families due to its sole presence in Ichthyosporea.

Exploring the utility of cross-laboratory RAD-sequencing datasets for phylogenetic analysis

BACKGROUND

Restriction site-Associated DNA sequencing (RAD-Seq) is widely applied to generate genome-wide sequence and genetic marker datasets. RAD-Seq has been extensively utilised, both at the population level and across species, for example in the construction of phylogenetic trees. However, the consistency of RAD-Seq data generated in different laboratories, and the potential use of cross-species orthologous RAD loci in the estimation of genetic relationships, have not been widely investigated. This study describes the use of SbfI RAD-Seq data for the estimation of evolutionary relationships amongst ten teleost fish species, using previously established phylogeny as a benchmark.

RESULTS

The number of orthologous SbfI RAD loci identified decreased with increasing evolutionary distance between the species, with several thousand loci conserved across five salmonid species (divergence ~50 MY), and several hundred conserved across the more distantly related teleost species (divergence ~100-360 MY). The majority (>70%) of loci identified between the more distantly related species were genic in origin, suggesting that the bias of SbfI towards genic regions is useful for identifying distant orthologs. Interspecific single nucleotide variants at each orthologous RAD locus were identified. Evolutionary relationships estimated using concatenated sequences of interspecific variants were congruent with previously published phylogenies, even for distantly (divergence up to ~360 MY) related species.

CONCLUSION

Overall, this study has demonstrated that orthologous SbfI RAD loci can be identified across closely and distantly related species. This has positive implications for the repeatability of SbfI RAD-Seq and its potential to address research questions beyond the scope of the original studies. Furthermore, the concordance in tree topologies and relationships estimated in this study with published teleost phylogenies suggests that similar meta-datasets could be utilised in the prediction of evolutionary relationships across populations and species with readily available RAD-Seq datasets, but for which relationships remain uncharacterised.

Conserved Noncoding Sequences Regulate lhx5 Expression in the Zebrafish Forebrain

The LIM homeobox family protein Lhx5 plays important roles in forebrain development in the vertebrates. The lhx5 gene exhibits complex temporal and spatial expression patterns during early development but its transcriptional regulation mechanisms are not well understood. Here, we have used transgenesis in zebrafish in order to define regulatory elements that drive lhx5 expression in the forebrain. Through comparative genomic analysis we identified 10 non-coding sequences conserved in five teleost species. We next examined the enhancer activities of these conserved non-coding sequences with Tol2 transposon mediated transgenesis. We found a proximately located enhancer gave rise to robust reporter EGFP expression in the forebrain regions. In addition, we identified an enhancer located at approximately 50 kb upstream of lhx5 coding region that is responsible for reporter gene expression in the hypothalamus. We also identify an enhancer located approximately 40 kb upstream of the lhx5 coding region that is required for expression in the prethalamus (ventral thalamus). Together our results suggest discrete enhancer elements control lhx5 expression in different regions of the forebrain.

Comparative Evolution of Duplicated Ddx3 Genes in Teleosts: Insights from Japanese Flounder, Paralichthys olivaceus

Following the two rounds of whole-genome duplication that occurred during deuterostome evolution, a third genome duplication event occurred in the stem lineage of ray-finned fishes. This teleost-specific genome duplication is thought to be responsible for the biological diversification of ray-finned fishes. DEAD-box polypeptide 3 (DDX3) belongs to the DEAD-box RNA helicase family. Although their functions in humans have been well studied, limited information is available regarding their function in teleosts. In this study, two teleost Ddx3 genes were first identified in the transcriptome of Japanese flounder (Paralichthys olivaceus). We confirmed that the two genes originated from teleost-specific genome duplication through synteny and phylogenetic analysis. Additionally, comparative analysis of genome structure, molecular evolution rate, and expression pattern of the two genes in Japanese flounder revealed evidence of subfunctionalization of the duplicated Ddx3 genes in teleosts. Thus, the results of this study reveal novel insights into the evolution of the teleost Ddx3 genes and constitute important groundwork for further research on this gene family.

Genome-wide identification, evolution and expression analysis of nuclear receptor superfamily in Nile tilapia, Oreochromis niloticus

The nuclear receptor (NR) superfamily, which is divided into 7 subfamilies, constitutes one of the largest classes of transcription factors. In this study, through comprehensive database search, we identified all NRs (including 4 novel members) from the tilapia (75), common carp (137), zebrafish (73), fugu (73), tetraodon (72), stickleback (70), medaka (69), coelacanth (55), spotted gar (51) and elephant shark (50). For 21 NRs, two duplicates were found in teleosts, while only one in tetrapods. These duplicates, except those of DAX1, SHP and GCNF found in the elephant shark, were derived from 3R (third round of genome duplication). The linkage duplication of 5 syntenic blocks (comprising 14 duplicated NR couples) in teleosts further supported their 3R origin. Based on transcriptome data from adult tilapia, 53 NRs were found to be expressed in more than one tissue (brain, head kidney, heart, liver, kidney, muscle, ovary and testis), and 4 were tissue-specific, indicating their essential roles in the corresponding tissue. Based on the XX and XY gonadal transcriptome data from four developmental stages, 65 NRs were detected in gonads, with 21, 31, 11 and 29 expressed sexual dimorphically at 5, 30, 90 and 180days after hatching, respectively. The expression of four selected genes was examined by in situ hybridization (ISH) and quantitative PCR (qPCR) to validate the spatial and temporal expression profiles of NRs. Comparative analyses of the expression profiles of duplicated NRs revealed divergence in gene expression as well as gene function. Our results demonstrated that NRs may play important roles in sex determination and gonadal development in teleosts.

The evolutionary fate of alternatively spliced homologous exons after gene duplication

Alternative splicing and gene duplication are the two main processes responsible for expanding protein functional diversity. Although gene duplication can generate new genes and alternative splicing can introduce variation through alternative gene products, the interplay between the two processes is complex and poorly understood. Here, we have carried out a study of the evolution of alternatively spliced exons after gene duplication to better understand the interaction between the two processes. We created a manually curated set of 97 human genes with mutually exclusively spliced homologous exons and analyzed the evolution of these exons across five distantly related vertebrates (lamprey, spotted gar, zebrafish, fugu, and coelacanth). Most of these exons had an ancient origin (more than 400 Ma). We found examples supporting two extreme evolutionary models for the behaviour of homologous axons after gene duplication. We observed 11 events in which gene duplication was accompanied by splice isoform separation, that is, each paralog specifically conserved just one distinct ancestral homologous exon. At other extreme, we identified genes in which the homologous exons were always conserved within paralogs, suggesting that the alternative splicing event cannot easily be separated from the function in these genes. That many homologous exons fall in between these two extremes highlights the diversity of biological systems and suggests that the subtle balance between alternative splicing and gene duplication is adjusted to the specific cellular context of each gene.

Functional Pairing of Class B1 Ligand-GPCR in Cephalochordate Provides Evidence of the Origin of PTH and PACAP/Glucagon Receptor Family

Several hypotheses have been proposed regarding the origin and evolution of the secretin family of peptides and receptors. However, identification of homologous ligand-receptor pairs in invertebrates and vertebrates is difficult because of the low levels of sequence identity between orthologs of distant species. In this study, five receptors structurally related to the vertebrate class B1 G protein-coupled receptor (GPCR) family were characterized from amphioxus (Branchiostoma floridae). Phylogenetic analysis showed that they clustered with vertebrate parathyroid hormone receptors (PTHR) and pituitary adenylate cyclase-activating polypeptide (PACAP)/glucagon receptors. These PTHR-like receptors shared synteny with several PTH and PACAP/glucagon receptors identified in spotted gar, Xenopus, and human, indicating that amphioxus preserves the ancestral chordate genomic organization of these receptor subfamilies. According to recent data by Mirabeau and Joly, amphioxus also expresses putative peptide ligands including homologs of PTH (bfPTH1 and 2) and PACAP/GLUC-like peptides (bfPACAP/GLUCs) that may interact with these receptors. Functional analyses showed that bfPTH1 and bfPTH2 activated one of the amphioxus receptors (bf98C) whereas bfPACAP/GLUCs strongly interacted with bf95. In summary, our data confirm the presence of PTH and PACAP/GLUC ligand-receptor pairs in amphioxus, demonstrating that functional homologs of vertebrate PTH and PACAP/glucagon GPCR subfamilies arose before the cephalochordate divergence from the ancestor of tunicates and vertebrates.

Construction of a high density SNP linkage map of kelp (Saccharina japonica) by sequencing Taq I site associated DNA and mapping of a sex determining locus

BACKGROUND

Kelp (Saccharina japonica) has been intensively cultured in China for almost a century. Its genetic improvement is comparable with that of rice. However, the development of its molecular tools is extremely limited, thus its genes, genetics and genomics. Kelp performs an alternative life cycle during which sporophyte generation alternates with gametophyte generation. The gametophytes of kelp can be cloned and crossed. Due to these characteristics, kelp may serve as a reference for the biological and genetic studies of Volvox, mosses and ferns.

RESULTS

We constructed a high density single nucleotide polymorphism (SNP) linkage map for kelp by restriction site associated DNA (RAD) sequencing. In total, 4,994 SNP-containing physical (tag-defined) RAD loci were mapped on 31 linkage groups. The map expanded a total genetic distance of 1,782.75 cM, covering 98.66% of the expected (1,806.94 cM). The length of RAD tags (85 bp) was extended to 400-500 bp with Miseq method, offering us an easiness of developing SNP chips and shifting SNP genotyping to a high throughput track. The number of linkage groups was in accordance with the documented with cytological methods. In addition, we identified a set of microsatellites (99 in total) from the extended RAD tags. A gametophyte sex determining locus was mapped on linkage group 2 in a window about 9.0 cM in width, which was 2.66 cM up to marker_40567 and 6.42 cM down to marker_23595.

CONCLUSIONS

A high density SNP linkage map was constructed for kelp, an intensively cultured brown alga in China. The RAD tags were also extended so that a SNP chip could be developed. In addition, a set of microsatellites were identified among mapped loci, and a gametophyte sex determining locus was mapped. This map will facilitate the genetic studies of kelp including for example the evaluation of germplasm and the decipherment of the genetic bases of economic traits.

The environmental genomics of metazoan thermal adaptation

Continued and accelerating change in the thermal environment places an ever-greater priority on understanding how organisms are going to respond. The paradigm of ‘move, adapt or die', regarding ways in which organisms can respond to environmental stressors, stimulates intense efforts to predict the future of biodiversity. Assuming that extinction is an unpalatable outcome, researchers have focussed attention on how organisms can shift in their distribution to stay in the same thermal conditions or can stay in the same place by adapting to a changing thermal environment. How likely these respective outcomes might be depends on the answer to a fundamental evolutionary question, namely what genetic changes underpin adaptation to the thermal environment. The increasing access to and decreasing costs of next-generation sequencing (NGS) technologies, which can be applied to both model and non-model systems, provide a much-needed tool for understanding thermal adaptation. Here we consider broadly what is already known from non-NGS studies about thermal adaptation, then discuss the benefits and challenges of different NGS methodologies to add to this knowledge base. We then review published NGS genomics and transcriptomics studies of thermal adaptation to heat stress in metazoans and compare these results with previous non-NGS patterns. We conclude by summarising emerging patterns of genetic response and discussing future directions using these increasingly common techniques.

A second corticotropin-releasing hormone gene (CRH2) is conserved across vertebrate classes and expressed in the hindbrain of a basal neopterygian fish, the spotted gar (Lepisosteus oculatus)

To investigate the origins of the vertebrate stress-response system, we searched sequenced vertebrate genomes for genes resembling corticotropin-releasing hormone (CRH). We found that vertebrate genomes possess, in addition to CRH, another gene that resembles CRH in sequence and syntenic environment. This paralogous gene was previously identified only in the elephant shark (a holocephalan), but we find it also in marsupials, monotremes, lizards, turtles, birds, and fishes. We examined the relationship of this second vertebrate CRH gene, which we name CRH2, to CRH1 (previously known as CRH) and urocortin1/urotensin1 (UCN1/UTS1) in primitive fishes, teleosts, and tetrapods. The paralogs CRH1 and CRH2 likely evolved via duplication of CRH during a whole-genome duplication early in the vertebrate lineage. CRH2 was subsequently lost in both teleost fishes and eutherian mammals but retained in other lineages. To determine where CRH2 is expressed relative to CRH1 and UTS1, we used in situ hybridization on brain tissue from spotted gar (Lepisosteus oculatus), a neopterygian fish closely related to teleosts. In situ hybridization revealed widespread distribution of both crh1 and uts1 in the brain. Expression of crh2 was restricted to the putative secondary gustatory/secondary visceral nucleus, which also expressed calcitonin-related polypeptide alpha (calca), a marker of parabrachial nucleus in mammals. Thus, the evolutionary history of CRH2 includes restricted expression in the brain, sequence changes, and gene loss, likely reflecting release of selective constraints following whole-genome duplication. The discovery of CRH2 opens many new possibilities for understanding the diverse functions of the CRH family of peptides across vertebrates.

Keeping track of the growing number of biological functions of chitin and its interaction partners in biomedical research

Chitin is a vital polysaccharide component of protective structures in many eukaryotic organisms but seems absent in vertebrates. Chitin or chitin oligomers are therefore prime candidates for non-self-molecules, which are recognized and degraded by the vertebrate immune system. Despite the absence of polymeric chitin in vertebrates, chitinases and chitinase-like proteins (CLPs) are well conserved in vertebrate species. In many studies, these proteins have been found to be involved in immune regulation and in mediating the degradation of chitinous external protective structures of invading pathogens. Several important aspects of chitin immunostimulation have recently been uncovered, advancing our understanding of the complex regulatory mechanisms that chitin mediates. Likewise, the last few years have seen large advances in our understanding of the mechanisms and molecular interactions of chitinases and CLPs in relation to immune response regulation. It is becoming increasingly clear that their function in this context is not exclusive to chitin producing pathogens, but includes bacterial infections and cancer signaling as well. Here we provide an overview of the immune signaling properties of chitin and other closely related biomolecules. We also review the latest literature on chitinases and CLPs of the GH18 family. Finally, we examine the existing literature on zebrafish chitinases, and propose the use of zebrafish as a versatile model to complement the existing murine models. This could especially be of benefit to the exploration of the function of chitinases in infectious diseases using high-throughput approaches and pharmaceutical interventions.

European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation

The European sea bass (Dicentrarchus labrax) is a temperate-zone euryhaline teleost of prime importance for aquaculture and fisheries. This species is subdivided into two naturally hybridizing lineages, one inhabiting the north-eastern Atlantic Ocean and the other the Mediterranean and Black seas. Here we provide a high-quality chromosome-scale assembly of its genome that shows a high degree of synteny with the more highly derived teleosts. We find expansions of gene families specifically associated with ion and water regulation, highlighting adaptation to variation in salinity. We further generate a genome-wide variation map through RAD-sequencing of Atlantic and Mediterranean populations. We show that variation in local recombination rates strongly influences the genomic landscape of diversity within and differentiation between lineages. Comparing predictions of alternative demographic models to the joint allele-frequency spectrum indicates that genomic islands of differentiation between sea bass lineages were generated by varying rates of introgression across the genome following a period of geographical isolation.

The European sea bass is an economically important fish species, which is subject to intense selective breeding. Here, the authors sequence the genome of the European sea bass and highlight gene family expansions underlying adaptation to salinity change, as well as the genomic architecture of speciation between two divergent sea bass lineages.

Deep conservation of wrist and digit enhancers in fish

There is no obvious morphological counterpart of the autopod (wrist/ankle and digits) in living fishes. Comparative molecular data may provide insight into understanding both the homology of elements and the evolutionary developmental mechanisms behind the fin to limb transition. In mouse limbs the autopod is built by a "late" phase of Hoxd and Hoxa gene expression, orchestrated by a set of enhancers located at the 5' end of each cluster. Despite a detailed mechanistic understanding of mouse limb development, interpretation of Hox expression patterns and their regulation in fish has spawned multiple hypotheses as to the origin and function of "autopod" enhancers throughout evolution. Using phylogenetic footprinting, epigenetic profiling, and transgenic reporters, we have identified and functionally characterized hoxD and hoxA enhancers in the genomes of zebrafish and the spotted gar, Lepisosteus oculatus, a fish lacking the whole genome duplication of teleosts. Gar and zebrafish "autopod" enhancers drive expression in the distal portion of developing zebrafish pectoral fins, and respond to the same functional cues as their murine orthologs. Moreover, gar enhancers drive reporter gene expression in both the wrist and digits of mouse embryos in patterns that are nearly indistinguishable from their murine counterparts. These functional genomic data support the hypothesis that the distal radials of bony fish are homologous to the wrist and/or digits of tetrapods.

The circadian clock of teleost fish: a comparative analysis reveals distinct fates for duplicated genes

The circadian clock is a central oscillator that coordinates endogenous rhythms. Members of six gene families underlie the metabolic machinery of this system. Although this machinery appears to correspond to a highly conserved genetic system in metazoans, it has been recognized that vertebrates possess a more diverse gene inventory than that of non-vertebrates. This difference could have originated in the two successive rounds of whole-genome duplications that took place in the common ancestor of the group. Teleost fish underwent an extra event of whole-genome duplication, which is thought to have provided an abundance of raw genetic material for the biological innovations that facilitated the radiation of the group. In this study, we assessed the relative contributions of whole-genome duplication and small-scale gene duplication to generate the repertoire of genes associated with the circadian clock of teleost fish. To achieve this goal, we annotated genes from six gene families associated with the circadian clock in eight teleost fish species, and we reconstructed their evolutionary history by inferring phylogenetic relationships. Our comparative analysis indicated that teleost species possess a variable repertoire of genes related to the circadian clock gene families and that the actual diversity of these genes has been shaped by a variety of phenomena, such as the complete deletion of ohnologs, the differential retention of genes, and lineage-specific gene duplications. From a functional perspective, the subfunctionalization of two ohnolog genes (PER1a and PER1b) in zebrafish highlights the power of whole-genome duplications to generate biological diversity.

New insights into the evolution of vertebrate CRH (corticotropin-releasing hormone) and invertebrate DH44 (diuretic hormone 44) receptors in metazoans

The corticotropin releasing hormone receptors (CRHR) and the arthropod diuretic hormone 44 receptors (DH44R) are structurally and functionally related members of the G protein-coupled receptors (GPCR) of the secretin-like receptor superfamily. We show here that they derive from a bilaterian predecessor. In protostomes, the receptor became DH44R that has been identified and functionally characterised in several arthropods but the gene seems to be absent from nematode genomes. Duplicate DH44R genes (DH44 R1 and DH44R2) have been described in some arthropods resulting from lineage-specific duplications. Recently, CRHR-DH44R-like receptors have been identified in the genomes of some lophotrochozoans (molluscs, which have a lineage-specific gene duplication, and annelids) as well as representatives of early diverging deuterostomes. Vertebrates have previously been reported to have two CRHR receptors that were named CRHR1 and CRHR2. To resolve their origin we have analysed recently assembled genomes from representatives of early vertebrate divergencies including elephant shark, spotted gar and coelacanth. We show here by analysis of synteny conservation that the two CRHR genes arose from a common ancestral gene in the early vertebrate tetraploidizations (2R) approximately 500 million years ago. Subsequently, the teleost-specific tetraploidization (3R) resulted in a duplicate of CRHR1 that has been lost in some teleost lineages. These results help distinguish orthology and paralogy relationships and will allow studies of functional conservation and changes during evolution of the individual members of the receptor family and their multiple native peptide agonists.

Evolution of endothelin receptors in vertebrates

Endothelin receptors are G protein coupled receptors (GPCRs) of the β-group of rhodopsin receptors that bind to endothelin ligands, which are 21 amino acid long peptides derived from longer prepro-endothelin precursors. The most basal Ednr-like GPCR is found outside vertebrates in the cephalochordate amphioxus, but endothelin ligands are only present among vertebrates, including the lineages of jawless vertebrates (lampreys and hagfishes), cartilaginous vertebrates (sharks, rays, and chimaeras), and bony vertebrates (ray-finned fishes and lobe-finned vertebrates including tetrapods). A bona fide endothelin system is thus a vertebrate-specific innovation with important roles for regulating the cardiovascular system, renal and pulmonary processes, as well as for the development of the vertebrate-specific neural crest cell population and its derivatives. Expectedly, dysregulation of endothelin receptors and the endothelin system leads to a multitude of human diseases. Despite the importance of different types of endothelin receptors for vertebrate development and physiology, current knowledge on endothelin ligand-receptor interactions, on the expression of endothelin receptors and their ligands, and on the functional roles of the endothelin system for embryonic development and in adult vertebrates is very much biased towards amniote vertebrates. Recent analyses from a variety of vertebrate lineages, however, have shown that the endothelin system in lineages such as teleost fish and lampreys is more diverse and is divergent from the mammalian endothelin system. This diversity is mainly based on differential evolution of numerous endothelin system components among vertebrate lineages generated by two rounds of whole genome duplication (three in teleosts) during vertebrate evolution. Here we review current understanding of the evolutionary history of the endothelin receptor family in vertebrates supplemented with surveys on the endothelin receptor gene complement of newly available genome assemblies from phylogenetically informative taxa. Our assessment further highlights the diversity of the vertebrate endothelin system and calls for detailed functional and pharmacological analyses of the endothelin system beyond tetrapods.

Evolution of gonadotropin-inhibitory hormone receptor and its ligand

Gonadotropin-inhibitory hormone (GnIH) is a neuropeptide inhibitor of gonadotropin secretion, which was first identified in the Japanese quail hypothalamus. GnIH peptides share a C-terminal LPXRFamide (X=L or Q) motif in most vertebrates. The receptor for GnIH (GnIHR) is the seven-transmembrane G protein-coupled receptor 147 (GPR147) that inhibits cAMP production. GPR147 is also named neuropeptide FF (NPFF) receptor 1 (NPFFR1), because it also binds NPFF that has a C-terminal PQRFamide motif. To understand the evolutionary history of the GnIH system in the animal kingdom, we searched for receptors structurally similar to GnIHR in the genome of six mammals (human, mouse, rat, cattle, cat, and rabbit), five birds (pigeon, chicken, turkey, budgerigar, and zebra finch), one reptile (green anole), one amphibian (Western clawed flog), six fishes (zebrafish, Nile tilapia, Fugu, coelacanth, spotted gar, and lamprey), one hemichordate (acorn worm), one echinoderm (purple sea urchin), one mollusk (California sea hare), seven insects (pea aphid, African malaria mosquito, honey bee, buff-tailed bumblebee, fruit fly, jewel wasp, and red flour beetle), one cnidarian (hydra), and constructed phylogenetic trees by neighbor joining (NJ) and maximum likelihood (ML) methods. A multiple sequence alignment of the receptors showed highly conserved seven-transmembrane domains as well as disulfide bridge sites between the first and second extracellular loops, including the receptor of hydra. Both NJ and ML analyses grouped the receptors of vertebrates into NPFFR1 and NPFFR2 (GPR74), and the receptors of insects into the receptor for SIFamide peptides that share a C-terminal YRKPPFNGSIFamide motif. Although human, quail and zebrafish GnIHR (NPFFR1) were most structurally similar to SIFamide receptor of fruit fly in the Famide peptide (FMRFamide, neuropeptide F, short neuropeptide F, drosulfakinin, myosuppressin, SIFamide) receptor families, the amino acid sequences and the peptide coding regions of GnIH precursors were most similar to FMRFamide precursor of fruit fly in the precursors of Famide peptide families. Chromosome synteny analysis of the precursor genes of human, quail and zebrafish GnIH and fruit fly Famide peptides further identified conserved synteny in vertebrate GnIH and fruit fly FMRFa precursor genes as well as other Famide peptide precursor genes. These results suggest that GnIH and its receptor pair and SIFamide and its receptor pair may have diverged and co-evolved independently in vertebrates and insects, respectively, from their ancestral Famide peptide and its receptor pair, during diversification and evolution of deuterostomian and protostomian species.

The lineage-specific evolution of aquaporin gene clusters facilitated tetrapod terrestrial adaptation

A major physiological barrier for aquatic organisms adapting to terrestrial life is dessication in the aerial environment. This barrier was nevertheless overcome by the Devonian ancestors of extant Tetrapoda, but the origin of specific molecular mechanisms that solved this water problem remains largely unknown. Here we show that an ancient aquaporin gene cluster evolved specifically in the sarcopterygian lineage, and subsequently diverged into paralogous forms of AQP2, -5, or -6 to mediate water conservation in extant Tetrapoda. To determine the origin of these apomorphic genomic traits, we combined aquaporin sequencing from jawless and jawed vertebrates with broad taxon assembly of >2,000 transcripts amongst 131 deuterostome genomes and developed a model based upon Bayesian inference that traces their convergent roots to stem subfamilies in basal Metazoa and Prokaryota. This approach uncovered an unexpected diversity of aquaporins in every lineage investigated, and revealed that the vertebrate superfamily consists of 17 classes of aquaporins (Aqp0 - Aqp16). The oldest orthologs associated with water conservation in modern Tetrapoda are traced to a cluster of three aqp2-like genes in Actinistia that likely arose >500 Ma through duplication of an aqp0-like gene present in a jawless ancestor. In sea lamprey, we show that aqp0 first arose in a protocluster comprised of a novel aqp14 paralog and a fused aqp01 gene. To corroborate these findings, we conducted phylogenetic analyses of five syntenic nuclear receptor subfamilies, which, together with observations of extensive genome rearrangements, support the coincident loss of ancestral aqp2-like orthologs in Actinopterygii. We thus conclude that the divergence of sarcopterygian-specific aquaporin gene clusters was permissive for the evolution of water conservation mechanisms that facilitated tetrapod terrestrial adaptation.

First-generation linkage map for the European tree frog (Hyla arborea) with utility in congeneric species

Background

Western Palearctic tree frogs (Hyla arborea group) represent a strong potential for evolutionary and conservation genetic research, so far underexploited due to limited molecular resources. New microsatellite markers have recently been developed for Hyla arborea, with high cross-species utility across the entire circum-Mediterranean radiation. Here we conduct sibship analyses to map available markers for use in future population genetic applications.

Findings

We characterized eight linkage groups, including one sex-linked, all showing drastically reduced recombination in males compared to females, as previously documented in this species. Mapping of the new 15 markers to the ~200 My diverged Xenopus tropicalis genome suggests a generally conserved synteny with only one confirmed major chromosome rearrangement.

Conclusions

The new microsatellites are representative of several chromosomes of H. arborea that are likely to be conserved across closely-related species. Our linkage map provides an important resource for genetic research in European Hylids, notably for studies of speciation, genome evolution and conservation.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-850) contains supplementary material, which is available to authorized users.

Phylogeny of zebrafish, a "model species," within Danio, a "model genus"

Zebrafish (Danio rerio) is an important model for vertebrate development, genomics, physiology, behavior, toxicology, and disease. Additionally, work on numerous Danio species is elucidating evolutionary mechanisms for morphological development. Yet, the relationships of zebrafish and its closest relatives remain unclear possibly due to incomplete lineage sorting, speciation with gene flow, and interspecies hybridization. To clarify these relationships, we first constructed phylogenomic data sets from 30,801 restriction-associated DNA (RAD)-tag loci (483,026 variable positions) with clear orthology to a single location in the sequenced zebrafish genome. We then inferred a well-supported species tree for Danio and tested for gene flow during the diversification of the genus. An approach independent of the sequenced zebrafish genome verified all inferred relationships. Although identification of the sister taxon to zebrafish has been contentious, multiple RAD-tag data sets and several analytical methods provided strong evidence for Danio aesculapii as the most closely related extant zebrafish relative studied to date. Data also displayed patterns consistent with gene flow during speciation and postspeciation introgression in the lineage leading to zebrafish. The incorporation of biogeographic data with phylogenomic analyses put these relationships in a phylogeographic context and supplied additional support for D. aesculapii as the sister species to D. rerio. The clear resolution of this study establishes a framework for investigating the evolutionary biology of Danio and the heterogeneity of genome evolution in the recent history of a model organism within an emerging model genus for genetics, development, and evolution.

Multiple thyrotropin β-subunit and thyrotropin receptor-related genes arose during vertebrate evolution

Thyroid-stimulating hormone (TSH) is composed of a specific β subunit and an α subunit that is shared with the two pituitary gonadotropins. The three β subunits derive from a common ancestral gene through two genome duplications (1R and 2R) that took place before the radiation of vertebrates. Analysis of genomic data from phylogenetically relevant species allowed us to identify an additional Tshβ subunit-related gene that was generated through 2R. This gene, named Tshβ2, present in cartilaginous fish, little skate and elephant shark, and in early lobe-finned fish, coelacanth and lungfish, was lost in ray-finned fish and tetrapods. The absence of a second type of TSH receptor (Tshr) gene in these species suggests that both TSHs act through the same receptor. A novel Tshβ sister gene, named Tshβ3, was generated through the third genomic duplication (3R) that occurred early in the teleost lineage. Tshβ3 is present in most teleost groups but was lostin tedraodontiforms. The 3R also generated a second Tshr, named Tshrb. Interestingly, the new Tshrb was translocated from its original chromosomic position after the emergence of eels and was then maintained in its new position. Tshrb was lost in tetraodontiforms and in ostariophysians including zebrafish although the latter species have two TSHs, suggesting that TSHRb may be dispensable. The tissue distribution of duplicated Tshβs and Tshrs was studied in the European eel. The endocrine thyrotropic function in the eel would be essentially mediated by the classical Tshβ and Tshra, which are mainly expressed in the pituitary and thyroid, respectively. Tshβ3 and Tshrb showed a similar distribution pattern in the brain, pituitary, ovary and adipose tissue, suggesting a possible paracrine/autocrine mode of action in these non-thyroidal tissues. Further studies will be needed to determine the binding specificity of the two receptors and how these two TSH systems are interrelated.

Daily and seasonal expression of clock genes in the pituitary of the European sea bass (Dicentrarchus labrax)

The expression of select clock genes (clock, bmal, per1, per2, cry1, cry2) was investigated throughout the day and across the four seasons for two consecutive years in the pituitary of adult sea bass (Dicentrarchus labrax). A rhythmic pattern of daily expression was consistently observed in summer and autumn, while arrhythmicity was observed for some clock genes during spring and winter, concomitant with low water temperatures. The expression of clock and bmal showed highest values at the end of the day and during the night, while that of per and cry was mostly antiphasic, with high values during the day. Melatonin affects clock-gene expression in the pituitary of mammals. We therefore sought to test the effect of melatonin on clock-gene expression in the pituitary of sea bass both in vivo and in vitro. Melatonin modestly affected the expression of some clock genes (in particular cry genes) when added to the fish diet or the culture medium of pituitary glands. Our data show that clock genes display rhythmic daily expression in the pituitary of adult sea bass, which are profoundly modified according to the season. We suggest that the effect of photoperiod on clock gene expression may be mediated, at least in part, by melatonin, and that temperature may have a key role adjusting seasonal variations.

Dynamic evolution of the GnRH receptor gene family in vertebrates

BACKGROUND

Elucidating the mechanisms underlying coevolution of ligands and receptors is an important challenge in molecular evolutionary biology. Peptide hormones and their receptors are excellent models for such efforts, given the relative ease of examining evolutionary changes in genes encoding for both molecules. Most vertebrates possess multiple genes for both the decapeptide gonadotropin releasing hormone (GnRH) and for the GnRH receptor. The evolutionary history of the receptor family, including ancestral copy number and timing of duplications and deletions, has been the subject of controversy.

RESULTS

We report here for the first time sequences of three distinct GnRH receptor genes in salamanders (axolotls, Ambystoma mexicanum), which are orthologous to three GnRH receptors from ranid frogs. To understand the origin of these genes within the larger evolutionary context of the gene family, we performed phylogenetic analyses and probabilistic protein homology searches of GnRH receptor genes in vertebrates and their near relatives. Our analyses revealed four points that alter previous views about the evolution of the GnRH receptor gene family. First, the "mammalian" pituitary type GnRH receptor, which is the sole GnRH receptor in humans and previously presumed to be highly derived because it lacks the cytoplasmic C-terminal domain typical of most G-protein coupled receptors, is actually an ancient gene that originated in the common ancestor of jawed vertebrates (Gnathostomata). Second, unlike previous studies, we classify vertebrate GnRH receptors into five subfamilies. Third, the order of subfamily origins is the inverse of previous proposed models. Fourth, the number of GnRH receptor genes has been dynamic in vertebrates and their ancestors, with multiple duplications and losses.

CONCLUSION

Our results provide a novel evolutionary framework for generating hypotheses concerning the functional importance of structural characteristics of vertebrate GnRH receptors. We show that five subfamilies of vertebrate GnRH receptors evolved early in the vertebrate phylogeny, followed by several independent instances of gene loss. Chief among cases of gene loss are humans, best described as degenerate with respect to GnRH receptors because we retain only a single, ancient gene.

Unexpected multiplicity of QRFP receptors in early vertebrate evolution

The neuropeptide QRFP, also called 26RFa, and its G protein-coupled receptor GPR103 have been identified in all vertebrates investigated. In mammals, this peptide-receptor pair has been found to have several effects including stimulation of appetite. Recently, we reported that a QRFP peptide is present in amphioxus, Branchiostoma floridae, and we also identified a QRFP receptor (QRFPR) that mediates a functional response to sub-nanomolar concentrations of the amphioxus peptide as well as short and long human QRFP (Xu et al., submitted). Because the ancestral vertebrate underwent two tetraploidizations, it might be expected that duplicates of the QRFP gene and its receptor gene may exist. Indeed, we report here the identification of multiple vertebrate QRFPR genes. Three QRFPR genes are present in the coelacanth Latimeria chalumnae, representing an early diverging sarcopterygian lineage. Three QRFPR genes are present in the basal actinopterygian fish, the spotted gar. Phylogenetic and chromosomal analyses show that only two of these receptor genes are orthologous between the two species, thus demonstrating a total of four distinct vertebrate genes. Three of the QRFPR genes resulted from the early vertebrate tetraploidizations and were copied along with syntenic neuropeptide Y receptor genes. The fourth QRFPR gene may be an even older and distinct lineage. Because mammals and birds have only a single QRFPR gene, this means that three genes have been lost in these lineages, and at least one of these was lost independently in mammals and birds because it is still present in a turtle. In conclusion, these results show that the QRFP system gained considerable complexity in the early stages of vertebrate evolution and still maintains much of this in some lineages, and that it has been secondarily reduced in mammals.