Phylogenomics of pike cichlids (Cichlidae: Crenicichla): the rapid ecological speciation of an incipient species flock

Phylogenomics of trans-Andean tetras of the genus Hyphessobrycon Durbin 1908 (Stethaprioninae: Characidae) and colonization patterns of Middle America

Hyphessobrycon is one of the most species rich and widely distributed genera in the family Characidae, with more than 160 species ranging from Veracruz, Mexico to Mar Chiquita Lagoon in Buenos Aires, Argentina. The majority of Hyphessobrycon diversity shows a cis-Andean distribution; only nine species are trans-Andean including H. compressus (Meek 1908). It is well established that Hyphessobrycon is not monophyletic but it has been suggested that natural groups can be identified within the larger Hyphessobrycon species group. In this study, we tested the monophyly of trans-Andean species of Hyphessobrycon and investigated the placement of H. compressus. We inferred the first phylogenomic hypothesis of trans-Andean Hyphessobrycon that includes nearly complete taxonomic sampling (eight of nine valid species) using ultraconserved elements (UCEs). We analyzed 75% (1682 UCEs), 90% (1258 UCEs), and 95% (838 UCEs) complete data matrices, and inferred phylogenomic hypotheses under concatenation and coalescent approaches. In all cases, we recovered the monophyly of trans-Andean Hyphessobrycon inclusive of H. compressus, strong support for three species groups, and evidence of cryptic diversity within the widespread H. compressus and H. condotensis. We used our phylogenomic hypothesis to investigate the biogeographic history of Hyphessobrycon in Middle America. Our ancestral range estimation analysis suggests a single event of cis- to trans-Andean colonization followed by stepwise colonization from the Pacific slope of northwestern South America (Chocó block) to northern Middle America (Maya block). Our work supports the recognition of the trans-Andean species as Hyphessobrycon sensu stricto and provides an evolutionary template to examine morphological characters that will allow us to better understand the diversity of Hyphessobrycon in Middle America.

Genome-wide species delimitation analyses of a silverside fish species complex in central Mexico indicate taxonomic over-splitting

BACKGROUND

Delimiting species across a speciation continuum is a complex task, as the process of species origin is not generally instantaneous. The use of genome-wide data provides unprecedented resolution to address convoluted species delimitation cases, often unraveling cryptic diversity. However, because genome-wide approaches based on the multispecies coalescent model are known to confound population structure with species boundaries, often resulting in taxonomic over-splitting, it has become increasingly evident that species delimitation research must consider multiple lines of evidence. In this study, we used phylogenomic, population genomic, and coalescent-based species delimitation approaches, and examined those in light of morphological and ecological information, to investigate species numbers and boundaries comprising the Chirostoma "humboltianum group" (family Atherinidae). The humboltianum group is a taxonomically controversial species complex where previous morphological and mitochondrial studies produced conflicting species delimitation outcomes. We generated ddRADseq data for 77 individuals representing the nine nominal species in the group, spanning their distribution range in the central Mexican plateau.

RESULTS

Our results conflict with the morphospecies and ecological delimitation hypotheses, identifying four independently evolving lineages organized in three geographically cohesive clades: (i) chapalae and sphyraena groups in Lake Chapala, (ii) estor group in Lakes Pátzcuaro and Zirahuén, and (iii) humboltianum sensu stricto group in Lake Zacapu and Lerma river system.

CONCLUSIONS

Overall, our study provides an atypical example where genome-wide analyses delineate fewer species than previously recognized on the basis of morphology. It also highlights the influence of the geological history of the Chapala-Lerma hydrological system in driving allopatric speciation in the humboltianum group.

Functional Trade-offs Asymmetrically Promote Phenotypic Evolution

