The tree of life and a new classification of bony fishes

The timing of Timezyme diversification in vertebrates

All biological functions in vertebrates are synchronized with daily and seasonal changes in the environment by the time keeping hormone melatonin. Its nocturnal surge is primarily due to the rhythmic activity of the arylalkylamine N-acetyl transferase AANAT, which thus became the focus of many investigations regarding its evolution and function. Various vertebrate isoforms have been reported from cartilaginous fish to mammals but their origin has not been clearly established. Using phylogeny and synteny, we took advantage of the increasing number of available genomes in order to test whether the various rounds of vertebrate whole genome duplications were responsible for the diversification of AANAT. We highlight a gene secondary loss of the AANAT2 in the Sarcopterygii, revealing for the first time that the AAANAT1/2 duplication occurred before the divergence between Actinopterygii (bony fish) and Sarcopterygii (tetrapods, lobe-finned fish, and lungfish). We hypothesize the teleost-specific whole genome duplication (WDG) generated the appearance of the AANAT1a/1b and the AANAT2/2′paralogs, the 2′ isoform being rapidly lost in the teleost common ancestor (ray-finned fish). We also demonstrate the secondary loss of the AANAT1a in a Paracantopterygii (Atlantic cod) and of the 1b in some Ostariophysi (zebrafish and cave fish). Salmonids present an even more diverse set of AANATs that may be due to their specific WGD followed by secondary losses. We propose that vertebrate AANAT diversity resulted from 3 rounds of WGD followed by previously uncharacterized secondary losses. Extant isoforms show subfunctionalized localizations, enzyme activities and affinities that have increased with time since their emergence.

Permian-Triassic Osteichthyes (bony fishes): diversity dynamics and body size evolution

The Permian and Triassic were key time intervals in the history of life on Earth. Both periods are marked by a series of biotic crises including the most catastrophic of such events, the end-Permian mass extinction, which eventually led to a major turnover from typical Palaeozoic faunas and floras to those that are emblematic for the Mesozoic and Cenozoic. Here we review patterns in Permian-Triassic bony fishes, a group whose evolutionary dynamics are understudied. Based on data from primary literature, we analyse changes in their taxonomic diversity and body size (as a proxy for trophic position) and explore their response to Permian-Triassic events. Diversity and body size are investigated separately for different groups of Osteichthyes (Dipnoi, Actinistia, 'Palaeopterygii', 'Subholostei', Holostei, Teleosteomorpha), within the marine and freshwater realms and on a global scale (total diversity) as well as across palaeolatitudinal belts. Diversity is also measured for different palaeogeographical provinces. Our results suggest a general trend from low osteichthyan diversity in the Permian to higher levels in the Triassic. Diversity dynamics in the Permian are marked by a decline in freshwater taxa during the Cisuralian. An extinction event during the end-Guadalupian crisis is not evident from our data, but 'palaeopterygians' experienced a significant body size increase across the Guadalupian-Lopingian boundary and these fishes upheld their position as large, top predators from the Late Permian to the Late Triassic. Elevated turnover rates are documented at the Permian-Triassic boundary, and two distinct diversification events are noted in the wake of this biotic crisis, a first one during the Early Triassic (dipnoans, actinistians, 'palaeopterygians', 'subholosteans') and a second one during the Middle Triassic ('subholosteans', neopterygians). The origination of new, small taxa predominantly among these groups during the Middle Triassic event caused a significant reduction in osteichthyan body size. Neopterygii, the clade that encompasses the vast majority of extant fishes, underwent another diversification phase in the Late Triassic. The Triassic radiation of Osteichthyes, predominantly of Actinopterygii, which only occurred after severe extinctions among Chondrichthyes during the Middle-Late Permian, resulted in a profound change within global fish communities, from chondrichthyan-rich faunas of the Permo-Carboniferous to typical Mesozoic and Cenozoic associations dominated by actinopterygians. This turnover was not sudden but followed a stepwise pattern, with leaps during extinction events.

The first virtual cranial endocast of a lungfish (sarcopterygii: dipnoi)

Lungfish, or dipnoans, have a history spanning over 400 million years and are the closest living sister taxon to the tetrapods. Most Devonian lungfish had heavily ossified endoskeletons, whereas most Mesozoic and Cenozoic lungfish had largely cartilaginous endoskeletons and are usually known only from isolated tooth plates or disarticulated bone fragments. There is thus a substantial temporal and evolutionary gap in our understanding of lungfish endoskeletal morphology, between the diverse and highly variable Devonian forms on the one hand and the three extant genera on the other. Here we present a virtual cranial endocast of Rhinodipterus kimberleyensis, from the Late Devonian Gogo Formation of Australia, one of the most derived fossil dipnoans with a well-ossified braincase. This endocast, generated from a Computed Microtomography (µCT) scan of the skull, is the first virtual endocast of any lungfish published, and only the third fossil dipnoan endocast to be illustrated in its entirety. Key features include long olfactory canals, a telencephalic cavity with a moderate degree of ventral expansion, large suparaotic cavities, and moderately enlarged utricular recesses. It has numerous similarities to the endocasts of Chirodipterus wildungensis and Griphognathus whitei, and to a lesser degree to 'Chirodipterus' australis and Dipnorhynchus sussmilchi. Among extant lungfish, it consistently resembles Neoceratodus more closely than Lepidosiren and Protopterus. Several trends in the evolution of the brains and labyrinth regions in dipnoans, such as the expansions of the utricular recess and telencephalic regions over time, are identified and discussed.

A morphospace for reef fishes: elongation is the dominant axis of body shape evolution

Tropical reef fishes are widely regarded as being perhaps the most morphologically diverse vertebrate assemblage on earth, yet much remains to be discovered about the scope and patterns of this diversity. We created a morphospace of 2,939 species spanning 56 families of tropical Indo-Pacific reef fishes and established the primary axes of body shape variation, the phylogenetic consistency of these patterns, and whether dominant patterns of shape change can be accomplished by diverse underlying changes. Principal component analysis showed a major axis of shape variation that contrasts deep-bodied species with slender, elongate forms. Furthermore, using custom methods to compare the elongation vector (axis that maximizes elongation deformation) and the main vector of shape variation (first principal component) for each family in the morphospace, we showed that two thirds of the families diversify along an axis of body elongation. Finally, a comparative analysis using a principal coordinate analysis based on the angles among first principal component vectors of each family shape showed that families accomplish changes in elongation with a wide range of underlying modifications. Some groups such as Pomacentridae and Lethrinidae undergo decreases in body depth with proportional increases in all body regions, while other families show disproportionate changes in the length of the head (e.g., Labridae), the trunk or caudal region in all combinations (e.g., Pempheridae and Pinguipedidae). In conclusion, we found that evolutionary changes in body shape along an axis of elongation dominates diversification in reef fishes. Changes in shape on this axis are thought to have immediate implications for swimming performance, defense from gape limited predators, suction feeding performance and access to some highly specialized habitats. The morphological modifications that underlie changes in elongation are highly diverse, suggesting a role for a range of developmental processes and functional consequences.

