Considering decoupled phenotypic diversification between ontogenetic phases in macroevolution: An example using Triggerfishes (Balistidae).

Across the Tree of Life, most studies of phenotypic disparity and diversification have been restricted to adult organisms. However, many lineages have distinct ontogenetic phases that differ from their adult forms in morphology and ecology. Focusing disproportionately on the evolution of adult forms unnecessarily hinders our understanding of the pressures shaping evolution over time. Non-adult disparity patterns are particularly important to consider for coastal ray-finned fishes, which often have juvenile phases with distinct phenotypes. These juvenile forms are often associated with sheltered nursery environments, with phenotypic shifts between adults and juvenile stages that are readily apparent in locomotor morphology. Whether this ontogenetic variation in locomotor morphology reflects a decoupling of diversification dynamics between life stages remains unknown. Here we investigate the evolutionary dynamics of locomotor morphology between adult and juvenile triggerfishes. We integrate a time-calibrated phylogenetic framework with geometric morphometric approaches and measurement data of fin aspect ratio and incidence, and reveal a mismatch between morphospace occupancy, the evolution of morphological disparity, and the tempo of trait evolution between life stages. Collectively, our results illuminate how the heterogeneity of morpho-functional adaptations can decouple the mode and tempo of morphological diversification between ontogenetic stages.

An inverse latitudinal gradient in speciation rate for marine fishes.

Far more species of organisms are found in the tropics than in temperate and polar regions, but the evolutionary and ecological causes of this pattern remain controversial1,2. Tropical marine fish communities are much more diverse than cold-water fish communities found at higher latitudes3,4, and several explanations for this latitudinal diversity gradient propose that warm reef environments serve as evolutionary 'hotspots' for species formation5-8. Here we test the relationship between latitude, species richness and speciation rate across marine fishes. We assembled a time-calibrated phylogeny of all ray-finned fishes (31,526 tips, of which 11,638 had genetic data) and used this framework to describe the spatial dynamics of speciation in the marine realm. We show that the fastest rates of speciation occur in species-poor regions outside the tropics, and that high-latitude fish lineages form new species at much faster rates than their tropical counterparts. High rates of speciation occur in geographical regions that are characterized by low surface temperatures and high endemism. Our results reject a broad class of mechanisms under which the tropics serve as an evolutionary cradle for marine fish diversity and raise new questions about why the coldest oceans on Earth are present-day hotspots of species formation.

A highly contiguous genome assembly for the Yellow Warbler (Setophaga petechia).

The Yellow Warbler (Setophaga petechia) is a small songbird in the wood-warbler family (Parulidae) that exhibits phenotypic and ecological differences across a widespread distribution and is important to California's riparian habitat conservation. Here, we present a high-quality de novo genome assembly of a vouchered female Yellow Warbler from southern California. Using HiFi long-read and Omni-C proximity sequencing technologies, we generated a 1.22 Gb assembly including 687 scaffolds with a contig N50 of 6.80 Mb, scaffold N50 of 21.18 Mb, and a BUSCO completeness score of 96.0%. This highly contiguous genome assembly provides an essential resource for understanding the history of gene flow, divergence, and local adaptation in Yellow Warblers and can inform conservation management of this charismatic bird species.

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 1288 ultraconserved element loci 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 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 2150 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. [Biodiversity; Characiformes; macroevolution; Neotropics; phylogenomics; ultraconserved elements.].

Using ultraconserved elements to track the influence of sea-level change on leafy seadragon populations.