Trade-offs are thought to bias evolution and are core features of many anatomical systems. Therefore, trade-offs may have far-reaching macroevolutionary consequences, including patterns of morphological, functional, and ecological diversity. Jaws, like many complex anatomical systems, are comprised of elements involved in biomechanical trade-offs. We test the impact of a core mechanical trade-off, transmission of velocity versus force (i.e., mechanical advantage), on rates of jaw evolution in Neotropical cichlids. Across 130 species representing a wide array of feeding ecologies, we find that the velocity-force trade-off impacts evolution of the surrounding jaw system. Specifically, rates of jaw evolution are faster at functional extremes than in more functionally intermediate or unspecialized jaws. Yet, surprisingly, the effect on jaw evolution is uneven across the extremes of the velocity-force continuum. Rates of jaw evolution are 4 to 10-fold faster in velocity-modified jaws, whereas force-modified jaws are 7 to 18-fold faster, compared to unspecialized jaws, depending on the extent of specialization. Further, we find that a more extreme mechanical trade-off resulted in faster rates of jaw evolution. The velocity-force trade-off reflects a gradient from specialization on capture-intensive (e.g., evasive or buried) to processing-intensive prey (e.g., attached or shelled), respectively. The velocity extreme of the trade-off is characterized by large magnitudes of trait change leading to functionally divergent specialists and ecological stasis. By contrast, the force extreme of the trade-off is characterized by enhanced ecological lability made possible by phenotypes more readily co-opted for different feeding ecologies. This asymmetry of macroevolutionary outcomes along each extreme is likely the result of an enhanced utility of the pharyngeal jaw system as force-modified oral jaws are adapted for prey that require intensive processing (e.g., algae, detritus, and molluscs). The velocity-force trade-off, a fundamental feature of many anatomical systems, promotes rapid phenotypic evolution of the surrounding jaw system in a canonical continental adaptive radiation. Considering that the velocity-force trade-off is an inherent feature of all jaw systems that involve a lower element that rotates at a joint, spanning the vast majority of vertebrates, our results may be widely applicable across the tree of life. [adaptive radiation; constraint; decoupling; jaws; macroevolution; specialization].

Riverscape genomics of cichlid fishes in the lower Congo: Uncovering mechanisms of diversification in an extreme hydrological regime

Freshwater fishes are notably diverse, given that freshwater habitat represents a tiny fraction of the earth's surface, but the mechanisms generating this diversity remain poorly understood. Rivers provide excellent models to understand how freshwater diversity is generated and maintained across heterogeneous habitats. In particular, the lower Congo River (LCR) consists of a dynamic hydroscape exhibiting extraordinary aquatic biodiversity, endemicity, morphological and ecological specialization. Previous studies have suggested that the numerous high-energy rapids throughout the LCR form physical barriers to gene flow, thus facilitating diversification and speciation, generating ichthyofaunal diversity. However, this hypothesis has not been fully explored using genome-wide SNPs for fish species distributed across the LCR. Here, we examined four lamprologine cichlids endemic to the LCR that are distributed along the river without range overlap. Using genome-wide SNP data, we tested the hypotheses that high-energy rapids serve as physical barriers to gene flow that generate genetic divergence at inter- and intraspecific levels, and that gene flow occurs primarily in a downstream direction. Our results are consistent with the prediction that powerful rapids sometimes act as a barrier to gene flow but also suggest that, at certain temporal and spatial scales, they may provide multidirectional dispersal opportunities for riverine rheophilic cichlid fishes. These results highlight the complexity of diversification processes in rivers and the importance of assessing such processes across different riverscapes.

Rapid Parallel Morphological and Mechanical Diversification of South American Pike Cichlids (Crenicichla)

