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

Biogeography of Greater Antillean freshwater fishes, with a review of competing hypotheses

In biogeography, vicariance and long-distance dispersal are often characterised as competing scenarios. However, they are related concepts, both relying on collective geological, ecological, and phylogenetic evidence. This is illustrated by freshwater fishes, which may immigrate to islands either when freshwater connections are temporarily present and later severed (vicariance), or by unusual means when ocean gaps are crossed (long-distance dispersal). Marine barriers have a strong filtering effect on freshwater fishes, limiting immigrants to those most capable of oceanic dispersal. The roles of vicariance and dispersal are debated for freshwater fishes of the Greater Antilles. We review three active hypotheses [Cretaceous vicariance, Greater Antilles-Aves Ridge (GAARlandia), long-distance dispersal] and propose long-distance dispersal to be an appropriate model due to limited support for freshwater fish use of landspans. Greater Antillean freshwater fishes have six potential source bioregions (defined from faunal similarity): Northern Gulf of México, Western Gulf of México, Maya Terrane, Chortís Block, Eastern Panamá, and Northern South America. Faunas of the Greater Antilles are composed of taxa immigrating from many of these bioregions, but there is strong compositional disharmony between island and mainland fish faunas (>90% of Antillean species are cyprinodontiforms, compared to <10% in Northern Gulf of México and Northern South America, and ≤50% elsewhere), consistent with a hypothesis of long-distance dispersal. Ancestral-area reconstruction analysis indicates there were 16 or 17 immigration events over the last 51 million years, 14 or 15 of these by cyprinodontiforms. Published divergence estimates and evidence available for each immigration event suggests they occurred at different times and by different pathways, possibly with rafts of vegetation discharged from rivers or washed to sea during storms. If so, ocean currents likely provide critical pathways for immigration when flowing from one landmass to another. On the other hand, currents create dispersal barriers when flowing perpendicularly between landmasses. In addition to high salinity tolerance, cyprinodontiforms collectively display a variety of adaptations that could enhance their ability to live with rafts (small body size, viviparity, low metabolism, amphibiousness, diapause, self-fertilisation). These adaptations likely also helped immigrants establish island populations after arrival and to persist long term thereafter. Cichlids may have used a pseudo bridge (Nicaragua Rise) to reach the Greater Antilles. Gars (Lepisosteidae) may have crossed the Straits of Florida to Cuba, a relatively short crossing that is not a barrier to gene flow for several cyprinodontiform immigrants. Indeed, widespread distributions of Quaternary migrants (Cyprinodon, Gambusia, Kryptolebias), within the Greater Antilles and among neighbouring bioregions, imply that long-distance dispersal is not necessarily inhibitory for well-adapted species, even though it appears to be virtually impossible for all other freshwater fishes.

Disentangling historical relationships within Poeciliidae (Teleostei: Cyprinodontiformes) using ultraconserved elements

Poeciliids (Cyprinodontiformes: Poeciliidae), commonly known as livebearers, are popular fishes in the aquarium trade (e.g., guppies, mollies, swordtails) that are widely distributed in the Americas, with 274 valid species in 27 genera. This group has undergone various taxonomic changes recently, spurred by investigations using traditional genetic markers. Here we used over 1,000 ultraconserved loci to infer the relationships within Poeciliidae in the first attempt at understanding their diversification based on genome-scale data. We explore gene tree discordance and investigate potential incongruence between concatenation and coalescent inference methods. Our aim is to examine the influence of incomplete lineage sorting and reticulate evolution on the poeciliids' evolutionary history and how these factors contribute to the observed gene tree discordace. Our concatenated and coalescent phylogenomic inferences recovered four major clades within Poeciliidae. Most supra-generic level relationships we inferred were congruent with previous molecular studies, but we found some disagreements; the Middle American taxa Phallichthys and Poecilia (Mollienesia) were recovered as non-monophyletic, and unlike other recent molecular studies, we recovered Brachyrhaphis as monophyletic. Our study is the first to provide signatures of reticulate evolution in Poeciliidae at the family level; however, continued finer-scale investigations are needed to understand the complex evolutionary history of the family along with a much-needed taxonomic re-evaluation.

