Diversity and disparity through time in the adaptive radiation of Antarctic notothenioid fishes

Convergent evolution in shape in European lineages of gobies

During their radiation, certain groups of animals evolved significant phenotypic disparity (morphological diversity), enabling them to thrive in diverse environments. Adaptations to the same type of environment can lead to convergent evolution in function and morphology. However, well-documented examples in repeated adaptations of teleost fishes to different habitats, which are not primarily related to trophic specialization, are still scarce. Gobies are a remarkable fish group, exhibiting a great species diversity, morphological variability, and extraordinary ability to colonize very different environments. A variety of lifestyles and body forms evolved also in European lineages of gobies. We conducted two-dimensional geometric morphometric and phylomorphospace analyses in European lineages of gobies and evaluated the extent of convergent evolution in shape associated with adaptation to various habitats. Our analyses revealed the change in shape along the nektonic-cryptobenthic axis, from very slender head and body to stout body and wide head. We showed convergent evolution related to mode of locomotion in the given habitat in four ecological groups: nektonic, hyperbenthic, cryptobenthic, and freshwater gobies. Gobies, therefore, emerge as a highly diversified lineage with unique lifestyle variations, offering invaluable insights into filling of ecomorphological space and mechanisms of adaptation to various aquatic environments with distinct locomotion requirements.

The complete mitochondrial genome of Trematomus hansoni Boulenger, 1902 (Perciformes, Nototheniidae)

The striped notothen Trematomus hansoni is an Antarctic fish species belonging to the family Nototheniidae (cod icefishes) that is distributed throughout the Southern Ocean. In this study, the complete mitochondrial genome of T. hansoni was sequenced using an Illumina MiSeq platform. The circular mitochondrial genome is 19,218 bp long and contains 13 protein-coding genes, 23 tRNA genes, two rRNA genes, and one control region. Notably, there are two trnG-UCC genes and the second gene, located between trnE-UUC and trnI-GAU, has no D-arm structure. The base composition is 56.18% of A + T and 43.82% of G + C. The phylogenetic analysis supports that T. hansoni is grouped into a single clade with T. bernacchii. This study will be a valuable resource for further research on the phylogeny and evolution of the genus Trematomus.

Morphology and diet are decoupled in nearshore notothenoids from King George Island, West Antarctica

Antarctic notothenioid fishes show wide adaptive morphological radiation, linked to habitat preferences and food composition. However, direct comparisons of phenotypic variability and feeding habits are still lacking, particularly in stages inhabiting nearshore areas. To assess these relationships, we collected juveniles and adults of the most common benthic species inhabiting shallow waters off the South Shetland Islands within a similar size range, the plunderfish Harpagifer antarcticus, the black rockcod Notothenia coriiceps, and the marbled rockcod Notothenia rossii. Individual size ranges varied from 44.0 to 98.9 mm standard length (LS) (H. antarcticus), from 95.8 to 109.3 mm LS (N. coriiceps), and from 63.0 to 113.0 mm LS (N. rossii). Notothenioid fish showed different morphospace variability, being larger for H. antarcticus than the other Notothenia species and associated with the position of the posterior end of the operculum, along with the location and relative size of the eye. The evolutionary allometry was low, but the static allometry was much higher, especially for H. antarcticus and N. rossii. The diet was mainly carnivorous, consisting of amphipods and euphausiids. Macroalgae were scarce or totally absent in the gut contents of all species. Only H. antarcticus showed an increase in the prey number and ingested prey volume with fish size. Finally, there was a significant covariation between shape changes and LS in all species (allometric effects), however, not with prey composition, probably due to the small size range or ontogenetic stage and the relative similarity (or lack of contrast) in the benthic environment that they utilized.

The hierarchical radiation of phyllostomid bats as revealed by adaptive molar morphology

Adaptive radiations are bursts in biodiversity that generate new evolutionary lineages and phenotypes. However, because they typically occur over millions of years, it is unclear how their macroevolutionary dynamics vary through time and among groups of organisms. Phyllostomid bats radiated extensively for diverse diets-from insects to vertebrates, fruit, nectar, and blood-and we use their molars as a model system to examine the dynamics of adaptive radiations. Three-dimensional shape analyses of lower molars of Noctilionoidea (Phyllostomidae and close relatives) indicate that different diet groups exhibit distinct morphotypes. Comparative analyses further reveal that phyllostomids are a striking example of a hierarchical radiation; phyllostomids' initial, higher-level diversification involved an "early burst" in molar morphological disparity as lineages invaded new diet-affiliated adaptive zones, followed by subsequent lower-level diversifications within adaptive zones involving less dramatic morphological changes. We posit that strong selective pressures related to initial shifts to derived diets may have freed molars from morpho-functional constraints associated with the ancestral molar morphotype. Then, lineages with derived diets (frugivores and nectarivores) diversified within broad adaptive zones, likely reflecting finer-scale niche partitioning. Importantly, the observed early burst pattern is only evident when examining molar traits that are strongly linked to diet, highlighting the value of ecomorphological traits in comparative studies. Our results support the hypothesis that adaptive radiations are commonly hierarchical and involve different tempos and modes at different phylogenetic levels, with early bursts being more common at higher levels.

De novo assembly and annotation of the Patagonian toothfish (Dissostichus eleginoides) genome

BACKGROUND

Patagonian toothfish (Dissostichus eleginoides) is an economically and ecologically important fish species in the family Nototheniidae. Juveniles occupy progressively deeper waters as they mature and grow, and adults have been caught as deep as 2500 m, living on or in just above the southern shelves and slopes around the sub-Antarctic islands of the Southern Ocean. As apex predators, they are a key part of the food web, feeding on a variety of prey, including krill, squid, and other fish. Despite its importance, genomic sequence data, which could be used for more accurate dating of the divergence between Patagonian and Antarctic toothfish, or establish whether it shares adaptations to temperature with fish living in more polar or equatorial climes, has so far been limited.

RESULTS

A high-quality D. eleginoides genome was generated using a combination of Illumina, PacBio and Omni-C sequencing technologies. To aid the genome annotation, the transcriptome derived from a variety of toothfish tissues was also generated using both short and long read sequencing methods. The final genome assembly was 797.8 Mb with a N50 scaffold length of 3.5 Mb. Approximately 31.7% of the genome consisted of repetitive elements. A total of 35,543 putative protein-coding regions were identified, of which 50% have been functionally annotated. Transcriptomics analysis showed that approximately 64% of the predicted genes (22,617 genes) were found to be expressed in the tissues sampled. Comparative genomics analysis revealed that the anti-freeze glycoprotein (AFGP) locus of D. eleginoides does not contain any AFGP proteins compared to the same locus in the Antarctic toothfish (Dissostichus mawsoni). This is in agreement with previously published results looking at hybridization signals and confirms that Patagonian toothfish do not possess AFGP coding sequences in their genome.

CONCLUSIONS

We have assembled and annotated the Patagonian toothfish genome, which will provide a valuable genetic resource for ecological and evolutionary studies on this and other closely related species.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

Genomics of Secondarily Temperate Adaptation in the Only Non-Antarctic Icefish

