The "fish-specific" Hox cluster duplication is coincident with the origin of teleosts

Fish Evo-Devo: Moving Toward Species-Specific and Knowledge-Based Interactome

A knowledge-based interactome maps interactions among proteins and molecules within a cell using experimental data, computational predictions, and literature mining. These interactomes are vital for understanding cellular functions, pathways, and the evolutionary conservation of protein interactions. They reveal how interactions regulate growth, differentiation, and development. Transitioning to functionally validated interactomes is crucial in evolutionary developmental biology (Evo-Devo), especially for non-model species, to uncover unique regulatory networks, evolutionary novelties, and reliable gene interaction models. This enhances our understanding of complex trait evolution across species. The European Evo-Devo 2024 conference in Helsinki hosted the first fish-specific Evo-Devo symposium, highlighting the growing interest in fish models. Advances in genome annotation, genome editing, imaging, and molecular screening are expanding fish Evo-Devo research. High-throughput molecular data have enabled the deduction of gene regulatory networks. The next steps involve creating species-specific interactomes, validating them functionally, and integrating additional molecular data to deepen the understanding of complex regulatory interactions in fish Evo-Devo. This short review aims to address the logical steps for this transition, as well as the necessities and limitations of this journey.

Zebrafish: unraveling genetic complexity through duplicated genes

The zebrafish is an invaluable model organism for genetic, developmental, and disease research. Although its high conservation with humans is often cited as justification for its use, the zebrafish harbors oft-ignored genetic characteristics that may provide unique insights into gene structure and function. Zebrafish, along with other teleost fish, underwent an additional round of whole genome duplication after their split from tetrapods-resulting in an abundance of duplicated genes when compared to other vertebrates. These duplicated genes have evolved in distinct ways over the ensuing 350 million years. Thus, each gene within a duplicated gene pair has nuanced differences that create a unique identity. By investigating both members of the gene pair together, we can elucidate the mechanisms that underly protein structure and function and drive the complex interplay within biological systems, such as signal transduction cascades, genetic regulatory networks, and evolution of tissue and organ function. It is crucial to leverage such studies to explore these molecular dynamics, which could have far-reaching implications for both basic science and therapeutic development. Here, we will review the role of gene duplications and the existing models for gene divergence and retention following these events. We will also highlight examples within each of these models where studies comparing duplicated genes in the zebrafish have yielded key insights into protein structure, function, and regulation.

Divergence time estimates for the hypoxia-inducible factor-1 alpha (HIF1α) reveal an ancient emergence of animals in low-oxygen environments

Unveiling the tempo and mode of animal evolution is necessary to understand the links between environmental changes and biological innovation. Although the earliest unambiguous metazoan fossils date to the late Ediacaran period, molecular clock estimates agree that the last common ancestor (LCA) of all extant animals emerged ~850 Ma, in the Tonian period, before the oldest evidence for widespread ocean oxygenation at ~635-560 Ma in the Ediacaran period. Metazoans are aerobic organisms, that is, they are dependent on oxygen to survive. In low-oxygen conditions, most animals have an evolutionarily conserved pathway for maintaining oxygen homeostasis that triggers physiological changes in gene expression via the hypoxia-inducible factor (HIFa). However, here we confirm the absence of the characteristic HIFa protein domain responsible for the oxygen sensing of HIFa in sponges and ctenophores, indicating the LCA of metazoans lacked the functional protein domain as well, and so could have maintained their transcription levels unaltered under the very low-oxygen concentrations of their environments. Using Bayesian relaxed molecular clock dating, we inferred that the ancestral gene lineage responsible for HIFa arose in the Mesoproterozoic Era, ~1273 Ma (Credibility Interval 957-1621 Ma), consistent with the idea that important genetic machinery associated with animals evolved much earlier than the LCA of animals. Our data suggest at least two duplication events in the evolutionary history of HIFa, which generated three vertebrate paralogs, products of the two successive whole-genome duplications that occurred in the vertebrate LCA. Overall, our results support the hypothesis of a pre-Tonian emergence of metazoans under low-oxygen conditions, and an increase in oxygen response elements during animal evolution.

The Impact of Whole Genome Duplication on the Evolution of the Arachnids

The proliferation of genomic resources for Chelicerata in the past 10 years has revealed that the evolution of chelicerate genomes is more dynamic than previously thought, with multiple waves of ancient whole genome duplications affecting separate lineages. Such duplication events are fascinating from the perspective of evolutionary history because the burst of new gene copies associated with genome duplications facilitates the acquisition of new gene functions (neofunctionalization), which may in turn lead to morphological novelties and spur net diversification. While neofunctionalization has been invoked in several contexts with respect to the success and diversity of spiders, the overall impact of whole genome duplications on chelicerate evolution and development remains imperfectly understood. The purpose of this review is to examine critically the role of whole genome duplication on the diversification of the extant arachnid orders, as well as assess functional datasets for evidence of subfunctionalization or neofunctionalization in chelicerates. This examination focuses on functional data from two focal model taxa: the spider Parasteatoda tepidariorum, which exhibits evidence for an ancient duplication, and the harvestman Phalangium opilio, which exhibits an unduplicated genome. I show that there is no evidence that taxa with genome duplications are more successful than taxa with unduplicated genomes. I contend that evidence for sub- or neofunctionalization of duplicated developmental patterning genes in spiders is indirect or fragmentary at present, despite the appeal of this postulate for explaining the success of groups like spiders. Available expression data suggest that the condition of duplicated Hox modules may have played a role in promoting body plan disparity in the posterior tagma of some orders, such as spiders and scorpions, but functional data substantiating this postulate are critically missing. Spatiotemporal dynamics of duplicated transcription factors in spiders may represent cases of developmental system drift, rather than neofunctionalization. Developmental system drift may represent an important, but overlooked, null hypothesis for studies of paralogs in chelicerate developmental biology. To distinguish between subfunctionalization, neofunctionalization, and developmental system drift, concomitant establishment of comparative functional datasets from taxa exhibiting the genome duplication, as well as those that lack the paralogy, is sorely needed.

Interlocus Gene Conversion, Natural Selection, and Paralog Homogenization

Following a duplication, the resulting paralogs tend to diverge. While mutation and natural selection can accelerate this process, they can also slow it. Here, we quantify the paralog homogenization that is caused by point mutations and interlocus gene conversion (IGC). Among 164 duplicated teleost genes, the median percentage of postduplication codon substitutions that arise from IGC rather than point mutation is estimated to be between 7% and 8%. By differentiating between the nonsynonymous codon substitutions that homogenize the protein sequences of paralogs and the nonhomogenizing nonsynonymous substitutions, we estimate the homogenizing nonsynonymous rates to be higher for 163 of the 164 teleost data sets as well as for all 14 data sets of duplicated yeast ribosomal protein-coding genes that we consider. For all 14 yeast data sets, the estimated homogenizing nonsynonymous rates exceed the synonymous rates.

Enrichment of intrinsically disordered residues in ohnologs facilitates abiotic stress resilience in Brassica rapa