During the Last Glacial Maximum (LGM), global sea levels were 120-130 m lower than today, resulting in the emergence of most continental shelves and extirpation of subtidal organisms from these areas. During the interglacial periods, rapid inundation of shelf regions created a dynamic environment for coastal organisms, such as the charismatic leafy seadragon (Phycodurus eques, Syngnathidae), a brooder with low dispersal ability inhabiting kelp beds in temperate Australia. Reconstructions of the palaeoshoreline revealed that the increase of shallow areas since the LGM was not uniform across the species' range and we investigated the effects of these asymmetries on genetic diversity and structuring. Using targeted capture of 857 variable ultraconserved elements (UCEs, 2,845 single nucleotide polymorphisms) in 68 individuals, we found that the regionally different shelf topographies were paralleled by contrasting population genetic patterns. In the west, populations may not have persisted through sea-level lows because shallow seabed was very limited. Shallow genetic structure, weak expansion signals and a westward cline in genetic diversity indicate a postglacial recolonization of the western part of the range from a more eastern location following sea-level rise. In the east, shallow seabed persisted during the LGM and increased considerably after the flooding of large bays, which resulted in strong demographic expansions, deeper genetic structure and higher genetic diversity. This study suggests that postglacial flooding with rising sea levels produced locally variable signatures in colonizing populations.

Estimating Diversification Rates on Incompletely Sampled Phylogenies: Theoretical Concerns and Practical Solutions.

Molecular phylogenies are a key source of information about the tempo and mode of species diversification. However, most empirical phylogenies do not contain representatives of all species, such that diversification rates are typically estimated from incompletely sampled data. Most researchers recognize that incomplete sampling can lead to biased rate estimates, but the statistical properties of methods for accommodating incomplete sampling remain poorly known. In this point of view, we demonstrate theoretical concerns with the widespread use of analytical sampling corrections for sparsely sampled phylogenies of higher taxonomic groups. In particular, corrections based on "sampling fractions" can lead to low statistical power to infer rate variation when it is present, depending on the likelihood function used for inference. In the extreme, the sampling fraction correction can lead to spurious patterns of diversification that are driven solely by unbalanced sampling across the tree in concert with low overall power to infer shifts. Stochastic polytomy resolution provides an alternative to sampling fraction approaches that avoids some of these biases. We show that stochastic polytomy resolvers can greatly improve the power of common analyses to estimate shifts in diversification rates. We introduce a new stochastic polytomy resolution method (Taxonomic Addition for Complete Trees [TACT]) that uses birth-death-sampling estimators across an ultrametric phylogeny to estimate branching times for unsampled taxa, with taxonomic information to compatibly place new taxa onto a backbone phylogeny. We close with practical recommendations for diversification inference under several common scenarios of incomplete sampling. [Birth-death process; diversification; incomplete sampling; phylogenetic uncertainty; rate heterogeneity; rate shifts; stochastic polytomy resolution.].

Supermatrix phylogeny resolves goby lineages and reveals unstable root of Gobiaria.

Gobies, sleepers, and cardinalfishes represent major clades of a species rich radiation of small bodied, ecologically diverse percomorphs (Gobiaria). Molecular phylogenetics has been crucial to resolving broad relationships of sleepers and gobies (Gobioidei), but the phylogenetic placements of cardinalfishes and nurseryfishes, as reciprocal or sequential sister clades to Gobioidei, are uncertain. In order to evaluate relationships among and within families we used a phylogenetic data mining approach to generate densely sampled trees inclusive of all higher taxa. We utilized conspecific amino acid homology to improve alignment accuracy, included ambiguously identified taxa to increase taxon sampling density, and resampled individual gene alignments to filter rogue sequences before concatenation. This approach yielded the most comprehensive tree yet of Gobiaria, inferred from a sparse (17 percent-complete) supermatrix of one ribosomal and 22 protein coding loci (18,065 characters), comprised of 50 outgroup and 777 ingroup taxa, representing 32 percent of species and 68 percent of genera. Our analyses confirmed the lineage-based classification of gobies with strong support, identified sleeper clades with unforeseen levels of systematic uncertainty, and quantified competing phylogenetic signals that confound resolution of the root topology. We also discovered that multilocus data completeness was related to maximum likelihood branch support, and verified that the phylogenetic uncertainty of shallow relationships observed within goby lineages could largely be explained by supermatrix sparseness. These results demonstrate the potential and limits of publicly available sequence data for producing densely-sampled phylogenetic trees of exceptionally biodiverse groups.