White-blooded Antarctic icefishes, a family within the adaptive radiation of Antarctic notothenioid fishes, are an example of extreme biological specialization to both the chronic cold of the Southern Ocean and life without hemoglobin. As a result, icefishes display derived physiology that limits them to the cold and highly oxygenated Antarctic waters. Against these constraints, remarkably one species, the pike icefish Champsocephalus esox, successfully colonized temperate South American waters. To study the genetic mechanisms underlying secondarily temperate adaptation in icefishes, we generated chromosome-level genome assemblies of both C. esox and its Antarctic sister species, Champsocephalus gunnari. The C. esox genome is similar in structure and organization to that of its Antarctic congener; however, we observe evidence of chromosomal rearrangements coinciding with regions of elevated genetic divergence in pike icefish populations. We also find several key biological pathways under selection, including genes related to mitochondria and vision, highlighting candidates behind temperate adaptation in C. esox. Substantial antifreeze glycoprotein (AFGP) pseudogenization has occurred in the pike icefish, likely due to relaxed selection following ancestral escape from Antarctica. The canonical AFGP locus organization is conserved in C. esox and C. gunnari, but both show a translocation of two AFGP copies to a separate locus, previously unobserved in cryonotothenioids. Altogether, the study of this secondarily temperate species provides an insight into the mechanisms underlying adaptation to ecologically disparate environments in this otherwise highly specialized group.

The mitochondrial genome of the red icefish (Channichthys rugosus) casts doubt on its species status

Antarctic notothenioid fishes are recognised as one of the rare examples of adaptive radiation in the marine system. Withstanding the freezing temperatures of Antarctic waters, these fishes have diversified into over 100 species within no more than 10–20 million years. However, the exact species richness of the radiation remains contested. In the genus Channichthys, between one and nine species are recognised by different authors. To resolve the number of Channichthys species, genetic information would be highly valuable; however, so far, only sequences of a single species, C. rhinoceratus, are available. Here, we present the nearly complete sequence of the mitochondrial genome of C. rugosus, obtained from a formalin-fixed museum specimen sampled in 1974. This sequence differs from the mitochondrial genome of C. rhinoceratus in no more than 27 positions, suggesting that the two species may be synonymous.

Environmental correlates of phenotypic evolution in ecologically diverse Liolaemus lizards

Evolutionary correlations between phenotypic and environmental traits characterize adaptive radiations. However, the lizard genus Liolaemus, one of the most ecologically diverse terrestrial vertebrate radiations on earth, has so far shown limited or mixed evidence of adaptive diversification in phenotype. Restricted use of comprehensive environmental data, incomplete taxonomic representation and not considering phylogenetic uncertainty may have led to contradictory evidence. We compiled a 26-taxon dataset for the Liolaemus gracilis species group, representing much of the ecological diversity represented within Liolaemus and used environmental data to characterize how environments occupied by species' relate to phenotypic evolution. Our analyses, explicitly accounting for phylogenetic uncertainty, suggest diversification in phenotypic traits toward the present, with body shape evolution rapidly evolving in this group. Body shape evolution correlates with the occupation of different structural habitats indicated by vegetation axes suggesting species have adapted for maximal locomotory performance in these habitats. Our results also imply that the effects of phylogenetic uncertainty and model misspecification may be more extensive on univariate, relative to multivariate analyses of evolutionary correlations, which is an important consideration in analyzing data from rapidly radiating adaptive radiations.

Adaptations and Diversity of Antarctic Fishes: A Genomic Perspective

Antarctic notothenioid fishes are the classic example of vertebrate adaptive radiation in a marine environment. Notothenioids diversified from a single common ancestor ∼25 Mya to more than 140 species today, and they represent ∼90% of fish biomass on the continental shelf of Antarctica. As they diversified in the cold Southern Ocean, notothenioids evolved numerous traits, including osteopenia, anemia, cardiomegaly, dyslipidemia, and aglomerular kidneys, that are beneficial or tolerated in their environment but are pathological in humans. Thus, notothenioids are models for understanding adaptive radiations, physiological and biochemical adaptations to extreme environments, and genetic mechanisms of human disease. Since 2014, 16 notothenioid genomes have been published, which enable a first-pass holistic analysis of the notothenioid radiation and the genetic underpinnings of novel notothenioid traits. Here, we review the notothenioid radiation from a genomic perspective and integrate our insights with recent observations from other fish radiations. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The evolutionary puzzle solution for the origins of the partial loss of the Cτ2 exon in notothenioid fishes

Cryonotothenioidea is the main group of fishes that thrive in the extremely cold Antarctic environment, thanks to the acquisition of peculiar morphological, physiological and molecular adaptations. We have previously disclosed that IgM, the main immunoglobulin isotype in teleosts, display typical cold-adapted features. Recently, we have analyzed the gene encoding the heavy chain constant region (CH) of the IgT isotype from the Antarctic teleost Trematomus bernacchii (family Nototheniidae), characterized by the near-complete deletion of the CH2 domain. Here, we aimed to track the loss of the CH2 domain along notothenioid phylogeny and to identify its ancestral origins. To this end, we obtained the IgT gene sequences from several species belonging to the Antarctic families Nototheniidae, Bathydraconidae and Artedidraconidae. All species display a CH2 remnant of variable size, encoded by a short Cτ2 exon, which retains functional splicing sites and therefore is included in the mature transcript. We also considered representative species from the three non-Antarctic families: Eleginopsioidea (Eleginops maclovinus), Pseudaphritioidea (Pseudaphritis urvillii) and Bovichtidae (Bovichtus diacanthus and Cottoperca gobio). Even though only E. maclovinus, the sister taxa of Cryonotothenioidea, shared the partial loss of Cτ2, the other non-Antarctic notothenioid species displayed early molecular signatures of this event. These results shed light on the evolutionary path that underlies the origins of this remarkable gene structural modification.

Phylogenomic species delimitation dramatically reduces species diversity in an Antarctic adaptive radiation

Application of genetic data to species delimitation often builds confidence in delimitations previously hypothesized using morphological, ecological, and geographic data and frequently yields recognition of previously-undescribed cryptic diversity. However, a recent critique of genomic data-based species delimitation approaches is that they have the potential to conflate population structure with species diversity, resulting in taxonomic oversplitting. The need for an integrative approach to species delimitation, in which molecular, morphological, ecological, and geographic lines of evidence are evaluated together, is becoming increasingly apparent. Here, we integrate phylogenetic, population genetic, and coalescent analyses of genome-wide sequence data with investigation of variation in multiple morphological traits to delimit species within the Antarctic barbeled plunderfishes (Artedidraconidae: Pogonophryne). Pogonophryne currently comprises 29 valid species, most of which are distinguished solely by variation in ornamentation of the mental barbel that projects from the lower jaw, a structure previously shown to vary widely within a single species. However, our genomic and phenotypic analyses result in a dramatic reduction in the number of distinct species recognized within the clade, providing evidence to support the recognition of no more than six species. We propose to synonymize 24 of the currently recognized species with five species of Pogonophryne. We find genomic and phenotypic evidence for a new species of Pogonophryne from specimens collected in the Ross Sea. Our findings represent a rare example in which application of molecular data provides evidence of taxonomic oversplitting on the basis of morphology, clearly demonstrating the utility of an integrative species delimitation framework.

Species Radiations in the Sea: What the Flock?

