Phylogenomics, Lineage Diversification Rates, and the Evolution of Diadromy in Clupeiformes (Anchovies, Herrings, Sardines, and Relatives).

Migration independently evolved numerous times in animals, with a myriad of ecological and evolutionary implications. In fishes, perhaps the most extreme form of migration is diadromy, the migration between marine and freshwater environments. Key and longstanding questions are: how many times has diadromy evolved in fishes, how frequently do diadromous clades give rise to non-diadromous species, and does diadromy influence lineage diversification rates? Many diadromous fishes have large geographic ranges with constituent populations that use isolated freshwater habitats. This may limit gene flow among some populations, increasing the likelihood of speciation in diadromous lineages relative to non-diadromous lineages. Alternatively, diadromy may reduce lineage diversification rates if migration is associated with enhanced dispersal capacity that facilitates gene flow within and between populations. Clupeiformes (herrings, sardines, shads and anchovies) is a model clade for testing hypotheses about the evolution of diadromy because it includes an exceptionally high proportion of diadromous species and several independent evolutionary origins of diadromy. However, relationships among major clupeiform lineages remain unresolved and existing phylogenies sparsely sampled diadromous species, limiting the resolution of phylogenetically-informed statistical analyses. We assembled a phylogenomic dataset and used multi-species coalescent and concatenation-based approaches to generate the most comprehensive, highly-resolved clupeiform phylogeny to date, clarifying associations among several major clades and identifying recalcitrant relationships needing further examination. We determined that variation in rates of sequence evolution (heterotachy) and base-composition (non-stationarity) had little impact on our results. Using this phylogeny, we characterized evolutionary patterns of diadromy and tested for differences in lineage diversification rates between diadromous, marine, and freshwater lineages. We identified thirteen transitions to diadromy, all during the Cenozoic Era (ten origins of anadromy, two origins of catadromy, and one origin of amphidromy), and seven losses of diadromy. Two diadromous lineages rapidly generated non-diadromous species, demonstrating that diadromy is not an evolutionary dead-end. We discovered considerably faster transition rates out of diadromy than to diadromy. The largest lineage diversification rate increase in Clupeiformes was associated with a transition to diadromy, but we uncovered little statistical support for categorically faster lineage diversification rates in diadromous versus non-diadromous fishes. We propose that diadromy may increase the potential for accelerated lineage diversification, particularly in species that migrate long distances. However, this potential may only be realized in certain biogeographic contexts, such as when diadromy allows access to ecosystems in which there is limited competition from incumbent species.

Conservation of Affinity Rather Than Sequence Underlies a Dynamic Evolution of the Motif-Mediated p53/MDM2 Interaction in Ray-Finned Fishes.

The transcription factor and cell cycle regulator p53 is marked for degradation by the ubiquitin ligase MDM2. The interaction between these 2 proteins is mediated by a conserved binding motif in the disordered p53 transactivation domain (p53TAD) and the folded SWIB domain in MDM2. The conserved motif in p53TAD from zebrafish displays a 20-fold weaker interaction with MDM2, compared to the interaction in human and chicken. To investigate this apparent difference, we tracked the molecular evolution of the p53TAD/MDM2 interaction among ray-finned fishes (Actinopterygii), the largest vertebrate clade. Intriguingly, phylogenetic analyses, ancestral sequence reconstructions, and binding experiments showed that different loss-of-affinity changes in the canonical binding motif within p53TAD have occurred repeatedly and convergently in different fish lineages, resulting in relatively low extant affinities (KD = 0.5 to 5 µM). However, for 11 different fish p53TAD/MDM2 interactions, nonconserved regions flanking the canonical motif increased the affinity 4- to 73-fold to be on par with the human interaction. Our findings suggest that compensating changes at conserved and nonconserved positions within the motif, as well as in flanking regions of low conservation, underlie a stabilizing selection of "functional affinity" in the p53TAD/MDM2 interaction. Such interplay complicates bioinformatic prediction of binding and calls for experimental validation. Motif-mediated protein-protein interactions involving short binding motifs and folded interaction domains are very common across multicellular life. It is likely that the evolution of affinity in motif-mediated interactions often involves an interplay between specific interactions made by conserved motif residues and nonspecific interactions by nonconserved disordered regions.

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.

The impact of paleoclimatic changes on body size evolution in marine fishes.