How Well Can We Estimate Diversity Dynamics for Clades in Diversity Decline?

The fossil record shows that the vast majority of all species that ever existed are extinct and that most lineages go through an expansion and decline in diversity. However, macroevolutionary analyses based upon molecular phylogenies have difficulty inferring extinction dynamics, raising questions about whether the neontological record can contribute to an understanding of the decline phenomenon. Two recently developed diversification methods for molecular phylogenies (RPANDA and BAMM) incorporate models that theoretically have the capacity to capture decline dynamics by allowing extinction to be higher than speciation. However, the performance of these frameworks over a wide range of decline scenarios has not been studied. Here, we investigate the behavior of these methods under decline scenarios caused by decreasing speciation and increasing extinction through time on simulated trees at fixed intervals over diversity trajectories with expansion and decline phases. We also compared method performance over a comprehensive data set of 214 empirical trees. Our results show that both methods perform equally well when varying speciation rates control decline. When decline was only caused by an increase in extinction rates both methods wrongly assign the variation in net diversification to a drop in speciation, even though the positive gamma values of those trees would suggest otherwise. We also found a tendency for RPANDA to favor increasing extinction and BAMM to favor decreasing speciation as the most common cause of decline in empirical trees. Overall our results shed light on the limitations of both methods, encouraging researchers to carefully interpret the results from diversification studies.

Lack of Signal for the Impact of Conotoxin Gene Diversity on Speciation Rates in Cone Snails.

Understanding why some groups of organisms are more diverse than others is a central goal in macroevolution. Evolvability, or the intrinsic capacity of lineages for evolutionary change, is thought to influence disparities in species diversity across taxa. Over macroevolutionary time scales, clades that exhibit high evolvability are expected to have higher speciation rates. Cone snails (family: Conidae, $>$900 spp.) provide a unique opportunity to test this prediction because their toxin genes can be used to characterize differences in evolvability between clades. Cone snails are carnivorous, use prey-specific venom (conotoxins) to capture prey, and the genes that encode venom are known and diversify through gene duplication. Theory predicts that higher gene diversity confers a greater potential to generate novel phenotypes for specialization and adaptation. Therefore, if conotoxin gene diversity gives rise to varying levels of evolvability, conotoxin gene diversity should be coupled with macroevolutionary speciation rates. We applied exon capture techniques to recover phylogenetic markers and conotoxin loci across 314 species, the largest venom discovery effort in a single study. We paired a reconstructed timetree using 12 fossil calibrations with species-specific estimates of conotoxin gene diversity and used trait-dependent diversification methods to test the impact of evolvability on diversification patterns. Surprisingly, we did not detect any signal for the relationship between conotoxin gene diversity and speciation rates, suggesting that venom evolution may not be the rate-limiting factor controlling diversification dynamics in Conidae. Comparative analyses showed some signal for the impact of diet and larval dispersal strategy on diversification patterns, though detection of a signal depended on the dataset and the method. If our results remain true with increased taxonomic sampling in future studies, they suggest that the rapid evolution of conid venom may cause other factors to become more critical to diversification, such as ecological opportunity or traits that promote isolation among lineages.

Resolving Deep Nodes in an Ancient Radiation of Neotropical Fishes in the Presence of Conflicting Signals from Incomplete Lineage Sorting.