Species flocks are proliferations of closely-related species, usually after colonization of depauperate habitat. These radiations are abundant on oceanic islands and in ancient freshwater lakes, but rare in marine habitats. This contrast is well documented in the Hawaiian Archipelago, where terrestrial examples include the speciose silverswords (sunflower family Asteraceae), Drosophila fruit flies, and honeycreepers (passerine birds), all derived from one or a few ancestral lineages. The marine fauna of Hawai'i is also the product of rare colonization events, but these colonizations usually yield only one species. Dispersal ability is key to understanding this evolutionary inequity. While terrestrial fauna rarely colonize between oceanic islands, marine fauna with pelagic larvae can make this leap in every generation. An informative exception is the marine fauna that lack a pelagic larval stage. These low-dispersal species emulate a "terrestrial" mode of reproduction (brooding, viviparity, crawl-away larvae), yielding marine species flocks in scattered locations around the world. Elsewhere, aquatic species flocks are concentrated in specific geographic settings, including the ancient lakes of Baikal (Siberia) and Tanganyika (eastern Africa), and Antarctica. These locations host multiple species flocks across a broad taxonomic spectrum, indicating a unifying evolutionary phenomenon. Hence marine species flocks can be singular cases that arise due to restricted dispersal or other intrinsic features, or they can be geographically clustered, promoted by extrinsic ecological circumstances. Here, we review and contrast intrinsic cases of species flocks in individual taxa, and extrinsic cases of geological/ecological opportunity, to elucidate the processes of species radiations.

Exploring adaptive landscapes across deep time: A case study using echinoid body size

Adaptive landscapes are a common way of conceptualizing the phenotypic evolution of lineages across deep time. Although multiple approaches exist to implement this concept into operational models of trait evolution, inferring adaptive landscapes from comparative datasets remains challenging. Here, I explore the macroevolutionary dynamics of echinoid body size using data from over 5000 specimens and a phylogenetic framework incorporating a dense fossil sampling and spanning approximately 270 million years. Furthermore, I implement a novel approach of exploring alternative parameterizations of adaptive landscapes that succeeds in finding simpler, yet better-fitting models. Echinoid body size has been constrained to evolve within a single adaptive optimum for much of the clade's history. However, most of the morphological disparity of echinoids was generated by multiple regime shifts that drove the repeated evolution of miniaturized and gigantic forms. Events of body size innovation occurred predominantly in the Late Cretaceous and were followed by a drastic slowdown following the Cretaceous-Paleogene mass extinction. The discovery of these patterns is contingent upon directly sampling fossil taxa. The macroevolution of echinoid body size is therefore characterized by a late increase in disparity (likely linked to an expansion of ecospace), generated by active processes driving lineages toward extreme morphologies.

Delayed adaptive radiation among New Zealand stream fishes: joint estimation of divergence time and trait evolution in a newly delineated island species flock

Adaptive radiations are generally thought to occur soon after a lineage invades a region offering high levels of ecological opportunity. However, few adaptive radiations beyond a handful of exceptional examples are known, so a comprehensive understanding of their dynamics is still lacking. Here, we present a novel case of an island species flock of freshwater fishes with a radically different tempo of adaptive history than that found in many popular evolutionary model systems. Using a phylogenomic data set combined with simultaneous Bayesian estimation of divergence times and trait-based speciation and extinction models, we show that the New Zealand Gobiomorphus gudgeons comprise a monophyletic assemblage, but surprisingly, the radiation did not fully occupy freshwater habitats and explosively speciate until more than 10 myr after the lineage invaded the islands. This shift in speciation rate was not accompanied by an acceleration in the rate of morphological evolution in the freshwater crown clade relative to the other species, but is correlated with a reduction in head pores and scales as well as an increase in egg size. Our results challenge the notion that clades always rapidly exploit ecological opportunities in the absence of competing lineages. Instead, we demonstrate that adaptive radiation can experience a slow start before undergoing accelerated diversification and that lineage and phenotypic diversification may be uncoupled in young radiations. [Adaptive radiation; Eleotridae; freshwater; Gobiomorphus; New Zealand.]

Toxin expression in snake venom evolves rapidly with constant shifts in evolutionary rates

Key innovations provide ecological opportunity by enabling access to new resources, colonization of new environments, and are associated with adaptive radiation. The most well-known pattern associated with adaptive radiation is an early burst of phenotypic diversification. Venoms facilitate prey capture and are widely believed to be key innovations leading to adaptive radiation. However, few studies have estimated their evolutionary rate dynamics. Here, we test for patterns of adaptive evolution in venom gene expression data from 52 venomous snake species. By identifying shifts in tempo and mode of evolution along with models of phenotypic evolution, we show that snake venom exhibits the macroevolutionary dynamics expected of key innovations. Namely, all toxin families undergo shifts in their rates of evolution, likely in response to changes in adaptive optima. Furthermore, we show that rapid-pulsed evolution modelled as a Lévy process better fits snake venom evolution than conventional early burst or Ornstein-Uhlenbeck models. While our results support the idea of snake venom being a key innovation, the innovation of venom chemistry lacks clear mechanisms that would lead to reproductive isolation and thus adaptive radiation. Therefore, the extent to which venom directly influences the diversification process is still a matter of contention.

Multiple independent chromosomal fusions accompanied the radiation of the Antarctic teleost genus Trematomus (Notothenioidei:Nototheniidae)

Background

Chromosomal rearrangements are thought to be an important driving force underlying lineage diversification, but their link to speciation continues to be debated. Antarctic teleost fish of the family Nototheniidae (Notothenioidei) diversified in a changing environmental context, which led to ecological, morphological, and genetic differentiation among populations. In addition, extensive chromosomal repatterning accompanied species divergence in several clades. The most striking karyotypic changes involved the recent species radiation (about 10 My) of the genus Trematomus, with chromosomal pair numbers ranging between 29 and 12. These dramatic reductions in chromosome number resulted mostly from large-scale chromosome fusions. Multiple centric and/or tandem fusions have been hypothesized in at least seven of the twelve recognized Trematomus species. To reconstruct their evolutionary history, we employed comparative cytogenomics (BAC-FISH and chromosome painting) to reveal patterns of interspecific chromosomal orthologies across several notothenioid clades.

Results

We defined orthologous chromosomal segments of reference, termed Structural Units (SUs). SUs were identified in a total of 18 notothenioid species. We demonstrated for the first time that SUs were strongly conserved across every specimen examined, with chromosomal syntenies highlighting a paucity of intrachromosomal macro-rearrangements. Multiple independent fusions of these SUs were inferred in the Trematomus species, in contrast to the shared SU fusions in species of the sister lineage Notothenia.

Conclusions

The SU segments were defined units of chromosomal rearrangement in the entire family Nototheiidae, which diverged from the other notothenioid families 20 My ago. Some of the identified chromosomal syntenies within the SUs were even conserved in their closest relatives, the family Eleginopsidae. Comparing the timing of acquisition of the fusions in the closely related genera Notothenia and Trematomus of the nototheniid species family, we conclude that they exhibit distinct chromosomal evolutionary histories, which may be relevant to different speciation scenarios.

Cytogenetic characterization of the Antarctic silverfish Pleuragramma antarctica (Boulenger 1902) through analysis of mitotic chromosomes from early larvae