Explosive bouts of diversification are one of the most conspicuous features of the tree of life. When such bursts are repeated in similar environments it suggests some degree of predictability in the evolutionary process. We assess parallel adaptive radiation of South American pike cichlids (Crenicichla) using phylogenomics and phylogenetic comparative methods. We find that species flocks in the Uruguay and Iguazú River basins rapidly diversified into the same set of ecomorphs that reflect feeding ecology. Both adaptive radiations involve expansion of functional morphology, resulting in unique jaw phenotypes. Yet, form and function were decoupled such that most ecomorphs share similar mechanical properties of the jaws (i.e., jaw motion during a feeding strike). Prey mobility explained six to nine-fold differences in the rate of morphological evolution, but had no effect on the rate of mechanical evolution. We find no evidence of gene flow between species flocks or with surrounding coastal lineages that may explain their rapid diversification. When compared to cichlids of the East African Great Lakes and other prominent adaptive radiations, pike cichlids share many themes, including rapid expansion of phenotypic diversity, specialization along the benthic-to-pelagic habitat and soft-to-hard prey axes, and the evolution of conspicuous functional innovations. Yet, decoupled evolution of form and function and the absence of hybridization as a catalyzing force are departures from patterns observed in other adaptive radiations. Many-to-one mapping of morphology to mechanical properties is a mechanism by which pike cichlids exhibit a diversity of feeding ecologies while avoiding exacerbating underlying mechanical trade-offs.

MtDNA species-level phylogeny and delimitation support significantly underestimated diversity and endemism in the largest Neotropical cichlid genus (Cichlidae: Crenicichla)

Crenicichla is the largest and most widely distributed genus of Neotropical cichlids. Here, we analyze a mtDNA dataset comprising 681 specimens (including Teleocichla, a putative ingroup of Crenicichla) and 77 out of 105 presently recognized valid species (plus 10 out of 36 nominal synonyms plus over 50 putatively new species) from 129 locations in 31 major river drainages throughout the whole distribution of the genus in South America. Based on these data we make an inventory of diversity and highlight taxa and biogeographic areas worthy of further sampling effort and conservation protection. Using three methods of molecular species delimitation, we find between 126 and 168 species-like clusters, i.e., an average increase of species diversity of 65–121% with a range of increase between species groups. The increase ranges from 0% in the Missioneira and Macrophthama groups, through 25–40% (Lacustris group), 50–87% (Reticulata group, Teleocichla), 68–168% (Saxatilis group), 125–200% (Wallacii group), and 158–241% in the Lugubris group. We found a high degree of congruence between clusters derived from the three used methods of species delimitation. Overall, our results recognize substantially underestimated diversity in Crenicichla including Teleocichla. Most of the newly delimited putative species are from the Amazon-Orinoco-Guiana (AOG) core area (Greater Amazonia) of the Neotropical region, especially from the Brazilian and Guiana shield areas of which the former is under the largest threat and largest degree of environmental degradation of all the Amazon.

Sexual dimorphism in an adaptive radiation: Does intersexual niche differentiation result in ecological character displacement?

Evolutionary radiations are one plausible explanation for the rich biodiversity on Earth. Adaptive radiations are the most studied form of evolutionary radiations, and ecological opportunity has been identified as one factor permitting them. Competition among individuals is supposedly highest in populations of conspecifics. Divergent modes of resource use might minimize trophic overlap, and thus intersexual competition, resulting in ecological character displacement between sexes. However, the role of intersexual differentiation in speciation processes is insufficiently studied. The few studies available suggest that intersexual niche differentiation exists in adaptive radiations, but their role within the radiation, and the extent of differentiation within the organism itself, remains largely unexplored. Here, we test the hypothesis that multiple morphological structures are affected by intersexual niche differentiation in "roundfin" Telmatherina, the first case where intersexual niche differentiation was demonstrated in an adaptive fish radiation. We show that sexes of two of the three morphospecies differ in several structural components of the head, all of these are likely adaptive. Sexual dimorphism is linked to the respective morphospecies-specific ecology and affects several axes of variation. Trait variation translates into different feeding modes, processing types, and habitat usages that add to interspecific variation in all three morphospecies. Intrasexual selection, that is, male-male competition, may contribute to variation in some of the traits, but appears unlikely in internal structures, which are invisible to other individuals. We conclude that intersexual variation adds to the adaptive diversity of roundfins and might play a key role in minimizing intersexual competition in emerging radiations.

Comparing Ultraconserved Elements and Exons for Phylogenomic Analyses of Middle American Cichlids: When Data Agree to Disagree