Historical factors that have shaped the evolution of tropical reef fishes: a review of phylogenies, biogeography, and remaining questions

Biodiversity patterns across the marine tropics have intrigued evolutionary biologists and ecologists alike. Tropical coral reefs host 1/3 of all marine species of fish on 0.1% of the ocean's surface. Yet our understanding of how mechanistic processes have underpinned the generation of this diversity is limited. However, it has become clear that the biogeographic history of the marine tropics has played an important role in shaping the diversity of tropical reef fishes we see today. In the last decade, molecular phylogenies and age estimation techniques have provided a temporal framework in which the ancestral biogeographic origins of reef fish lineages have been inferred, but few have included fully sampled phylogenies or made inferences at a global scale. We are currently at a point where new sequencing technologies are accelerating the reconstruction and the resolution of the Fish Tree of Life. How will a complete phylogeny of fishes benefit the study of biodiversity in the tropics? Here, I review the literature concerning the evolutionary history of reef-associated fishes from a biogeographic perspective. I summarize the major biogeographic and climatic events over the last 65 million years that have regionalized the tropical marine belt and what effect they have had on the molecular record of fishes and global biodiversity patterns. By examining recent phylogenetic trees of major reef associated groups, I identify gaps to be filled in order to obtain a clearer picture of the origins of coral reef fish assemblages. Finally, I discuss questions that remain to be answered and new approaches to uncover the mechanistic processes that underpin the evolution of biodiversity on coral reefs.

Drivers of parasite sharing among Neotropical freshwater fishes

Because host-parasite interactions are so ubiquitous, it is of primary interest for ecologists to understand the factors that generate, maintain and constrain these associations. Phylogenetic comparative studies have found abundant evidence for host-switching to relatively unrelated hosts, sometimes related to diversification events, in a variety of host-parasite systems. For Monogenoidea (Platyhelminthes) parasites, it has been suggested that the co-speciation model alone cannot explain host occurrences, hence host-switching and/or non-vicariant modes of speciation should be associated with the origins and diversification of several monogenoid taxa. The factors that shape broad patterns of parasite sharing were investigated using path analysis as a way to generate hypotheses about the origins of host-parasite interactions between monogenoid gill parasites and Neotropical freshwater fishes. Parasite sharing was assessed from an interaction matrix, and explanatory variables included phylogenetic relationships, environmental preferences, biological traits and geographic distribution for each host species. Although geographic distribution of hosts and host ecology are important factors to understand host-parasite interactions, especially within host lineages that share a relatively recent evolutionary history, phylogeny had the strongest overall direct effect on parasite sharing. Phylogenetic contiguity of host communities may allow a 'stepping-stone' mode of host-switching, which increases parasite sharing. Our results reinforce the importance of including evolutionary history in the study of ecological associations, including emerging infectious diseases risk assessment.

An evaluation of fossil tip-dating versus node-age calibrations in tetraodontiform fishes (Teleostei: Percomorphaceae)

Time-calibrated phylogenies based on molecular data provide a framework for comparative studies. Calibration methods to combine fossil information with molecular phylogenies are, however, under active development, often generating disagreement about the best way to incorporate paleontological data into these analyses. This study provides an empirical comparison of the most widely used approach based on node-dating priors for relaxed clocks implemented in the programs BEAST and MrBayes, with two recently proposed improvements: one using a new fossilized birth-death process model for node dating (implemented in the program DPPDiv), and the other using a total-evidence or tip-dating method (implemented in MrBayes and BEAST). These methods are applied herein to tetraodontiform fishes, a diverse group of living and extinct taxa that features one of the most extensive fossil records among teleosts. Previous estimates of time-calibrated phylogenies of tetraodontiforms using node-dating methods reported disparate estimates for their age of origin, ranging from the late Jurassic to the early Paleocene (ca. 150-59Ma). We analyzed a comprehensive dataset with 16 loci and 210 morphological characters, including 131 taxa (95 extant and 36 fossil species) representing all families of fossil and extant tetraodontiforms, under different molecular clock calibration approaches. Results from node-dating methods produced consistently younger ages than the tip-dating approaches. The older ages inferred by tip dating imply an unlikely early-late Jurassic (ca. 185-119Ma) origin for this order and the existence of extended ghost lineages in their fossil record. Node-based methods, by contrast, produce time estimates that are more consistent with the stratigraphic record, suggesting a late Cretaceous (ca. 86-96Ma) origin. We show that the precision of clade age estimates using tip dating increases with the number of fossils analyzed and with the proximity of fossil taxa to the node under assessment. This study suggests that current implementations of tip dating may overestimate ages of divergence in calibrated phylogenies. It also provides a comprehensive phylogenetic framework for tetraodontiform systematics and future comparative studies.

The infrabranchial musculature and its bearing on the phylogeny of percomorph fishes (Osteichthyes: Teleostei)

The muscles serving the ventral portion of the gill arches ( = infrabranchial musculature) are poorly known in bony fishes. A comparative analysis of the infrabranchial muscles in the major percomorph lineages reveals a large amount of phylogenetically-relevant information. Characters derived from this anatomical system are identified and discussed in light of current hypotheses of phylogenetic relationships among percomorphs. New evidence supports a sister-group relationship between the Batrachoidiformes and Lophiiformes and between the Callionymoidei and Gobiesocoidei. Investigated data also corroborate the existence of two monophyletic groups, one including the Pristolepididae, Badidae, and Nandidae, and a second clade consisting of all non-amarsipid stromateiforms. New synapomorphies are proposed for the Atherinomorphae, Blenniiformes, Lophiiformes, Scombroidei (including Sphyraenidae), and Gobiiformes. Within the latter order, the Rhyacichthyidae and Odontobutidae are supported as the successive sister families of all remaining gobiiforms. The present analysis further confirms the validity of infrabranchial musculature characters previously proposed to support the grouping of the Mugiliformes with the Atherinomorphae and the monophyly of the Labriformes with the possible inclusion of the Pholidichthyiformes. Interestingly, most hypotheses of relationships supported by the infrabranchial musculature have been advanced by preceding anatomists on the basis of distinct data sources, but were never recovered in recent molecular phylogenies. These conflicts clearly indicate the current unsatisfactory resolution of the higher-level phylogeny of percomorphs.