Arabidopsis thaliana and Brassica rapa are in the same evolutionary lineage, although the latter experienced an additional whole genome triplication event. Therefore, it would be intriguing to investigate the traits that gene duplication imposes to mediate plant stress tolerance. Here, we noticed that B. rapa abiotic stress resistance (ASR) genes which code at least one stress responsive domain have a significantly higher number of paralogs than A. thaliana. Analysing the disordered content of the ASR genes in both species, we found that intrinsically disordered residues (IDR) are specifically enriched in whole genome duplication (WGD) derived paralogs. Subsequently, domain similarity analysis between WGD pairs of both species has revealed that majority of WGD pairs in B. rapa did not share domains with each other. Furthermore, domain enrichment analysis has shown that B. rapa paralogs contain 36 distinct stress responsive enriched domains, significantly higher than A. thaliana paralogs. Next, we performed MSA to investigate the domain conservation between orthologs and ohnologs pairs, we found that 80.13% of B. rapa ohnologs acquire new domains, depicting the fact that ohnologs play a significant role in stress-related behaviours. The average IDR content of the ohnologs enriching new domains after gene duplication in B. rapa (0.19), is also significantly higher than A. thaliana (0.04). Interestingly, we also found that all of these attributes i.e., exhibiting higher number of WGD paralogs and enhancement of IDR in ASR genes of B. rapa compared to A. thaliana is exclusive for ASR genes only. No such significant differences were observed in randomly selected non-ASR genes between the two species. Together these results provide strong support for the hypothesis that augmentation of IDR content followed by a whole genome duplication event imposes the stress resistance potentiality in B. rapa. This research will shed light on the mechanism of how B. rapa is able to successfully adapt to stress over the evolutionary timescale.

Genomic Characterization of hox Genes in Senegalese Sole (Solea senegalensis, Kaup 1858): Clues to Evolutionary Path in Pleuronectiformes

The Senegalese sole (Solea senegalensis, Kaup 1858), a marine flatfish, belongs to the Pleuronectiformes order. It is a commercially important species for fisheries and aquaculture. However, in aquaculture, several production bottlenecks have still to be resolved, including skeletal deformities and high mortality during the larval and juvenile phase. The study aims to characterize the hox gene clusters in S. senegalensis to understand better the developmental and metamorphosis process in this species. Using a BAC library, the clones that contain hox genes were isolated, sequenced by NGS and used as BAC-FISH probes. Subsequently the hox clusters were studied by sequence analysis, comparative genomics, and cytogenetic and phylogenetic analysis. Cytogenetic analysis demonstrated the localization of four BAC clones on chromosome pairs 4, 12, 13, and 16 of the Senegalese sole cytogenomic map. Comparative and phylogenetic analysis showed a highly conserved organization in each cluster and different phylogenetic clustering in each hox cluster. Analysis of structural and repetitive sequences revealed accumulations of polymorphisms mediated by repetitive elements in the hoxba cluster, mainly retroelements. Therefore, a possible loss of the hoxb7a gene can be established in the Pleuronectiformes lineage. This work allows the organization and regulation of hox clusters to be understood, and is a good base for further studies of expression patterns.

Might Gene Duplication and Neofunctionalization Contribute to the Sexual Lability Observed in Fish?

Sex determination and differentiation varies widely across vertebrates, but is most dramatically diverse in fishes. Among fishes sex reversal and sex change are observed in 41 teleost families spanning 7 orders. These sex-changing fish perhaps highlight better than any other system that sex determination is not the narrow and fixed construct we once thought, but a plastic trait that is better viewed as a reaction norm. However, while this stunning transformation is increasingly understood, a fundamental question arises, which is why some fish species have retained this inherent plasticity in sexual fate, while others have not? Here, we explore our current understanding of sex change in fish, some of the factors that permit and constrain sex reversal, and posit that gene duplication and neofunctionalization contribute to the sexual lability observed in fish.

Genotyping single nucleotide polymorphisms and inferring ploidy by amplicon sequencing for polyploid, ploidy-variable organisms

Whole genome duplication is hypothesized to have played a critical role in the evolution of several major taxa, including vertebrates, and while many lineages have rediploidized, some retain polyploid genomes. Additionally, variation in ploidy can occur naturally or be artificially induced within select plant and animal species. Modern genetic techniques have not been widely applied to polyploid or ploidy-variable species, in part due to the difficulty of obtaining genotype data from polyploids. In this study, we demonstrate a strategy for developing an amplicon sequencing panel of single nucleotide polymorphisms for high-throughput genotyping of polyploid organisms. We then develop a method to infer ploidy of individuals from amplicon sequencing data that is generalized to apply to any ploidy and does not require prior identification of heterozygous genotypes. Combining these two techniques will allow researchers to both infer ploidy and generate ploidy-aware genotypes with the same amplicon sequencing panel. We demonstrate this approach with white sturgeon Acipenser transmontanus, a ploidy-variable (octoploid, decaploid and dodecaploid) imperiled species under conservation management in the Pacific Northwest and obtained a panel of 325 loci. These loci were validated by examining inheritance in known-cross families, and the ploidy inference method was validated with known ploidy samples. We provide scripts that adapt existing pipelines to genotype polyploids and an R package for application of the ploidy inference method. We expect that these techniques will empower studies of genetic variation and inheritance in polyploid organisms that vary in ploidy level, either naturally or as a result of artificial propagation practices.

Genomic Survey of Tyrosine Kinases Repertoire in Electrophorus electricus With an Emphasis on Evolutionary Conservation and Diversification

Tyrosine kinases (TKs) play key roles in the regulation of multicellularity in organisms and involved primarily in cell growth, differentiation, and cell-to-cell communication. Genome-wide characterization of TKs has been conducted in many metazoans; however, systematic information regarding this superfamily in Electrophorus electricus (electric eel) is still lacking. In this study, we identified 114 TK genes in the E electricus genome and investigated their evolution, molecular features, and domain architecture using phylogenetic profiling to gain a better understanding of their similarities and specificity. Our results suggested that the electric eel TK (EeTK) repertoire was shaped by whole-genome duplications (WGDs) and tandem duplication events. Compared with other vertebrate TKs, gene members in Jak, Src, and EGFR subfamily duplicated specifically, but with members lost in Eph, Axl, and Ack subfamily in electric eel. We also conducted an exhaustive survey of TK genes in genomic databases, identifying 1674 TK proteins in 31 representative species covering all the main metazoan lineages. Extensive evolutionary analysis indicated that TK repertoire in vertebrates tended to be remarkably conserved, but the gene members in each subfamily were very variable. Comparative expression profile analysis showed that electric organ tissues and muscle shared a similar pattern with specific highly expressed TKs (ie, epha7, musk, jak1, and pdgfra), suggesting that regulation of TKs might play an important role in specifying an electric organ identity from its muscle precursor. We further identified TK genes exhibiting tissue-specific expression patterns, indicating that members in TKs participated in subfunctionalization representing an evolutionary divergence required for the performance of different tissues. This work generates valuable information for further gene function analysis and identifying candidate TK genes reflecting their unique tissue-function specializations in electric eel.

The lasting after-effects of an ancient polyploidy on the genomes of teleosts