This paper describes the cytogenetic features of the Antarctic silverfish Pleuragramma antarctica (Boulenger 1902), a keystone species of the Antarctic coastal marine ecosystem. Conventional cytogenetic analyses and physical mapping of repetitive DNA sequences were performed on metaphase plates obtained through direct chromosome preparation from P. antarctica early larvae. The Antarctic silverfish have a diploid number (2n) = 48, and a karyotype made up of a majority of two-armed chromosomes (karyotype formula36m/sm + 10st + 2a, fundamental number = 94). Major ribosomal gene repeats were detected on three chromosome pairs (20, 21, and 23), in correspondence of dim DAPI stained regions. Long Interspersed Nuclear Elements (LINEs) were abundant and wide spread over all chromosomes. Overall, the cytogenetic data presented herein are consistent with a long independent cytogenetic and evolutionary history for the species. The large number of two-armed chromosomes, indicative of highly-rearranged karyotype, coupled with a diploid number of 48, a presumed primitive character for this fish group, and the spread of the major ribosomal genes on three chromosome pairs, make the Antarctic silverfish distinct from all other notothenioid species.

Selective Sampling of Species and Fossils Influences Age Estimates Under the Fossilized Birth-Death Model

The fossilized birth-death (FBD) model allows the estimation of species divergence times from molecular and fossil information in a coherent framework of diversification and fossil sampling. Some assumptions of the FBD model, however, are difficult to meet in phylogenetic analyses of highly diverse groups. Here, I use simulations to assess the impact of extreme model violations, including diversified sampling of species and the exclusive use of the oldest fossils per clade, on divergence times estimated with the FBD model. My results demonstrate that selective sampling of fossils can produce dramatically overestimated divergence times when the FBD model is used for inference, due to an interplay of underestimates for the model parameters net diversification rate, turnover, and fossil-sampling proportion. In contrast, divergence times estimated with CladeAge, a method that uses information about the oldest fossils per clade together with estimates of sampling and diversification rates, are accurate under these conditions. Practitioners of Bayesian divergence-time estimation should therefore ensure that the dataset conforms to the expectations of the FBD model, or estimates of sampling and diversification rates should be obtained a priori so that CladeAge can be used for the inference.

Evidence that DNA repair genes, a family of tumor suppressor genes, are associated with evolution rate and size of genomes

Adaptive radiation and evolutionary stasis are characterized by very different evolution rates. The main aim of this study was to investigate if any genes have a special role to a high or low evolution rate. The availability of animal genomes permitted comparison of gene content of genomes of 24 vertebrate species that evolved through adaptive radiation (representing high evolutionary rate) and of 20 vertebrate species that are considered as living fossils (representing a slow evolutionary rate or evolutionary stasis). Mammals, birds, reptiles, and bony fishes were included in the analysis. Pathway analysis was performed for genes found to be specific in adaptive radiation or evolutionary stasis respectively. Pathway analysis revealed that DNA repair and cellular response to DNA damage are important (false discovery rate = 8.35 × 10⁻⁵; 7.15 × 10⁻⁶, respectively) for species evolved through adaptive radiation. This was confirmed by further genetic in silico analysis (p = 5.30 × 10⁻³). Nucleotide excision repair and base excision repair were the most significant pathways. Additionally, the number of DNA repair genes was found to be linearly related to the genome size and the protein number (proteome) of the 44 animals analyzed (p < 1.00 × 10⁻⁴), this being compatible with Drake’s rule. This is the first study where radiated and living fossil species have been genetically compared. Evidence has been found that cancer-related genes have a special role in radiated species. Linear association of the number of DNA repair genes with the species genome size has also been revealed. These comparative genetics results can support the idea of punctuated equilibrium evolution.

Intergeneric hybrids inform reproductive isolating barriers in the Antarctic icefish radiation

Interspecific hybridization or barriers to hybridization may have contributed to the diversification of Antarctic icefishes (Channichthyidae), but data supporting these hypotheses is scarce. To understand the potential for hybridization and to investigate reproductive isolating mechanisms among icefish species, we performed in vitro fertilization experiments using eggs from a female blackfin icefish Chaenocephalus aceratus and sperm from a male of another genera, the ocellated icefish Chionodraco rastrospinosus. Sequencing of genomic and mitochondrial DNA confirmed the intergeneric hybrid nature of resulting embryos which successfully developed and hatched as active larvae at about four and a half months during the Antarctic winter. This result demonstrates the compatibility of gametes of these two species and the viability of resulting zygotes and larvae. Due to logistic constraints and the slow developmental rate of icefishes, we could not test for long-term hybrid viability, fertility, fitness, or hybrid breakdown. Analysis of our fishing records and available literature, however, suggests that the strongest barriers to hybridization among parapatric icefish species are likely to be behavioral and characterized by assortative mating and species-specific courtship and nesting behaviors. This conclusion suggests that, in long-lived fish species with late sexual maturity and high energetic investment in reproduction like icefishes, pre-mating barriers are energetically more efficient than post-mating barriers to prevent hybridization.

Insertion Hot Spots of DIRS1 Retrotransposon and Chromosomal Diversifications among the Antarctic Teleosts Nototheniidae

By their faculty to transpose, transposable elements are known to play a key role in eukaryote genomes, impacting both their structuration and remodeling. Their integration in targeted sites may lead to recombination mechanisms involved in chromosomal rearrangements. The Antarctic fish family Nototheniidae went through several waves of species radiations. It is a suitable model to study transposable element (TE)-mediated mechanisms associated to genome and chromosomal diversifications. After the characterization of Gypsy (GyNoto), Copia (CoNoto), and DIRS1 (YNoto) retrotransposons in the genomes of Nototheniidae (diversity, distribution, conservation), we focused on their chromosome location with an emphasis on the three identified nototheniid radiations (the Trematomus, the plunderfishes, and the icefishes). The strong intrafamily TE conservation and wide distribution across species of the whole family suggest an ancestral acquisition with potential secondary losses in some lineages. GyNoto and CoNoto (including Hydra and GalEa clades) mostly produced interspersed signals along chromosomal arms. On the contrary, insertion hot spots accumulating in localized regions (mainly next to centromeric and pericentromeric regions) highlighted the potential role of YNoto in chromosomal diversifications as facilitator of the fusions which occurred in many nototheniid lineages, but not of the fissions.

Phylogenomics of an extra-Antarctic notothenioid radiation reveals a previously unrecognized lineage and diffuse species boundaries

Background

The impressive adaptive radiation of notothenioid fishes in Antarctic waters is generally thought to have been facilitated by an evolutionary key innovation, antifreeze glycoproteins, permitting the rapid evolution of more than 120 species subsequent to the Antarctic glaciation. By way of contrast, the second-most species-rich notothenioid genus, Patagonotothen, which is nested within the Antarctic clade of Notothenioidei, is almost exclusively found in the non-Antarctic waters of Patagonia. While the drivers of the diversification of Patagonotothen are currently unknown, they are unlikely to be related to antifreeze glycoproteins, given that water temperatures in Patagonia are well above freezing point. Here we performed a phylogenetic analysis based on genome-wide single nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing (RADseq) in a total of twelve Patagonotothen species.

Results

We present a well-supported, time-calibrated phylogenetic hypothesis including closely and distantly related outgroups, confirming the monophyly of the genus Patagonotothen with an origin approximately 3 million years ago and the paraphyly of both the sister genus Lepidonotothen and the family Notothenidae. Our phylogenomic and population genetic analyses highlight a previously unrecognized linage and provide evidence for shared genetic variation between some closely related species. We also provide a mitochondrial phylogeny showing mitonuclear discordance.

Conclusions

