Genome-wide interrogation advances resolution of recalcitrant groups in the tree of life

Evolution across the adaptive landscape in a hyperdiverse beetle radiation

The extraordinary diversification of beetles on Earth is a textbook example of adaptive evolution. Yet, the tempo and drivers of this super-radiation remain largely unclear. Here, we address this problem by investigating macroevolutionary dynamics in darkling beetles (Coleoptera: Tenebrionidae), one of the most ecomorphologically diverse beetle families (with over 30,000 species). Using multiple genomic datasets and analytical approaches, we resolve the long-standing inconsistency over deep relationships in the family. In conjunction with a landmark-based dataset of body shape morphology, we show that the evolutionary history of darkling beetles is marked by ancient rapid radiations, frequent ecological transitions, and rapid bursts of morphological diversification. On a global scale, our analyses uncovered a significant pulse of phenotypic diversification proximal to the Cretaceous-Palaeogene (K/Pg) mass extinction and convergence of body shape associated with recurrent ecological specializations. On a regional scale, two major Australasian radiations, the Adeliini and the Heleine clade, exhibited contrasting patterns of ecomorphological diversification, representing phylogenetic niche conservatism versus adaptive radiation. Our findings align with the Simpsonian model of adaptive evolution across the macroevolutionary landscape and highlight a significant role of ecological opportunity in driving the immense ecomorphological diversity in a hyperdiverse beetle group.

The role of abiotic and biotic factors in the unequal body shape diversification of a Gondwanan fish radiation (Otophysi: Characiformes)

Understanding why some clades diversify greatly while others do not is a major goal of evolutionary biology. Both abiotic and biotic factors are important in driving unequal morphological diversity across the tree of life. However, few studies have quantified how abiotic habitat and community composition differences influence unequal morphological diversification in spatiotemporally diffuse radiations. We use geometric morphometrics, abiotic habitat data generated by Geographic Information Systems (GIS) analyses, evolutionary simulations, and phylogenetic comparative methods to determine whether random evolution, habitat variation, competition for niches, or a combination of factors influenced the unequal body shape diversity of a Gondwanan freshwater fish radiation. We find that Neotropical characiform lineages, which exhibit substantially more body shape diversity than their African counterparts, occupy significantly more slope and elevation habitats than African lineages. Differences in habitat occupation between the continental radiations occur through a combination of likely competition with cypriniform fishes in Africa restricting access to higher slope and elevation habitats and significantly more low-elevation and slope habitats available in the Neotropics. Our findings suggest that spatiotemporally widespread radiations, like Characiformes, do not diversify across homogenized habitats and biotic assemblages, with differences in community structure and physical habitat playing an important role in driving unequal morphological diversification.

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.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

Marine and freshwater fishes of Alabama: a revised checklist and discussion of taxonomic issues