The ancestor of most teleost fishes underwent a whole-genome duplication event three hundred million years ago. Despite its antiquity, the effects of this event are evident both in the structure of teleost genomes and in how the surviving duplicated genes still operate to drive form and function. I inferred a set of shared syntenic regions that survive from the teleost genome duplication (TGD) using eight teleost genomes and the outgroup gar genome (which lacks the TGD). I then phylogenetically modeled the TGD's resolution via shared and independent gene losses and applied a new simulation-based statistical test for the presence of bias toward the preservation of genes from one parental subgenome. On the basis of that test, I argue that the TGD was likely an allopolyploidy. I find that duplicate genes surviving from this duplication in zebrafish are less likely to function in early embryo development than are genes that have returned to single copy at some point in this species' history. The tissues these ohnologs are expressed in, as well as their biological functions, lend support to recent suggestions that the TGD was the source of a morphological innovation in the structure of the teleost retina. Surviving duplicates also appear less likely to be essential than singletons, despite the fact that their single-copy orthologs in mouse are no less essential than other genes.

Three-dimensional characterization of osteocyte volumes at multiple scales, and its relationship with bone biology and genome evolution in ray-finned fishes

Osteocytes, cells embedded within the bone mineral matrix, inform on key aspects of vertebrate biology. In particular, a relationship between volumes of the osteocytes and bone growth and/or genome size has been proposed for several tetrapod lineages. However, the variation in osteocyte volume across different scales is poorly characterized and mostly relies on incomplete, two-dimensional information. In this study, we characterize the variation of osteocyte volumes in ray-finned fishes (Actinopterygii), a clade including more than half of modern vertebrate species in which osteocyte biology is poorly known. We use X-ray synchrotron micro-computed tomography (SRµCT) to achieve a three-dimensional visualization of osteocyte lacunae and direct measurement of their size (volumes). Our specimen sample is designed to characterize variation in osteocyte lacuna morphology at three scales: within a bone, among the bones of one individual and among species. At the intra-bone scale, we find that osteocyte lacunae vary noticeably in size between zones of organized and woven bone (being up to six times larger in woven bone), and across cyclical bone deposition. This is probably explained by differences in bone deposition rate, with larger osteocyte lacunae contained in bone that deposits faster. Osteocyte lacuna volumes vary 3.5-fold among the bones of an individual, and this cannot readily be explained by variation in bone growth rate or other currently observable factors. Finally, we find that genome size provides the best explanation of variation in osteocyte lacuna volume among species: actinopterygian taxa with larger genomes (polyploid taxa in particular) have larger osteocyte lacunae (with a ninefold variation in median osteocyte volume being measured). Our findings corroborate previous two-dimensional studies in tetrapods that also observed similar patterns of intra-individual variation and found a correlation with genome size. This opens new perspectives for further studies on bone evolution, physiology and palaeogenomics in actinopterygians, and vertebrates as a whole.

Phylogenetic perspective on the relationships and evolutionary history of the Acipenseriformes

The Acipenseriformes, as one of the earliest extant vertebrates, plays an important role in the evolution of fishes and even the whole vertebrates. Here we collected and analyzed all complete mitochondrial genomes of Acipenseriformes species. Phylogenetic analyses demonstrated that the polytomous branch included Acipenseridae and Polyodontidae formed five clades. The Polyodontidae clade and the Scaphirhynchus clade both were monophyletic group, whereas the Acipenser species and the Huso species both were polyphyletic group. The Bayesian divergence times showed that the origin time for Acipenseriformes was at 318.0 Mya, which was similar to the some previous results of 312.1 Mya, 346.9 Mya and 389.7 Mya. The result was in good consistent with the paleontological data available and the split time of the Pacific and Atlantic Oceans from the Jurassic to the Cretaceous (Laurasia splits in North America and Eurasia). The dN/dS ratios showed the evolutionary rates gradually slow down in five major Acipenseriformes clades from the Clade A (the Pacific sturgeons species) to Clade C (the genus Scaphirhynchus), which was related to the process of geographical formation.

Rapid genomic DNA variation in newly hybridized carp lineages derived from Cyprinus carpio (♀) × Megalobrama amblycephala (♂)

Background

Distant hybridization can generate changes in phenotypes and genotypes that lead to the formation of new hybrid lineages with genetic variation. In this study, the establishment of two bisexual fertile carp lineages, including the improved diploid common carp (IDC) lineage and the improved diploid scattered mirror carp (IDMC) lineage, from the interspecific hybridization of common carp (Cyprinus carpio, 2n = 100) (♀) × blunt snout bream (Megalobrama amblycephala, 2n = 48) (♂), provided a good platform to investigate the genetic relationship between the parents and their hybrid progenies.

Result

In this study, we investigated the genetic variation of 12 Hox genes in the two types of improved carp lineages derived from common carp (♀) × blunt snout bream (♂). Hox gene clusters were abundant in the first generation of IDC, but most were not stably inherited in the second generation. In contrast, we did not find obvious mutations in Hox genes in the first generation of IDMC, and almost all the Hox gene clusters were stably inherited from the first generation to the second generation of IDMC. Interestingly, we found obvious recombinant clusters of Hox genes in both improved carp lineages, and partially recombinant clusters of Hox genes were stably inherited from the first generation to the second generation in both types of improved carp lineages. On the other hand, some Hox genes were gradually becoming pseudogenes, and some genes were completely pseudogenised in IDC or IDMC.

Conclusions

Our results provided important evidence that distant hybridization produces rapid genomic DNA changes that may or may not be stably inherited, providing novel insights into the function of hybridization in the establishment of improved lineages used as new fish resources for aquaculture.

Genome-wide identification of nonvisual opsin family reveals amplification of RPE-retinal G protein receptor gene (RGR) and offers novel insights into functions of RGR(s) in Paralichthys olivaceus (Paralichthyidae, Teleostei)

Opsins play important roles in the image-forming visual pathways and numerous biological systems such as the biological clock and circadian rhythm. However, the nonvisual opsins involved in nonimage forming process are not clear to date. The aim of this study was to characterize nonvisual opsins in Paralichthys olivaceus. A total of 24 nonvisual opsin genes were identified. Expressions of these genes in eye, brain, heart, testis, and fin were investigated by quantitative real-time polymerase chain reaction (qRT-PCR). Testis contained a surprisingly large number of nonvisual opsins including Opn4m2, Tmt2a, Tmt3b, Opn3, RRH, Opn7a, and Opn7b. Syntenic and phylogenetic analyses confirmed that the RGRa and RGRb originated from the teleost-specific genome duplication (TSGD). qRT-PCR results demonstrated high RGRa and RGRb expression in the eye, while the expression levels in the brain, heart, testis, and fin were relatively weak. In situ hybridization results presented here revealed the presence of both RGRa and RGRb mRNA-positive signals in the ganglion cell layer but absence in the intracellular compartment of retinal pigment epithelium (RPE) and Müller glial cells. Therefore, we hypothesized that RGRa and RGRb had undergone subfunctionalization in P. olivaceus after TSGD. In conclusion, this study provides novel insights into the evolutionary fates of the RGR genes, still, further studies need to be done to explore the mechanism about the lack of RGR genes' expression in RPE.

The Molecular Evolution of Circadian Clock Genes in Spotted Gar (Lepisosteus oculatus)