Resolving patterns of ancient and rapid diversifications is one of the most challenging tasks in evolutionary biology. These difficulties arise from confusing phylogenetic signals that are associated with the interplay of incomplete lineage sorting (ILS) and homoplasy. Phylogenomic analyses of hundreds, or even thousands, of loci offer the potential to resolve such contentious relationships. Yet, how much useful phylogenetic information these large data sets contain remains uncertain and often goes untested. Here, we assess the utility of different data filtering approaches to maximize phylogenetic information and minimize noise when reconstructing an ancient radiation of Neotropical electric knifefishes (Order Gymnotiformes) using ultraconserved elements. We found two contrasting hypotheses of gymnotiform evolutionary relationships depending on whether phylogenetic inferences were based on concatenation or coalescent methods. In the first case, all analyses inferred a previously-and commonly-proposed hypothesis, where the family Apteronotidae was found as the sister group to all other gymnotiform families. In contrast, coalescent-based analyses suggested a novel hypothesis where families producing pulse-type (viz., Gymnotidae, Hypopomidae, and Rhamphichthyidae) and wave-type electric signals (viz., Apteronotidae, Sternopygidae) were reciprocally monophyletic. Nodal support for this second hypothesis increased when analyzing loci with the highest phylogenetic information content and further increased when data were pruned using targeted filtering methods that maximized phylogenetic informativeness at the deepest nodes of the Gymnotiformes. Bayesian concordance analyses and topology tests of individual gene genealogies demonstrated that the difficulty of resolving this radiation was likely due to high gene-tree incongruences that resulted from ILS. We show that data filtering reduces gene-tree heterogeneity and increases nodal support and consistency of species trees using coalescent methods; however, we failed to observe the same effect when using concatenation methods. Furthermore, the targeted filtering strategies applied here support the use of "gene data interrogation" rather than "gene genealogy interrogation" approaches in phylogenomic analyses, to extract phylogenetic signal from intractable portions of the Tree of Life.

Pleiotropic jaw morphology links the evolution of mechanical modularity and functional feeding convergence in Lake Malawi cichlids.

Complexity in how mechanistic variation translates into ecological novelty could be critical to organismal diversification. For instance, when multiple distinct morphologies can generate the same mechanical or functional phenotype, this could mitigate trade-offs and/or provide alternative ways to meet the same ecological challenge. To investigate how this type of complexity shapes diversity in a classic adaptive radiation, we tested several evolutionary consequences of the anterior jaw four-bar linkage for Lake Malawi cichlid trophic diversification. Using a novel phylogenetic framework, we demonstrated that different mechanical outputs of the same four jaw elements are evolutionarily associated with both jaw protrusion distance and jaw protrusion angle. However, these two functional aspects of jaw protrusion have evolved independently. Additionally, although four-bar morphology showed little evidence for attraction to optima, there was substantial evidence of adaptive peaks for emergent four-bar linkage mechanics and jaw protrusion abilities among Malawi feeding guilds. Finally, we highlighted a clear case of two cichlid species that have -independently evolved to graze algae in less than 2 Myr and have converged on similar jaw protrusion abilities as well as four-bar linkage mechanics, but have evolved these similarities via non-convergent four-bar morphologies.

A phylogenomic framework for pelagiarian fishes (Acanthomorpha: Percomorpha) highlights mosaic radiation in the open ocean.

The fish clade Pelagiaria, which includes tunas as its most famous members, evolved remarkable morphological and ecological variety in a setting not generally considered conducive to diversification: the open ocean. Relationships within Pelagiaria have proven elusive due to short internodes subtending major lineages suggestive of rapid early divergences. Using a novel sequence dataset of over 1000 ultraconserved DNA elements (UCEs) for 94 of the 286 species of Pelagiaria (more than 70% of genera), we provide a time-calibrated phylogeny for this widely distributed clade. Some inferred relationships have clear precedents (e.g. the monophyly of 'core' Stromateoidei, and a clade comprising 'Gempylidae' and Trichiuridae), but others are unexpected despite strong support (e.g. Chiasmodontidae + Tetragonurus). Relaxed molecular clock analysis using node-based fossil calibrations estimates a latest Cretaceous origin for Pelagiaria, with crown-group families restricted to the Cenozoic. Estimated mean speciation rates decline from the origin of the group in the latest Cretaceous, although credible intervals for root and tip rates are broad and overlap in most cases, and there is higher-than-expected partitioning of body shape diversity (measured as fineness ratio) between clades concentrated during the Palaeocene-Eocene. By contrast, more direct measures of ecology show either no substantial deviation from a null model of diversification (diet) or patterns consistent with evolutionary constraint or high rates of recent change (depth habitat). Collectively, these results indicate a mosaic model of diversification. Pelagiarians show high morphological disparity and modest species richness compared to better-studied fish radiations in contrasting environments. However, this pattern is also apparent in other clades in open-ocean or deep-sea habitats, and suggests that comparative study of such groups might provide a more inclusive model of the evolution of diversity in fishes.