Choosing among types of genomic markers to be used in a phylogenomic study can have a major influence on the cost, design, and results of a study. Yet few attempts have been made to compare categories of next-generation sequence markers limiting our ability to compare the suitability of these different genomic fragment types. Here, we explore properties of different genomic markers to find if they vary in the accuracy of component phylogenetic trees and to clarify the causes of conflict obtained from different data sets or inference methods. As a test case, we explore the causes of discordance between phylogenetic hypotheses obtained using a novel data set of ultraconserved elements (UCEs) and a recently published exon data set of the cichlid tribe Heroini. Resolving relationships among heroine cichlids has historically been difficult, and the processes of colonization and diversification in Middle America and the Greater Antilles are not yet well understood. Despite differences in informativeness and levels of gene tree discordance between UCEs and exons, the resulting phylogenomic hypotheses generally agree on most relationships. The independent data sets disagreed in areas with low phylogenetic signal that were overwhelmed by incomplete lineage sorting and nonphylogenetic signals. For UCEs, high levels of incomplete lineage sorting were found to be the major cause of gene tree discordance, whereas, for exons, nonphylogenetic signal is most likely caused by a reduced number of highly informative loci. This paucity of informative loci in exons might be due to heterogeneous substitution rates that are problematic to model (i.e., computationally restrictive) resulting in systematic errors that UCEs (being less informative individually but more uniform) are less prone to. These results generally demonstrate the robustness of phylogenomic methods to accommodate genomic markers with different biological and phylogenetic properties. However, we identify common and unique pitfalls of different categories of genomic fragments when inferring enigmatic phylogenetic relationships.

Comparing diversification rates in lakes, rivers, and the sea

The diversity of species inhabiting freshwater relative to marine habitats is striking, given that freshwater habitats encompass <1% of Earth's water. The most commonly proposed explanation for this pattern is that freshwater habitats are more fragmented than marine habitats, allowing more opportunities for allopatric speciation and thus increased diversification rates in freshwater. However, speciation may be generally faster in sympatry than in allopatry, as illustrated by lacustrine radiations such as African cichlids. Such differences between rivers and lakes may be important to consider when comparing diversification broadly among freshwater and marine groups. Here I compared diversification rates of teleost fishes in marine, riverine and lacustrine habitats. I found that lakes had faster speciation and net diversification rates than other aquatic habitats. However, most freshwater diversity arose in rivers. Surprisingly, riverine and marine habitats had similar rates of net diversification on average. Biogeographic models suggest that lacustrine habitats are evolutionarily unstable, explaining the dearth of lacustrine species in spite of their rapid diversification. Collectively, these results suggest that strong diversification rate differences are unlikely to explain the freshwater paradox. Instead, this pattern may be attributable to the comparable amount of time spent in riverine and marine habitats over the 200-million-year history of teleosts.

Distribution of freshwater fish from the Southern Neotropics reveals three new areas of endemism and show diffuse limits among freshwater ecoregions

Abstract Ecoregions and areas of endemism are central concepts in biogeography. Based on collection records and the Endemic Analyses implemented with NDM/VNDM method we analyzed fish areas of endemism in the junction of three freshwater ecoregions related to the Rio de la Plata estuary (Lower Parana, Lower Uruguay, Laguna dos Patos) in Southern Neotropical region. Using two grid cell sizes, results obtained showed the same general patterns. Areas of endemism recovered were mainly associated either to Lower Uruguay or to Laguna dos Patos. In both ecoregions nested areas of endemism were identified within larger patterns of endemism. Noteworthy, one area recovered occurred across Lower Uruguay and Laguna dos Patos limits. Our results also suggest a revision of the Lower Uruguay and Lower Parana ecoregion limits, and highlight the relevance of the Rio de la Plata estuary as a barrier and corridor for freshwater fishes in the area.

Diversification of Neotropical Freshwater Fishes