Body size is an important species trait, correlating with life span, fecundity, and other ecological factors. Over Earth's geological history, climate shifts have occurred, potentially shaping body size evolution in many clades. General rules attempting to summarize body size evolution include Bergmann's rule, which states that species reach larger sizes in cooler environments and smaller sizes in warmer environments, and Cope's rule, which poses that lineages tend to increase in size over evolutionary time. Tetraodontiform fishes (including pufferfishes, boxfishes, and ocean sunfishes) provide an extraordinary clade to test these rules in ectotherms owing to their exemplary fossil record and the great disparity in body size observed among extant and fossil species. We examined Bergmann's and Cope's rules in this group by combining phylogenomic data (1,103 exon loci from 185 extant species) with 210 anatomical characters coded from both fossil and extant species. We aggregated data layers on paleoclimate and body size from the species examined, and inferred a set of time-calibrated phylogenies using tip-dating approaches for downstream comparative analyses of body size evolution by implementing models that incorporate paleoclimatic information. We found strong support for a temperature-driven model in which increasing body size over time is correlated with decreasing oceanic temperatures. On average, extant tetraodontiforms are two to three times larger than their fossil counterparts, which otherwise evolved during periods of warmer ocean temperatures. These results provide strong support for both Bergmann's and Cope's rules, trends that are less studied in marine fishes compared to terrestrial vertebrates and marine invertebrates.

Patterns of Phenotypic Evolution Associated with Marine/Freshwater Transitions in Fishes.

Evolutionary transitions between marine and freshwater ecosystems have occurred repeatedly throughout the phylogenetic history of fishes. The theory of ecological opportunity predicts that lineages that colonize species-poor regions will have greater potential for phenotypic diversification than lineages invading species-rich regions. Thus, transitions between marine and freshwaters may promote phenotypic diversification in trans-marine/freshwater fish clades. We used phylogenetic comparative methods to analyze body size data in nine major fish clades that have crossed the marine/freshwater boundary. We explored how habitat transitions, ecological opportunity, and community interactions influenced patterns of phenotypic diversity. Our analyses indicated that transitions between marine and freshwater habitats did not drive body size evolution, and there are few differences in body size between marine and freshwater lineages. We found that body size disparity in freshwater lineages is not correlated with the number of independent transitions to freshwaters. We found a positive correlation between body size disparity and overall species richness of a given area, and a negative correlation between body size disparity and diversity of closely related species. Our results indicate that the diversity of incumbent freshwater species does not restrict phenotypic diversification, but the diversity of closely related taxa can limit body size diversification. Ecological opportunity arising from colonization of novel habitats does not seem to have a major effect in the trajectory of body size evolution in trans-marine/freshwater clades. Moreover, competition with closely related taxa in freshwaters has a greater effect than competition with distantly related incumbent species.

Does genome size increase with water depth in marine fishes?

A growing body of research suggests that genome size in animals can be affected by ecological factors. Half a century ago, Ebeling et al. proposed that genome size increases with depth in some teleost fish groups and discussed a number of biological mechanisms that may explain this pattern (e.g., passive accumulation, adaptive acclimation). Using phylogenetic comparative approaches, we revisit this hypothesis based on genome size and ecological data from up to 708 marine fish species in combination with a set of large-scale phylogenies, including a newly inferred tree. We also conduct modeling approaches of trait evolution and implement a variety of regression analyses to assess the relationship between genome size and depth. Our reanalysis of Ebeling et al.'s dataset shows a weak association between these variables, but the overall pattern in their data is driven by a single clade. Although new analyses based on our "all-species" dataset resulted in positive correlations, providing some evidence that genome size evolves as a function of depth, only one subclade consistently yielded statistically significant correlations. By contrast, negative correlations are rare and nonsignificant. All in all, we find modest evidence for an increase in genome size along the depth axis in marine fishes. We discuss some mechanistic explanations for the observed trends.

Phylogenomics of Bony-Tongue Fishes (Osteoglossomorpha) Shed Light on the Craniofacial Evolution and Biogeography of the Weakly Electric Clade (Mormyridae).