Creatine metabolism differs between mammals and rainbow trout (Oncorhynchus mykiss)

Creatine plays an important role in the cell as an energy buffer. As the energy system is a basic element of the organism it may possibly contribute to differences between rainbow trout strains selected for the traits growth and robustness, respectively. The cDNA sequences of creatine-related genes encoding glycine amidinotransferase (GATM), guanidinoacetate N-methyltransferase (GAMT), creatine kinase muscle-type (CKM) and creatine transporter 1 (CT1, encoded by gene solute carrier family 6, member 8 (SLC6A8)) were characterized in rainbow trout. Transcripts of the respective genes were quantified in kidney, liver, brain and skeletal muscle in both trout strains that had been acclimated to different temperatures. Several differences between the compared trout strains were found as well as between temperatures indicating that the energy system may contribute to differences between both strains. In addition to that, the expression data showed clear differences between the creatine system in rainbow trout and mammals, as the spatial distribution of the enzyme-encoding gene expression was clearly different from the patterns described for mammals. In rainbow trout, creatine synthesis seems to take place to a big extent in the skeletal muscle.

Mitochondrial genomic investigation of flatfish monophyly

We present the first study to use whole mitochondrial genome sequences to examine phylogenetic affinities of the flatfishes (Pleuronectiformes). Flatfishes have attracted attention in evolutionary biology since the early history of the field because understanding the evolutionary history and patterns of diversification of the group will shed light on the evolution of novel body plans. Because recent molecular studies based primarily on DNA sequences from nuclear loci have yielded conflicting results, it is important to examine phylogenetic signal in different genomes and genome regions. We aligned and analyzed mitochondrial genome sequences from thirty-nine pleuronectiforms including nine that are newly reported here, and sixty-six non-pleuronectiforms (twenty additional clade L taxa [Carangimorpha or Carangimorpharia] and forty-six secondary outgroup taxa). The analyses yield strong support for clade L and weak support for the monophyly of Pleuronectiformes. The suborder Pleuronectoidei receives moderate support, and as with other molecular studies the putatively basal lineage of Pleuronectiformes, the Psettodoidei is frequently not most closely related to other pleuronectiforms. Within the Pleuronectoidei, the basal lineages in the group are poorly resolved, however several flatfish subclades receive consistent support. The affinities of Lepidoblepharon and Citharoides among pleuronectoids are particularly uncertain with these data.

Multiple groups of endogenous epsilon-like retroviruses conserved across primates

Several types of cancer in fish are caused by retroviruses, including those responsible for major outbreaks of disease, such as walleye dermal sarcoma virus and salmon swim bladder sarcoma virus. These viruses form a phylogenetic group often described as the epsilonretrovirus genus. Epsilon-like retroviruses have become endogenous retroviruses (ERVs) on several occasions, integrating into germ line cells to become part of the host genome, and sections of fish and amphibian genomes are derived from epsilon-like retroviruses. However, epsilon-like ERVs have been identified in very few mammals. We have developed a pipeline to screen full genomes for ERVs, and using this pipeline, we have located over 800 endogenous epsilon-like ERV fragments in primate genomes. Genomes from 32 species of mammals and birds were screened, and epsilon-like ERV fragments were found in all primate and tree shrew genomes but no others. These viruses appear to have entered the genome of a common ancestor of Old and New World monkeys between 42 million and 65 million years ago. Based on these results, there is an ancient evolutionary relationship between epsilon-like retroviruses and primates. Clearly, these viruses had the potential to infect the ancestors of primates and were at some point a common pathogen in these hosts. Therefore, this result raises questions about the potential of epsilonretroviruses to infect humans and other primates and about the evolutionary history of these retroviruses.

IMPORTANCE

Epsilonretroviruses are a group of retroviruses that cause several important diseases in fish. Retroviruses have the ability to become a permanent part of the DNA of their host by entering the germ line as endogenous retroviruses (ERVs), where they lose their infectivity over time but can be recognized as retroviruses for millions of years. Very few mammals are known to have epsilon-like ERVs; however, we have identified over 800 fragments of endogenous epsilon-like ERVs in the genomes of all major groups of primates, including humans. These viruses seem to have circulated and infected primate ancestors 42 to 65 million years ago. We are now interested in how these viruses have evolved and whether they have the potential to infect modern humans or other primates.

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.

Phylogenetic informativeness reconciles ray-finned fish molecular divergence times

Background

Discordance among individual molecular age estimates, or between molecular age estimates and the fossil record, is observed in many clades across the Tree of Life. This discordance is attributed to a variety of variables including calibration age uncertainty, calibration placement, nucleotide substitution rate heterogeneity, or the specified molecular clock model. However, the impact of changes in phylogenetic informativeness of individual genes over time on phylogenetic inferences is rarely analyzed. Using nuclear and mitochondrial sequence data for ray-finned fishes (Actinopterygii) as an example, we extend the utility of phylogenetic informativeness profiles to predict the time intervals when nucleotide substitution saturation results in discordance among molecular ages estimated.

Results

We demonstrate that even with identical calibration regimes and molecular clock methods, mitochondrial based molecular age estimates are systematically older than those estimated from nuclear sequences. This discordance is most severe for highly nested nodes corresponding to more recent (i.e., Jurassic-Recent) divergences. By removing data deemed saturated, we reconcile the competing age estimates and highlight that the older mtDNA based ages were driven by nucleotide saturation.

Conclusions

Homoplasious site patterns in a DNA sequence alignment can systematically bias molecular divergence time estimates. Our study demonstrates that PI profiles can provide a non-arbitrary criterion for data exclusion to mitigate the influence of homoplasy on time calibrated branch length estimates. Analyses of actinopterygian molecular clocks demonstrate that scrutiny of the time scale on which sequence data is informative is a fundamental, but generally overlooked, step in molecular divergence time estimation.

The genome and linkage map of the northern pike (Esox lucius): conserved synteny revealed between the salmonid sister group and the Neoteleostei