Neotropical freshwater fishes (NFFs) constitute the most diverse continental vertebrate fauna on Earth, with more than 6,200 named species compressed into an aquatic footprint <0.5% of the total regional land-surface area and representing the greatest phenotypic disparity and functional diversity of any continental ichthyofauna. Data from the fossil record and time-calibrated molecular phylogenies indicate that most higher taxa (e.g., genera, families) diversified relatively continuously through the Cenozoic, across broad geographic ranges of the South American platform. Biodiversity data for most NFF clades support a model of continental radiation rather than adaptive radiation, in which speciation occurs mainly in allopatry, and speciation and adaptation are largely decoupled. These radiations occurred under the perennial influence of river capture and sea-level oscillations, which episodically fragmented and merged portions of adjacent river networks. The future of the NFF fauna into the Anthropocene is uncertain, facing numerous threats at local, regional, and continental scales.

A Cautionary Note on the Use of Genotype Callers in Phylogenomics

Next-generation-sequencing genotype callers are commonly used in studies to call variants from newly sequenced species. However, due to the current availability of genomic resources, it is still common practice to use only one reference genome for a given genus, or even one reference for an entire clade of a higher taxon. The problem with traditional genotype callers, such as the one from GATK, is that they are optimized for variant calling at the population level. However, when these callers are used at the phylogenetic level, the consequences for downstream analyses can be substantial. Here, we performed simulations to compare the performance between the genotype callers of GATK and ATLAS, and present their differences at various phylogenetic scales. We show that the genotype caller of GATK substantially underestimates the number of variants at the phylogenetic level, but not at the population level. We also found that the accuracy of heterozygote calls declines with increasing distance to the reference genome. We quantified this decline and found that it is very sharp in GATK, while ATLAS maintains high accuracy even at moderately divergent species from the reference. We further suggest that efforts should be taken towards acquiring more reference genomes per species, before pursuing high-scale phylogenomic studies. [ATLAS; efficiency of SNP calling; GATK; heterozygote calling; next-generation sequencing; reference genome; variant calling.]

The role of head shape and trophic variation in the diversification of the genus Herichthys in sympatry and allopatry

In the present study we evaluated the putative cases of sympatric speciation in the genus Herichthys by studying the variation in head shape using principal component analysis, phylomorphospace and reconstructions of the ancestral states of feeding preferences. Herichthys includes both allopatric and sympatric sister species, as well as sympatric unrelated species and thus offers great potential for evolutionary studies of putatively sympatric speciation. Herichthys is the northernmost group of cichlids in America and one of the most ecologically disparate genera within Middle American cichlids. Fifteen anatomical points were recorded on the heads of 293 specimens of the 11 species recognized within the genus. The results show that in spite of having wide variation in consumed diets, most species of Herichthys are close in morphospace. However, morphological variation was great among the two pairs of sympatric sister species in agreement with the suggested sympatric model of speciation.

Unraveling the systematics and evolution of the 'Geophagus' brasiliensis (Cichliformes: Cichlidae) species complex

The 'Geophagus' brasiliensis complex is one of the most abundant groups of cichlids from eastern coastal basins in South America. Traditionally, this fish group has been recognized as incertae sedis because of phylogenetic uncertainties and unclear taxonomy. In addition, the remarkable morphological, chromosomal, and DNA variation reported over recent years in several populations of these cichlids has increased the debate about their species richness and their distributional range. Here, we tested the presence of independent evolutionary lineages within the 'G.' brasiliensis complex, addressing their taxonomic status and evolutionary relationships, including a comparative analysis of genetic and morphological patterns, based on an extensive dataset, comprising 172 sampling sites along most of their known range using a mitochondrial marker, RADseq data and geometric morphometrics. The number of putative species in the present study varied from 9 to 11 depending on the molecular species delimitation methods used. Our results revealed at least two putative new taxa ('Geophagus' sp. Doce and 'Geophagus' sp. Upper Contas). Morphometric analyses, particularly those based on Canonical Variate Analysis (CVA), revealed significant morphological differentiation between species within the main clades. On the other hand, analyses of morphological phylogenetic signal and phylomorphospace provided no evidence of adaptive differentiation among these species. Thus, diversification in the 'G.' brasiliensis complex seems to have been influenced by hydrogeological events that promoted allopatry, such as the presence of paleodrainages and distributional reconfiguration through river captures. We propose major changes in the known distribution of some species within the complex and conservatively suggest the recognition of 10 species within the 'Geophagus' brasiliensis complex, with the potential for further dividing 'G.' rufomarginatus after additional taxonomic evaluation.