Bonytongues (Osteoglossomorpha) constitute an ancient clade of teleost fishes distributed in freshwater habitats throughout the world. The group includes well-known species such as arowanas, featherbacks, pirarucus, and the weakly electric fishes in the family Mormyridae. Their disjunct distribution, extreme morphologies, and electrolocating capabilities (Gymnarchidae and Mormyridae) have attracted much scientific interest, but a comprehensive phylogenetic framework for comparative analysis is missing, especially for the species-rich family Mormyridae. Of particular interest are disparate craniofacial morphologies among mormyrids which might constitute an exceptional model system to study convergent evolution. We present a phylogenomic analysis based on 546 exons of 179 species (out of 260), 28 out of 29 genera, and all six families of extant bonytongues. Based on a recent reassessment of the fossil record of osteoglossomorphs, we inferred dates of divergence among transcontinental clades and the major groups. The estimated ages of divergence among extant taxa (e.g., Osteoglossomorpha, Osteoglossiformes, and Mormyroidea) are older than previous reports, but most of the divergence dates obtained for clades on separate continents are too young to be explained by simple vicariance hypotheses. Biogeographic analysis of mormyrids indicates that their high species diversity in the Congo Basin is a consequence of range reductions of previously widespread ancestors and that the highest diversity of craniofacial morphologies among mormyrids originated in this basin. Special emphasis on a taxon-rich representation for mormyrids revealed pervasive misalignment between our phylogenomic results and mormyrid taxonomy due to repeated instances of convergence for extreme craniofacial morphologies. Estimation of ancestral phenotypes revealed contingent evolution of snout elongation and unique projections from the lower jaw to form the distinctive Schnauzenorgan. Synthesis of comparative analyses suggests that the remarkable craniofacial morphologies of mormyrids evolved convergently due to niche partitioning, likely enabled by interactions between their exclusive morphological and electrosensory adaptations. [Africa; ancestral state estimation; diversity; exon capture; freshwater fishes; Phylogenomics.].

Exon-capture data and locus screening provide new insights into the phylogeny of flatfishes (Pleuronectoidei).

There is an extensive collection of literature on the taxonomy and phylogenetics of flatfishes (Pleuronectiformes) that extends over two centuries, but consensus on many of their evolutionary relationships remains elusive. Phylogenetic uncertainty stems from highly divergent results derived from morphological and genetic characters, and between various molecular datasets. Deciphering relationships is complicated by rapid diversification early in the Pleuronectiformes tree and an abundance of studies that incompletely and inconsistently sample taxa and genetic markers. We present phylogenies based on a genome-wide dataset (4,434 nuclear markers via exon-capture) and wide taxon sampling (86 species spanning 12 of 16 families) of the largest flatfish suborder (Pleuronectoidei). Nine different subsets of the data and two tree construction approaches (eighteen phylogenies in total) are remarkably consistent with other recent molecular phylogenies, and show strong support for the monophyly of all families included except Pleuronectidae. Analyses resolved a novel phylogenetic hypothesis for the family Rhombosoleidae as being within the Pleuronectoidea rather than the Soleoidea, and failed to support the subfamily Hippoglossinae as a monophyletic group. Our results were corroborated with evidence from previous phylogenetic studies to outline regions of persistent phylogenetic uncertainty and identify groups in need of further phylogenetic inference.

Deep reefs are not refugium for shallow-water fish communities in the southwestern Atlantic.

The deep reef refugia hypothesis (DRRH) predicts that deep reef ecosystems may act as refugium for the biota of disturbed shallow waters. Because deep reefs are among the most understudied habitats on Earth, formal tests of the DRRH remain scarce. If the DRRH is valid at the community level, the diversity of species, functions, and lineages of fish communities of shallow reefs should be encapsulated in deep reefs.We tested the DRRH by assessing the taxonomic, functional, and phylogenetic diversity of 22 Brazilian fish communities between 2 and 62 m depth. We partitioned the gamma diversity of shallow (<30 m) and deep reefs (>30 m) into independent alpha and beta components, accounted for species' abundance, and assessed whether beta patterns were mostly driven by spatial turnover or nestedness.We recorded 3,821 fishes belonging to 85 species and 36 families. Contrary to DRRH expectations, only 48% of the species occurred in both shallow and deep reefs. Alpha diversity of rare species was higher in deep reefs as expected, but alpha diversity of typical and dominant species did not vary with depth. Alpha functional diversity was higher in deep reefs only for rare and typical species, but not for dominant species. Alpha phylogenetic diversity was consistently higher in deep reefs, supporting DRRH expectations.Profiles of taxonomic, functional, and phylogenetic beta diversity indicated that deep reefs were not more heterogeneous than shallow reefs, contradicting expectations of biotic homogenization near sea surface. Furthermore, pairwise beta-diversity analyses revealed that the patterns were mostly driven by spatial turnover rather than nestedness at any depth. Conclusions. Although some results support the DRRH, most indicate that the shallow-water reef fish diversity is not fully encapsulated in deep reefs. Every reef contributes significantly to the regional diversity and must be managed and protected accordingly.