The northern pike is the most frequently studied member of the Esociformes, the closest order to the diverse and economically important Salmoniformes. The ancestor of all salmonids purportedly experienced a whole-genome duplication (WGD) event, making salmonid species ideal for studying the early impacts of genome duplication while complicating their use in wider analyses of teleost evolution. Studies suggest that the Esociformes diverged from the salmonid lineage prior to the WGD, supporting the use of northern pike as a pre-duplication outgroup. Here we present the first genome assembly, reference transcriptome and linkage map for northern pike, and evaluate the suitability of this species to provide a representative pre-duplication genome for future studies of salmonid and teleost evolution. The northern pike genome sequence is composed of 94,267 contigs (N50 = 16,909 bp) contained in 5,688 scaffolds (N50 = 700,535 bp); the total scaffolded genome size is 878 million bases. Multiple lines of evidence suggest that over 96% of the protein-coding genome is present in the genome assembly. The reference transcriptome was constructed from 13 tissues and contains 38,696 transcripts, which are accompanied by normalized expression data in all tissues. Gene-prediction analysis produced a total of 19,601 northern pike-specific gene models. The first-generation linkage map identifies 25 linkage groups, in agreement with northern pike's diploid karyotype of 2N = 50, and facilitates the placement of 46% of assembled bases onto linkage groups. Analyses reveal a high degree of conserved synteny between northern pike and other model teleost genomes. While conservation of gene order is limited to smaller syntenic blocks, the wider conservation of genome organization implies the northern pike exhibits a suitable approximation of a non-duplicated Protacanthopterygiian genome. This dataset will facilitate future studies of esocid biology and empower ongoing examinations of the Atlantic salmon and rainbow trout genomes by facilitating their comparison with other major teleost groups.

Mitogenomic phylogeny of the Percichthyidae and Centrarchiformes (Percomorphaceae): comparison with recent nuclear gene-based studies and simultaneous analysis

Delineation of the fish family Percichthyidae (Percomorphaceae) has a long and convoluted history, with recent morphological-based studies restricting species members to South American and Australian freshwater and catadromous temperate perches. Four recent nuclear gene-based phylogenetic studies, however, found that the Percichthyidae was not monophyletic and was nested within a newly discovered inter-familial clade of Percomorphaceae, the Centrarchiformes, which comprises the Centrarchidae and 12 other families. Here, we reexamined the systematics of the Percichthyidae and Centrarchiformes based on new mitogenomic information. Our mitogenomic results are globally congruent with the recent nuclear gene-based studies although the overall amount of phylogenetic signal of the mitogenome is lower. They do not support the monophyly of the Percichthyidae, because the catadromous genus Percalates is not exclusively related to the freshwater percichthyids. The Percichthyidae (minus Percalates) and Percalates belong to a larger clade, equivalent to the Centrarchiformes, but their respective sister groups are unresolved. Because all recent analyses recover a monophyletic Centrarchiformes but with substantially different intra-relationships, we performed a simultaneous analysis for a character set combining the mitogenome and 19 nuclear genes previously published, for 22 centrarchiform taxa. This analysis furthermore indicates that the Centrarchiformes are divided into three lineages and the superfamily Cirrhitoidea is monophyletic as well as the temperate and freshwater centrarchiform perch-like fishes. It also clarifies some of the relationships within the freshwater Percichthyidae.

A new Liopropoma sea bass (Serranidae, Epinephelinae, Liopropomini) from deep reefs off Curaçao, southern Caribbean, with comments on depth distributions of western Atlantic liopropomins

Collecting reef-fish specimens using a manned submersible diving to 300 m off Curaçao, southern Caribbean, is resulting in the discovery of numerous new fish species. The new Liopropoma sea bass described here differs from other western Atlantic members of the genus in having VIII, 13 dorsal-fin rays; a moderately indented dorsal-fin margin; a yellow-orange stripe along the entire upper lip; a series of approximately 13 white, chevron-shaped markings on the ventral portion of the trunk; and a reddish-black blotch on the tip of the lower caudal-fin lobe. The new species, with predominantly yellow body and fins, closely resembles the other two “golden basses” found together with it at Curaçao: L. aberrans and L. olneyi. It also shares morphological features with the other western Atlantic liopropomin genus, Bathyanthias. Preliminary phylogenetic data suggest that western Atlantic liopropomins, including Bathyanthias, are monophyletic with respect to Indo-Pacific Liopropoma, and that Bathyanthias is nested within Liopropoma, indicating a need for further study of the generic limits of Liopropoma. The phylogenetic data also suggest that western Atlantic liopropomins comprise three monophyletic clades that have overlapping depth distributions but different depth maxima (3–135 m, 30–150 m, 133–411 m). The new species has the deepest depth range (182–241 m) of any known western Atlantic Liopropoma species. Both allopatric and depth-mediated ecological speciation may have contributed to the evolution of western Atlantic Liopropomini.

Connectivity of vertebrate genomes: Paired-related homeobox (Prrx) genes in spotted gar, basal teleosts, and tetrapods

Teleost fish are important models for human biology, health, and disease. Because genome duplication in a teleost ancestor (TGD) impacts the evolution of teleost genome structure and gene repertoires, we must discriminate gene functions that are shared and ancestral from those that are lineage-specific in teleosts or tetrapods to accurately apply inferences from teleost disease models to human health. Generalizations must account both for the TGD and for divergent evolution between teleosts and tetrapods after the likely two rounds of genome duplication shared by all vertebrates. Progress in sequencing techniques provides new opportunities to generate genomic and transcriptomic information from a broad range of phylogenetically informative taxa that facilitate detailed understanding of gene family and gene function evolution. We illustrate here the use of new sequence resources from spotted gar (Lepisosteus oculatus), a rayfin fish that diverged from teleosts before the TGD, as well as RNA-Seq data from gar and multiple teleost lineages to reconstruct the evolution of the Paired-related homeobox (Prrx) transcription factor gene family, which is involved in the development of mesoderm and neural crest-derived mesenchyme. We show that for Prrx genes, the spotted gar genome and gene expression patterns mimic mammals better than teleosts do. Analyses force the seemingly paradoxical conclusion that regulatory mechanisms for the limb expression domains of Prrx genes existed before the evolution of paired appendages. Detailed evolutionary analyses like those reported here are required to identify fish species most similar to the human genome to optimally connect fish models to human gene functions in health and disease.

Dispersal capacity predicts both population genetic structure and species richness in reef fishes