Molecular clocks, biogeography and species diversity in Herichthys with evaluation of the role of Punta del Morro as a vicariant brake along the Mexican Transition Zone in the context of local and global time frame of cichlid diversification

Using molecular dated phylogenies and biogeographic reconstructions, the species diversity, biogeography and time frame of evolution of the genus Herichthys were evaluated. In particular, we test the role of Punta del Morro (PdM) as a vicariant brake along the Mexican Transition Zone in the context of local and global time frame of cichlid diversification using several sets of calibrations. Species diversity in Herichthys is complex and the here employed dating methods suggest young age and rapid divergence for many species while species delimitation methods did not resolve these young species including both sympatric species pairs. Based on our molecular clock dating analyses, Herichthys has colonized its present distribution area significantly prior to the suggested vicariance by PdM (10-17.1 Ma vs. 5 to 7.5 Ma). The PdM constraint is in conflict with all other paleogeographic and fossil constraints including novel ones introduced in this study that are, however, congruent among each other. Our study demonstrates that any cichlid datings significantly older or younger than the bounds presented by our analyses and discussion have to be taken as highly questionable from the point of view of Middle American paleogeography and cichlid biogeography unless we allow the option that cichlid biogeography is completely independent from ecological and geological constraints.

Shared ecological traits influence shape of the skeleton in flatfishes (Pleuronectiformes)

In the age of phylogenetic comparative methods, evolutionary biologists have been able to explore evolutionary trends in form in unique and extraordinarily diverse groups of animals. Pleuronectiformes, commonly known as flatfishes, is a diverse and specialized order of fishes that have remarkable asymmetry induced by ocular migration and a benthic life style. Although flatfishes are unique from other fishes, species within the group are morphologically diverse. The origin of ocular migration has been a primary focus of research; however, little is known about overall shape diversification among the flatfishes. In this study, we use integrative methods to examine how body shape evolved within the flatfishes. Shape was quantified from X-rays using geometric morphometrics for 389 individuals across 145 species. The most recent and robust phylogeny was overlaid onto the morphospace and phylogenetic signal was calculated to ascertain convergence in the morphospace. In addition, phylogenetic linear models were employed to determine if ecological traits were correlated with shape and if size had an effect on overall body shape. Results revealed that the majority of variation evolved recently, within the past 15–10-million-years in the middle Miocene, and is highly variable within the flatfishes. These changes are best summarized by body depth, jaw length and medial fin length. Dorsal and anal fin length are correlated, which may be due to the unique mode of locomotion used by flatfishes. A phylogenetic linear model and phylomorphospace analysis suggested that several ecological traits are correlated with shape, which indicates an ecological role in the diversification of flatfishes.

Fine-scale species delimitation: speciation in process and periodic patterns in nudibranch diversity

Using the nudibranch genus Amphorina as a model, ongoing speciation is demonstrated, as well as how periodic-like patterns in colouration can be included in an integrated method of fine-scale species delimitation. By combining several methods, including BPP analysis and the study of molecular, morphological, and ecological data from a large number of specimens within a broad geographic range from northern Europe to the Mediterranean, five species are recognised within the genus Amphorina, reviewed here for the first time. Two new species from the southwestern coast of Sweden are described, A. viriola sp. nov. and A. andra sp. nov. Evidence is provided of a recent speciation process between the two closely related, yet separate, species which inhabit the same geographic localities but demonstrate strict water depth differentiation, with one species inhabiting the shallow brackish top layer above the halocline and the other species inhabiting the underlying saltier water. The results presented here are of relevance for currently debated issues such as conservation in relation to speciation, fine species delimitation, and integration of molecular, morphological and ecological information in biodiversity studies. The periodic approach to biological taxonomy has considerable practical potential for various organismal groups.

Macroevolutionary analyses indicate that repeated adaptive shifts towards predatory diets affect functional diversity in Neotropical cichlids