Circadian rhythms are biological rhythms with a period of approximately 24 h. While canonical circadian clock genes and their regulatory mechanisms appear highly conserved, the evolution of clock gene families is still unclear due to several rounds of whole genome duplication in vertebrates. The spotted gar (Lepisosteus oculatus), as a non-teleost ray-finned fish, represents a fish lineage that diverged before the teleost genome duplication (TGD), providing an outgroup for exploring the evolutionary mechanisms of circadian clocks after whole-genome duplication. In this study, we interrogated the spotted gar draft genome sequences and found that spotted gar contains 26 circadian clock genes from 11 families. Phylogenetic analysis showed that 9 of these 11 spotted gar circadian clock gene families have the same number of genes as humans, while the members of the nfil3 and cry families are different between spotted gar and humans. Using phylogenetic and syntenic analyses, we found that nfil3-1 is conserved in vertebrates, while nfil3-2 and nfil3-3 are maintained in spotted gar, teleost fish, amphibians, and reptiles, but not in mammals. Following the two-round vertebrate genome duplication (VGD), spotted gar retained cry1a, cry1b, and cry2, and cry3 is retained in spotted gar, teleost fish, turtles, and birds, but not in mammals. We hypothesize that duplication of core clock genes, such as (nfil3 and cry), likely facilitated diversification of circadian regulatory mechanisms in teleost fish. We also found that the transcription factor binding element (Ahr::Arnt) is retained only in one of the per1 or per2 duplicated paralogs derived from the TGD in the teleost fish, implicating possible subfuctionalization cases. Together, these findings help decipher the repertoires of the spotted gar's circadian system and shed light on how the vertebrate circadian clock systems have evolved.

Deciphering The Potential Role of Hox Genes in Pancreatic Cancer

The Hox gene family plays an important role in organogenesis and animal development. Currently, 39 Hox genes that are clustered in four chromosome regions have been identified in humans. Emerging evidence suggests that Hox genes are involved in the development of the pancreas. However, the expression of Hox genes in pancreatic tumor tissues has been investigated in only a few studies. In addition, whether specific Hox genes can promote or suppress cancer metastasis is not clear. In this article, we first review the recent progress in studies on the role of Hox genes in pancreatic cancer. By comparing the expression profiles of pancreatic cancer cells isolated from genetically engineered mice established in our laboratory with three different proliferative and metastatic abilities, we identified novel Hox genes that exhibited tumor-promoting activity in pancreatic cancer. Finally, a potential oncogenic mechanism of the Hox genes was hypothesized.

Gonadal soma controls ovarian follicle proliferation through Gsdf in zebrafish

BACKGROUND

Aberrant signaling between germ cells and somatic cells can lead to reproductive disease and depends on diffusible signals, including transforming growth factor-beta (TGFB) -family proteins. The TGFB-family protein Gsdf (gonadal soma derived factor) controls sex determination in some fish and is a candidate for mediating germ cell/soma signaling.

RESULTS

Zebrafish expressed gsdf in somatic cells of bipotential gonads and expression continued in ovarian granulosa cells and testicular Sertoli cells. Homozygous gsdf knockout mutants delayed leaving the bipotential gonad state, but then became a male or a female. Mutant females ovulated a few oocytes, then became sterile, accumulating immature follicles. Female mutants stored excess lipid and down-regulated aromatase, gata4, insulin receptor, estrogen receptor, and genes for lipid metabolism, vitellogenin, and steroid biosynthesis. Mutant females contained less estrogen and more androgen than wild-types. Mutant males were fertile. Genomic analysis suggests that Gsdf, Bmp15, and Gdf9, originated as paralogs in vertebrate genome duplication events.

CONCLUSIONS

In zebrafish, gsdf regulates ovarian follicle maturation and expression of genes for steroid biosynthesis, obesity, diabetes, and female fertility, leading to ovarian and extra-ovarian phenotypes that mimic human polycystic ovarian syndrome (PCOS), suggesting a role for a related TGFB signaling molecule in the etiology of PCOS. Developmental Dynamics 246:925-945, 2017. © 2017 Wiley Periodicals, Inc.

Support for Lungfish as the Closest Relative of Tetrapods by Using Slowly Evolving Ray-Finned Fish as the Outgroup

In a previous analysis of the phylogenetic relationships of coelacanths, lungfishes and tetrapods, using cartilaginous fish (CF) as the outgroup, the sister relationship of lungfishes and tetrapods was constructed with high statistical support. However, using as the outgroup ray-finned fish (RF), which are more taxonomically closely related to the three lineages than CF, the sister relationship of coelacanths and tetrapods was most often constructed depending on the methods and the data sets, but the statistical support was generally low except in the cases in which the data set including a small number of species was analyzed. In this study, instead of the fast evolving ray-finned fish, teleost fish (TF), in the previous data sets, by using two slowly evolving RF, gar and bowfin, as the outgroup, we showed that the sister relationship of lungfishes and tetrapods was reconstructed with high statistical support. In our analysis the evolutionary rates of gar and bowfin were similar to each other and one third to one half of TF. The difference of the amino acid frequencies of the two species with other lineages was larger than those of TF. This study provides a strong support for lungfishes as the closest relative of tetrapods and indicates the importance of using an appropriate outgroup with small divergence in phylogenetic construction.

Evolutionary conservation and expression of miR-10a-3p in olive flounder and rock bream

MicroRNAs (miRNAs) are small non-coding RNAs (ncRNAs) that mainly bind to the seed sequences located within the 3' untranslated region (3' UTR) of target genes. They perform an important biological function as regulators of gene expression. Different genes can be regulated by the same miRNA, whilst different miRNAs can be regulated by the same genes. Here, the evolutionary conservation and expression pattern of miR-10a-3p in olive flounder and rock bream was examined. Binding sites (AAAUUC) to seed region of the 3' UTR of target genes were highly conserved in various species. The expression pattern of miR-10a-3p was ubiquitous in the examined tissues, whilst its expression level was decreased in gill tissues infected by viral hemorrhagic septicemia virus (VHSV) compared to the normal control. In the case of rock bream, the spleen, kidney, and liver tissues showed dominant expression levels of miR-10a-3p. Only the liver tissues in the rock bream samples infected by the iridovirus indicated a dominant miR-10a-3p expression. The gene ontology (GO) analysis of predicted target genes for miR-10a-3p revealed that multiple genes are related to binding activity, catalytic activity, cell components as well as cellular and metabolic process. Overall the results imply that the miR-10a-3p could be used as a biomarker to detect VHSV infection in olive flounder and iridovirus infection in rock bream. In addition, the data provides fundamental information for further study of the complex interaction between miR-10a-3p and gene expression.

Lineage-specific rediploidization is a mechanism to explain time-lags between genome duplication and evolutionary diversification

BACKGROUND

The functional divergence of duplicate genes (ohnologues) retained from whole genome duplication (WGD) is thought to promote evolutionary diversification. However, species radiation and phenotypic diversification are often temporally separated from WGD. Salmonid fish, whose ancestor underwent WGD by autotetraploidization ~95 million years ago, fit such a 'time-lag' model of post-WGD radiation, which occurred alongside a major delay in the rediploidization process. Here we propose a model, 'lineage-specific ohnologue resolution' (LORe), to address the consequences of delayed rediploidization. Under LORe, speciation precedes rediploidization, allowing independent ohnologue divergence in sister lineages sharing an ancestral WGD event.

RESULTS

Using cross-species sequence capture, phylogenomics and genome-wide analyses of ohnologue expression divergence, we demonstrate the major impact of LORe on salmonid evolution. One-quarter of each salmonid genome, harbouring at least 4550 ohnologues, has evolved under LORe, with rediploidization and functional divergence occurring on multiple independent occasions >50 million years post-WGD. We demonstrate the existence and regulatory divergence of many LORe ohnologues with functions in lineage-specific physiological adaptations that potentially facilitated salmonid species radiation. We show that LORe ohnologues are enriched for different functions than 'older' ohnologues that began diverging in the salmonid ancestor.