Based on a combination of phylogenomic and population genomic approaches, we provide evidence for the existence of a new, potentially cryptic, Patagonotothen species, and demonstrate that genetic boundaries between some closely related species are diffuse, likely due to recent introgression and/or incomplete linage sorting. The detected mitonuclear discordance highlights the limitations of relying on a single locus for species barcoding. In addition, our time-calibrated phylogenetic hypothesis shows that the early burst of diversification roughly coincides with the onset of the intensification of Quaternary glacial cycles and that the rate of species accumulation may have been stepwise rather than constant. Our phylogenetic framework not only advances our understanding of the origin of a high-latitude marine radiation, but also provides the basis for the study of the ecology and life history of the genus Patagonotothen, as well as for their conservation and commercial management.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1345-z) contains supplementary material, which is available to authorized users.

Information Criteria for Comparing Partition Schemes

When inferring phylogenies, one important decision is whether and how nucleotide substitution parameters should be shared across different subsets or partitions of the data. One sort of partitioning error occurs when heterogeneous subsets are mistakenly lumped together and treated as if they share parameter values. The opposite kind of error is mistakenly treating homogeneous subsets as if they result from distinct sets of parameters. Lumping and splitting errors are not equally bad. Lumping errors can yield parameter estimates that do not accurately reflect any of the subsets that were combined whereas splitting errors yield estimates that did not benefit from sharing information across partitions. Phylogenetic partitioning decisions are often made by applying information criteria such as the Akaike information criterion (AIC). As with other information criteria, the AIC evaluates a model or partition scheme by combining the maximum log-likelihood value with a penalty that depends on the number of parameters being estimated. For the purpose of selecting an optimal partitioning scheme, we derive an adjustment to the AIC that we refer to as the AIC$^{(p)}$ and that is motivated by the idea that splitting errors are less serious than lumping errors. We also introduce a similar adjustment to the Bayesian information criterion (BIC) that we refer to as the BIC$^{(p)}$. Via simulation and empirical data analysis, we contrast AIC and BIC behavior to our suggested adjustments. We discuss these results and also emphasize why we expect the probability of lumping errors with the AIC$^{(p)}$ and the BIC$^{(p)}$ to be relatively robust to model parameterization.

Phylogenetic analysis of Antarctic notothenioids illuminates the utility of RADseq for resolving Cenozoic adaptive radiations

Notothenioids are a clade of ∼120 species of marine fishes distributed in extreme southern hemisphere temperate near-shore habitats and in the Southern Ocean surrounding Antarctica. Over the past 25 years, molecular and morphological approaches have redefined hypotheses of relationships among notothenioid lineages as well as their relationships among major lineages of percomorph teleosts. These phylogenies provide a basis for investigation of mechanisms of evolutionary diversification within the clade and have enhanced our understanding of the notothenioid adaptive radiation. Despite extensive efforts, there remain several questions concerning the phylogeny of notothenioids. In this study, we deploy DNA sequences of ∼100,000 loci obtained using RADseq to investigate the phylogenetic relationships of notothenioids and to assess the utility of RADseq loci for lineages that exhibit divergence times ranging from the Paleogene to the Quaternary. The notothenioid phylogenies inferred from the RADseq loci provide unparalleled resolution and node support for several long-standing problems including, (1) relationships among species of Trematomus, (2) resolution of Indonotothenia cyanobrancha as the sister lineage of Trematomus, (3) the deep paraphyly of Nototheniidae, (4) the paraphyly of Lepidonotothen s.l., (5) paraphyly of Artedidraco, and 6) the monophyly of the Bathydraconidae. Assessment of site rates demonstrates that RADseq loci are similar to mtDNA protein coding genes and exhibit peak phylogenetic informativeness at the time interval during which the major Antarctic notothenioid lineages originated and diversified. In addition to providing a well-resolved phylogenetic hypothesis for notothenioids, our analyses quantify the predicted utility of RADseq loci for Cenozoic phylogenetic inferences.

Morphospace saturation in the stem-gnathostomes pteraspidiformes heterostracans: an early radiation of a 'bottom' heavy clade

Ostracoderms (fossil armoured jawless fishes) shed light on early vertebrate evolution by revealing the step-wise acquisition of jawed vertebrate characters, and were important constituents of Middle Palaeozoic vertebrate faunas. A wide variety of head shield shapes are observed within and between the ostracoderm groups, but the timing of these diversifications and the consistency between different measures of their morphospace are unclear. Here, we present the first disparity (explored morphospace) versus diversity (number of taxa) analysis of Pteraspidiformes heterostracans using continuous and discrete characters. Patterns of taxic diversity and morphological disparity are in accordance: they both show a rise to a peak in the Lochkovian followed by a gradual decline in the Middle-Late Devonian. Patterns are largely consistent for disparity measures using sum of ranges or total variance, and when using continuous or discrete characters. Pteraspidiformes heterostracans can be classified as a "bottom-heavy clade", i.e., a group where a high initial disparity decreasing over time is detected. In fact, the group explored morphospace early in its evolutionary history, with much of the subsequent variation in dermal armour occurring as variation in the proportions of already evolved anatomical features. This Early Devonian radiation is also in agreement with the paleobiogeographic distribution of the group, with a maximum of dispersal and explored morphospace during the Lochkovian and Pragian time bins.

Mobilization of retrotransposons as a cause of chromosomal diversification and rapid speciation: the case for the Antarctic teleost genus Trematomus

Background

The importance of transposable elements (TEs) in the genomic remodeling and chromosomal rearrangements that accompany lineage diversification in vertebrates remains the subject of debate. The major impediment to understanding the roles of TEs in genome evolution is the lack of comparative and integrative analyses on complete taxonomic groups. To help overcome this problem, we have focused on the Antarctic teleost genus Trematomus (Notothenioidei: Nototheniidae), as they experienced rapid speciation accompanied by dramatic chromosomal diversity. Here we apply a multi-strategy approach to determine the role of large-scale TE mobilization in chromosomal diversification within Trematomus species.

Results

Despite the extensive chromosomal rearrangements observed in Trematomus species, our measurements revealed strong interspecific genome size conservation. After identifying the DIRS1, Gypsy and Copia retrotransposon superfamilies in genomes of 13 nototheniid species, we evaluated their diversity, abundance (copy numbers) and chromosomal distribution. Four families of DIRS1, nine of Gypsy, and two of Copia were highly conserved in these genomes; DIRS1 being the most represented within Trematomus genomes. Fluorescence in situ hybridization mapping showed preferential accumulation of DIRS1 in centromeric and pericentromeric regions, both in Trematomus and other nototheniid species, but not in outgroups: species of the Sub-Antarctic notothenioid families Bovichtidae and Eleginopsidae, and the non-notothenioid family Percidae.

Conclusions

In contrast to the outgroups, High-Antarctic notothenioid species, including the genus Trematomus, were subjected to strong environmental stresses involving repeated bouts of warming above the freezing point of seawater and cooling to sub-zero temperatures on the Antarctic continental shelf during the past 40 millions of years (My). As a consequence of these repetitive environmental changes, including thermal shocks; a breakdown of epigenetic regulation that normally represses TE activity may have led to sequential waves of TE activation within their genomes. The predominance of DIRS1 in Trematomus species, their transposition mechanism, and their strategic location in “hot spots” of insertion on chromosomes are likely to have facilitated nonhomologous recombination, thereby increasing genomic rearrangements. The resulting centric and tandem fusions and fissions would favor the rapid lineage diversification, characteristic of the nototheniid adaptive radiation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4714-x) contains supplementary material, which is available to authorized users.

De Novo Gene Evolution of Antifreeze Glycoproteins in Codfishes Revealed by Whole Genome Sequence Data