Phylogenomic reappraisal of the Neotropical catfish family Loricariidae (Teleostei: Siluriformes) using ultraconserved elements.

Neotropical freshwaters host more than 6000 fish species, of which 983 are suckermouth armored catfishes of the family Loricariidae - the most-diverse catfish family and fifth most species-rich vertebrate family on Earth. Given their diversity and ubiquitous distribution across many habitat types, loricariids are an excellent system in which to investigate factors that create and maintain Neotropical fish diversity, yet robust phylogenies needed to support such ecological and evolutionary studies are lacking. We sought to buttress the systematic understanding of loricariid catfishes by generating a genome-scale data set (1041 loci, 328,330 bp) for 140 species spanning 75 genera and five of six previously proposed subfamilies. Both maximum likelihood and Bayesian analyses strongly supported the monophyly of Loricariidae. Our results also reinforced the established backbone of loricariid interrelationships: Delturinae as sister to all other analyzed loricariids, with subfamily Rhinelepinae diverging next, followed by Loricariinae sister to Hypostominae + Hypoptopomatinae. Previous DNA-based relationships within Hypostominae and Loricariinae were strongly supported. However, we evaluated for the first time DNA-based relationships among many Hypoptopomatinae genera and found significant differences with this subfamily's current genus-level classification, prompting several taxonomic changes. Finally, we placed our topological results within a fossil-calibrated temporal context indicating that early Loricariidae diversification occurred across the Cretaceous-Paleogene boundary ∼65 million years ago (Ma). Our study lays a strong foundation for future research to focus on relationships among species and the macroevolutionary processes affecting loricariid diversification rates and patterns.

Phylogenetic Factor Analysis.

Phylogenetic comparative methods explore the relationships between quantitative traits adjusting for shared evolutionary history. This adjustment often occurs through a Brownian diffusion process along the branches of the phylogeny that generates model residuals or the traits themselves. For high-dimensional traits, inferring all pair-wise correlations within the multivariate diffusion is limiting. To circumvent this problem, we propose phylogenetic factor analysis (PFA) that assumes a small unknown number of independent evolutionary factors arise along the phylogeny and these factors generate clusters of dependent traits. Set in a Bayesian framework, PFA provides measures of uncertainty on the factor number and groupings, combines both continuous and discrete traits, integrates over missing measurements and incorporates phylogenetic uncertainty with the help of molecular sequences. We develop Gibbs samplers based on dynamic programming to estimate the PFA posterior distribution, over 3-fold faster than for multivariate diffusion and a further order-of-magnitude more efficiently in the presence of latent traits. We further propose a novel marginal likelihood estimator for previously impractical models with discrete data and find that PFA also provides a better fit than multivariate diffusion in evolutionary questions in columbine flower development, placental reproduction transitions and triggerfish fin morphometry.

Phylogenomics of a putatively convergent novelty: did hypertrophied lips evolve once or repeatedly in Lake Malawi cichlid fishes?