Testing the Utility of Alternative Metrics of Branch Support to Address the Ancient Evolutionary Radiation of Tunas, Stromateoids, and Allies (Teleostei: Pelagiaria).

The use of high-throughput sequencing technologies to produce genome-scale data sets was expected to settle some long-standing controversies across the Tree of Life, particularly in areas where short branches occur at deep timescales. Instead, these data sets have often yielded many well-supported but conflicting topologies, and highly variable gene-tree distributions. A variety of branch-support metrics beyond the nonparametric bootstrap are now available to assess how robust a phylogenetic hypothesis may be, as well as new methods to quantify gene-tree discordance. We applied multiple branch-support metrics to a study of an ancient group of marine fishes (Teleostei: Pelagiaria) whose interfamilial relationships have proven difficult to resolve due to a rapid accumulation of lineages very early in its history. We analyzed hundreds of loci including published ultraconserved elements and newly generated exonic data along with their flanking regions to represent all 16 extant families for more than 150 out of 284 valid species in the group. Branch support was typically lower at inter- than intra-familial relationships regardless of the type of marker used. Several nodes that were highly supported with bootstrap had a very low site and gene-tree concordance, revealing underlying conflict. Despite this conflict, we were able to identify four consistent interfamilial clades, each comprised of two or three families. Combining exons with their flanking regions also produced increased branch lengths at the deep branches of the pelagiarian tree. Our results demonstrate the limitations of employing current metrics of branch support and species-tree estimation when assessing the confidence of ancient evolutionary radiations and emphasize the necessity to embrace alternative measurements to explore phylogenetic uncertainty and discordance in phylogenomic data sets.[Concatenation; exons; introns; phylogenomics; species-tree methods; target capture.].

Exon probe sets and bioinformatics pipelines for all levels of fish phylogenomics.

Exon markers have a long history of use in phylogenetics of ray-finned fishes, the most diverse clade of vertebrates with more than 35,000 species. As the number of published genomes increases, it has become easier to test exons and other genetic markers for signals of ancient duplication events and filter out paralogues that can mislead phylogenetic analysis. We present seven new probe sets for current target-capture phylogenomic protocols that capture 1,104 exons explicitly filtered for paralogues using gene trees. These seven probe sets span the diversity of teleost fishes, including four sets that target five hyperdiverse percomorph clades which together comprise ca. 17,000 species (Carangaria, Ovalentaria, Eupercaria, and Syngnatharia + Pelagiaria combined). We additionally included probes to capture legacy nuclear exons and mitochondrial markers that have been commonly used in fish phylogenetics (despite some exons being flagged for paralogues) to facilitate integration of old and new molecular phylogenetic matrices. We tested these probes experimentally for 56 fish species (eight species per probe set) and merged new exon-capture sequence data into an existing data matrix of 1,104 exons and 300 ray-finned fish species. We provide an optimized bioinformatics pipeline to assemble exon capture data from raw reads to alignments for downstream analysis. We show that legacy loci with known paralogues are at risk of assembling duplicated sequences with target-capture, but we also assembled many useful orthologous sequences that can be integrated with many PCR-generated matrices. These probe sets are a valuable resource for advancing fish phylogenomics because targeted exons can easily be extracted from increasingly available whole genome and transcriptome data sets, and also may be integrated with existing PCR-based exon and mitochondrial data.

Evolutionary determinism and convergence associated with water-column transitions in marine fishes.

Repeatable, convergent outcomes are prima facie evidence for determinism in evolutionary processes. Among fishes, well-known examples include microevolutionary habitat transitions into the water column, where freshwater populations (e.g., sticklebacks, cichlids, and whitefishes) recurrently diverge toward slender-bodied pelagic forms and deep-bodied benthic forms. However, the consequences of such processes at deeper macroevolutionary scales in the marine environment are less clear. We applied a phylogenomics-based integrative, comparative approach to test hypotheses about the scope and strength of convergence in a marine fish clade with a worldwide distribution (snappers and fusiliers, family Lutjanidae) featuring multiple water-column transitions over the past 45 million years. We collected genome-wide exon data for 110 (∼80%) species in the group and aggregated data layers for body shape, habitat occupancy, geographic distribution, and paleontological and geological information. We also implemented approaches using genomic subsets to account for phylogenetic uncertainty in comparative analyses. Our results show independent incursions into the water column by ancestral benthic lineages in all major oceanic basins. These evolutionary transitions are persistently associated with convergent phenotypes, where deep-bodied benthic forms with truncate caudal fins repeatedly evolve into slender midwater species with furcate caudal fins. Lineage diversification and transition dynamics vary asymmetrically between habitats, with benthic lineages diversifying faster and colonizing midwater habitats more often than the reverse. Convergent ecological and functional phenotypes along the benthic-pelagic axis are pervasive among different lineages and across vastly different evolutionary scales, achieving predictable high-fitness solutions for similar environmental challenges, ultimately demonstrating strong determinism in fish body-shape evolution.