New genes can arise through duplication of a pre-existing gene or de novo from non-coding DNA, providing raw material for evolution of new functions in response to a changing environment. A prime example is the independent evolution of antifreeze glycoprotein genes (afgps) in the Arctic codfishes and Antarctic notothenioids to prevent freezing. However, the highly repetitive nature of these genes complicates studies of their organization. In notothenioids, afgps evolved from an extant gene, yet the evolutionary origin of afgps in codfishes is unknown. Here, we demonstrate that afgps in codfishes have evolved de novo from non-coding DNA 13-18 Ma, coinciding with the cooling of the Northern Hemisphere. Using whole-genome sequence data from several codfishes and notothenioids, we find higher copy number of afgp in species exposed to more severe freezing suggesting a gene dosage effect. Notably, antifreeze function is lost in one lineage of codfishes analogous to the afgp losses in non-Antarctic notothenioids. This indicates that selection can eliminate the antifreeze function when freezing is no longer imminent. In addition, we show that evolution of afgp-assisting antifreeze potentiating protein genes (afpps) in notothenioids coincides with origin and lineage-specific losses of afgp. The origin of afgps in codfishes is one of the first examples of an essential gene born from non-coding DNA in a non-model species. Our study underlines the power of comparative genomics to uncover past molecular signatures of genome evolution, and further highlights the impact of de novo gene origin in response to a changing selection regime.

Lineage Diversity and Size Disparity in Musteloidea: Testing Patterns of Adaptive Radiation Using Molecular and Fossil-Based Methods

Adaptive radiation is hypothesized to be a primary mechanism that drives the remarkable species diversity and morphological disparity across the Tree of Life. Tests for adaptive radiation in extant taxa are traditionally estimated from calibrated molecular phylogenies with little input from extinct taxa. With 85 putative species in 33 genera and over 400 described extinct species, the carnivoran superfamily Musteloidea is a prime candidate to investigate patterns of adaptive radiation using both extant- and fossil-based macroevolutionary methods. The species diversity and equally impressive ecological and phenotypic diversity found across Musteloidea is often attributed to two adaptive radiations coinciding with two major climate events, the Eocene-Oligocene transition and the Mid-Miocene Climate Transition. Here, we compiled a novel time-scaled phylogeny for 88% of extant musteloids and used it as a framework for testing the predictions of adaptive radiation hypotheses with respect to rates of lineage diversification and phenotypic evolution. Contrary to expectations, we found no evidence for rapid bursts of lineage diversification at the origin of Musteloidea, and further analyses of lineage diversification rates using molecular and fossil-based methods did not find associations between rates of lineage diversification and the Eocene-Oligocene transition or Mid-Miocene Climate Transition as previously hypothesized. Rather, we found support for decoupled diversification dynamics driven by increased clade carrying capacity in the branches leading to a subclade of elongate mustelids. Supporting decoupled diversification dynamics between the subclade of elongate mustelids and the ancestral musteloid regime is our finding of increased rates of body length evolution, but not body mass evolution, within the decoupled mustelid subclade. The lack of correspondence in rates of body mass and length evolution suggest that phenotypic evolutionary rates under a single morphological metric, even one as influential as mass, may not capture the evolution of diversity in clades that exhibit elongate body shapes. The discordance in evolutionary rates between body length and body mass along with evidence of decoupled diversification dynamics suggests that body elongation might be an innovation for the exploitation of novel Mid-Miocene resources, resulting in the radiation of some musteloids.

Bayesian Phylogenetic Estimation of Clade Ages Supports Trans-Atlantic Dispersal of Cichlid Fishes

Divergence-time estimation based on molecular phylogenies and the fossil record has provided insights into fundamental questions of evolutionary biology. In Bayesian node dating, phylogenies are commonly time calibrated through the specification of calibration densities on nodes representing clades with known fossil occurrences. Unfortunately, the optimal shape of these calibration densities is usually unknown and they are therefore often chosen arbitrarily, which directly impacts the reliability of the resulting age estimates. As possible solutions to this problem, two nonexclusive alternative approaches have recently been developed, the “fossilized birth–death” (FBD) model and “total-evidence dating.” While these approaches have been shown to perform well under certain conditions, they require including all (or a random subset) of the fossils of each clade in the analysis, rather than just relying on the oldest fossils of clades. In addition, both approaches assume that fossil records of different clades in the phylogeny are all the product of the same underlying fossil sampling rate, even though this rate has been shown to differ strongly between higher level taxa. We here develop a flexible new approach to Bayesian age estimation that combines advantages of node dating and the FBD model. In our new approach, calibration densities are defined on the basis of first fossil occurrences and sampling rate estimates that can be specified separately for all clades. We verify our approach with a large number of simulated data sets, and compare its performance to that of the FBD model. We find that our approach produces reliable age estimates that are robust to model violation, on par with the FBD model. By applying our approach to a large data set including sequence data from over 1000 species of teleost fishes as well as 147 carefully selected fossil constraints, we recover a timeline of teleost diversification that is incompatible with previously assumed vicariant divergences of freshwater fishes. Our results instead provide strong evidence for transoceanic dispersal of cichlids and other groups of teleost fishes.

Disentangling Incomplete Lineage Sorting and Introgression to Refine Species-Tree Estimates for Lake Tanganyika Cichlid Fishes

Adaptive radiation is thought to be responsible for the evolution of a great portion of the past and present diversity of life. Instances of adaptive radiation, characterized by the rapid emergence of an array of species as a consequence to their adaptation to distinct ecological niches, are important study systems in evolutionary biology. However, because of the rapid lineage formation in these groups, and occasional gene flow between the participating species, it is often difficult to reconstruct the phylogenetic history of species that underwent an adaptive radiation. In this study, we present a novel approach for species-tree estimation in rapidly diversifying lineages, where introgression is known to occur, and apply it to a multimarker data set containing up to 16 specimens per species for a set of 45 species of East African cichlid fishes (522 individuals in total), with a main focus on the cichlid species flock of Lake Tanganyika. We first identified, using age distributions of most recent common ancestors in individual gene trees, those lineages in our data set that show strong signatures of past introgression. This led us to formulate three hypotheses of introgression between different lineages of Tanganyika cichlids: the ancestor of Boulengerochromini (or of Boulengerochromini and Bathybatini) received genomic material from the derived H-lineage; the common ancestor of Cyprichromini and Perissodini experienced, in turn, introgression from Boulengerochromini and/or Bathybatini; and the Lake Tanganyika Haplochromini and closely related riverine lineages received genetic material from Cyphotilapiini. We then applied the multispecies coalescent model to estimate the species tree of Lake Tanganyika cichlids, but excluded the lineages involved in these introgression events, as the multispecies coalescent model does not incorporate introgression. This resulted in a robust species tree, in which the Lamprologini were placed as sister lineage to the H-lineage (including the Eretmodini), and we identify a series of rapid splitting events at the base of the H-lineage. Divergence ages estimated with the multispecies coalescent model were substantially younger than age estimates based on concatenation, and agree with the geological history of the Great Lakes of East Africa. Finally, we formally tested the three hypotheses of introgression using a likelihood framework, and find strong support for introgression between some of the cichlid tribes of Lake Tanganyika.

Discrete and continuous character-based disparity analyses converge to the same macroevolutionary signal: a case study from captorhinids