During adaptive radiation, diversification within clades is limited by adaptation to the available ecological niches, and this may drive patterns of both trait and species diversity. However, adaptation to disparate niches may result in varied impacts on the timing, pattern and rate of morphological evolution. In this study, we examined the relationship between feeding ecology and functional diversification across a diverse clade of freshwater fishes, the Neotropical cichlids. Species dietary niches were ordinated via multivariate analysis of stomach content data. We investigated changes in the rate and pattern of morphological diversification associated with feeding, including dietary niche and degree of dietary specialization. A major division in dietary niche space was observed between predators that consume fish and macroinvertebrates vs. other groups with diets dominated by small invertebrates, detritus or vegetation. These trophic niches were strongly associated with groupings defined by functional morphospace. Clades within the piscivore/macroinvertivore group rarely transitioned to other dietary niches. Comparatively, high dietary specialization enhanced functional diversification, driving the evolution of more extreme morphologies. Divergent patterns of trophic diversification among Neotropical cichlids appear to derive from different performance demands in regional abiotic and biotic environments associated with biogeographical history.

Head Shape Modulates Diversification of a Classic Cichlid Pharyngeal Jaw Innovation

Functional innovations are often invoked to explain the uneven distribution of ecological diversity. Innovations may provide access to new adaptive zones by expanding available ecological opportunities and may serve as catalysts of adaptive radiation. However, diversity is often unevenly distributed within clades that share a key innovation, highlighting the possibility that the impact of the innovation is mediated by other traits. Pharyngognathy is a widely recognized innovation of the pharyngeal jaws that enhances the ability to process hard and tough prey in several major radiations of fishes, including marine wrasses and freshwater cichlids. We explored diversification of lower pharyngeal jaw shape, a key feature of pharyngognathy, and the extent to which it is influenced by head shape in Neotropical cichlids. While pharyngeal jaw shape was unaffected by either head length or head depth, its disparity declined dramatically with increasing head width. Head width also predicted the rate of pharyngeal jaw evolution such that higher rates were associated with narrow heads. Wide heads are associated with exploiting prey that require intense processing by pharyngeal jaws that have expanded surfaces for the attachment of enlarged muscles. However, we show that a wide head constrains access to adaptive peaks associated with several trophic roles. A constraint on the independent evolution of pharyngeal jaw and head shape may explain the uneven distribution of ecological diversity within a clade that shares a major functional innovation.

Comparative Functional and Phylogenomic Analyses of Host Association in the Remoras (Echeneidae), a Family of Hitchhiking Fishes

The family Echeneidae consists of eight species of marine fishes that hitchhike by adhering to a wide variety of vertebrate hosts via a sucking disc. While several studies have focused on the interrelationships of the echeneids and the adhesion performance of a single species, no clear phylogenetic hypothesis has emerged and the morphological basis of adhesion remains largely unknown. We first set out to resolve the interrelationships of the Echeneidae by taking a phylogenomic approach using ultraconserved elements. Then, within this framework, we characterized disc morphology through µ-CT analysis, evaluated host specificity through an analysis of host phylogenetic distance, and determined which axes of disc morphological variation are associated with host diversity, skin surface properties, mean pairwise phylogenetic distance (MPD obs.), and swimming regime. We recovered an extremely well-supported topology, found that the specificity of host choice is more variable in a pelagic group and less variable in a reef-generalist group than previously proposed, and that axes of disc morphospace are best explained by models that include host skin surface roughness, host MPD obs., and maximum host Reynolds number. This suggests that ecomorphological diversification was driven by the selection pressures of host skin surface roughness, host specialization, and hydrodynamic regime.

Genomic, ecological, and morphological approaches to investigating species limits: A case study in modern taxonomy from Tropical Eastern Pacific surgeonfishes