New insights about species delimitation in red snappers (Lutjanus purpureus and L. campechanus) using multilocus data.

Lutjanus campechanus and Lutjanus purpureus are two commercially important lutjanid fishes (snappers) with non-sympatric distribution throughout Western Atlantic. Even though both taxa have traditionally been regarded as valid species, their taxonomic status remains under debate. In the present study, we used phylogeographic approaches and molecular methods of species delimitation to elucidate the taxonomic issues between both species, based on 1478 base pairs from four genomic regions. We found haplotypes shared between the two species, particularly in relation to nuclear DNA (nuDNA) sequences. The molecular delimitation of species supported the discrimination of L. purpureus and L. campechanus as distinct evolutionary units. Nonetheless, a unidirectional gene flow was found from L. campechanus towards L. purpureus. Therefore, it seems plausible to infer that L. campechanus and L. purpureus are two evolutionary units in which the apparent sharing of haplotypes should be driven by introgression.

Biogeography, habitat transitions and hybridization in a radiation of South American silverside fishes revealed by mitochondrial and genomic RAD data.

Rivers and lake systems in the southern cone of South America have been widely influenced by historical glaciations, carrying important implications for the evolution of aquatic organisms, including prompting transitions between marine and freshwater habitats and by triggering hybridization among incipient species via waterway connectivity and stream capture events. Silverside fishes (Odontesthes) in the region comprise a radiation of 19 marine and freshwater species that have been hypothesized on the basis of morphological or mitochondrial DNA data to have either transitioned repeatedly into continental waters from the sea or colonized marine habitats following freshwater diversification. New double digest restriction-site associated DNA data presented here provide a robust framework to investigate the biogeographical history of and habitat transitions in Odontesthes. We show that Odontesthes silversides originally diversified in the Pacific but independently colonized the Atlantic three times, producing three independent marine-to-freshwater transitions. Our results also indicate recent introgression of marine mitochondrial haplotypes into two freshwater clades, with more recurring instances of hybridization among Atlantic- versus Pacific-slope species. In Pacific freshwater drainages, hybridization with a marine species appears to be geographically isolated and may be related to glaciation events. Substantial structural differences of estuarine gradients between these two geographical areas may have influenced the frequency, intensity and evolutionary effects of hybridization events.

On trends and patterns in macroevolution: Williston's law and the branchiostegal series of extant and extinct osteichthyans.

BACKGROUND: The branchiostegal series consists of an alignment of bony elements in the posterior portion of the skull of osteichthyan vertebrates. We trace the evolution of the number of elements in a comprehensive survey that includes 440 extant and 66 extinct species. Using a newly updated actinopterygian tree in combination with phylogenetic comparative analyses, we test whether osteichthyan branchiostegals follow an evolutionary trend under 'Williston's law', which postulates that osteichthyan lineages experienced a reduction of bony elements over time. RESULTS: We detected no overall macroevolutionary trend in branchiostegal numbers, providing no support for 'Williston's law'. This result is robust to the subsampling of palaeontological data, but the estimation of the model parameters is much more ambiguous. CONCLUSIONS: We find substantial evidence for a macroevolutionary dynamic favouring an 'early burst' of trait evolution over alternative models. Our study highlights the challenges of accurately reconstructing macroevolutionary dynamics even with large amounts of data about extant and extinct taxa.

Genomics overrules mitochondrial DNA, siding with morphology on a controversial case of species delimitation.