Congruence between ultraconserved element-based matrices and phylotranscriptomic datasets in the scorpion Tree of Life

Scorpions are ancient and historically renowned for their potent venom. Traditionally, the systematics of this group of arthropods was supported by morphological characters, until recent phylogenomic analyses (using RNAseq data) revealed most of the higher-level taxa to be non-monophyletic. While these phylogenomic hypotheses are stable for almost all lineages, some nodes have been hard to resolve due to minimal taxonomic sampling (e.g. family Chactidae). In the same line, it has been shown that some nodes in the Arachnid Tree of Life show disagreement between hypotheses generated using transcritptomes and other genomic sources such as the ultraconserved elements (UCEs). Here, we compared the phylogenetic signal of transcriptomes vs. UCEs by retrieving UCEs from new and previously published scorpion transcriptomes and genomes, and reconstructed phylogenies using both datasets independently. We reexamined the monophyly and phylogenetic placement of Chactidae, sampling an additional chactid species using both datasets. Our results showed that both sets of genome-scale datasets recovered highly similar topologies, with Chactidae rendered paraphyletic owing to the placement of Nullibrotheas allenii. As a first step toward redressing the systematics of Chactidae, we establish the family Anuroctonidae (new family) to accommodate the genus Anuroctonus.

The genomic landscape, causes, and consequences of extensive phylogenomic discordance in Old World mice and rats

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is often widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents (Old World mice and rats) - a large and ecologically diverse group that gave rise to the mouse and rat model systems. Combining new linked-read genome assemblies for seven murine species with eleven published rodent genomes, we first used ultra-conserved elements (UCEs) to infer a robust species tree. We then used whole genomes to examine finer-scale patterns of discordance and found that phylogenies built from proximate chromosomal regions had similar phylogenies. However, there was no relationship between tree similarity and local recombination rates in house mice, suggesting that genetic linkage influences phylogenetic patterns over deeper timescales. This signal may be independent of contemporary recombination landscapes. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to less discordance, while genes experiencing positive selection showed more discordant and variable phylogenetic signals. Finally, we show that assuming a single species tree can result in high error rates when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

Comparative analysis of cranial morphology in Middle-American heroine cichlids (Actinopterygii: Cichliformes)

Heroine cichlids are characterized by high morphological diversity, mainly in structures related to the capture and processing of food. The existence of ecomorphological groups has been proposed based on feeding behavior, where it is common for some phylogenetically unrelated species to show evolutionary convergence. Using geometric morphometrics and comparative phylogenetic methods, the variation in cranial morphology was evaluated for 17 species of heroine cichlids representing 5 ecomorphs. Cranial ecomorphs were recovered and significant differences were determined. Morphological variation of the ecomorphs was mainly explained by two axes: (1) the position of the mouth determined by the shape of the bones of the oral jaw and (2) the height of the head, defined by the size and position of the supraoccipital crest and the distance to the interopercle-subopercle junction. Cranial variation among species was related to phylogeny. To better understand the evolution of cranial morphology, it is necessary to evaluate the morphofunctional relationship of other anatomical structures related to feeding, as well as to increase the number of study species in each ecomorph by including other lineages.

Simultaneously collecting coding and non-coding phylogenomic data using homemade full-length cDNA probes, tested by resolving the high-level relationships of Colubridae

Resolving intractable phylogenetic relationships often requires simultaneously analyzing a large number of coding and non-coding orthologous loci. To gather both coding and non-coding data, traditional sequence capture methods require custom-designed commercial probes. Here, we present a cost-effective sequence capture method based on homemade probes, to capture thousands of coding and non-coding orthologous loci simultaneously, suitable for all organisms. This approach, called “FLc-Capture,” synthesizes biotinylated full-length cDNAs from mRNA as capture probes, eliminates the need for costly commercial probe design and synthesis. To demonstrate the utility of FLc-Capture, we prepared full-length cDNA probes from mRNA extracted from a common colubrid snake. We performed capture experiments with these homemade cDNA probes and successfully obtained thousands of coding and non-coding genomic loci from 24 Colubridae species and 12 distantly related snake species of other families. The average capture specificity of FLc-Capture across all tested snake species is 35%, similar to the previously published EecSeq method. We constructed two phylogenomic data sets, one including 1,075 coding loci (∼817,000 bp) and the other including 1,948 non-coding loci (∼1,114,000 bp), to study the phylogeny of Colubridae. Both data sets yielded highly similar and well-resolved trees, with 85% of nodes having &gt;95% bootstrap support. Our experimental tests show that FLc-Capture is a flexible, fast, and cost-effective sequence capture approach for simultaneously gathering coding and non-coding phylogenomic data sets to study intractable phylogenetic questions. We hope that this method will serve as a new data collection tool for evolutionary biologists working in the era of phylogenomics.