BACKGROUND: Phylogenies provide critical information about convergence during adaptive radiation. To test whether there have been multiple origins of a distinctive trophic phenotype in one of the most rapidly radiating groups known, we used ultra-conserved elements (UCEs) to examine the evolutionary affinities of Lake Malawi cichlids lineages exhibiting greatly hypertrophied lips. RESULTS: The hypertrophied lip cichlids Cheilochromis euchilus, Eclectochromis ornatus, Placidochromis "Mbenji fatlip", and Placidochromis milomo are all nested within the non-mbuna clade of Malawi cichlids based on both concatenated sequence and single nucleotide polymorphism (SNP) inferred phylogenies. Lichnochromis acuticeps that exhibits slightly hypertrophied lips also appears to have evolutionary affinities to this group. However, Chilotilapia rhoadesii that lacks hypertrophied lips was recovered as nested within the species Cheilochromis euchilus. Species tree reconstructions and analyses of introgression provided largely ambiguous patterns of Malawi cichlid evolution. CONCLUSIONS: Contrary to mitochondrial DNA phylogenies, bifurcating trees based on our 1024 UCE loci supported close affinities of Lake Malawi lineages with hypertrophied lips. However, incomplete lineage sorting in Malawi tends to render these inferences more tenuous. Phylogenomic analyses will continue to provide powerful inferences about whether phenotypic novelties arose once or multiple times during adaptive radiation.

Filtering nucleotide sites by phylogenetic signal to noise ratio increases confidence in the Neoaves phylogeny generated from ultraconserved elements.

Despite genome scale analyses, high-level relationships among Neoaves birds remain contentious. The placements of the Neoaves superorders are notoriously difficult to resolve because they involve deep splits followed by short internodes. Using our approach, we investigate whether filtering UCE loci on their phylogenetic signal to noise ratio helps to resolve key nodes in the Neoaves tree of life. We find that our analysis of data sets filtered for high signal to noise ratio results in topologies that are inconsistent with unfiltered results but that are congruent with whole-genome analyses. These relationships include the Columbea + Passerea sister relationship and the Phaethontimorphae + Aequornithia sister relationship. We also find increased statistical support for more recent nodes (i.e. the Pelecanidae + Ardeidae sister relationship, the Eucavitaves clade, and the Otidiformes + Musophagiformes sister relationship). We also find instances where support is reduced for well-established clades, possibly due to the removal of sites with moderate signal-to-noise ratio. Our results suggest that filtering on the basis of signal to noise ratio is a useful tool for resolving problematic splits in phylogenomic data sets.

Phylogenetics and geography of speciation in New World Halichoeres wrasses.

The New World Halichoeres comprises about 30 small to medium sized wrasse species that are prominent members of reef communities throughout the tropical Western Atlantic and Eastern Pacific. We conducted a phylogenetic analysis of this group and related lineages using new and previously published sequence data. We estimated divergence times, evaluated the monophyly of this group, their relationship to other labrids, as well as the time-course and geography of speciation. These analyses show that all members of New World Halichoeres form a monophyletic group that includes Oxyjulis and Sagittalarva. New World Halichoeres is one of numerous labrid groups that appear to have radiated rapidly about 32 Ma and form a large polytomy within the julidine wrasses. We reconstruct the tropical Western Atlantic to be the ancestral area of New World Halichoeres, with four invasions of the Eastern Pacific and one reversal from East Pacific to Western Atlantic. These five speciation events were spread across the history of the group, with none corresponding closely to the time of the closure of the Isthmus of Panama. Three speciation events within the Atlantic occurred across the Orinoco-Amazon outflow and within the Pacific, five involve splits between lineages that occupy coastal reef systems and offshore islands. Of eight sister species pairs, seven show complete allopatry and one is fully sympatric.

A phylogenomic perspective on the robust capuchin monkey (Sapajus) radiation: First evidence for extensive population admixture across South America.