The relationship between diversity and disparity during the evolutionary history of a clade provides unique insights into evolutionary radiations and the biological response to bottlenecks and to extinctions. Here we present the first comprehensive comparison of diversity and disparity of captorhinids, a group of basal amniotes that is important for understanding the early evolution of high-fiber herbivory. A new fully resolved phylogeny is presented, obtained by the inclusion of 31 morphometric characters. The new dataset is used to calculate diversity and disparity through the evolutionary history of the clade, using both discrete and continuous characters. Captorhinids do not show a decoupling between diversity and disparity, and are characterized by a rather symmetric disparity distribution, with a peak in occupied morphospace at about the midpoint of the clade's evolutionary history (Kungurian). This peak represents a delayed adaptive radiation, identified by the first appearance of several high-fiber herbivores in the clade, along with numerous omnivorous taxa. The discrete characters and continuous morphometric characters indicate the same disparity trends. Therefore, we argue that in the absence of one of these two possible proxies, the disparity obtained from just one source can be considered robust and representative of a general disparity pattern.

Phenotypic disparity in Iberian short-horned grasshoppers (Acrididae): the role of ecology and phylogeny

BACKGROUND

The combination of model-based comparative techniques, disparity analyses and ecomorphological correlations constitutes a powerful method to gain insight into the evolutionary mechanisms that shape morphological variation and speciation processes. In this study, we used a time-calibrated phylogeny of 70 Iberian species of short-horned grasshoppers (Acrididae) to test for patterns of morphological disparity in relation to their ecology and phylogenetic history. Specifically, we examined the role of substrate type and level of ecological specialization in driving different aspects of morphological evolution (locomotory traits, chemosensitive organs and cranial morphology) in this recent radiation.

RESULTS

We found a bimodal distribution of locomotory attributes corresponding to the two main substrate type guilds (plant vs. ground); plant-perching species tend to exhibit larger wings and thicker femora than those that remain on the ground. This suggests that life form (i.e., substrate type) is an important driving force in the evolution of morphological traits in short-horned grasshoppers, irrespective of ancestry. Substrate type and ecological specialization had no significant influence on head shape, a trait that showed a strong phylogenetic conservatism. Finally, we also found a marginal significant association between the length of antennae and the level of ecological specialization, suggesting that the development of sensory organs may be favored in specialist species.

CONCLUSIONS

Our results provide evidence that even in taxonomic groups showing limited morphological and ecological disparity, natural selection seems to play a more important role than genetic drift in driving the speciation process. Overall, this study suggests that morphostatic radiations should not necessarily be considered as "non-adaptive" and that the speciation process can bind both adaptive divergence mechanisms and neutral speciation processes related with allopatric and/or reproductive isolation.

Studying the evolutionary significance of thermal adaptation in ectotherms: The diversification of amphibians' energetics

A fundamental problem in evolutionary biology is the understanding of the factors that promote or constrain adaptive evolution, and assessing the role of natural selection in this process. Here, comparative phylogenetics, that is, using phylogenetic information and traits to infer evolutionary processes has been a major paradigm . In this study, we discuss Ornstein-Uhlenbeck models (OU) in the context of thermal adaptation in ectotherms. We specifically applied this approach to study amphibians's evolution and energy metabolism. It has been hypothesized that amphibians exploit adaptive zones characterized by low energy expenditure, which generate specific predictions in terms of the patterns of diversification in standard metabolic rate (SMR). We complied whole-animal metabolic rates for 122 species of amphibians, and adjusted several models of diversification. According to the adaptive zone hypothesis, we expected: (1) to find "accelerated evolution" in SMR (i.e., diversification above Brownian Motion expectations, BM), (2) that a model assuming evolutionary optima (i.e., an OU model) fits better than a white-noise model and (3) that a model assuming multiple optima (according to the three amphibians's orders) fits better than a model assuming a single optimum. As predicted, we found that the diversification of SMR occurred most of the time, above BM expectations. Also, we found that a model assuming an optimum explained the data in a better way than a white-noise model. However, we did not find evidence that an OU model with multiple optima fits the data better, suggesting a single optimum in SMR for Anura, Caudata and Gymnophiona. These results show how comparative phylogenetics could be applied for testing adaptive hypotheses regarding history and physiological performance in ectotherms.

Phylogenetic Analysis Supports the Aerobic-Capacity Model for the Evolution of Endothermy

The evolution of endothermy is a controversial topic in evolutionary biology, although several hypotheses have been proposed to explain it. To a great extent, the debate has centered on the aerobic-capacity model (AC model), an adaptive hypothesis involving maximum and resting rates of metabolism (MMR and RMR, respectively; hereafter "metabolic traits"). The AC model posits that MMR, a proxy of aerobic capacity and sustained activity, is the target of directional selection and that RMR is also influenced as a correlated response. Associated with this reasoning are the assumptions that (1) factorial aerobic scope (FAS; MMR/RMR) and net aerobic scope (NAS; MMR - RMR), two commonly used indexes of aerobic capacity, show different evolutionary optima and (2) the functional link between MMR and RMR is a basic design feature of vertebrates. To test these assumptions, we performed a comparative phylogenetic analysis in 176 vertebrate species, ranging from fish and amphibians to birds and mammals. Using disparity-through-time analysis, we also explored trait diversification and fitted different evolutionary models to study the evolution of metabolic traits. As predicted, we found (1) a positive phylogenetic correlation between RMR and MMR, (2) diversification of metabolic traits exceeding that of random-walk expectations, (3) that a model assuming selection fits the data better than alternative models, and (4) that a single evolutionary optimum best fits FAS data, whereas a model involving two optima (one for ectotherms and another for endotherms) is the best explanatory model for NAS. These results support the AC model and give novel information concerning the mode and tempo of physiological evolution of vertebrates.

Continental cichlid radiations: functional diversity reveals the role of changing ecological opportunity in the Neotropics

Adaptive radiations have been hypothesized to contribute broadly to the diversity of organisms. Models of adaptive radiation predict that ecological opportunity and ecological release, the availability of empty ecological niches and the response by adapting lineages to occupy them, respectively, drive patterns of phenotypic and lineage diversification. Adaptive radiations driven by 'ecological opportunity' are well established in island systems; it is less clear if ecological opportunity influences continent-wide diversification. We use Neotropical cichlid fishes to test if variation in rates of functional evolution is consistent with changing ecological opportunity. Across a functional morphological axis associated with ram-suction feeding traits, evolutionary rates declined through time as lineages diversified in South America. Evolutionary rates of ram-suction functional morphology also appear to have accelerated as cichlids colonized Central America and encountered renewed opportunity. Our results suggest that ecological opportunity may play an important role in shaping patterns of morphological diversity of even broadly distributed lineages like Neotropical cichlids.

Evolution of opercle bone shape along a macrohabitat gradient: species identification using mtDNA and geometric morphometric analyses in neotropical sea catfishes (Ariidae)