Dispersal is a fundamental species characteristic that should directly affect both rates of gene flow among spatially distributed populations and opportunities for speciation. Yet no single trait associated with dispersal has been demonstrated to affect both micro- and macroevolutionary patterns of diversity across a diverse biological assemblage. Here, we examine patterns of genetic differentiation and species richness in reef fishes, an assemblage of over 7,000 species comprising approximately one-third of the extant bony fishes and over one-tenth of living vertebrates. In reef fishes, dispersal occurs primarily during a planktonic larval stage. There are two major reproductive and parental investment syndromes among reef fishes, and the differences between them have implications for dispersal: (1) benthic guarding fishes lay negatively buoyant eggs, typically guarded by the male parent, and from these eggs hatch large, strongly swimming larvae; in contrast, (2) pelagic spawning fishes release small floating eggs directly into the water column, which drift unprotected before small weakly swimming larvae hatch. Using phylogenetic comparative methods, we show that benthic guarders have significantly greater population structure than pelagic spawners and additionally that taxonomic families of benthic guarders are more species rich than families of pelagic spawners. Our findings provide a compelling case for the continuity between micro- and macroevolutionary processes of biological diversification and underscore the importance of dispersal-related traits in influencing the mode and tempo of evolution.

Mitochondrial-nuclear interactions: compensatory evolution or variable functional constraint among vertebrate oxidative phosphorylation genes?

Oxidative phosphorylation (OXPHOS), the major energy-producing pathway in aerobic organisms, includes protein subunits encoded by both mitochondrial (mt) and nuclear (nu) genomes. How these independent genomes have coevolved is a long-standing question in evolutionary biology. Although mt genes evolve faster than most nu genes, maintenance of OXPHOS structural stability and functional efficiency may involve correlated evolution of mt and nu OXPHOS genes. The nu OXPHOS genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt OXPHOS subunits with which they interact. Evolutionary rates of nu OXPHOS genes should, therefore, be higher than that of nu genes that are not involved in OXPHOS (nu non-OXPHOS). We tested the compensatory evolution hypothesis by comparing the evolutionary rates (synonymous substitution rate dS and nonsynonymous substitution rate dN) among 13 mt OXPHOS genes, 60 nu OXPHOS genes, and 77 nu non-OXPHOS genes in vertebrates (7 fish and 40 mammal species). The results from a combined analysis of all OXPHOS subunits fit the predictions of the hypothesis. However, results from two OXPHOS complexes did not fit this pattern when analyzed separately. We found that the d(N) of nu OXPHOS genes for "core" subunits (those involved in the major catalytic activity) was lower than that of "noncore" subunits, whereas there was no significant difference in d(N) between genes for nu non-OXPHOS and core subunits. This latter finding suggests that compensatory changes play a minor role in the evolution of OXPHOS genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits.

Diversification of substrate specificities in teleostei Fads2: characterization of Δ4 and Δ6Δ5 desaturases of Chirostoma estor

Currently existing data show that the capability for long-chain PUFA (LC-PUFA) biosynthesis in teleost fish is more diverse than in other vertebrates. Such diversity has been primarily linked to the subfunctionalization that teleostei fatty acyl desaturase (Fads)2 desaturases have undergone during evolution. We previously showed that Chirostoma estor, one of the few representatives of freshwater atherinopsids, had the ability for LC-PUFA biosynthesis from C18 PUFA precursors, in agreement with this species having unusually high contents of DHA. The particular ancestry and pattern of LC-PUFA biosynthesis activity of C. estor make this species an excellent model for study to gain further insight into LC-PUFA biosynthetic abilities among teleosts. The present study aimed to characterize cDNA sequences encoding fatty acyl elongases and desaturases, key genes involved in the LC-PUFA biosynthesis. Results show that C. estor expresses an elongase of very long-chain FA (Elovl)5 elongase and two Fads2 desaturases displaying Δ4 and Δ6/Δ5 specificities, thus allowing us to conclude that these three genes cover all the enzymatic abilities required for LC-PUFA biosynthesis from C18 PUFA. In addition, the specificities of the C. estor Fads2 enabled us to propose potential evolutionary patterns and mechanisms for subfunctionalization of Fads2 among fish lineages.

Expression evolution facilitated the convergent neofunctionalization of a sodium channel gene

Ion channels have played a substantial role in the evolution of novel traits across all of the domains of life. A fascinating example of a novel adaptation is the convergent evolution of electric organs in the Mormyroid and Gymnotiform electric fishes. The regulated currents that flow through ion channels directly generate the electrical signals which have evolved in these fish. Here, we investigated how the expression evolution of two sodium channel paralogs (Scn4aa and Scn4ab) influenced their convergent molecular evolution following the teleost-specific whole-genome duplication. We developed a reliable assay to accurately measure the expression stoichiometry of these genes and used this technique to analyze relative expression of the duplicate genes in a phylogenetic context. We found that before a major shift in expression from skeletal muscle and neofunctionalization in the muscle-derived electric organ, Scn4aa was first downregulated in the ancestors of both electric lineages. This indicates that underlying the convergent evolution of this gene, there was a greater propensity toward neofunctionalization due to its decreased expression relative to its paralog Scn4ab. We investigated another derived muscle tissue, the sonic organ of Porichthys notatus, and show that, as in the electric fishes, Scn4aa again shows a radical shift in expression away from the ancestral muscle cells into the evolutionarily novel muscle-derived tissue. This study presents evidence that expression downregulation facilitates neofunctionalization after gene duplication, a pattern that may often set the stage for novel trait evolution after gene duplication.

Deepwater Horizon crude oil impacts the developing hearts of large predatory pelagic fish

The Deepwater Horizon disaster released more than 636 million L of crude oil into the northern Gulf of Mexico. The spill oiled upper surface water spawning habitats for many commercially and ecologically important pelagic fish species. Consequently, the developing spawn (embryos and larvae) of tunas, swordfish, and other large predators were potentially exposed to crude oil-derived polycyclic aromatic hydrocarbons (PAHs). Fish embryos are generally very sensitive to PAH-induced cardiotoxicity, and adverse changes in heart physiology and morphology can cause both acute and delayed mortality. Cardiac function is particularly important for fast-swimming pelagic predators with high aerobic demand. Offspring for these species develop rapidly at relatively high temperatures, and their vulnerability to crude oil toxicity is unknown. We assessed the impacts of field-collected Deepwater Horizon (MC252) oil samples on embryos of three pelagic fish: bluefin tuna, yellowfin tuna, and an amberjack. We show that environmentally realistic exposures (1-15 µg/L total PAH) cause specific dose-dependent defects in cardiac function in all three species, with circulatory disruption culminating in pericardial edema and other secondary malformations. Each species displayed an irregular atrial arrhythmia following oil exposure, indicating a highly conserved response to oil toxicity. A considerable portion of Gulf water samples collected during the spill had PAH concentrations exceeding toxicity thresholds observed here, indicating the potential for losses of pelagic fish larvae. Vulnerability assessments in other ocean habitats, including the Arctic, should focus on the developing heart of resident fish species as an exceptionally sensitive and consistent indicator of crude oil impacts.