CONCLUSIONS

LORe has unappreciated significance as a nested component of post-WGD divergence that impacts the functional properties of genes, whilst providing ohnologues available solely for lineage-specific adaptation. Under LORe, which is predicted following many WGD events, the functional outcomes of WGD need not appear 'explosively', but can arise gradually over tens of millions of years, promoting lineage-specific diversification regimes under prevailing ecological pressures.

"Holostei versus Halecostomi" Problem: Insight from Cytogenetics of Ancient Nonteleost Actinopterygian Fish, Bowfin Amia calva

Bowfin belongs to an ancient lineage of nonteleost ray-finned fishes (actinopterygians) and is the only extant survivor of a once diverged group, the Halecomorphi or Amiiformes. Owing to the scarcity of extant nonteleost ray-finned lineages, also referred as "living fossils," their phylogenetic interrelationships have been the target of multiple hypotheses concerning their sister group relationships. Molecular and morphological data sets have produced controversial results; bowfin is considered as either the sister group to genome-duplicated teleosts (together forming the group of Halecostomi) or to gars (Lepisosteiformes; together forming the group of Holostei). However, any detailed cytogenetic analysis of bowfin chromosomes has never been performed to address this issue. Here we examined bowfin chromosomes by conventional (Giemsa-staining, C-banding, base-specific fluorescence and silver staining) and molecular (FISH with rDNA probes) cytogenetic protocols. We identified diploid chromosome number 2n = 46 with a middle-sized submetacentric chromosome pair as the major ribosomal DNA-bearing (45S rDNA), GC-positive and silver-positive element. The minor rDNA (5S rDNA) sites were localized in the pericentromeric region of one middle-sized acrocentric chromosome pair. Comparison with available cytogenetic data of other nonteleost actinopterygians (bichirs, sturgeons, gars) and teleost species including representative of basally branching lineages showed bowfin chromosomal characteristics more similar to the teleost type than to any other nonteleosts. Particularly striking differences were identified between bowfin and gars, the latter of which were found to mimic mammalian AT/GC genomic organisation. Such conclusion however contradicts the most recent phylogenomic results and raises the question what states are ancestral and what are derived.

Identification and conservation of gene loss events of Hox gene clusters in the marine medaka (Oryzias melastigma)

In this study, the identification of the whole Hox gene clusters (46 Hox genes) in the marine medaka Oryzias melastigma was investigated using genome assembly and RNA-seq information. Moreover, the gene loss events of Hox gene clusters, which may occur during fish evolution, were examined for a better understanding of the evolutionary status of the gene lost events of the Hox gene cluster across fish species, particularly in the genus Oryzias.

Molecular phylogeny and patterns of diversification in syngnathid fishes

The family Syngnathidae is a large and diverse clade of morphologically unique bony fishes, with 57 genera and 300 described species of seahorses, pipefishes, pipehorses, and seadragons. They primarily inhabit shallow coastal waters in temperate and tropical oceans, and are characterized by a fused jaw, male brooding, and extraordinary crypsis. Phylogenetic relationships within the Syngnathidae remain poorly resolved due to lack of generic taxon sampling, few diagnostic morphological characters, and limited molecular data. The phylogenetic placement of the threatened, commercially exploited seahorses remains a topic of intense interest, with conflicting topologies based on morphology and predominantly mitochondrial genetic data. In this study, we integrate eight nuclear and mitochondrial markers and 17 morphological characters to investigate the phylogenetic structure of the family Syngnathidae at the generic level. We include 91 syngnathid species representing 48 of the 57 recognized genera, all major ocean basins, and a broad array of temperate and tropical habitats including rocky and coral reefs, sand and silt, mangroves, seagrass beds, estuaries, and rivers. Maximum likelihood and Bayesian analyses of 5160bp from eight loci produced high congruence among alternate topologies, defining well-supported and sometimes novel clades. We present a hypothesis that confirms a deep phylogenetic split between lineages with trunk- or tail-brood pouch placement, and provides significant new insights into the morphological evolution and biogeography of this highly derived fish clade. Based on the fundamental division between lineages - the tail brooding "Urophori" and the trunk brooding "Gastrophori" - we propose a revision of Syngnathidae classification into only two subfamilies: the Nerophinae and the Syngnathinae. We find support for distinct principal clades within the trunk-brooders and tail-brooders, the latter of which include seahorses, seadragons, independent lineages of pipehorses, and clades that originated in southern Australia and the Western Atlantic. We suggest the seahorse genus Hippocampus is of Indo-Pacific origin and its sister clade is an unexpected grouping of several morphologically disparate Indo-Pacific genera, including the Pacific pygmy pipehorses. Taxonomic revision is required for multiple genera, particularly to reflect deep evolutionary splits in nominal lineages from the Atlantic versus the Indo-Pacific.

Little evidence for enhanced phenotypic evolution in early teleosts relative to their living fossil sister group

Since Darwin, biologists have been struck by the extraordinary diversity of teleost fishes, particularly in contrast to their closest "living fossil" holostean relatives. Hypothesized drivers of teleost success include innovations in jaw mechanics, reproductive biology and, particularly at present, genomic architecture, yet all scenarios presuppose enhanced phenotypic diversification in teleosts. We test this key assumption by quantifying evolutionary rate and capacity for innovation in size and shape for the first 160 million y (Permian-Early Cretaceous) of evolution in neopterygian fishes (the more extensive clade containing teleosts and holosteans). We find that early teleosts do not show enhanced phenotypic evolution relative to holosteans. Instead, holostean rates and innovation often match or can even exceed those of stem-, crown-, and total-group teleosts, belying the living fossil reputation of their extant representatives. In addition, we find some evidence for heterogeneity within the teleost lineage. Although stem teleosts excel at discovering new body shapes, early crown-group taxa commonly display higher rates of shape evolution. However, the latter reflects low rates of shape evolution in stem teleosts relative to all other neopterygian taxa, rather than an exceptional feature of early crown teleosts. These results complement those emerging from studies of both extant teleosts as a whole and their sublineages, which generally fail to detect an association between genome duplication and significant shifts in rates of lineage diversification.

Identification of a novel leptin receptor duplicate in Atlantic salmon: Expression analyses in different life stages and in response to feeding status

In recent years rapidly growing research has led to identification of several fish leptin orthologs and numerous duplicated paralogs possibly arisen from the third and fourth round whole genome duplication (3R and 4R WGD) events. In this study we identify in Atlantic salmon a duplicated LepRA gene, named LepRA2, that further extend possible evolutionary scenarios of the leptin and leptin receptor system. The 1121 amino acid sequence of the novel LepRA2 shares 80% sequence identity with the LepRA1 paralog, and contains the protein motifs typical of the functional (long form) leptin receptor in vertebrates. In silico predictions showed similar electrostatic properties of LepRA1 and LepRA2 and high sequence conservation at the leptin interaction surfaces within the CHR/leptin-binding and FNIII domains, suggesting conserved functional specificity between the two duplicates. Analysis of temporal expression profiles during pre-hatching stages indicate that both transcripts are involved in modulating leptin developmental functions, although the LepRA1 paralog may play a major role as the embryo complexity increases. There is ubiquitous distribution of LepRs underlying pleiotropism of leptin in all tissues investigated. LepRA1 and LepRA2 are differentially expressed with LepRA1 more abundant than LepRA2 in most of the tissues investigated, with the only exception of liver. Analysis of constitutive LepRA1 and LepRA2 expression in brain and liver at parr, post-smolt and adult stages reveal striking spatial divergence between the duplicates at all stages investigated. This suggests that, beside increased metabolic requirements, leptin sensitivity in the salmon brain might be linked to important variables such as habitat, ecology and life cycle. Furthermore, leptins and LepRs mRNAs in the brain showed gene-specific variability in response to long term fasting, suggesting that leptin's roles as modulator of nutritional status in Atlantic salmon might be governed by distinct genetic evolutionary processes and distinct functions between the paralogs.