Paleogene emergence and evolutionary history of the Amazonian fossorial fish genus Tarumania (Teleostei: Tarumaniidae)

Tarumania walkerae is a rare fossorial freshwater fish species from the lower Rio Negro, Central Amazonia, composing the monotypic and recently described family Tarumaniidae. The family has been proposed as the sister group of Erythrinidae by both morphological and molecular studies despite distinct arrangements of the superfamily Erythrinoidea within Characiformes. Recent phylogenomic studies and time-calibrated analyses of characoid fishes have not included specimens of Tarumania in their analyses. We obtained genomic data for T. walkerae and constructed a phylogeny based on 1795 nuclear loci with 488,434 characters of ultraconserved elements (UCEs) for 108 terminals including specimens of all 22 characiform families. The phylogeny confirms the placement of Tarumaniidae as sister to Erythrinidae but differs from the morphological hypothesis in the placement of the two latter families as sister to the clade with Hemiodontidae, Cynodontidae, Serrasalmidae, Parodontidae, Anostomidae, Prochilodontidae, Chilodontidae, and Curimatidae. The phylogeny calibrated with five characoid fossils indicates that Erythrinoidea diverged from their relatives during the Late Cretaceous circa 90 Ma (108–72 Ma), and that Tarumania diverged from the most recent common ancestor of Erythrinidae during the Paleogene circa 48 Ma (66–32 Ma). The occurrence of the erythrinoid-like †Tiupampichthys in the Late Cretaceous–Paleogene formations of the El Molino Basin of Bolivia supports our hypothesis for the emergence of the modern Erythrinidae and Tarumaniidae during the Paleogene.

Prolonged morphological expansion of spiny-rayed fishes following the end-Cretaceous

Spiny-rayed fishes (Acanthomorpha) dominate modern marine habitats and account for more than a quarter of all living vertebrate species. Previous time-calibrated phylogenies and patterns from the fossil record explain this dominance by correlating the origin of major acanthomorph lineages with the Cretaceous-Palaeogene mass extinction. Here we infer a time-calibrated phylogeny using ultraconserved elements that samples 91.4% of all acanthomorph families and investigate patterns of body shape disparity. Our results show that acanthomorph lineages steadily accumulated throughout the Cenozoic and underwent a significant expansion of among-clade morphological disparity several million years after the end-Cretaceous. These acanthomorph lineages radiated into and diversified within distinct regions of morphospace that characterize iconic lineages, including fast-swimming open-ocean predators, laterally compressed reef fishes, bottom-dwelling flatfishes, seahorses and pufferfishes. The evolutionary success of spiny-rayed fishes is the culmination of multiple species-rich and phenotypically disparate lineages independently diversifying across the globe under a wide range of ecological conditions.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the firemouth cichlid, Thorichthys meeki (Perciformes: Cichlidae)

Here, we report information on the complete mitochondrial genome of the firemouth cichlid, Thorichthys meeki (Brind 1918). Illumina HiSeq genome sequencing produced the assembly of a circular mitogenome of 16,527 base pairs (bp) from T. meeki consisting of 46.8% GC nucleotides, 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a putative control region as shown in the typical teleost gene composition. The gene order of the T. meeki mitogenome was identical to that of other cichlid species. A maximum likelihood phylogenetic tree based on mitochondrial PCGs showed a close relationship of T. meeki with Thorichthys aureus (Gunther 1862) within Heroini tribe.