Two waves of colonization straddling the K-Pg boundary formed the modern reef fish fauna

Living reef fishes are one of the most diverse vertebrate assemblages on Earth. Despite its prominence and ecological importance, the origins and assembly of the reef fish fauna is poorly described. A patchy fossil record suggests that the major colonization of reef habitats must have occurred in the Late Cretaceous and early Palaeogene, with the earliest known modern fossil coral reef fish assemblage dated to 50 Ma. Using a phylogenetic approach, we analysed the early evolutionary dynamics of modern reef fishes. We find that reef lineages successively colonized reef habitats throughout the Late Cretaceous and early Palaeogene. Two waves of invasion were accompanied by increasing morphological convergence: one in the Late Cretaceous from 90 to 72 Ma and the other immediately following the end-Cretaceous mass extinction. The surge in reef invasions after the Cretaceous–Palaeogene boundary continued for 10 Myr, after which the pace of transitions to reef habitats slowed. Combined, these patterns match a classic niche-filling scenario: early transitions to reefs were made rapidly by morphologically distinct lineages and were followed by a decrease in the rate of invasions and eventual saturation of morphospace. Major alterations in reef composition, distribution and abundance, along with shifts in climate and oceanic currents, occurred during the Late Cretaceous and early Palaeogene interval. A causal mechanism between these changes and concurrent episodes of reef invasion remains obscure, but what is clear is that the broad framework of the modern reef fish fauna was in place within 10 Myr of the end-Cretaceous extinction.

Body size and geographic range do not explain long term variation in fish populations: a Bayesian phylogenetic approach to testing assembly processes in stream fish assemblages

We combine evolutionary biology and community ecology to test whether two species traits, body size and geographic range, explain long term variation in local scale freshwater stream fish assemblages. Body size and geographic range are expected to influence several aspects of fish ecology, via relationships with niche breadth, dispersal, and abundance. These traits are expected to scale inversely with niche breadth or current abundance, and to scale directly with dispersal potential. However, their utility to explain long term temporal patterns in local scale abundance is not known. Comparative methods employing an existing molecular phylogeny were used to incorporate evolutionary relatedness in a test for covariation of body size and geographic range with long term (1983 - 2010) local scale population variation of fishes in West Fork White River (Indiana, USA). The Bayesian model incorporating phylogenetic uncertainty and correlated predictors indicated that neither body size nor geographic range explained significant variation in population fluctuations over a 28 year period. Phylogenetic signal data indicated that body size and geographic range were less similar among taxa than expected if trait evolution followed a purely random walk. We interpret this as evidence that local scale population variation may be influenced less by species-level traits such as body size or geographic range, and instead may be influenced more strongly by a taxon's local scale habitat and biotic assemblages.

Deepwater Horizon crude oil impacts the developing hearts of large predatory pelagic fish

The Deepwater Horizon disaster released more than 636 million L of crude oil into the northern Gulf of Mexico. The spill oiled upper surface water spawning habitats for many commercially and ecologically important pelagic fish species. Consequently, the developing spawn (embryos and larvae) of tunas, swordfish, and other large predators were potentially exposed to crude oil-derived polycyclic aromatic hydrocarbons (PAHs). Fish embryos are generally very sensitive to PAH-induced cardiotoxicity, and adverse changes in heart physiology and morphology can cause both acute and delayed mortality. Cardiac function is particularly important for fast-swimming pelagic predators with high aerobic demand. Offspring for these species develop rapidly at relatively high temperatures, and their vulnerability to crude oil toxicity is unknown. We assessed the impacts of field-collected Deepwater Horizon (MC252) oil samples on embryos of three pelagic fish: bluefin tuna, yellowfin tuna, and an amberjack. We show that environmentally realistic exposures (1-15 µg/L total PAH) cause specific dose-dependent defects in cardiac function in all three species, with circulatory disruption culminating in pericardial edema and other secondary malformations. Each species displayed an irregular atrial arrhythmia following oil exposure, indicating a highly conserved response to oil toxicity. A considerable portion of Gulf water samples collected during the spill had PAH concentrations exceeding toxicity thresholds observed here, indicating the potential for losses of pelagic fish larvae. Vulnerability assessments in other ocean habitats, including the Arctic, should focus on the developing heart of resident fish species as an exceptionally sensitive and consistent indicator of crude oil impacts.

The origins of adipose fins: an analysis of homoplasy and the serial homology of vertebrate appendages

Adipose fins are appendages found on the dorsal midline between the dorsal and caudal fins in more than 6000 living species of teleost fishes. It has been consistently argued that adipose fins evolved once and have been lost repeatedly across teleosts owing to limited function. Here, we demonstrate that adipose fins originated repeatedly by using phylogenetic and anatomical evidence. This suggests that adipose fins are adaptive, although their function remains undetermined. To test for generalities in the evolution of form in de novo vertebrate fins, we studied the skeletal anatomy of adipose fins across 620 species belonging to 186 genera and 55 families. Adipose fins have repeatedly evolved endoskeletal plates, anterior dermal spines and fin rays. The repeated evolution of fin rays in adipose fins suggests that these fins can evolve new tissue types and increased structural complexity by expressing fin-associated developmental modules in these new territories. Patterns of skeletal elaboration differ between the various occurrences of adipose fins and challenge prevailing hypotheses for vertebrate fin origin. Adipose fins represent a powerful and, thus far, barely studied model for exploring the evolution of vertebrate limbs and the roles of adaptation and generative biases in morphological evolution.

Divergence in skeletal mass and bone morphology in antarctic notothenioid fishes