Whole Genome Duplications Shaped the Receptor Tyrosine Kinase Repertoire of Jawed Vertebrates

The receptor tyrosine kinase (RTK) gene family, involved primarily in cell growth and differentiation, comprises proteins with a common enzymatic tyrosine kinase intracellular domain adjacent to a transmembrane region. The amino-terminal portion of RTKs is extracellular and made of different domains, the combination of which characterizes each of the 20 RTK subfamilies among mammals. We analyzed a total of 7,376 RTK sequences among 143 vertebrate species to provide here the first comprehensive census of the jawed vertebrate repertoire. We ascertained the 58 genes previously described in the human and mouse genomes and established their phylogenetic relationships. We also identified five additional RTKs amounting to a total of 63 genes in jawed vertebrates. We found that the vertebrate RTK gene family has been shaped by the two successive rounds of whole genome duplications (WGD) called 1R and 2R (1R/2R) that occurred at the base of the vertebrates. In addition, the Vegfr and Ephrin receptor subfamilies were expanded by single gene duplications. In teleost fish, 23 additional RTK genes have been retained after another expansion through the fish-specific third round (3R) of WGD. Several lineage-specific gene losses were observed. For instance, birds have lost three RTKs, and different genes are missing in several fish sublineages. The RTK gene family presents an unusual high gene retention rate from the vertebrate WGDs (58.75% after 1R/2R, 64.4% after 3R), resulting in an expansion that might be correlated with the evolution of complexity of vertebrate cellular communication and intracellular signaling.

Resolving the Phylogenetic Position of Coelacanth: The Closest Relative Is Not Always the Most Appropriate Outgroup

Determining the phylogenetic relationship of two extant lineages of lobe-finned fish, coelacanths and lungfishes, and tetrapods is important for understanding the origin of tetrapods. We analyzed data sets from two previous studies along with a newly collected data set, each of which had varying numbers of species and genes and varying extent of missing sites. We found that in all the data sets the sister relationship of lungfish and tetrapods was constructed with the use of cartilaginous fish as the outgroup with a high degree of statistical support. In contrast, when ray-finned fish were used as the outgroup, which is taxonomically an immediate outgroup of lobe-finned fish and tetrapods, the sister relationship of coelacanth and tetrapods was supported most strongly, although the statistical support was weaker. Even though it is generally accepted that the closest relative is an appropriate outgroup, our analysis suggested that the large divergence of the ray-finned fish as indicated by their long branch lengths and different amino acid frequencies made them less suitable as an outgroup than cartilaginous fish.

Evolutionary significance and diversification of the phosphoglucose isomerase genes in vertebrates

Background

Phosphoglucose isomerase (PGI) genes are important multifunctional proteins whose evolution has, until now, not been well elucidated because of the limited number of completely sequenced genomes. Although the multifunctionality of this gene family has been considered as an original and innate characteristic, PGI genes may have acquired novel functions through changes in coding sequences and exon/intron structure, which are known to lead to functional divergence after gene duplication. A whole-genome comparative approach was used to estimate the rates of molecular evolution of this protein family.

Results

The results confirm the presence of two isoforms in teleost fishes and only one variant in all other vertebrates. Phylogenetic reconstructions grouped the PGI genes into five main groups: lungfishes/coelacanth/cartilaginous fishes, teleost fishes, amphibians, reptiles/birds and mammals, with the teleost group being subdivided into two subclades comprising PGI1 and PGI2. This PGI partitioning into groups is consistent with the synteny and molecular evolution results based on the estimation of the ratios of nonsynonymous to synonymous changes (Ka/Ks) and divergence rates between both PGI paralogs and orthologs. Teleost PGI2 shares more similarity with the variant found in all other vertebrates, suggesting that it has less evolved than PGI1 relative to the PGI of common vertebrate ancestor.

Conclusions

The diversification of PGI genes into PGI1 and PGI2 is consistent with a teleost-specific duplication before the radiation of this lineage, and after its split from the other infraclasses of ray-finned fishes. The low average Ka/Ks ratios within teleost and mammalian lineages suggest that both PGI1 and PGI2 are functionally constrained by purifying selection and may, therefore, have the same functions. By contrast, the high average Ka/Ks ratios and divergence rates within reptiles and birds indicate that PGI may be involved in different functions. The synteny analyses show that the genomic region harbouring PGI genes has independently undergone genomic rearrangements in mammals versus the reptile/bird lineage in particular, which may have contributed to the actual functional diversification of this gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1683-x) contains supplementary material, which is available to authorized users.

Diversification of non-visual photopigment parapinopsin in spectral sensitivity for diverse pineal functions

Background

Recent genome projects of various animals have uncovered an unexpectedly large number of opsin genes, which encode protein moieties of photoreceptor molecules, in most animals. In visual systems, the biological meanings of this diversification are clear; multiple types of visual opsins with different spectral sensitivities are responsible for color vision. However, the significance of the diversification of non-visual opsins remains uncertain, in spite of the importance of understanding the molecular mechanism and evolution of varied non-visual photoreceptions.

Results

Here, we investigated the diversification of the pineal photopigment parapinopsin, which serves as the UV-sensitive photopigment for the pineal wavelength discrimination in the lamprey, linking it with other pineal photoreception. Spectroscopic analyses of the recombinant pigments of the two teleost parapinopsins PP1 and PP2 revealed that PP1 is a UV-sensitive pigment, similar to lamprey parapinopsin, but PP2 is a blue-sensitive pigment, with an absorption maximum at 460–480 nm, showing the diversification of non-visual pigment with respect to spectral sensitivity. We also found that PP1 and PP2 exhibit mutually exclusive expressions in the pineal organs of three teleost species. By using transgenic zebrafish in which these parapinopsin-expressing cells are labeled, we found that PP1-expressing cells basically possess neuronal processes, which is consistent with their involvement in wavelength discrimination. Interestingly, however, PP2-expressing cells rarely possess neuronal processes, raising the possibility that PP2 could be involved in non-neural responses rather than neural responses. Furthermore, we found that PP2-expressing cells contain serotonin and aanat2, the key enzyme involved in melatonin synthesis from serotonin, whereas PP1-expressing cells do not contain either, suggesting that blue-sensitive PP2 is instead involved in light-regulation of melatonin secretion.

Conclusions

In this paper, we have clearly shown the different molecular properties of duplicated non-visual opsins by demonstrating the diversification of parapinopsin with respect to spectral sensitivity. Moreover, we have shown a plausible link between the diversification and its physiological impact by discovering a strong candidate for the underlying pigment in light-regulated melatonin secretion in zebrafish; the diversification could generate a new contribution of parapinopsin to pineal photoreception. Current findings could also provide an opportunity to understand the “color” preference of non-visual photoreception.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0174-9) contains supplementary material, which is available to authorized users.