Checklists are fundamental and important tools for organizing information about biodiversity that provide a basis for conservation and additional scientific research. While Alabama is recognized as an aquatic biodiversity hotspot with the highest native freshwater fish diversity in the contiguous United States, we currently lack an up-to-date list of the states fishes. In particular, much has changed over the past ~20 years regarding our knowledge of fishes from Alabama and the Mobile River Basin, rendering past comprehensive treatments by Mettee et al. (1996) and Boschung and Mayden (2004) out of date. Here, we provide a revised checklist of marine and freshwater fishes known from the coastal and inland waters of Alabama that includes 463 species (335 primarily freshwater fishes, and 128 marine or diadromous fishes) in 35 orders, 78 families, and 176 genera. Extant, extirpated, and extinct species are included, as are putative candidate species. The checklist is based on prior work, searches of the literature and online sources, as well as parsing a large compilation of >140,000 fish records for Alabama and the Mobile River Basin from 37 data providers in the global Fishnet2 database (www.fishnet2.net) and >4000 marine survey records from the SEAMAP database (https://www.gsmfc.org/seamap.php). After editing and quality control checks, the final combined database contained 144,215 collection records, ~95% of which were georeferenced. We discuss the species descriptions, nomenclatural changes, and updates to marine species that account for changes to the state list, and we close with a discussion of ~13 candidate species forms that remain undescribed, which represent outstanding taxonomic issues in need of further research attention.

Unravelling the genomic features, phylogeny and genetic basis of tooth ontogenesis in Characiformes through analysis of four genomes

Characiformes is a diverse and evolutionarily significant order of freshwater fish encompassing over 2,300 species. Despite its diversity, our understanding of Characiformes' evolutionary relationships and adaptive mechanisms is limited due to insufficient genome sequences. In this study, we sequenced and assembled the genomes of four Characiformes species, three of which were chromosome-level assemblies. Our analyses revealed dynamic changes in gene family evolution, repeat sequences and variations in chromosomal collinearity within these genomes. With the assembled genomes, we were not only able to elucidate the evolutionary relationship of the four main orders in Otophysi but also indicated Characiformes as the paraphyletic group. Comparative genomic analysis with other available fish genomes shed light on the evolution of genes related to tooth development in Characiformes. Notably, variations in the copy number of secretory calcium-binding phosphoproteins (SCPP) genes were observed among different orders of Otophysi, indicating their potential contribution to the diversity of tooth types. Our study offers invaluable genome sequences and novel insights into Characiformes' evolution, paving the way for further genomic and evolutionary research in fish.

Chromosome-Aware Phylogenomics of Assassin Bugs (Hemiptera: Reduvioidea) Elucidates Ancient Gene Conflict

Though the phylogenetic signal of loci on sex chromosomes can differ from those on autosomes, chromosomal-level genome assemblies for nonvertebrates are still relatively scarce and conservation of chromosomal gene content across deep phylogenetic scales has therefore remained largely unexplored. We here assemble a uniquely large and diverse set of samples (17 anchored hybrid enrichment, 24 RNA-seq, and 70 whole-genome sequencing samples of variable depth) for the medically important assassin bugs (Reduvioidea). We assess the performance of genes based on multiple features (e.g., nucleotide vs. amino acid, nuclear vs. mitochondrial, and autosomal vs. X chromosomal) and employ different methods (concatenation and coalescence analyses) to reconstruct the unresolved phylogeny of this diverse (∼7,000 spp.) and old (>180 Ma) group. Our results show that genes on the X chromosome are more likely to have discordant phylogenies than those on autosomes. We find that the X chromosome conflict is driven by high gene substitution rates that impact the accuracy of phylogenetic inference. However, gene tree clustering showed strong conflict even after discounting variable third codon positions. Alternative topologies were not particularly enriched for sex chromosome loci, but spread across the genome. We conclude that binning genes to autosomal or sex chromosomes may result in a more accurate picture of the complex evolutionary history of a clade.

Confusion will be my epitaph: genome-scale discordance stifles phylogenetic resolution of Holothuroidea

Sea cucumbers (Holothuroidea) are a diverse clade of echinoderms found from intertidal waters to the bottom of the deepest oceanic trenches. Their reduced skeletons and limited number of phylogenetically informative traits have long obfuscated morphological classifications. Sanger-sequenced molecular datasets have also failed to constrain the position of major lineages. Noteworthy, topological uncertainty has hindered a resolution for Neoholothuriida, a highly diverse clade of Permo-Triassic age. We perform the first phylogenomic analysis of Holothuroidea, combining existing datasets with 13 novel transcriptomes. Using a highly curated dataset of 1100 orthologues, our efforts recapitulate previous results, struggling to resolve interrelationships among neoholothuriid clades. Three approaches to phylogenetic reconstruction (concatenation under both site-homogeneous and site-heterogeneous models, and coalescent-aware inference) result in alternative resolutions, all of which are recovered with strong support and across a range of datasets filtered for phylogenetic usefulness. We explore this intriguing result using gene-wise log-likelihood scores and attempt to correlate these with a large set of gene properties. While presenting novel ways of exploring and visualizing support for alternative trees, we are unable to discover significant predictors of topological preference, and our efforts fail to favour one topology. Neoholothuriid genomes seem to retain an amalgam of signals derived from multiple phylogenetic histories.

Phylogenomics, evolution of trophic traits and divergence times of hemiodontid fishes (Ostariophysi: Characiformes)

The South American characiform family Hemiodontidae comprises five genera and 34 species. The family lacks comprehensive phylogenetic hypotheses resolving its species relationships. The studies that addressed these questions exhibited a narrow taxon sampling or used single-locus markers. Herein we surveyed hundreds of ultraconserved elements (UCEs) loci to provide the first molecular phylogenetic hypothesis and divergence time estimates for hemiodontids encompassing all its genera and most species (27 of the 34 valid species). We also tracked the history of the protractile upper jaw in the genera Argonectes and Bivibranchia across the recovered phylogenies through ancestral state reconstruction. Our results corroborate the monophyly of Hemiodontidae and the genera Argonectes and Bivibranchia in all phylogenetic methods with maximum clade support. The genera Anodus and Hemiodus were not monophyletic because Anodus elongatus was sister to the monotypic Micromischodus instead of A. orinocensis, and H. immaculatus did not form a clade with its other congeners, but instead was sister to the clade including Anodus and Micromischodus. All remaining species of Hemiodus were placed together into a monophyletic group, where they were arranged into four major subclades. The relationship in the family is summarised as: (Bivibranchia, (Argonectes, ((H. immaculatus, (Anodus, Micromischodus)), Hemiodus clade))), in discordance with the morphological phylogeny that placed all genera monophyletic and resolved the family as: ((Anodus, Micromischodus), (Hemiodus, (Argonectes, Bivibranchia))). The origin of Hemiodontidae was estimated from the Late Cretaceous to the Middle Paleogene, with the mean age in the Paleocene, while the origin of most hemiodontid genera except Bivibranchia occurred in the Miocene. Unordered parsimony and likelihood reconstruction indicates that Argonectes and Bivibranchia developed their protractile upper jaw independently.

Exploring Conflicts in Whole Genome Phylogenetics: A Case Study Within Manakins (Aves: Pipridae)

Some phylogenetic problems remain unresolved even when large amounts of sequence data are analyzed and methods that accommodate processes such as incomplete lineage sorting are employed. In addition to investigating biological sources of phylogenetic incongruence, it is also important to reduce noise in the phylogenomic dataset by using appropriate filtering approach that addresses gene tree estimation errors. We present the results of a case study in manakins, focusing on the very difficult clade comprising the genera Antilophia and Chiroxiphia. Previous studies suggest that Antilophia is nested within Chiroxiphia, though relationships among Antilophia+Chiroxiphia species have been highly unstable. We extracted more than 11,000 loci (ultra-conserved elements and introns) from whole genomes and conducted analyses using concatenation and multispecies coalescent methods. Topologies resulting from analyses using all loci differed depending on the data type and analytical method, with 2 clades (Antilophia+Chiroxiphia and Manacus+Pipra+Machaeopterus) in the manakin tree showing incongruent results. We hypothesized that gene trees that conflicted with a long coalescent branch (e.g., the branch uniting Antilophia+Chiroxiphia) might be enriched for cases of gene tree estimation error, so we conducted analyses that either constrained those gene trees to include monophyly of Antilophia+Chiroxiphia or excluded these loci. While constraining trees reduced some incongruence, excluding the trees led to completely congruent species trees, regardless of the data type or model of sequence evolution used. We found that a suite of gene metrics (most importantly the number of informative sites and likelihood of intralocus recombination) collectively explained the loci that resulted in non-monophyly of Antilophia+Chiroxiphia. We also found evidence for introgression that may have contributed to the discordant topologies we observe in Antilophia+Chiroxiphia and led to deviations from expectations given the multispecies coalescent model. Our study highlights the importance of identifying factors that can obscure phylogenetic signal when dealing with recalcitrant phylogenetic problems, such as gene tree estimation error, incomplete lineage sorting, and reticulation events. [Birds; c-gene; data type; gene estimation error; model fit; multispecies coalescent; phylogenomics; reticulation].

Multiple Outgroups Can Cause Random Rooting in Phylogenomics

Outgroup selection has been a major challenge since the rise of phylogenetics, and it has remained so in the phylogenomic era. Our goal here is to use large phylogenomic animal datasets to examine the impact of outgroup selection on the final topology. The results of our analyses further solidify the fact that distant outgroups can cause random rooting, and that this holds for concatenated and coalescent-based methods. The results also indicate that the standard practice of using multiple outgroups often causes random rooting. Most researchers go out of their way to get multiple outgroups, as this has been standard practice for decades. Based on our findings, this practice should stop. Instead, our results suggest that a single (most closely) related relative should be selected as the outgroup, unless all outgroups are roughly equally closely related to the ingroup.

Climate-influenced boreotropical survival and rampant introgressions explain the thriving of New World grapes in the north temperate zone

The north temperate region was characterized by a warm climate and a rich thermophilic flora before the Eocene, but early diversifications of the temperate biome under global climate change and biome shift remain uncertain. Moreover, it is becoming clear that hybridization/introgression is an important driving force of speciation in plant diversity. Here, we applied analyses from biogeography and phylogenetic networks to account for both introgression and incomplete lineage sorting based on genomic data from the New World Vitis, a charismatic component of the temperate North American flora with known and suspected gene flow among species. Biogeographic inference and fossil evidence suggest that the grapes were widely distributed from North America to Europe during the Paleocene to the Eocene, followed by widespread extinction and survival of relicts in the tropical New World. During the climate warming in the early Miocene, a Vitis ancestor migrated northward from the refugia with subsequent diversification in the North American region. We found strong evidence for widespread incongruence and reticulate evolution among nuclear genes within both recent and ancient lineages of the New World Vitis. Furthermore, the organellar genomes showed strong conflicts with the inferred species tree from the nuclear genomes. Our phylogenomic analyses provided an important assessment of the wide occurrence of reticulate introgression in the New World Vitis, which potentially represents one of the most important mechanisms for the diversification of Vitis species in temperate North America and even the entire temperate Northern Hemisphere. The scenario we report here may be a common model of temperate diversification of flowering plants adapted to the global climate cooling and fluctuation in the Neogene.

Integrated phylogenomic analyses unveil reticulate evolution in Parthenocissus (Vitaceae), highlighting speciation dynamics in the Himalayan-Hengduan Mountains

Hybridization caused by frequent environmental changes can lead both to species diversification (speciation) and to speciation reversal (despeciation), but the latter has rarely been demonstrated. Parthenocissus, a genus with its trifoliolate lineage in the Himalayan-Hengduan Mountains (HHM) region showing perplexing phylogenetic relationships, provides an opportunity for investigating speciation dynamics based on integrated evidence. We investigated phylogenetic discordance and reticulate evolution in Parthenocissus based on rigorous analyses of plastome and transcriptome data. We focused on reticulations in the trifoliolate lineage in the HHM region using a population-level genome resequencing dataset, incorporating evidence from morphology, distribution, and elevation. Comprehensive analyses confirmed multiple introgressions within Parthenocissus in a robust temporal-spatial framework. Around the HHM region, at least three hybridization hot spots were identified, one of which showed evidence of ongoing speciation reversal. We present a solid case study using an integrative methodological approach to investigate reticulate evolutionary history and its underlying mechanisms in plants. It demonstrates an example of speciation reversal through frequent hybridizations in the HHM region, which provides new perspectives on speciation dynamics in mountainous areas with strong topographic and environmental heterogeneity.

Monsoon boosted radiation of the endemic East Asian carps

Major historical events often trigger the rapid flourishing of a few lineages, which in turn shape established biodiversity patterns. How did this process occur and develop? This study provides a window into this issue. The endemic East Asian carps (EEAC) dominated the ichthyofauna of East Asia and exhibited a high degree of adaptation to monsoonal river-lake ecosystems. A series of evidence, including ecogeography, phylogenetics, and macroevolution, suggests that the EEAC is a lineage that arose with the East Asian monsoon and thrived intimately with subsequent monsoon activities. We further deduce the evolution of the EEAC and find that a range of historical events in the monsoon setting (e.g., marine transgression and regression and glacial-interglacial cycle) have further reshaped the distribution patterns of EEAC's members. Comparative genomics analyses reveal that introgressions during the initial period of EEAC radiation and innovations in the regulation of the brain and nervous system may have aided their adaptation to river-lake ecosystems in a monsoon setting, which boosted radiation. Overall, this study strengthens knowledge of the evolutionary patterns of freshwater fishes in East Asia and provides a model case for understanding the impact of major historical events on the evolution of biota.

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.

Integrative analysis reveals cryptic speciation linked to habitat differentiation within Albanian populations of the anomalous blues (Lepidoptera, Lycaenidae, Polyommatus Latreille, 1804)

The Balkan Peninsula is one of the greatest hotspots for biodiversity in Europe. While the region has been investigated thoroughly, some parts remain understudied and may still harbour undiscovered diversity, even in well-studied organisms such as Lepidoptera. Here we investigated the group of the so-called anomalous blue butterflies, also known as 'brown complex' of the subgenus Agrodiaetus Hübner, 1822 including the taxa of the entire Polyommatusaroaniensis (Brown, 1976) species complex. This species complex is distributed in the southern part of the Balkan Peninsula and known to be represented by three closely related allopatric species, differentiated by their chromosome numbers (n) and mitochondrial (mt) DNA. These are P.aroaniensis sensu stricto (Southern Greece, Peloponnese, n=47-48; mt haplogroup aroa1), P.timfristos Lukhtanov, Vishnevskaya et Shapoval, 2016 (Central Greece, Attika, n=38, aroa2) and P.orphicus Kolev, 2005 (North-Eastern Greece, Southern Bulgaria, n=41-42, orph1). Based on an analysis of chromosomal, molecular and morphological markers, we demonstrate that a fourth taxon of this species complex exists in Albania. This taxon possesses the mt haplogroup aroa3, which is the most differentiated within the entire P.aroaniensis species complex, and the karyotype (n=42-43), which differs by one fixed chromosome fission from P.orphicus. The Albanian taxon seems to be ecologically specialised (habitat on dark-coloured, ophiolitic substrate soils) and differs in colouration (wing reflectance) from the others taxa of the P.aroaniensis species group. Based on the evidence here presented and following the current view of the taxonomy of the group, we propose considering the Albanian taxon as a new species, here described as Polyommatuslurae sp. nov. At the contact zone between the new species and P.orphicus, in addition to typical ones, we detected specimens with haplogroup orph2, karyotype n=43 and intermediate morphology, which seem to represent P.lurae × P.orphicus hybrids.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

Phylogenomics and classification of Notropis and related shiners (Cypriniformes: Leuciscidae) and the utility of exon capture on lower taxonomic groups

North American minnows of the Shiner Clade, within the family Leuciscidae, represent one of the most taxonomically complex clades of the order Cypriniformes due to the large number of taxa coupled with conserved morphologies. Species within this clade were moved between genera and subgenera until the community decided to lump many of the unclassified taxa with similar morphologies into one genus, Notropis, which has held up to 325 species. Despite phylogentic studies that began to re-elevate some genera merged into Notropis, such as Cyprinella, Luxilus, Lythrurus, and Pteronotropis, the large genus Notropis remained as a taxonomic repository for many shiners of uncertain placement. Recent molecular advances in sequencing technologies have provided the opportunity to re-examine the Shiner Clade using phylogenomic markers. Using a fish probe kit, we sequenced 90 specimens in 87 species representing 16 genera included in the Shiner Clade, with a resulting dataset of 1,004 loci and 286,455 base pairs. Despite the large dataset, only 32,349 bp (11.29%) were phylogenetically informative. In our maximum likelihood tree, 78% of nodes are 100% bootstrap supported demonstrating the utility of the phylogenomic markers at lower taxonomic levels. Unsurprisingly, species within Notropis as well as Hudsonius, Luxilus, and Alburnops are not resolved as monophyletic groups. Cyprinella is monophyletic if Cyprinella callistia is excluded, and Pteronotropis is monophyletic if it includes Hudsonius cummingsae. Taxonomic changes we propose are: restriction of species included in Alburnops and Notropis, elevation of the subgenus Hydrophlox, expansion of species included in Miniellus, movement of Hudsonius cummingsae to Pteronotropis, and resurrection of the genera Coccotis and Paranotropis. We additionally had two specimens of three species, Notropis atherinoides, Ericymba amplamala, and Pimephales vigilax and found signficant differences between the localities (1,086, 1,424, and 845 nucleotides respectively).

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

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.

Phylogenomic analysis of the parrots of the world distinguishes artifactual from biological sources of gene tree discordance

Gene tree discordance is expected in phylogenomic trees and biological processes are often invoked to explain it. However, heterogeneous levels of phylogenetic signal among individuals within datasets may cause artifactual sources of topological discordance. We examined how the information content in tips and subclades impacts topological discordance in the parrots (Order: Psittaciformes), a diverse and highly threatened clade of nearly 400 species. Using ultraconserved elements from 96% of the clade's species-level diversity, we estimated concatenated and species trees for 382 ingroup taxa. We found that discordance among tree topologies was most common at nodes dating between the late Miocene and Pliocene, and often at the taxonomic level of genus. Accordingly, we used two metrics to characterize information content in tips and assess the degree to which conflict between trees was being driven by lower quality samples. Most instances of topological conflict and non-monophyletic genera in the species tree could be objectively identified using these metrics. For subclades still discordant after tip-based filtering, we used a machine learning approach to determine whether phylogenetic signal or noise was the more important predictor of metrics supporting the alternative topologies. We found that when signal favored one of the topologies, noise was the most important variable in poorly performing models that favored the alternative topology. In sum, we show that artifactual sources of gene tree discordance, which are likely a common phenomenon in many datasets, can be distinguished from biological sources by quantifying the information content in each tip and modeling which factors support each topology.

First cytogenetic data on Afrotropical lutefishes (Citharinidae) in the light of karyotype evolution in Characiformes

The Afrotropical lutefish family Citharinidae (Citharinoidei, Characiformes) comprises three genera with eight species in total. Although Citharinidae have been studied in terms of taxonomy and systematics, no cytogenetic information was available for any representative of the family. Furthermore, only one species out of 116 in Citharinoidei (Distichodusaffinis Günther, 1873) has been studied cytogenetically. Here, we report the karyotypes of Citharinuscitharus (Geoffroy St. Hilaire, 1809) from West Africa and Citharinuslatus Müller et Troschel, 1844 from Northeast Africa. The former has the diploid chromosome number 2n = 40 and the fundamental number FN = 80, while the latter has 2n = 44 and FN = 88. Hence, these karyotypes consist exclusively of bi-armed chromosomes. Such karyotypes were previously found in D.affinis and in many lineages of Neotropical species of another suborder of Characiformes, Characoidei. In contrast, the karyotypes dominated by uni-armed elements are typical for a number of phylogenetically basal lineages of Afrotropical and Neotropical Characoidei. We discuss the importance of our data on Citharinidae for the understanding of the karyotype evolution within the order Characiformes.

Divergent cis-regulatory evolution underlies the convergent loss of sodium channel expression in electric fish

South American and African weakly electric fish independently evolved electric organs from muscle. In both groups, a voltage-gated sodium channel gene independently lost expression from muscle and gained it in the electric organ, allowing the channel to become specialized for generating electric signals. It is unknown how this voltage-gated sodium channel gene is targeted to muscle in any vertebrate. We describe an enhancer that selectively targets sodium channel expression to muscle. Next, we demonstrate how the loss of this enhancer, but not trans-activating factors, drove the loss of sodium channel gene expression from muscle in South American electric fish. While this enhancer is also altered in African electric fish, key transcription factor binding sites and enhancer activity are retained, suggesting that the convergent loss of sodium channel expression from muscle in these two electric fish lineages occurred via different processes.

New phylogenetic insights into the African catfish families Mochokidae and Austroglanididae

Several hundred catfish species (order: Siluriformes) belonging to 11 families inhabit Africa, of which at least six families are endemic to the continent. Although four of those families are well-known to belong to the 'Big-Africa clade', no previous study has addressed the phylogenetic placement of the endemic African catfish family Austroglanididae in a comprehensive framework with molecular data. Furthermore, interrelationships within the 'Big-Africa clade', including the most diverse family Mochokidae, remain unclear. This study was therefore designed to help reconstruct inter- and intrarelationships of all currently valid mochokid genera, to infer their position within the 'Big Africa clade' and to establish a first molecular phylogenetic hypothesis of the relationships of the enigmatic Austroglanididae within the Siluriformes. We assembled a comprehensive mitogenomic dataset comprising all protein coding genes and representing almost all recognized catfish families (N = 33 of 39) with carefully selected species (N = 239). We recovered the monophyly of the previously identified multifamily clades 'Big Asia' and 'Big Africa' and determined Austroglanididae to be closely related to Pangasiidae, Ictaluroidea and Ariidae. Mochokidae was recovered as the sister group to a clade encompassing Auchenoglanididae, Claroteidae, Malapteruridae and the African Schilbeidae, albeit with low statistical support. The two mochokid subfamilies Mochokinae and Chiloglanidinae as well as the chiloglanid tribe Atopochilini were recovered as reciprocally monophyletic. The genus Acanthocleithron forms the sister group of all remaining Mochokinae, although with low support. The genus Atopodontus is the sister group of all remaining Atopochilini. In contrast to morphological reconstructions, the monophyly of the genus Chiloglanis was strongly supported in our analysis, with Chiloglanis macropterus nested within a Chiloglanis sublineage encompassing only other taxa from the Congo drainage. This is an important result because the phylogenetic relationships of C. macropterus have been controversial in the past, and because we and other researchers assumed that this species would be resolved as sister to most or all other members of Chiloglanis. The apparent paraphyly of Synodontis with respect to Microsynodontis provided an additional surprise, with Synodontis punu turning out to be the sister group of the latter genus.

Supergene origin and maintenance in Atlantic cod

Supergenes are sets of genes that are inherited as a single marker and encode complex phenotypes through their joint action. They are identified in an increasing number of organisms, yet their origins and evolution remain enigmatic. In Atlantic cod, four megabase-scale supergenes have been identified and linked to migratory lifestyle and environmental adaptations. Here we investigate the origin and maintenance of these four supergenes through analysis of whole-genome-sequencing data, including a new long-read-based genome assembly for a non-migratory Atlantic cod individual. We corroborate the finding that chromosomal inversions underlie all four supergenes, and we show that they originated at different times between 0.40 and 1.66 million years ago. We reveal gene flux between supergene haplotypes where migratory and stationary Atlantic cod co-occur and conclude that this gene flux is driven by gene conversion, on the basis of an increase in GC content in exchanged sites. Additionally, we find evidence for double crossover between supergene haplotypes, leading to the exchange of an ~275 kilobase fragment with genes potentially involved in adaptation to low salinity in the Baltic Sea. Our results suggest that supergenes can be maintained over long timescales in the same way as hybridizing species, through the selective purging of introduced genetic variation.

Atlantic cod carries four supergenes linked to migratory lifestyle and environmental adaptations. Using whole-genome sequencing, the authors show that the genome inversions that underlie the supergenes originated at different times and show gene flux between supergene haplotypes.

Tracking the evolutionary pathways among Brazilian Lebiasina species (Teleostei: Lebiasinidae): a chromosomal and genomic comparative investigation

Abstract Despite several difficulties in chromosomal analyses of small-sized fishes, the cytogenetics of the Lebiasinidae was largely improved in the last years, showing differential patterns in the chromosomal evolution inside the family. In this context, it has been shown that genus Lebiasina preserves its karyotypic macrostructure, composed of 2n = 36 chromosomes, whereas the other genera generally present higher 2n. This study focused on the comparative cytogenetics of three Lebiasina species, one of them analyzed here for the first time, using conventional and molecular procedures. The results reinforced the differentiated evolutionary path of the genus Lebiasina while, at the same time, highlighted the genomic particularities that have accompanied the evolution of each species. In this sense, the repetitive components of the genome played a significant role in the differentiation of each species. It is also notable that L. minuta and L. melanoguttata, the two species that occur exclusively in the Brazilian territory, show greater chromosomal similarities to each other than to the trans-Andean sister species, L. bimaculata.

OUP accepted manuscript

The scale of data sets used to infer phylogenies has grown dramatically in the last decades, providing researchers with an enormous amount of information with which to draw inferences about evolutionary history. However, standard approaches to assessing confidence in those inferences (e.g., nonparametric bootstrap proportions [BP] and Bayesian posterior probabilities [PPs]) are still deeply influenced by statistical procedures and frameworks that were developed when information was much more limited. These approaches largely quantify uncertainty caused by limited amounts of data, which is often vanishingly small with modern, genome-scale sequence data sets. As a consequence, today’s phylogenomic studies routinely report near-complete confidence in their inferences, even when different studies reach strongly conflicting conclusions and the sites and loci in a single data set contain much more heterogeneity than our methods assume or can accommodate. Therefore, we argue that BPs and marginal PPs of bipartitions have outlived their utility as the primary means of measuring phylogenetic support for modern phylogenomic data sets with large numbers of sites relative to the number of taxa. Continuing to rely on these measures will hinder progress towards understanding remaining sources of uncertainty in the most challenging portions of the Tree of Life. Instead, we encourage researchers to examine the ideas and methods presented in this special issue of Systematic Biology and to explore the area further in their own work. The papers in this special issue outline strategies for assessing confidence and uncertainty in phylogenomic data sets that move beyond stochastic error due to limited data and offer promise for more productive dialogue about the challenges that we face in reaching our shared goal of understanding the history of life on Earth.[Big data; gene tree variation; genomic era; statistical bias.]

Molecular diversity and historical phylogeography of the widespread genus Mastiglanis (Siluriformes: Heptapteridae) based on palaeogeographical events in South America

The catfish family Heptapteridae, endemic in the Neotropical region, is composed of small- to medium-sized species. The genus Mastiglanis was monotypic until very recently and is often misidentified as Imparfinis because of similarities in diagnostic characters, including the length of maxillary barbels and limited pigmentation. We provide the first molecular inference of diversity for samples identified as Mastiglanis. Partial 16S and COI sequences were produced for 84 samples identified morphologically as Mastiglanis from the Amazon, Orinoco and Essequibo river basins. Species delimitation and phylogenetic methods recovered the genus as monophyletic, but samples assigned to the species Mastiglanis asopos yielded 21 distinct operational taxonomic units, often in sympatry. The first cladogenesis event, at ~12 Mya, was associated with marine incursions and/or vicariance events between the northward-flowing systems in the west and the rivers that drain to the east. Subsequent diversification during the Miocene and Pliocene was inferred to have resulted from environmental changes associated with Andean orogeny, whereas recent diversification in the late Pliocene and Pleistocene was inferred to result from environmental changes during glacial cycles. Despite the many probable cryptic species found across the large geographical distribution of the genus, further sampling is expected to increase taxonomic richness in this genus.

Whole-Genome Phylogenetic Reconstruction as a Powerful Tool to Reveal Homoplasy and Ancient Rapid Radiation in Waterflea Evolution

Although phylogeny estimation is notoriously difficult in radiations that occurred several hundred million years ago, phylogenomic approaches offer new ways to examine relationships among ancient lineages and evaluate hypotheses that are key to evolutionary biology. Here, we reconstruct the deep-rooted relationships of one of the oldest living arthropod clades, the branchiopod crustaceans, using a kaleidoscopic approach. We use concatenation and coalescent tree-building methods to analyze a large multigene data set at the nucleotide and amino acid level and examine gene tree versus species tree discordance. We unequivocally resolve long-debated relationships among extant orders of the Cladocera, the waterfleas, an ecologically relevant zooplankton group in global aquatic and marine ecosystems that is famous for its model systems in ecology and evolution. To build the data set, we assembled eight de novo genomes of key taxa including representatives of all extant cladoceran orders and suborders. Our phylogenetic analysis focused on a BUSCO-based set of 823 conserved single-copy orthologs shared among 23 representative taxa spanning all living branchiopod orders, including 11 cladoceran families. Our analysis supports the monophyly of the Cladocera and reveals remarkable homoplasy in their body plans. We found large phylogenetic distances between lineages with similar ecological specializations, indicating independent evolution in major body plans, such as in the pelagic predatory orders Haplopoda and Onychopoda (the “Gymnomera”). In addition, we assessed rapid cladogenesis by estimating relative timings of divergence in major lineages using reliable fossil-calibrated priors on eight nodes in the branchiopod tree, suggesting a Paleozoic origin around 325 Ma for the cladoceran ancestor and an ancient rapid radiation around 252 Ma at the Perm/Triassic boundary. These findings raise new questions about the roles of homoplasy and rapid radiation in the diversification of the cladocerans and help examine trait evolution from a genomic perspective in a functionally well understood, ancient arthropod group. [Cladocera; Daphnia; evolution; homoplasy; molecular clock; phylogenomics; systematics; waterfleas.]

The Emerging Phylogenetic Perspective on the Evolution of Actinopterygian Fishes

The emergence of a new phylogeny of ray-finned fishes at the turn of the twenty-first century marked a paradigm shift in understanding the evolutionary history of half of living vertebrates. We review how the new ray-finned fish phylogeny radically departs from classical expectations based on morphology. We focus on evolutionary relationships that span the backbone of ray-finned fish phylogeny, from the earliest divergences among teleosts and nonteleosts to the resolution of major lineages of Percomorpha. Throughout, we feature advances gained by the new phylogeny toward a broader understanding of ray-finned fish evolutionary history and the implications for topics that span from the genetics of human health to reconsidering the concept of living fossils. Additionally, we discuss conceptual challenges that involve reconciling taxonomic classification with phylogenetic relationships and propose an alternate higher-level classification for Percomorpha. Our review highlights remaining areas of phylogenetic uncertainty and opportunities for comparative investigations empowered by this new phylogenetic perspective on ray-finned fishes.

Tracking the Evolutionary Trends Among Small-Size Fishes of the Genus Pyrrhulina (Characiforme, Lebiasinidae): New Insights From a Molecular Cytogenetic Perspective

Miniature fishes have always been a challenge for cytogenetic studies due to the difficulty in obtaining chromosomal preparations, making them virtually unexplored. An example of this scenario relies on members of the family Lebiasinidae which include miniature to medium-sized, poorly known species, until very recently. The present study is part of undergoing major cytogenetic advances seeking to elucidate the evolutionary history of lebiasinids. Aiming to examine the karyotype diversification more deeply in Pyrrhulina, here we combined classical and molecular cytogenetic analyses, including Giemsa staining, C-banding, repetitive DNA mapping, comparative genomic hybridization (CGH), and whole chromosome painting (WCP) to perform the first analyses in five Pyrrhulina species (Pyrrhulina aff. marilynae, Pyrrhulina sp., P. obermulleri, P. marilynae and Pyrrhulina cf. laeta). The diploid number (2n) ranged from 40 to 42 chromosomes among all analyzed species, but P. marilynae is strikingly differentiated by having 2n = 32 chromosomes and a karyotype composed of large meta/submetacentric chromosomes, whose plesiomorphic status is discussed. The distribution of microsatellites does not markedly differ among species, but the number and position of the rDNA sites underwent significant changes among them. Interspecific comparative genome hybridization (CGH) found a moderate divergence in the repetitive DNA content among the species’ genomes. Noteworthy, the WCP reinforced our previous hypothesis on the origin of the X₁X₂Y multiple sex chromosome system in P. semifasciata. In summary, our data suggest that the karyotype differentiation in Pyrrhulina has been driven by major structural rearrangements, accompanied by high dynamics of repetitive DNAs.

Interrogating Phylogenetic Discordance Resolves Deep Splits in the Rapid Radiation of Old World Fruit Bats (Chiroptera: Pteropodidae)

The family Pteropodidae (Old World fruit bats) comprises $>$200 species distributed across the Old World tropics and subtropics. Most pteropodids feed on fruit, suggesting an early origin of frugivory, although several lineages have shifted to nectar-based diets. Pteropodids are of exceptional conservation concern with $>$50% of species considered threatened, yet the systematics of this group has long been debated, with uncertainty surrounding early splits attributed to an ancient rapid diversification. Resolving the relationships among the main pteropodid lineages is essential if we are to fully understand their evolutionary distinctiveness, and the extent to which these bats have transitioned to nectar-feeding. Here we generated orthologous sequences for $>$1400 nuclear protein-coding genes (2.8 million base pairs) across 114 species from 43 genera of Old World fruit bats (57% and 96% of extant species- and genus-level diversity, respectively), and combined phylogenomic inference with filtering by information content to resolve systematic relationships among the major lineages. Concatenation and coalescent-based methods recovered three distinct backbone topologies that were not able to be reconciled by filtering via phylogenetic information content. Concordance analysis and gene genealogy interrogation show that one topology is consistently the best supported, and that observed phylogenetic conflicts arise from both gene tree error and deep incomplete lineage sorting. In addition to resolving long-standing inconsistencies in the reported relationships among major lineages, we show that Old World fruit bats have likely undergone at least seven independent dietary transitions from frugivory to nectarivory. Finally, we use this phylogeny to identify and describe one new genus. [Chiroptera; coalescence; concordance; incomplete lineage sorting; nectar feeder; species tree; target enrichment.].

Comparative cytogenetics of Serrasalmidae (Teleostei: Characiformes): The relationship between chromosomal evolution and molecular phylogenies

Serrasalmidae has high morphological and chromosomal diversity. Based on molecular hypotheses, the family is currently divided into two subfamilies, Colossomatinae and Serrasalminae, with Serrasalminae composed of two tribes: Myleini (comprising most of pacus species) and Serrasalmini (represented by Metynnis, Catoprion, and remaining piranha’s genera). This study aimed to analyze species of the tribes Myleini (Myloplus asterias, M. lobatus, M. rubripinnis, M. schomburgki, and Tometes camunani) and Serrasalmini (Metynnis cuiaba, M. hypsauchen, and M. longipinnis) using classical and molecular cytogenetic techniques in order to understand the chromosomal evolution of the family. The four species of the genus Myloplus and T. camunani presented 2n = 58 chromosomes, while the species of Metynnis presented 2n = 62 chromosomes. The distribution of heterochromatin occurred predominantly in pericentromeric regions in all species. Tometes camunani and Myloplus spp. presented only one site with 5S rDNA. Multiple markers of 18S rDNA were observed in T. camunani, M. asterias, M. lobatus, M. rubripinnis, and M. schomburgkii. For Metynnis, however, synteny of the 18S and 5S rDNA was observed in the three species, in addition to an additional 5S marker in M. longipinnis. These data, when superimposed on the phylogeny of the family, suggest a tendency to increase the diploid chromosome number from 54 to 62 chromosomes, which occurred in a nonlinear manner and is the result of several chromosomal rearrangements. In addition, the different karyotype formulas and locations of ribosomal sequences can be used as cytotaxonomic markers and assist in the identification of species.

Ultra-Conserved Elements and morphology reciprocally illuminate conflicting phylogenetic hypotheses in Chalcididae (Hymenoptera, Chalcidoidea)

Recent technical advances combined with novel computational approaches have promised the acceleration of our understanding of the tree of life. However, when it comes to hyperdiverse and poorly known groups of invertebrates, studies are still scarce. As published phylogenies will be rarely challenged by future taxonomists, careful attention must be paid to potential analytical bias. We present the first molecular phylogenetic hypothesis for the family Chalcididae, a group of parasitoid wasps, with a representative sampling (144 ingroups and seven outgroups) that covers all described subfamilies and tribes, and 82% of the known genera. Analyses of 538 Ultra-Conserved Elements (UCEs) with supermatrix (RAxML and IQTREE) and gene tree reconciliation approaches (ASTRAL, ASTRID) resulted in highly supported topologies in overall agreement with morphology but reveal conflicting topologies for some of the deepest nodes. To resolve these conflicts, we explored the phylogenetic tree space with clustering and gene genealogy interrogation methods, analyzed marker and taxon properties that could bias inferences and performed a thorough morphological analysis (130 characters encoded for 40 taxa representative of the diversity). This joint analysis reveals that UCEs enable attainment of resolution between ancestry and convergent/divergent evolution when morphology is not informative enough, but also shows that a systematic exploration of bias with different analytical methods and a careful analysis of morphological features is required to prevent publication of artifactual results. We highlight a GC content bias for maximum-likelihood approaches, an artifactual mid-point rooting of the ASTRAL tree and a deleterious effect of high percentage of missing data (>85% missing UCEs) on gene tree reconciliation methods. Based on the results we propose a new classification of the family into eight subfamilies and ten tribes that lay the foundation for future studies on the evolutionary history of Chalcididae.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Integrating Sequence Capture and Restriction Site-Associated DNA Sequencing to Resolve Recent Radiations of Pelagic Seabirds

The diversification of modern birds has been shaped by a number of radiations. Rapid diversification events make reconstructing the evolutionary relationships among taxa challenging due to the convoluted effects of incomplete lineage sorting (ILS) and introgression. Phylogenomic data sets have the potential to detect patterns of phylogenetic incongruence, and to address their causes. However, the footprints of ILS and introgression on sequence data can vary between different phylogenomic markers at different phylogenetic scales depending on factors such as their evolutionary rates or their selection pressures. We show that combining phylogenomic markers that evolve at different rates, such as paired-end double-digest restriction site-associated DNA (PE-ddRAD) and ultraconserved elements (UCEs), allows a comprehensive exploration of the causes of phylogenetic discordance associated with short internodes at different timescales. We used thousands of UCE and PE-ddRAD markers to produce the first well-resolved phylogeny of shearwaters, a group of medium-sized pelagic seabirds that are among the most phylogenetically controversial and endangered bird groups. We found that phylogenomic conflict was mainly derived from high levels of ILS due to rapid speciation events. We also documented a case of introgression, despite the high philopatry of shearwaters to their breeding sites, which typically limits gene flow. We integrated state-of-the-art concatenated and coalescent-based approaches to expand on previous comparisons of UCE and RAD-Seq data sets for phylogenetics, divergence time estimation, and inference of introgression, and we propose a strategy to optimize RAD-Seq data for phylogenetic analyses. Our results highlight the usefulness of combining phylogenomic markers evolving at different rates to understand the causes of phylogenetic discordance at different timescales. [Aves; incomplete lineage sorting; introgression; PE-ddRAD-Seq; phylogenomics; radiations; shearwaters; UCEs.].

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.

Uneven Missing Data Skew Phylogenomic Relationships within the Lories and Lorikeets

The resolution of the Tree of Life has accelerated with advances in DNA sequencing technology. To achieve dense taxon sampling, it is often necessary to obtain DNA from historical museum specimens to supplement modern genetic samples. However, DNA from historical material is generally degraded, which presents various challenges. In this study, we evaluated how the coverage at variant sites and missing data among historical and modern samples impacts phylogenomic inference. We explored these patterns in the brush-tongued parrots (lories and lorikeets) of Australasia by sampling ultraconserved elements in 105 taxa. Trees estimated with low coverage characters had several clades where relationships appeared to be influenced by whether the sample came from historical or modern specimens, which were not observed when more stringent filtering was applied. To assess if the topologies were affected by missing data, we performed an outlier analysis of sites and loci, and a data reduction approach where we excluded sites based on data completeness. Depending on the outlier test, 0.15% of total sites or 38% of loci were driving the topological differences among trees, and at these sites, historical samples had 10.9× more missing data than modern ones. In contrast, 70% data completeness was necessary to avoid spurious relationships. Predictive modeling found that outlier analysis scores were correlated with parsimony informative sites in the clades whose topologies changed the most by filtering. After accounting for biased loci and understanding the stability of relationships, we inferred a more robust phylogenetic hypothesis for lories and lorikeets.

Concordance-based approaches for the inference of relationships and molecular rates with phylogenomic datasets

Gene tree conflict is common and finding methods to analyze and alleviate the negative effects that conflict has on species tree analysis is a crucial part of phylogenomics. This study aims to expand the discussion of inferring species trees and molecular branch lengths when conflict is present. Conflict is typically examined in two ways: inferring its prevalence, and inferring the influence of the individual genes (how strongly one gene supports any given topology compared to an alternative topology). Here, we examine a procedure for incorporating both conflict and the influence of genes in order to infer evolutionary relationships. All supported relationships in the gene trees are analyzed and the likelihood of the genes constrained to these relationships is summed to provide a likelihood for the relationship. Consensus tree assembly is conducted based on the sum of likelihoods for a given relationship and choosing relationships based on the most likely relationship assuming it does not conflict with a relationship that has a higher likelihood score. If it is not possible for all most likely relationships to be combined into a single bifurcating tree then multiple trees are produced and a consensus tree with a polytomy is created. This procedure allows for more influential genes to have greater influence on an inferred relationship, does not assume conflict has arisen from any one source, and does not force the dataset to produce a single bifurcating tree. Using this approach on three empirical datasets, we examine and discuss the relationship between influence and prevalence of gene tree conflict. We find that in one of the datasets, assembling a bifurcating consensus tree solely composed of the most likely relationships is impossible. To account for conflict in molecular rate analysis we also introduce a concordance-based approach to the summary and estimation of branch lengths suitable for downstream comparative analyses. We demonstrate through simulation that even under high levels of stochastic conflict, the mean and median of the concordant rates recapitulate the true molecular rate better than using a supermatrix approach. Using a large phylogenomic dataset, we examine rate heterogeneity across concordant genes with a focus on the branch subtending crown angiosperms. Notably, we find highly variable rates of evolution along the branch subtending crown angiosperms. The approaches outlined here have several limitations, but they also represent some alternative methods for harnessing the complexity of phylogenomic datasets and enrich our inferences of both species' relationships and evolutionary processes.

Accelerated Diversification Explains the Exceptional Species Richness of Tropical Characoid Fishes

The Neotropics harbor the most species-rich freshwater fish fauna on the planet, but the timing of that exceptional diversification remains unclear. Did the Neotropics accumulate species steadily throughout their long history, or attain their remarkable diversity recently? Biologists have long debated the relative support for these museum and cradle hypotheses, but few phylogenies of megadiverse tropical clades have included sufficient taxa to distinguish between them. We used 1,288 ultraconserved element loci (UCE) spanning 293 species, 211 genera and 21 families of characoid fishes to reconstruct a new, fossil-calibrated phylogeny and infer the most likely diversification scenario for a clade that includes a third of Neotropical fish diversity. This phylogeny implies paraphyly of the traditional delimitation of Characiformes because it resolves the largely Neotropical Characoidei as the sister lineage of Siluriformes (catfishes), rather than the African Citharinodei. Time-calibrated phylogenies indicate an ancient origin of major characoid lineages and reveal a much more recent emergence of most characoid species. Diversification rate analyses infer increased speciation and decreased extinction rates during the Oligocene at around 30 million years ago (Ma) during a period of mega-wetland formation in the proto-Orinoco-Amazonas. Three species-rich and ecomorphologically diverse lineages (Anostomidae, Serrasalmidae, and Characidae) that originated more than 60 Ma in the Paleocene experienced particularly notable bursts of Oligocene diversification and now account collectively for 68% of the approximately 2,150 species of Characoidei. In addition to paleogeographic changes, we discuss potential accelerants of diversification in these three lineages. While the Neotropics accumulated a museum of ecomorphologically diverse characoid lineages long ago, this geologically dynamic region also cradled a much more recent birth of remarkable species-level diversity.

Development and phylogenetic utilities of a new set of single-/low-copy nuclear genes in Senecioneae (Asteraceae), with new insights into the tribal position and the relationships within subtribe Tussilagininae

The tribe Senecioneae is one of the largest tribes in Asteraceae, with a nearly cosmopolitan distribution. Despite great efforts devoted to elucidate the evolution of Senecioneae, many questions still remain concerning the systematics of this group, from the tribal circumscription and position to species relationships in many genera. The hybridization-based target enrichment method of next-generation sequencing has been accepted as a promising approach to resolve phylogenetic problems. We herein develop a set of single-/low-copy genes for Senecioneae, and test their phylogenetic utilities. Our results demonstrate that these genes work highly efficiently for Senecioneae, with a high average gene recovery of 98.8% across the tribe and recovering robust phylogenetic hypotheses at different levels. In particular, the delimitation of the Senecioneae has been confirmed to include Abrotanella and exclude Doronicum, with the former sister to core Senecioneae and the latter shown to be more closely related to Calenduleae. Moreover, Doronicum and Calenduleae are inferred to be the closest relatives of Senecioneae, which is a new hypothesis well supported by statistical topology tests, morphological evidence, and the profile of pyrrolizidine alkaloids, a special kind of chemical characters generally used to define Senecioneae. Furthermore, this study suggests a complex reticulation history in the diversification of Senecioneae, accounting for the prevalence of polyploid groups in the tribe. With subtribe Tussilagininae s.str. as a case study showing a more evident pattern of gene duplication, we further explored reconstructing the phylogeny in the groups with high ploidy levels. Our results also demonstrate that tree topologies based on sorted paralogous copies are stable across different methods of phylogenetic inference, and more congruent with the morphological evidence and the results of previous phylogenetic studies.

Multi-locus phylogeny with dense Guiana Shield sampling supports new suprageneric classification of the neotropical three-barbeled catfishes (Siluriformes: Heptapteridae)

The catfish family Heptapteridae is ubiquitous across a range of freshwater habitats from southern Mexico to northern Argentina and contains 23 genera and 228 valid species. After a century of mostly morphology-based systematic analyses of these fishes, we provide the first molecular phylogenetic hypothesis spanning most valid Heptapteridae genera (16 of 23). We examined eight of 14 valid genera in the Nemuroglanis-subclade (Heptapterini), all valid genera in the Brachyglanis-subclade (Brachyglaniini) and most valid Brachyglaniini species (11 of 15). Maximum likelihood and Bayesian analyses of a 4156-base alignment of five gene regions (three mitochondrial: COI, Cyt b, and ND2; two nuclear: RAG2, Glyt) yielded thoroughly resolved and statistically robust phylogenies that were largely congruent with each other and with previous morphology-based hypotheses. We propose a revised phylogenetic classification consisting of two subfamilies (Rhamdiinae, Heptapterinae) each with two tribes. Dense taxonomic sampling of Brachyglaniini, including type species of Brachyglanis, Gladioglanis, Leptorhamdia, and Myoglanis, revealed widespread paraphyly. Newly recovered clades within Brachyglaniini are closely associated with either the upper Orinoco or the Essequibo suggesting repeated dispersals and/or range expansions/contractions across the western Guiana Shield highlands and from there to the upper Amazon and Brazilian Shield. These biogeographical processes appear to have been an important driver of allopatric diversification in the clade.