Species delimitation is a major quest in biology and is essential for adequate management of the organismal diversity. A challenging example comprises the fish species of red snappers in the Western Atlantic. Red snappers have been traditionally recognized as two separate species based on morphology: Lutjanus campechanus (northern red snapper) and L. purpureus (southern red snapper). Recent genetic studies using mitochondrial markers, however, failed to delineate these nominal species, leading to the current lumping of the northern and southern populations into a single species ( L. campechanus). This decision carries broad implications for conservation and management as red snappers have been commercially over-exploited across the Western Atlantic and are currently listed as vulnerable. To address this conflict, we examine genome-wide data collected throughout the range of the two species. Population genomics, phylogenetic and coalescent analyses favour the existence of two independent evolutionary lineages, a result that confirms the morphology-based delimitation scenario in agreement with conventional taxonomy. Despite finding evidence of introgression in geographically neighbouring populations in northern South America, our genomic analyses strongly support isolation and differentiation of these species, suggesting that the northern and southern red snappers should be treated as distinct taxonomic entities.

Chinchaysuyoa, a new genus of the fish family Ariidae (Siluriformes), with a redescription of Chinchaysuyoa labiata from Ecuador and a new species description from Peru.

In recent years, morphological and molecular studies have improved our understanding about the relationships and classification schemes of the marine catfishes of the family Ariidae. A taxonomic issue that is still contentious concerns the limits and status of the freshwater Neotropical ariid diversity, in particular the species in the genus Potamarius. The delimitation of Potamarius is currently uncertain given the disjunct distribution of the species in Mesoamerica (Potamarius izabalensis, P. nelsoni and P. usumacintae, from Lake Izabal and Usumacinta River basins in Mexico to Guatemala) and Brazil (P. grandoculis, from coastal lakes in southeastern Brazil). The freshwater Arius labiatus and Hexanematichthys henni from the Peripa and Daule rivers in Ecuador that drain to the Eastern Pacific (EP), have also at times been listed as species inquirenda in Potamarius. Here, we redescribe Arius labiatus, redefine the taxonomic status of Hexanematichthys henni, as junior synonym of Arius labiatus, and describe a new species from Peru that is closely related to Arius labiatus. Based on morphological and molecular phylogenetic evidence, we also describe a new genus (Chinchaysuyoa) for the two South American species.

Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers.

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome-wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans-Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter-riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.

Phylogenomic incongruence, hypothesis testing, and taxonomic sampling: The monophyly of characiform fishes.

Phylogenomic studies using genome-wide datasets are quickly becoming the state of the art for systematics and comparative studies, but in many cases, they result in strongly supported incongruent results. The extent to which this conflict is real depends on different sources of error potentially affecting big datasets (assembly, stochastic, and systematic error). Here, we apply a recently developed methodology (GGI or gene genealogy interrogation) and data curation to new and published datasets with more than 1000 exons, 500 ultraconserved element (UCE) loci, and transcriptomic sequences that support incongruent hypotheses. The contentious non-monophyly of the order Characiformes proposed by two studies is shown to be a spurious outcome induced by sample contamination in the transcriptomic dataset and an ambiguous result due to poor taxonomic sampling in the UCE dataset. By exploring the effects of number of taxa and loci used for analysis, we show that the power of GGI to discriminate among competing hypotheses is diminished by limited taxonomic sampling, but not equally sensitive to gene sampling. Taken together, our results reinforce the notion that merely increasing the number of genetic loci for a few representative taxa is not a robust strategy to advance phylogenetic knowledge of recalcitrant groups. We leverage the expanded exon capture dataset generated here for Characiformes (206 species in 23 out of 24 families) to produce a comprehensive phylogeny and a revised classification of the order.

Comprehensive phylogeny of ray-finned fishes (Actinopterygii) based on transcriptomic and genomic data.

Our understanding of phylogenetic relationships among bony fishes has been transformed by analysis of a small number of genes, but uncertainty remains around critical nodes. Genome-scale inferences so far have sampled a limited number of taxa and genes. Here we leveraged 144 genomes and 159 transcriptomes to investigate fish evolution with an unparalleled scale of data: >0.5 Mb from 1,105 orthologous exon sequences from 303 species, representing 66 out of 72 ray-finned fish orders. We apply phylogenetic tests designed to trace the effect of whole-genome duplication events on gene trees and find paralogy-free loci using a bioinformatics approach. Genome-wide data support the structure of the fish phylogeny, and hypothesis-testing procedures appropriate for phylogenomic datasets using explicit gene genealogy interrogation settle some long-standing uncertainties, such as the branching order at the base of the teleosts and among early euteleosts, and the sister lineage to the acanthomorph and percomorph radiations. Comprehensive fossil calibrations date the origin of all major fish lineages before the end of the Cretaceous.