Although notothenioid fishes lack swim bladders, some species live temporarily or permanently in the water column. Given its relatively high density, skeletal mass is a key determinant of buoyancy. Notothenioids have reduced skeletal ossification, but there is little quantitative data on the phylogenetic distribution of this trait. We obtained dry skeletal masses for 54 specimens representing 20 species from six notothenioid families. Although comparative data are sparse, notothenioid skeletons comprise a smaller percentage of body mass, <3.5%, than those of three non-notothenioid perciforms. With relatively high skeletal mass, the non-Antarctic Bovichtus diacanthus is similar in skeletal mass to some non-notothenioids. Eleginops maclovinus, the non-Antarctic sister group of the Antarctic clade, has a relatively light skeleton (<2% of body mass) similar to many species in the Antarctic clade. Low skeletal mass is therefore a synapomorphy shared by Eleginops plus the Antarctic clade. We provide gross, histological, and micro-CT documentation of the structure and location of bone and cartilage in skulls, pectoral girdles, and vertebrae, with emphasis on the bovichtid B. diacanthus, the eleginopsid E. maclovinus, and the channichthyid Chaenodraco wilsoni. In Eleginops and the Antarctic clade, most bone is spongy and most species have persisting cartilage in the skull and appendicular skeleton. We also measured the relative size of the notochordal canal in adult vertebral centra of 38 species representing all eight families. There is considerable interspecific variation in this pedomorphic trait and all species show an ontogenetic reduction in the relative size of the canal. However, large persisting canals are present in adults of the Antarctic clade, especially in the nototheniids Pleuragramma and Aethotaxis and in a number of bathydraconid and channichthyid genera.

Phylogeography and demographic history of two widespread Indo-Pacific mudskippers (Gobiidae: Periophthalmus)

This study provides a first description of the phylogeographic patterns and evolutionary history of two species of the mudskipper genus Periophthalmus. These amphibious gobies are distributed throughout the whole Indo-Pacific region and Atlantic coast of Africa, in peritidal habitats of soft-bottom coastal ecosystems. Three sequence datasets of two widely distributed species, Periophthalmus argentilineatus and P. kalolo, were obtained by amplifying and sequencing two mtDNA markers (D-loop and 16S rDNA) and the nDNA rag1 region. The three datasets were then used to perform phylogeographic, demographic and population genetic analyses. Our results indicate that tectonic events and past climatic oscillations strongly contributed to shape present genetic differentiation, phylogeographic and demographic patterns. We found support for the monophyly of P. kalolo, and only shallow genetic differentiation between East-African and Indo-Malayan populations of this species. However, our collections of the morphospecies P. argentilineatus include three molecularly distinct lineages, one of them more closely related to P. kalolo. The presence of Miocenic timings for the most recent common ancestors of some of these morphologically similar clades, suggests the presence of strong stabilising selection in mudskippers' habitats. At population level, demographic analyses and palaeoecological records of mangrove ecosystems suggest that Pleistocene bottlenecks and expansion plus secondary contact events of the studied species were associated with recurrent sea transgressions during interglacials, and sea regressions or stable regimes during glacials, respectively.

Differential expression and ligand binding indicate alternative functions for zebrafish polymeric immunoglobulin receptor (pIgR) and a family of pIgR-like (PIGRL) proteins

The polymeric immunoglobulin (Ig) receptor (pIgR) is an integral transmembrane glycoprotein that plays an important role in the mammalian immune response by transporting soluble polymeric Igs across mucosal epithelial cells. Single pIgR genes, which are expressed in lymphoid organs including mucosal tissues, have been identified in several teleost species. A single pigr gene has been identified on zebrafish chromosome 2 along with a large multigene family consisting of 29 pigr-like (PIGRL) genes. Full-length transcripts from ten different PIGRL genes that encode secreted and putative inhibitory membrane-bound receptors have been characterized. Although PIGRL and pigr transcripts are detected in immune tissues, only PIGRL transcripts can be detected in lymphoid and myeloid cells. In contrast to pIgR which binds Igs, certain PIGRL proteins bind phospholipids. PIGRL transcript levels are increased after infection with Streptococcus iniae, suggesting a role for PIGRL genes during bacterial challenge. Transcript levels of PIGRL genes are decreased after infection with Snakehead rhabdovirus, suggesting that viral infection may suppress PIGRL function.

Boom and bust: ancient and recent diversification in bichirs (Polypteridae: Actinopterygii), a relictual lineage of ray-finned fishes

Understanding the history that underlies patterns of species richness across the Tree of Life requires an investigation of the mechanisms that not only generate young species-rich clades, but also those that maintain species-poor lineages over long stretches of evolutionary time. However, diversification dynamics that underlie ancient species-poor lineages are often hidden due to a lack of fossil evidence. Using information from the fossil record and time calibrated molecular phylogenies, we investigate the history of lineage diversification in Polypteridae, which is the sister lineage of all other ray-finned fishes (Actinopterygii). Despite originating at least 390 million years (Myr) ago, molecular timetrees support a Neogene origin for the living polypterid species. Our analyses demonstrate polypterids are exceptionally species depauperate with a stem lineage duration that exceeds 380 million years (Ma) and is significantly longer than the stem lineage durations observed in other ray-finned fish lineages. Analyses of the fossil record show an early Late Cretaceous (100.5-83.6 Ma) peak in polypterid genus richness, followed by 60 Ma of low richness. The Neogene species radiation and evidence for high-diversity intervals in the geological past suggest a "boom and bust" pattern of diversification that contrasts with common perceptions of relative evolutionary stasis in so-called "living fossils."

Are flatfishes (Pleuronectiformes) monophyletic?

All extant species of flatfish (order Pleuronectiformes) are thought to descend from a common ancestor, and therefore to represent a monophyletic group. This hypothesis is based largely on the dramatic bilateral asymmetry and associated ocular migration characteristics of all flatfish. Yet, molecular-based phylogenetic studies have been inconclusive on this premise. Support for flatfish monophyly has varied with differences in taxonomic and gene region sampling schemes. Notably, the genus Psettodes has been found to be more related to non-flatfishes than to other flatfishes in many recent studies. The polyphyletic nature of the Pleuronectiformes is often inferred to be the result of weak historical signal and/or artifact of phylogenetic inference due to a bias in the data. In this study, we address the question of pleuronectiform monophyly with a broad set of markers (from six phylogenetically informative nuclear loci) and inference methods designed to limit the influence of phylogenetic artifacts. Concomitant with a character-rich analytical strategy, an extensive taxonomic sampling of flatfish and potential close relatives is used to increase power and resolution. Results of our analyses are most consistent with a non-monophyletic Pleuronectiformes with Psettodes always being excluded. A fossil-calibrated Bayesian relaxed clock analysis estimates the age of Pleuronectoidei to be 73 Ma, and the time to most recent common ancestor of Pleuronectoidei, Psettodes, and other relative taxa to be 77 Ma. The ages are much older than the records of any fossil pleuronectiform currently recognized. We discuss our findings in the context of the available morphological evidence and discuss the compatibility of our molecular hypothesis with morphological data regarding extinct and extant flatfish forms.