A wide variety of species are distinguished by slight color variations. However, molecular analyses have repeatedly demonstrated that coloration does not always correspond to distinct evolutionary histories between closely related groups, suggesting that this trait is labile and can be misleading for species identification. In the present study, we analyze the evolutionary history of sister species of Prionurus surgeonfishes in the Tropical Eastern Pacific (TEP), which are distinguished by the presence or absence of dark spots on their body. We examined the species limits in this system using comparative specimen-based approaches, a mitochondrial gene (COI), more than 800 nuclear loci (Ultraconserved Elements), and abiotic niche comparisons. The results indicate there is a complete overlap of meristic counts and morphometric measurements between the two species. Further, we detected multiple individuals with intermediate spotting patterns suggesting that coloration is not diagnostic. Mitochondrial data recovered a single main haplotype shared between the species and all locations resulting in a complete lack of structure (ΦST = 0). Genomic analyses also suggest low levels of genetic differentiation (F ST = 0.013), and no alternatively fixed SNPs were detected between the two phenotypes. Furthermore, niche comparisons could not reject niche equivalency or similarity between the species. These results suggest that these two phenotypes are conspecific and widely distributed in the TEP. Here, we recognize Prionurus punctatus Gill 1862 as a junior subjective synonym of P. laticlavius (Valenciennes 1846). The underlying causes of phenotypic variation in this species are unknown. However, this system gives insight into general evolutionary dynamics within the TEP.

Island- and lake-like parallel adaptive radiations replicated in rivers

Parallel adaptive radiations have arisen following the colonization of islands by lizards and lakes by fishes. In these classic examples, parallel adaptive radiation is a response to the ecological opportunities afforded by the colonization of novel ecosystems and similar adaptive landscapes that favour the evolution of similar suites of ecomorphs, despite independent evolutionary histories. Here, we demonstrate that parallel adaptive radiations of cichlid fishes arose in South American rivers. Speciation-assembled communities of pike cichlids (Crenicichla) have independently diversified into similar suites of novel ecomorphs in the Uruguay and Paraná Rivers, including crevice feeders, periphyton grazers and molluscivores. There were bursts in phenotypic evolution associated with the colonization of each river and the subsequent expansion of morphospace following the evolution of the ecomorphs. These riverine clades demonstrate that characteristics emblematic of textbook parallel adaptive radiations of island- and lake-dwelling assemblages are feasible evolutionary outcomes even in labile ecosystems such as rivers.

Exploring data processing strategies in NGS target enrichment to disentangle radiations in the tribe Cardueae (Compositae)

Target enrichment is a cost-effective sequencing technique that holds promise for elucidating evolutionary relationships in fast-evolving lineages. However, potential biases and impact of bioinformatic sequence treatments in phylogenetic inference have not been thoroughly explored yet. Here, we investigate this issue with an ultimate goal to shed light into a highly diversified group of Compositae (Asteraceae) constituted by four main genera: Arctium, Cousinia, Saussurea, and Jurinea. Specifically, we compared sequence data extraction methods implemented in two easy-to-use workflows, PHYLUCE and HybPiper, and assessed the impact of two filtering practices intended to reduce phylogenetic noise. In addition, we compared two phylogenetic inference methods: (1) the concatenation approach, in which all loci were concatenated in a supermatrix; and (2) the coalescence approach, in which gene trees were produced independently and then used to construct a species tree under coalescence assumptions. Here we confirm the usefulness of the set of 1061 COS targets (a nuclear conserved orthology loci set developed for the Compositae) across a variety of taxonomic levels. Intergeneric relationships were completely resolved: there are two sister groups, Arctium-Cousinia and Saussurea-Jurinea, which are in agreement with a morphological hypothesis. Intrageneric relationships among species of Arctium, Cousinia, and Saussurea are also well defined. Conversely, conflicting species relationships remain for Jurinea. Methodological choices significantly affected phylogenies in terms of topology, branch length, and support. Across all analyses, the phylogeny obtained using HybPiper and the strictest scheme of removing fast-evolving sites was estimated as the optimal. Regarding methodological choices, we conclude that: (1) trees obtained under the coalescence approach are topologically more congruent between them than those inferred using the concatenation approach; (2) refining treatments only improved support values under the concatenation approach; and (3) branch support values are maximized when fast-evolving sites are removed in the concatenation approach, and when a higher number of loci is analyzed in the coalescence approach.