Genome-wide analysis of homeobox genes from Mesobuthus martensii reveals Hox gene duplication in scorpions

Homeobox genes belong to a large gene group, which encodes the famous DNA-binding homeodomain that plays a key role in development and cellular differentiation during embryogenesis in animals. Here, one hundred forty-nine homeobox genes were identified from the Asian scorpion, Mesobuthus martensii (Chelicerata: Arachnida: Scorpiones: Buthidae) based on our newly assembled genome sequence with approximately 248 × coverage. The identified homeobox genes were categorized into eight classes including 82 families: 67 ANTP class genes, 33 PRD genes, 11 LIM genes, five POU genes, six SINE genes, 14 TALE genes, five CUT genes, two ZF genes and six unclassified genes. Transcriptome data confirmed that more than half of the genes were expressed in adults. The homeobox gene diversity of the eight classes is similar to the previously analyzed Mandibulata arthropods. Interestingly, it is hypothesized that the scorpion M. martensii may have two Hox clusters. The first complete genome-wide analysis of homeobox genes in Chelicerata not only reveals the repertoire of scorpion, arachnid and chelicerate homeobox genes, but also shows some insights into the evolution of arthropod homeobox genes.

Molecular cloning, functional characterization, and evolutionary analysis of vitamin D receptors isolated from basal vertebrates

The vertebrate genome is a result of two rapid and successive rounds of whole genome duplication, referred to as 1R and 2R. Furthermore, teleost fish have undergone a third whole genome duplication (3R) specific to their lineage, resulting in the retention of multiple gene paralogs. The more recent 3R event in teleosts provides a unique opportunity to gain insight into how genes evolve through specific evolutionary processes. In this study we compare molecular activities of vitamin D receptors (VDR) from basal species that diverged at key points in vertebrate evolution in order to infer derived and ancestral VDR functions of teleost paralogs. Species include the sea lamprey (Petromyzon marinus), a 1R jawless fish; the little skate (Leucoraja erinacea), a cartilaginous fish that diverged after the 2R event; and the Senegal bichir (Polypterus senegalus), a primitive 2R ray-finned fish. Saturation binding assays and gel mobility shift assays demonstrate high affinity ligand binding and classic DNA binding characteristics of VDR has been conserved across vertebrate evolution. Concentration response curves in transient transfection assays reveal EC50 values in the low nanomolar range, however maximum transactivational efficacy varies significantly between receptor orthologs. Protein-protein interactions were investigated using co-transfection, mammalian 2-hybrid assays, and mutations of coregulator activation domains. We then combined these results with our previous study of VDR paralogs from 3R teleosts into a bioinformatics analysis. Our results suggest that 1, 25D3 acts as a partial agonist in basal species. Furthermore, our bioinformatics analysis suggests that functional differences between VDR orthologs and paralogs are influenced by differential protein interactions with essential coregulator proteins. We speculate that we may be observing a change in the pharmacodynamics relationship between VDR and 1, 25D3 throughout vertebrate evolution that may have been driven by changes in protein-protein interactions between VDR and essential coregulators.

Functional diversification of vitamin D receptor paralogs in teleost fish after a whole genome duplication event

The diversity and success of teleost fishes (Actinopterygii) has been attributed to three successive rounds of whole-genome duplication (WGD). WGDs provide a source of raw genetic material for evolutionary forces to act upon, resulting in the divergence of genes with altered or novel functions. The retention of multiple gene pairs (paralogs) in teleosts provides a unique opportunity to study how genes diversify and evolve after a WGD. This study examines the hypothesis that vitamin D receptor (VDR) paralogs (VDRα and VDRβ) from two distantly related teleost orders have undergone functional divergence subsequent to the teleost-specific WGD. VDRα and VDRβ paralogs were cloned from the Japanese medaka (Beloniformes) and the zebrafish (Cypriniformes). Initial transactivation studies using 1α, 25-dihydroxyvitamin D3 revealed that although VDRα and VDRβ maintain similar ligand potency, the maximum efficacy of VDRβ was significantly attenuated compared with VDRα in both species. Subsequent analyses revealed that VDRα and VDRβ maintain highly similar ligand affinities; however, VDRα demonstrated preferential DNA binding compared with VDRβ. Protein-protein interactions between the VDR paralogs and essential nuclear receptor coactivators were investigated using transactivation and mammalian two-hybrid assays. Our results imply that functional differences between VDRα and VDRβ occurred early in teleost evolution because they are conserved between distantly related species. Our results further suggest that the observed differences may be associated with differential protein-protein interactions between the VDR paralogs and coactivators. We speculate that the observed functional differences are due to subtle ligand-induced conformational differences between the two paralogs, leading to divergent downstream functions.

Molecular cloning and function analysis of insulin-like growth factor-binding protein 1a in blunt snout bream (Megalobrama amblycephala)

Insulin-like growth factor-binding protein 1 (IGFBP-1), a hypoxia-induced protein, is a member of the IGFBP family that regulates vertebrate growth and development. In this study, full-length IGFBP-1a cDNA was cloned from a hypoxia-sensitive Cyprinidae fish species, the blunt snout bream (Megalobrama amblycephala). IGFBP-1a was expressed in various organs of adult blunt snout bream, including strongly in the liver and weakly in the gonads. Under hypoxia, IGFBP-1a mRNA levels increased sharply in the skin, liver, kidney, spleen, intestine and heart tissues of juvenile blunt snout bream, but recovered to normal levels after 24-hour exposure to normal dissolved oxygen. In blunt snout bream embryos, IGFBP-1a mRNA was expressed at very low levels at both four and eight hours post-fertilization, and strongly at later stages. Embryonic growth and development rates decreased significantly in embryos injected with IGFBP-1a mRNA. The average body length of IGFBP-1a-overexpressed embryos was 82.4% of that of the control group, and somite numbers decreased to 85.2%. These findings suggest that hypoxia-induced IGFBP-1a may inhibit growth in this species under hypoxic conditions.

Enigmatic orthology relationships between Hox clusters of the African butterfly fish and other teleosts following ancient whole-genome duplication

Numerous ancient whole-genome duplications (WGD) have occurred during eukaryote evolution. In vertebrates, duplicated developmental genes and their functional divergence have had important consequences for morphological evolution. Although two vertebrate WGD events (1R/2R) occurred over 525 Ma, we have focused on the more recent 3R or TGD (teleost genome duplication) event which occurred approximately 350 Ma in a common ancestor of over 26,000 species of teleost fishes. Through a combination of whole genome and bacterial artificial chromosome clone sequencing we characterized all Hox gene clusters of Pantodon buchholzi, a member of the early branching teleost subdivision Osteoglossomorpha. We find 45 Hox genes organized in only five clusters indicating that Pantodon has suffered more Hox cluster loss than other known species. Despite strong evidence for homology of the five Pantodon clusters to the four canonical pre-TGD vertebrate clusters (one HoxA, two HoxB, one HoxC, and one HoxD), we were unable to confidently resolve 1:1 orthology relationships between four of the Pantodon clusters and the eight post-TGD clusters of other teleosts. Phylogenetic analysis revealed that many Pantodon genes segregate outside the conventional "a" and "b" post-TGD orthology groups, that extensive topological incongruence exists between genes physically linked on a single cluster, and that signal divergence causes ambivalence in assigning 1:1 orthology in concatenated Hox cluster analyses. Out of several possible explanations for this phenomenon we favor a model which keeps with the prevailing view of a single TGD prior to teleost radiation, but which also considers the timing of diploidization after duplication, relative to speciation events. We suggest that although the duplicated hoxa clusters diploidized prior to divergence of osteoglossomorphs, the duplicated hoxb, hoxc, and hoxd clusters concluded diploidization independently in osteoglossomorphs and other teleosts. We use the term "tetralogy" to describe the homology relationship which exists between duplicated sequences which originate through a shared WGD, but which diploidize into distinct paralogs from a common allelic pool independently in two lineages following speciation.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.