Phylogenetic relationships amongst the robust capuchin monkeys (genus Sapajus) are poorly understood. Morphology-based taxonomies have recognized anywhere from one to twelve different species. The current IUCN (2017) classification lists eight robust capuchins: S. xanthosternos, S. nigritus, S. robustus, S. flavius, S. libidinosus, S. cay, S. apella and S. macrocephalus. Here, we assembled the first phylogenomic data set for Sapajus using ultra-conserved elements (UCEs) to reconstruct a capuchin phylogeny. All phylogenomic analyses strongly supported a deep divergence of Sapajus and Cebus clades within the capuchin monkeys, and provided support for Sapajus nigritus, S. robustus and S. xanthosternos as distinct species. However, the UCE phylogeny lumped the putative species S. cay, S. libidinosus, S. apella, S. macrocephalus, and S. flavius together as a single widespread lineage. A SNP phylogeny constructed from the UCE data was better resolved and recovered S. flavius and S. libidinosus as sister species; however, S. apella, S. macrocephalus, and S. cay individuals were recovered in two geographic clades, from northeastern and southwestern Amazon, rather than clustering by currently defined morphospecies. STRUCTURE analysis of population clustering revealed widespread admixture among Sapajus populations within the Amazon and even into the Cerrado and Atlantic Forest. Difficulty in assigning species by morphology may be a result of widespread population admixture facilitated through frequent movement across major rivers and even ecosystems by robust capuchin monkeys.

Phylogenomic Systematics of Ostariophysan Fishes: Ultraconserved Elements Support the Surprising Non-Monophyly of Characiformes.

Ostariophysi is a superorder of bony fishes including more than 10,300 species in 1100 genera and 70 families. This superorder is traditionally divided into five major groups (orders): Gonorynchiformes (milkfishes and sandfishes), Cypriniformes (carps and minnows), Characiformes (tetras and their allies), Siluriformes (catfishes), and Gymnotiformes (electric knifefishes). Unambiguous resolution of the relationships among these lineages remains elusive, with previous molecular and morphological analyses failing to produce a consensus phylogeny. In this study, we use over 350 ultraconserved element (UCEs) loci comprising 5 million base pairs collected across 35 representative ostariophysan species to compile one of the most data-rich phylogenies of fishes to date. We use these data to infer higher level (interordinal) relationships among ostariophysan fishes, focusing on the monophyly of the Characiformes-one of the most contentiously debated groups in fish systematics. As with most previous molecular studies, we recover a non-monophyletic Characiformes with the two monophyletic suborders, Citharinoidei and Characoidei, more closely related to other ostariophysan clades than to each other. We also explore incongruence between results from different UCE data sets, issues of orthology, and the use of morphological characters in combination with our molecular data. [Conserved sequence; ichthyology; massively parallel sequencing; morphology; next-generation sequencing; UCEs.].

Phylogenomic analysis of Lake Malawi cichlid fishes: Further evidence that the three-stage model of diversification does not fit.

Adaptive radiations could often occur in discrete stages. For instance, the species flock of ∼1000 species of Lake Malawi cichlid fishes might have only diverged once between rocky and sandy environments during the initial stage of their diversification. All further diversification within the rock-dwelling (mbuna) or sand-dwelling (utaka) cichlids would have occurred during a subsequent second stage of extensive trophic evolution that was followed by a third stage of sexual trait divergence. We provide an improved phylogenetic framework for Malawi cichlids to test this three-stage hypothesis based on newly reconstructed phylogenetic relationships among 32 taxonomically disparate Malawi cichlids species. Using several reconstruction methods and 1037 ultra-conserved element (UCE) markers, we recovered a molecular phylogeny that confidently resolved relationships among most of the Malawi lineages sampled when a bifurcating framework was enforced. These bifurcating reconstructions also indicated that the sand-dwelling species Cyathochromis obliquidens was well-nested within the primarily rock-dwelling radiation known as the mbuna. In contrast to predictions from the three-stage model of vertebrate diversification, the recovered phylogeny reveals an initial colonization of rocky reefs, followed by substantial diversification of rock-dwelling lineages, and then at least one instance of subsequent evolution back into sandy habitats. This repeated evolution into major habitat types provides further evidence that the three-stage model of Malawi cichlid diversification has numerous exceptions.