Transitions between the marine and freshwater macrohabitat have occurred repeatedly in the evolution of teleost fishes. For example, ariid catfishes have moved from freshwater to marine environments, and vice versa. Opercles, a skeletal feature that has been shown to change during such transitions, were subjected to 2D geometric morphometric analyses in order to investigate evolutionary shape changes during habitat transition in ariid catfishes and to test the influence of habitat on shape changes. A mtDNA marker, which proved useful in previous studies, was used to verify species identities. It greatly improved the assignment of specimens to a species, which are difficult to assign by morphology alone. The application of a mtDNA marker confirmed the occurrence of Notarius biffi in Central America, South of El Salvador. Molecular identification together with principal component analysis (PCA) and further morphological inspection of neurocrania indicated the existence of a cryptic species within Bagre pinnimaculatus. Principal component (PC) scores of individual specimens clustered in morphospace by genus rather than by habitat. Strong phylogenetic structure was detected using a permutation test of PC scores of species means on a phylogenetic tree. Calculation of Pagel's λ suggested that opercle shape evolved according to a Brownian model of evolution. Yet canonical variate analysis (CVA) conducted on the habitat groups showed significant differences in opercle shapes among freshwater and marine species. Overall, opercle shape in tropical American Ariidae appears to be phylogenetically constrained. This verifies the application of opercle shape as a taxonomic tool for species identification in fossil ariid catfishes. At the same time, adaptation to freshwater habitats shows characteristic opercle shape trajectories in ariid catfishes, which might be used to detect habitat preferences in fossils.

Evolution in an extreme environment: developmental biases and phenotypic integration in the adaptive radiation of antarctic notothenioids

Background

Over the past 40 million years water temperatures have dramatically dropped in the Southern Ocean, which has led to the local extinction of most nearshore fish lineages. The evolution of antifreeze glycoproteins in notothenioids, however, enabled these ancestrally benthic fishes to survive and adapt as temperatures reached the freezing point of seawater (−1.86 °C). Antarctic notothenioids now represent the primary teleost lineage in the Southern Ocean and are of fundamental importance to the local ecosystem. The radiation of notothenioids has been fostered by the evolution of “secondary pelagicism”, the invasion of pelagic habitats, as the group diversified to fill newly available foraging niches in the water column. While elaborate craniofacial modifications have accompanied this adaptive radiation, little is known about how these morphological changes have contributed to the evolutionary success of notothenioids.

Results

We used a 3D-morphometrics approach to investigate patterns of morphological variation in the craniofacial skeleton among notothenioids, and show that variation in head shape is best explained by divergent selection with respect to foraging niche. We document further an accelerated rate of morphological evolution in the icefish family Channichthyidae, and show that their rapid diversification was accompanied by the evolution of relatively high levels of morphological integration. Whereas most studies suggest that extensive integration should constrain phenotypic evolution, icefish stand out as a rare example of increased integration possibly facilitating evolutionary potential. Finally, we show that the unique feeding apparatus in notothenioids in general, and icefish in particular, can be traced to shifts in early developmental patterning mechanisms and ongoing growth of the pharyngeal skeleton.

Conclusion

Our work suggests that ecological opportunity is a major factor driving craniofacial variation in this group. Further, the observation that closely related lineages can differ dramatically in integration suggests that this trait can evolve quickly. We propose that the evolution of high levels of phenotypic integration in icefishes may be considered a key innovation that facilitated their morphological evolution and subsequent ecological expansion.

Electronic supplementary material

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

A phylogenetic analysis of macroevolutionary patterns in fermentative yeasts

When novel sources of ecological opportunity are available, physiological innovations can trigger adaptive radiations. This could be the case of yeasts (Saccharomycotina), in which an evolutionary novelty is represented by the capacity to exploit simple sugars from fruits (fermentation). During adaptive radiations, diversification and morphological evolution are predicted to slow‐down after early bursts of diversification. Here, we performed the first comparative phylogenetic analysis in yeasts, testing the “early burst” prediction on species diversification and also on traits of putative ecological relevance (cell‐size and fermentation versatility). We found that speciation rates are constant during the time‐range we considered (ca., 150 millions of years). Phylogenetic signal of both traits was significant (but lower for cell‐size), suggesting that lineages resemble each other in trait‐values. Disparity analysis suggested accelerated evolution (diversification in trait values above Brownian Motion expectations) in cell‐size. We also found a significant phylogenetic regression between cell‐size and fermentation versatility (R² = 0.10), which suggests correlated evolution between both traits. Overall, our results do not support the early burst prediction both in species and traits, but suggest a number of interesting evolutionary patterns, that warrant further exploration. For instance, we show that the Whole Genomic Duplication that affected a whole clade of yeasts, does not seems to have a statistically detectable phenotypic effect at our level of analysis. In this regard, further studies of fermentation under common‐garden conditions combined with comparative analyses are warranted.

Evolution of opercle shape in cichlid fishes from Lake Tanganyika - adaptive trait interactions in extant and extinct species flocks

Phenotype-environment correlations and the evolution of trait interactions in adaptive radiations have been widely studied to gain insight into the dynamics underpinning rapid species diversification. In this study we explore the phenotype-environment correlation and evolution of operculum shape in cichlid fishes using an outline-based geometric morphometric approach combined with stable isotope indicators of macrohabitat and trophic niche. We then apply our method to a sample of extinct saurichthyid fishes, a highly diverse and near globally distributed group of actinopterygians occurring throughout the Triassic, to assess the utility of extant data to inform our understanding of ecomorphological evolution in extinct species flocks. A series of comparative methods were used to analyze shape data for 54 extant species of cichlids (N = 416), and 6 extinct species of saurichthyids (N = 44). Results provide evidence for a relationship between operculum shape and feeding ecology, a concentration in shape evolution towards present along with evidence for convergence in form, and significant correlation between the major axes of shape change and measures of gut length and body elongation. The operculum is one of few features that can be compared in extant and extinct groups, enabling reconstruction of phenotype-environment interactions and modes of evolutionary diversification in deep time.

Selection towards different adaptive optima drove the early diversification of locomotor phenotypes in the radiation of Neotropical geophagine cichlids

Background

Simpson envisaged a conceptual model of adaptive radiation in which lineages diversify into “adaptive zones” within a macroevolutionary adaptive landscape. However, only a handful of studies have empirically investigated this adaptive landscape and its consequences for our interpretation of the underlying mechanisms of phenotypic evolution. In fish radiations the evolution of locomotor phenotypes may represent an important dimension of ecomorphological diversification given the implications of locomotion for feeding and habitat use. Neotropical geophagine cichlids represent a newly identified adaptive radiation and provide a useful system for studying patterns of locomotor diversification and the implications of selective constraints on phenotypic divergence in general.

Results

We use multivariate ordination, models of phenotypic evolution and posterior predictive approaches to investigate the macroevolutionary adaptive landscape and test for evidence of early divergence of locomotor phenotypes in Geophagini. The evolution of locomotor phenotypes was characterized by selection towards at least two distinct adaptive peaks and the early divergence of modern morphological disparity. One adaptive peak included the benthic and epibenthic invertivores and was characterized by fishes with deep, laterally compressed bodies that optimize precise, slow-swimming manoeuvres. The second adaptive peak resulted from a shift in adaptive optima in the species-rich ram-feeding/rheophilic Crenicichla-Teleocichla clade and was characterized by species with streamlined bodies that optimize fast starts and rapid manoeuvres. Evolutionary models and posterior predictive approaches favoured an early shift to a new adaptive peak over decreasing rates of evolution as the underlying process driving the early divergence of locomotor phenotypes.

Conclusions

The influence of multiple adaptive peaks on the divergence of locomotor phenotypes in Geophagini is compatible with the expectations of an ecologically driven adaptive radiation. This study confirms that the diversification of locomotor phenotypes represents an important dimension of phenotypic evolution in the geophagine adaptive radiation. It also suggests that the commonly observed early burst of phenotypic evolution during adaptive radiations may be better explained by the concentration of shifts to new adaptive peaks deep in the phylogeny rather than overall decreasing rates of evolution.

Electronic supplementary material

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