Molecular and fossil evidence place the origin of cichlid fishes long after Gondwanan rifting

Cichlid fishes are a key model system in the study of adaptive radiation, speciation and evolutionary developmental biology. More than 1600 cichlid species inhabit freshwater and marginal marine environments across several southern landmasses. This distributional pattern, combined with parallels between cichlid phylogeny and sequences of Mesozoic continental rifting, has led to the widely accepted hypothesis that cichlids are an ancient group whose major biogeographic patterns arose from Gondwanan vicariance. Although the Early Cretaceous (ca 135 Ma) divergence of living cichlids demanded by the vicariance model now represents a key calibration for teleost molecular clocks, this putative split pre-dates the oldest cichlid fossils by nearly 90 Myr. Here, we provide independent palaeontological and relaxed-molecular-clock estimates for the time of cichlid origin that collectively reject the antiquity of the group required by the Gondwanan vicariance scenario. The distribution of cichlid fossil horizons, the age of stratigraphically consistent outgroup lineages to cichlids and relaxed-clock analysis of a DNA sequence dataset consisting of 10 nuclear genes all deliver overlapping estimates for crown cichlid origin centred on the Palaeocene (ca 65-57 Ma), substantially post-dating the tectonic fragmentation of Gondwana. Our results provide a revised macroevolutionary time scale for cichlids, imply a role for dispersal in generating the observed geographical distribution of this important model clade and add to a growing debate that questions the dominance of the vicariance paradigm of historical biogeography.

Speciation in fishes

The field of speciation has seen much renewed interest in the past few years, with theoretical and empirical advances that have moved it from a descriptive field to a predictive and testable one. The goal of this review is to provide a general background on research on speciation as it pertains to fishes. Three major components to the question are first discussed: the spatial, ecological and sexual factors that influence speciation mechanisms. We then move to the latest developments in the field of speciation genomics. Affordable and rapidly available, massively parallel sequencing data allow speciation studies to converge into a single comprehensive line of investigation, where the focus has shifted to the search for speciation genes and genomic islands of speciation. We argue that fish present a very diverse array of scenarios, making them an ideal model to study speciation processes.

Are characiform fishes Gondwanan in origin? Insights from a time-scaled molecular phylogeny of the Citharinoidei (Ostariophysi: Characiformes)

Fishes of the order Characiformes are a diverse and economically important teleost clade whose extant members are found exclusively in African and Neotropical freshwaters. Although their transatlantic distribution has been primarily attributed to the Early Cretaceous fragmentation of western Gondwana, vicariance has not been tested with temporal information beyond that contained in their fragmentary fossil record and a recent time-scaled phylogeny focused on the African family Alestidae. Because members of the suborder Citharinoidei constitute the sister lineage to the entire remaining Afro-Neotropical characiform radiation, we inferred a time-calibrated molecular phylogeny of citharinoids using a popular Bayesian approach to molecular dating in order to assess the adequacy of current vicariance hypotheses and shed light on the early biogeographic history of characiform fishes. Given that the only comprehensive phylogenetic treatment of the Citharinoidei has been a morphology-based analysis published over three decades ago, the present study also provided an opportunity to further investigate citharinoid relationships and update the evolutionary framework that has laid the foundations for the current classification of the group. The inferred chronogram is robust to changes in calibration priors and suggests that the origins of citharinoids date back to the Turonian (ca 90 Ma) of the Late Cretaceous. Most modern citharinoid genera, however, appear to have originated and diversified much more recently, mainly during the Miocene. By reconciling molecular-clock- with fossil-based estimates for the origins of the Characiformes, our results provide further support for the hypothesis that attributes the disjunct distribution of the order to the opening of the South Atlantic Ocean. The striking overlap in tempo of diversification and biogeographic patterns between citharinoids and the African-endemic family Alestidae suggests that their evolutionary histories could have been strongly and similarly influenced by Miocene geotectonic events that modified the landscape and produced the drainage pattern of Central Africa seen today.

Phylogeny and biogeography of a shallow water fish clade (Teleostei: Blenniiformes)

BACKGROUND

The Blenniiformes comprises six families, 151 genera and nearly 900 species of small teleost fishes closely associated with coastal benthic habitats. They provide an unparalleled opportunity for studying marine biogeography because they include the globally distributed families Tripterygiidae (triplefin blennies) and Blenniidae (combtooth blennies), the temperate Clinidae (kelp blennies), and three largely Neotropical families (Labrisomidae, Chaenopsidae, and Dactyloscopidae). However, interpretation of these distributional patterns has been hindered by largely unresolved inter-familial relationships and the lack of evidence of monophyly of the Labrisomidae.

RESULTS

We explored the phylogenetic relationships of the Blenniiformes based on one mitochondrial (COI) and four nuclear (TMO-4C4, RAG1, Rhodopsin, and Histone H3) loci for 150 blenniiform species, and representative outgroups (Gobiesocidae, Opistognathidae and Grammatidae). According to the consensus of Bayesian Inference, Maximum Likelihood, and Maximum Parsimony analyses, the monophyly of the Blenniiformes and the Tripterygiidae, Blenniidae, Clinidae, and Dactyloscopidae is supported. The Tripterygiidae is the sister group of all other blennies, and the Blenniidae is the sister group of the remaining blennies. The monophyly of the Labrisomidae is supported with the exclusion of the Cryptotremini and inclusion of Stathmonotus, and we elevate two subgenera of Labrisomus to establish a monophyletic classification within the family. The monophyly of the Chaenopsidae is supported with the exclusion of Stathmonotus (placed in the Stathmonotini) and Neoclinus and Mccoskerichthys (placed in the Neoclinini). The origin of the Blenniiformes was estimated in the present-day IndoPacific region, corresponding to the Tethys Sea approximately 60.3 mya. A largely Neotropical lineage including the Labrisomidae, Chaenopsidae and Dactyloscopidae (node IV) evolved around 37.6 mya when the Neotropics were increasingly separated from the IndoPacific, but well before the closure of the Tethys Sea.

CONCLUSIONS

Relationships recovered in this study are similar to those of earlier analyses within the Clinidae and Chaenopsidae, and partially similar within the Blenniidae, but tripterygiid relationships remain poorly resolved. We present the first comprehensive phylogenetic hypothesis for a monophyletic Labrisomidae with five tribes (Labrisomini, Mnierpini, Paraclinini, Stathmonotini and Starksiini). Global distributions of blenny genera included in our analysis support the evolution of a largely Neotropical clade whose closest relatives (clinids and cryptotremines) are temperate in distribution.