Ploidy-dependent survival of progeny arising from crosses between natural allotriploid Cobitis females and diploid C. taenia males (Pisces, Cobitidae)

Crosses between 21 triploid hybrid Cobitis females and 19 C. taenia (2n = 48) males led to viable progeny; whereas no embryonic development was observed in crosses with tetraploid males (4n = 98). The ploidy status of 491 progenies randomly selected with flow cytometry (316) or chromosome analysis (175) revealed an average of 55.2 % triploids and 44.8 % tetraploids, but the ratio of 3n versus 4n fish did change during development. In the first 2 days after hatching, approximately 65.1 % of tetraploid larvae were observed. Their number decreased significantly to 30.8 and 6.2 % on average during 2-5 and 10-15 months of life, respectively. The karyotype of tetraploid progeny (4n = 98) included 3n = 74 chromosomes of the parental female and n = 24 of C. taenia male. The number of tetraploid progeny indicated indirectly that about 66 % of eggs from 3n females were fertilized with C. taenia. The rest of the eggs developed clonally via gynogenesis or hemiclonally via hybridogenesis into triploids of the same karyotype structure as parental females. We have documented for the first time that (at least under experimental conditions) tetraploids are commonly formed, but are less viable than triploids, and a ratio similar to what is found under natural conditions is finally attained. The current explanation concerning the ploidy and karyotype structure of the progeny confirms that the eggs of 3n Cobitis females are not only capable of maintaining all chromosomes but are also capable of incorporating the sperm genome, thus creating the potential to produce tetraploids.

MicroRNA in teleost fish

MicroRNAs (miRNAs) are transcriptional and posttranscriptional regulators involved in nearly all known biological processes in distant eukaryotic clades. Their discovery and functional characterization have broadened our understanding of biological regulatory mechanisms in animals and plants. They show both evolutionary conserved and unique features across Metazoa. Here, we present the current status of the knowledge about the role of miRNA in development, growth, and physiology of teleost fishes, in comparison to other vertebrates. Infraclass Teleostei is the most abundant group among vertebrate lineage. Fish are an important component of aquatic ecosystems and human life, being the prolific source of animal proteins worldwide and a vertebrate model for biomedical research. We review miRNA biogenesis, regulation, modifications, and mechanisms of action. Specific sections are devoted to the role of miRNA in teleost development, organogenesis, tissue differentiation, growth, regeneration, reproduction, endocrine system, and responses to environmental stimuli. Each section discusses gaps in the current knowledge and pinpoints the future directions of research on miRNA in teleosts.

Evolution of the vertebrate Pax4/6 class of genes with focus on its novel member, the Pax10 gene

The members of the paired box (Pax) family regulate key developmental pathways in many metazoans as tissue-specific transcription factors. Vertebrate genomes typically possess nine Pax genes (Pax1-9), which are derived from four proto-Pax genes in the vertebrate ancestor that were later expanded through the so-called two-round (2R) whole-genome duplication. A recent study proposed that pax6a genes of a subset of teleost fishes (namely, acanthopterygians) are remnants of a paralog generated in the 2R genome duplication, to be renamed pax6.3, and reported one more group of vertebrate Pax genes (Pax6.2), most closely related to the Pax4/6 class. We propose to designate this new member Pax10 instead and reconstruct the evolutionary history of the Pax4/6/10 class with solid phylogenetic evidence. Our synteny analysis showed that Pax4, -6, and -10 originated in the 2R genome duplications early in vertebrate evolution. The phylogenetic analyses of relationships between teleost pax6a and other Pax4, -6, and -10 genes, however, do not support the proposed hypothesis of an ancient origin of the acanthopterygian pax6a genes in the 2R genome duplication. Instead, we confirmed the traditional scenario that the acanthopterygian pax6a is derived from the more recent teleost-specific genome duplication. Notably, Pax6 is present in all vertebrates surveyed to date, whereas Pax4 and -10 were lost multiple times in independent vertebrate lineages, likely because of their restricted expression patterns: Among Pax6-positive domains, Pax10 has retained expression in the adult retina alone, which we documented through in situ hybridization and quantitative reverse transcription polymerase chain reaction experiments on zebrafish, Xenopus, and anole lizard.

How Hox genes can shed light on the place of echinoderms among the deuterostomes

Background

The Hox gene cluster ranks among the greatest of biological discoveries of the past 30 years. Morphogenetic patterning genes are remarkable for the systems they regulate during major ontogenetic events, and for their expressions of molecular, temporal, and spatial colinearity. Recent descriptions of exceptions to these colinearities are suggesting deep phylogenetic signal that can be used to explore origins of entire deuterostome phyla. Among the most enigmatic of these deuterostomes in terms of unique body patterning are the echinoderms. However, there remains no overall synthesis of the correlation between this signal and the variations observable in the presence/absence and expression patterns of Hox genes.

Results

Recent data from Hox cluster analyses shed light on how the bizarre shift from bilateral larvae to radial adults during echinoderm ontogeny can be accomplished by equally radical modifications within the Hox cluster. In order to explore this more fully, a compilation of observations on the genetic patterns among deuterostomes is integrated with the body patterning trajectories seen across the deuterostome clade.

Conclusions

Synthesis of available data helps to explain morphogenesis along the anterior/posterior axis of echinoderms, delineating the origins and fate of that axis during ontogeny. From this, it is easy to distinguish between ‘seriality’ along echinoderm rays and true A/P axis phenomena such as colinearity within the somatocoels, and the ontogenetic outcomes of the unique translocation and inversion of the anterior Hox class found within the Echinodermata. An up-to-date summary and integration of the disparate lines of research so far produced on the relationship between Hox genes and pattern formation for all deuterostomes allows for development of a phylogeny and scenario for the evolution of deuterostomes in general, and the Echinodermata in particular.

Why do a wide variety of animals retain multiple isoforms of cyclooxygenase?

Cyclooxygenase (COX) has been cloned from the phyla Cnidaria, Mollusca, Arthropoda, and Chordata of the animal kingdom. Many organisms have multiple COX isoforms that have arisen from gene duplication. It is not well understood why there are multiple COX isoforms in the same organism, or when duplication of the COX gene occurred. Here, we summarize the current knowledge of the evolutionary history of COX in the animal kingdom and discuss the reasons why the multiple COX system has been retained so widely. The phylogenetic analysis suggests that all COX genes in animals may descend from a common ancestor and that the duplication of an ancestral COX gene might occur within each lineage after the divergence of the animal. In most instances, the expressions of multiple COX isoforms are separately regulated and these isoforms play different and important pathophysiological roles in each organism. This may be the reason why multiple COX isoforms are widely retained.