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.

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.

Mechanisms of equality and inequality in mammalian societies

The extent of (in)equality is highly diverse across species of social mammals, but we have a poor understanding of the factors that produce or inhibit equitable social organizations. Here, we adopt a comparative evolutionary perspective to test whether the evolution of social dominance hierarchies, a measure of social inequality in animals, exhibits phylogenetic conservatism and whether interspecific variation in these traits can be explained by sex, age or captivity. We find that hierarchy steepness and directional consistency evolve rapidly without any apparent constraint from evolutionary history. Given this extraordinary variability, we next consider multiple factors that have evolved to mitigate social inequality. Social networks, coalitionary support and knowledge transfer advantage to privilege some individuals over others. Nutritional access and prenatal stressors can impact the development of offspring, generating health disparities with intergenerational consequences. Intergenerational transfer of material resources (e.g. stone tools, food stashes, territories) advantage those who receive. Nonetheless, many of the same social species that experience unequal access to food (survival) and mates (reproduction) engage in levelling mechanisms such as food sharing, adoption, revolutionary coalitions, forgiveness and inequity aversion. Taken together, mammals rely upon a suite of mechanisms of (in)equality to balance the costs and benefits of group living. This article is part of the theme issue 'Evolutionary ecology of inequality'.

Keels of boxfish carapaces strongly improve stabilization against roll

Boxfish (Ostraciidae) have peculiar body shapes, with conspicuous keels formed by their bony carapaces. Previous studies have proposed various hydrodynamic roles for these keels, including reducing drag during swimming, contributing to passive stabilization of the swimming course, or providing resistance against roll rotations. Here, we tested these hypotheses using computational fluid dynamics simulations of five species of Ostraciidae with a range of carapace shapes. The hydrodynamic performance of the original carapace surface models, obtained from laser scanning of museum specimens, was compared with models where the keels had been digitally reduced. The original carapaces showed no reduced drag or increased passive stability against pitch and yaw compared to the reduced-keel carapaces. However, consistently for all studied species, a strong increase in roll drag and roll-added mass was observed for the original carapaces compared to the reduced-keel carapaces, despite the relatively small differences in keel height. In particular, the damping of roll movement by resistive drag torques increased considerably by the presence of keels. Our results suggest that the shape of the boxfish carapace is important in enabling the observed roll-free forward swimming of boxfish and may facilitate the control of manoeuvres.

The nature of privilege: intergenerational wealth in animal societies

Wealth inequality is widespread across human societies, from pastoral and small-scale agricultural groups to large modern social structures. The intergenerational transfer of wealth privileges some individuals over others through the transmission of resources external to an individual organism. Privileged access to household wealth (e.g., land, shelter, silver) positively influences the destinies of some (and their descendants) over others in human societies. Strikingly parallel phenomena exist in animal societies. Inheritance of nongenetic commodities (e.g., a nest, territory, tool) external to an individual also contributes greatly to direct fitness in animals. Here, we illustrate the evolutionary diversity of privilege and its disparity-generating effects on the evolutionary trajectories of lineages across the Tree of Life. We propose that integration of approaches used to study these patterns in humans may offer new insights into a core principle from behavioral ecology—differential access to inherited resources—and help to establish a broad, comparative framework for studying inequality in animals.

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.

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.].

Addressing incomplete lineage sorting and paralogy in the inference of uncertain salmonid phylogenetic relationships

Recent and continued progress in the scale and sophistication of phylogenetic research has yielded substantial advances in knowledge of the tree of life; however, segments of that tree remain unresolved and continue to produce contradicting or unstable results. These poorly resolved relationships may be the product of methodological shortcomings or of an evolutionary history that did not generate the signal traits needed for its eventual reconstruction. Relationships within the euteleost fish family Salmonidae have proven challenging to resolve in molecular phylogenetics studies in part due to ancestral autopolyploidy contributing to conflicting gene trees. We examine a sequence capture dataset from salmonids and use alternative strategies to accommodate the effects of gene tree conflict based on aspects of salmonid genome history and the multispecies coalescent. We investigate in detail three uncertain relationships: (1) subfamily branching, (2) monophyly of Coregonus and (3) placement of Parahucho. Coregoninae and Thymallinae are resolved as sister taxa, although conflicting topologies are found across analytical strategies. We find inconsistent and generally low support for the monophyly of Coregonus, including in results of analyses with the most extensive dataset and complex model. The most consistent placement of Parahucho is as sister lineage of Salmo.

A quantitative workflow for modeling diversification in material culture

Questions about the evolution of material culture are widespread in the humanities and social sciences. Statistical modeling of long-term changes in material culture is less common due to a lack of appropriate frameworks. Our goal is to close this gap and provide robust statistical methods for examining changes in the diversity of material culture. We provide an open-source and quantitative workflow for estimating rates of origination, extinction, and preservation, as well as identifying key shift points in the diversification histories of material culture. We demonstrate our approach using two distinct kinds of data: age ranges for the production of American car models, and radiocarbon dates associated with archaeological cultures of the European Neolithic. Our approach improves on existing frameworks by disentangling the relative contributions of origination and extinction to diversification. Our method also permits rigorous statistical testing of competing hypotheses to explain changes in diversity. Finally, we stress the value of a flexible approach that can be applied to data in various forms; this flexibility allows scholars to explore commonalities between forms of material culture and ask questions about the general properties of cultural change.

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.

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.

The Evolution of Color Pattern in Butterflyfishes (Chaetodontidae)

Coral reef fishes constitute one of the most diverse assemblages of vertebrates on the planet. Color patterns are known to serve a number of functions including intra- and inter-specific signaling, camouflage, mimicry, and defense. However, the relative importance of these and other factors in shaping color pattern evolution is poorly understood. Here we conduct a comparative phylogenetic analysis of color pattern evolution in the butterflyfishes (Chaetodontidae). Using recently developed tools for quantifying color pattern geometry as well as machine learning approaches, we investigate the tempo of evolution of color pattern elements and test whether ecological variables relating to defense, depth, and social behavior predict color pattern evolution. Butterflyfishes exhibit high diversity in measures of chromatic conspicuousness and the degrees of fine versus gross scale color patterning. Surprisingly, most diversity in color pattern was not predicted by any of the measures of ecology in our study, although we did find a significant but weak relationship between the level of fine scale patterning and some aspects of defensive morphology. We find that the tempo of color pattern diversification in butterflyfishes has increased toward the present and suggest that rapid evolution, presumably in response to evolutionary pressures surrounding speciation and lineage divergence, has effectively decoupled color pattern geometry from some aspects of ecology. Machine learning classification of color pattern appears to rely on a set of features that are weakly correlated with current color pattern geometry descriptors, but that may be better suited for the detection of discrete components of color pattern. A key challenge for future studies lies in determining whether rapid evolution has generally decoupled color patterns from ecology, or whether convergence in function produces convergence in color pattern at phylogenetic scales.

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.

Contemporaneous radiations of fungi and plants linked to symbiosis

Interactions between fungi and plants, including parasitism, mutualism, and saprotrophy, have been invoked as key to their respective macroevolutionary success. Here we evaluate the origins of plant-fungal symbioses and saprotrophy using a time-calibrated phylogenetic framework that reveals linked and drastic shifts in diversification rates of each kingdom. Fungal colonization of land was associated with at least two origins of terrestrial green algae and preceded embryophytes (as evidenced by losses of fungal flagellum, ca. 720 Ma), likely facilitating terrestriality through endomycorrhizal and possibly endophytic symbioses. The largest radiation of fungi (Leotiomyceta), the origin of arbuscular mycorrhizae, and the diversification of extant embryophytes occurred ca. 480 Ma. This was followed by the origin of extant lichens. Saprotrophic mushrooms diversified in the Late Paleozoic as forests of seed plants started to dominate the landscape. The subsequent diversification and explosive radiation of Agaricomycetes, and eventually of ectomycorrhizal mushrooms, were associated with the evolution of Pinaceae in the Mesozoic, and establishment of angiosperm-dominated biomes in the Cretaceous.

Plants and fungi interact widely and in diverse ways, from mutualism to parasitism and decomposition. Here, Lutzoni et al. analyse the timing of plant and fungal evolutionary radiations and identify four major periods in which plant-fungal interactions likely drove lineage diversification.

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.

Ice ages and butterflyfishes: Phylogenomics elucidates the ecological and evolutionary history of reef fishes in an endemism hotspot

For tropical marine species, hotspots of endemism occur in peripheral areas furthest from the center of diversity, but the evolutionary processes that lead to their origin remain elusive. We test several hypotheses related to the evolution of peripheral endemics by sequencing ultraconserved element (UCE) loci to produce a genome-scale phylogeny of 47 butterflyfish species (family Chaetodontidae) that includes all shallow water butterflyfish from the coastal waters of the Arabian Peninsula (i.e., Red Sea to Arabian Gulf) and their close relatives. Bayesian tree building methods produced a well-resolved phylogeny that elucidated the origins of butterflyfishes in this hotspots of endemism. We show that UCEs, often used to resolve deep evolutionary relationships, represent an important tool to assess the mechanisms underlying recently diverged taxa. Our analyses indicate that unique environmental conditions in the coastal waters of the Arabian Peninsula probably contributed to the formation of endemic butterflyfishes. Older endemic species are also associated with narrow versus broad depth ranges, suggesting that adaptation to deeper coral reefs in this region occurred only recently (<1.75 Ma). Even though deep reef environments were drastically reduced during the extreme low sea level stands of glacial ages, shallow reefs persisted, and as such there was no evidence supporting mass extirpation of fauna in this region.

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.

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.

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.].

Early-branching euteleost relationships: areas of congruence between concatenation and coalescent model inferences

Phylogenetic inference based on evidence from DNA sequences has led to significant strides in the development of a stable and robustly supported framework for the vertebrate tree of life. To date, the bulk of those advances have relied on sequence data from a small number of genome regions that have proven unable to produce satisfactory answers to consistently recalcitrant phylogenetic questions. Here, we re-examine phylogenetic relationships among early-branching euteleostean fish lineages classically grouped in the Protacanthopterygii using DNA sequence data surrounding ultraconserved elements. We report and examine a dataset of thirty-four OTUs with 17,957 aligned characters from fifty-three nuclear loci. Phylogenetic analysis is conducted in concatenated, joint gene trees and species tree estimation and summary coalescent frameworks. All analytical frameworks yield supporting evidence for existing hypotheses of relationship for the placement of Lepidogalaxias salamandroides, monophyly of the Stomiatii and the presence of an esociform + salmonid clade. Lepidogalaxias salamandroides and the Esociformes + Salmoniformes are successive sister lineages to all other euteleosts in the majority of analyses. The concatenated and joint gene trees and species tree analysis types produce high support values for this arrangement. However, inter-relationships of Argentiniformes, Stomiatii and Neoteleostei remain uncertain as they varied by analysis type while receiving strong and contradictory indices of support. Topological differences between analysis types are also apparent within the otomorph and the percomorph taxa in the data set. Our results identify concordant areas with strong support for relationships within and between early-branching euteleost lineages but they also reveal limitations in the ability of larger datasets to conclusively resolve other aspects of that phylogeny.

Body shape convergence driven by small size optimum in marine angelfishes

Convergent evolution of small body size occurs across many vertebrate clades and may reflect an evolutionary response to shared selective pressures. However it remains unclear if other aspects of phenotype undergo convergent evolution in miniaturized lineages. Here we present a comparative analysis of body size and shape evolution in marine angelfishes (Pomacanthidae), a reef fish family characterized by repeated transitions to small body size. We ask if lineages that evolve small sizes show convergent evolution in body shape. Our results reveal that angelfish lineages evolved three different stable size optima with one corresponding to the group of pygmy angelfishes (Centropyge). Then, we test if the observed shifts in body size are associated with changes to new adaptive peaks in shape. Our data suggest that independent evolution to small size optima have induced repeated convergence upon deeper body and steeper head profile in Centropyge These traits may favour manoeuvrability and visual awareness in these cryptic species living among corals, illustrating that functional demands on small size may be related to habitat specialization and predator avoidance. The absence of shape convergence in large marine angelfishes also suggests that more severe requirements exist for small than for large size optima.

Phylogenomic analysis of carangimorph fishes reveals flatfish asymmetry arose in a blink of the evolutionary eye

Background

Flatfish cranial asymmetry represents one of the most remarkable morphological innovations among vertebrates, and has fueled vigorous debate on the manner and rate at which strikingly divergent phenotypes evolve. A surprising result of many recent molecular phylogenetic studies is the lack of support for flatfish monophyly, where increasingly larger DNA datasets of up to 23 loci have either yielded a weakly supported flatfish clade or indicated the group is polyphyletic. Lack of resolution for flatfish relationships has been attributed to analytical limitations for dealing with processes such as nucleotide non-stationarity and incomplete lineage sorting (ILS). We tackle this phylogenetic problem using a sequence dataset comprising more than 1,000 ultraconserved DNA element (UCE) loci covering 45 carangimorphs, the broader clade containing flatfishes and several other specialized lineages such as remoras, billfishes, and archerfishes.

Results

We present a phylogeny based on UCE loci that unequivocally supports flatfish monophyly and a single origin of asymmetry. We document similar levels of discordance among UCE loci as in previous, smaller molecular datasets. However, relationships among flatfishes and carangimorphs recovered from multilocus concatenated and species tree analyses of our data are robust to the analytical framework applied and size of data matrix used. By integrating the UCE data with a rich fossil record, we find that the most distinctive carangimorph bodyplans arose rapidly during the Paleogene (66.0–23.03 Ma). Flatfish asymmetry, for example, likely evolved over an interval of no more than 2.97 million years.

Conclusions

The longstanding uncertainty in phylogenetic hypotheses for flatfishes and their carangimorph relatives highlights the limitations of smaller molecular datasets when applied to successive, rapid divergences. Here, we recovered significant support for flatfish monophyly and relationships among carangimorphs through analysis of over 1,000 UCE loci. The resulting time-calibrated phylogeny points to phenotypic divergence early within carangimorph history that broadly matches with the predictions of adaptive models of lineage diversification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0786-x) contains supplementary material, which is available to authorized users.

Replicated divergence in cichlid radiations mirrors a major vertebrate innovation

Decoupling of the upper jaw bones--jaw kinesis--is a distinctive feature of the ray-finned fishes, but it is not clear how the innovation is related to the extraordinary diversity of feeding behaviours and feeding ecology in this group. We address this issue in a lineage of ray-finned fishes that is well known for its ecological and functional diversity--African rift lake cichlids. We sequenced ultraconserved elements to generate a phylogenomic tree of the Lake Tanganyika and Lake Malawi cichlid radiations. We filmed a diverse array of over 50 cichlid species capturing live prey and quantified the extent of jaw kinesis in the premaxillary and maxillary bones. Our combination of phylogenomic and kinematic data reveals a strong association between biting modes of feeding and reduced jaw kinesis, suggesting that the contrasting demands of biting and suction feeding have strongly influenced cranial evolution in both cichlid radiations.

Evolution: Bioluminescent Courtship as an Engine of Diversity

A new study finds that the evolution of bioluminescent sexual displays drives high species richness across animal lineages, providing a crucial link between microevolutionary and macroevolutionary explanations of biodiversity.

Dietary breadth is positively correlated with venom complexity in cone snails

BACKGROUND

Although diet is believed to be a major factor underlying the evolution of venom, few comparative studies examine both venom composition and diet across a radiation of venomous species. Cone snails within the family, Conidae, comprise more than 700 species of carnivorous marine snails that capture their prey by using a cocktail of venomous neurotoxins (conotoxins or conopeptides). Venom composition across species has been previously hypothesized to be shaped by (a) prey taxonomic class (i.e., worms, molluscs, or fish) and (b) dietary breadth. We tested these hypotheses under a comparative phylogenetic framework using ecological data from past studies in conjunction with venom duct transcriptomes sequenced from 12 phylogenetically disparate cone snail species, including 10 vermivores (worm-eating), one molluscivore, and one generalist.

RESULTS

We discovered 2223 unique conotoxin precursor peptides that encoded 1864 unique mature toxins across all species, >90 % of which are new to this study. In addition, we identified two novel gene superfamilies and 16 novel cysteine frameworks. Each species exhibited unique venom profiles, with venom composition and expression patterns among species dominated by a restricted set of gene superfamilies and mature toxins. In contrast with the dominant paradigm for interpreting Conidae venom evolution, prey taxonomic class did not predict venom composition patterns among species. We also found a significant positive relationship between dietary breadth and measures of conotoxin complexity.

CONCLUSIONS

The poor performance of prey taxonomic class in predicting venom components suggests that cone snails have either evolved species-specific expression patterns likely as a consequence of the rapid evolution of conotoxin genes, or that traditional means of categorizing prey type (i.e., worms, mollusc, or fish) and conotoxins (i.e., by gene superfamily) do not accurately encapsulate evolutionary dynamics between diet and venom composition. We also show that species with more generalized diets tend to have more complex venoms and utilize a greater number of venom genes for prey capture. Whether this increased gene diversity confers an increased capacity for evolutionary change remains to be tested. Overall, our results corroborate the key role of diet in influencing patterns of venom evolution in cone snails and other venomous radiations.

Boxfish swimming paradox resolved: forces by the flow of water around the body promote manoeuvrability

The shape of the carapace protecting the body of boxfishes has been attributed an important hydrodynamic role in drag reduction and in providing automatic, flow-direction realignment and is therefore used in bioinspired design of cars. However, tight swimming-course stabilization is paradoxical given the frequent, high-performance manoeuvring that boxfishes display in their spatially complex, coral reef territories. Here, by performing flow-tank measurements of hydrodynamic drag and yaw moments together with computational fluid dynamics simulations, we reverse several assumptions about the hydrodynamic role of the boxfish carapace. Firstly, despite serving as a model system in aerodynamic design, drag-reduction performance was relatively low compared with more generalized fish morphologies. Secondly, the current theory of course stabilization owing to flow over the boxfish carapace was rejected, as destabilizing moments were found consistently. This solves the boxfish swimming paradox: destabilizing moments enhance manoeuvrability, which is in accordance with the ecological demands for efficient turning and tilting.

The evolution of capture myopathy in hooved mammals: a model for human stress cardiomyopathy?

Background and objectives: Capture myopathy (CM) syndromes in wildlife may be a model for human stress cardiomyopathy, including Takotsubo cardiomyopathy. Emotional stress or grief may trigger heart attack-like symptoms, and occasionally, sudden death in some humans. Similarly, wildlife exposed to predatory stresses, chase, or capture occasionally results in sudden death. To better understand the nature of vulnerability to stress-induced sudden death, we studied cases of CM in hooved mammals—ungulates—and hypothesized that CM would be associated with a syndrome of longevity-related traits.

Methodology: We reconstructed the evolution of CM in ungulates then determined how a set of life history traits explained variation in the likelihood that CM was reported.

Results: CM is broadly reported, but not in all genera, and phylogenetic analyses suggest that it is an evolutionarily labile trait. We found that the following traits were significantly associated with reports of CM: greater brain mass, faster maximum running speed, greater minimum group size and greater maximum longevity.

Conclusions and implications: CM may be an unavoidable consequence of adaptations to reduce predation risk that include increased running speed, sociality and having larger brains. Moreover, longer-lived species seem to be more likely to be susceptible to CM. Exploring variable susceptibility to CM highlights the evolutionary origins of the disorder, potential basic mechanisms that underlie vulnerability to the phenomenon, and the potential for reduction of risk through modification of life history trajectory.

Genome-wide ultraconserved elements exhibit higher phylogenetic informativeness than traditional gene markers in percomorph fishes

Ultraconserved elements (UCEs) have become popular markers in phylogenomic studies because of their cost effectiveness and their potential to resolve problematic phylogenetic relationships. Although UCE datasets typically contain a much larger number of loci and sites than more traditional datasets of PCR-amplified, single-copy, protein coding genes, a fraction of UCE sites are expected to be part of a nearly invariant core, and the relative performance of UCE datasets versus protein coding gene datasets is poorly understood. Here we use phylogenetic informativeness (PI) to compare the resolving power of multi-locus and UCE datasets in a sample of percomorph fishes with sequenced genomes (genome-enabled). We compare three data sets: UCE core regions, flanking sequence adjacent to the UCE core and a set of ten protein coding genes commonly used in fish systematics. We found the net informativeness of UCE core and flank regions to be roughly ten-fold and 100-fold more informative than that of the protein coding genes. On a per locus basis UCEs and protein coding genes exhibited similar levels of phylogenetic informativeness. Our results suggest that UCEs offer enormous potential for resolving relationships across the percomorph tree of life.

Trait decoupling promotes evolutionary diversification of the trophic and acoustic system of damselfishes

Trait decoupling, wherein evolutionary release of constraints permits specialization of formerly integrated structures, represents a major conceptual framework for interpreting patterns of organismal diversity. However, few empirical tests of this hypothesis exist. A central prediction, that the tempo of morphological evolution and ecological diversification should increase following decoupling events, remains inadequately tested. In damselfishes (Pomacentridae), a ceratomandibular ligament links the hyoid bar and lower jaws, coupling two main morphofunctional units directly involved in both feeding and sound production. Here, we test the decoupling hypothesis by examining the evolutionary consequences of the loss of the ceratomandibular ligament in multiple damselfish lineages. As predicted, we find that rates of morphological evolution of trophic structures increased following the loss of the ligament. However, this increase in evolutionary rate is not associated with an increase in trophic breadth, but rather with morphofunctional specialization for the capture of zooplanktonic prey. Lineages lacking the ceratomandibular ligament also shows different acoustic signals (i.e. higher variation of pulse periods) from others, resulting in an increase of the acoustic diversity across the family. Our results support the idea that trait decoupling can increase morphological and behavioural diversity through increased specialization rather than the generation of novel ecotypes.

Biogeography of the marmosets and tamarins (Callitrichidae)

The marmosets and tamarins, Family Callitrichidae, are Neotropical primates with over 60 species and subspecies that inhabit much of South America. Although callitrichids exhibit a remarkable widespread distribution, attempts to unravel their biogeographic history have been limited by taxonomic confusion and the lack of an appropriate statistical biogeographic framework. Here, we construct a time-calibrated multi-locus phylogeny from GenBank data and the callitrichid literature for 38 taxa. We use this framework to conduct statistical biogeographic analyses of callitrichids using BioGeoBEARS. The DIVAj model is the best supported reconstruction of biogeographic history among our analyses and suggests that the most recent common ancestor to the callitrichids was widespread across forested regions c. 14 Ma. There is also support for multiple colonizations of the Atlantic forest region from the Amazon basin, first by Leontopithecus c. 11 Ma and later by Callithrix c. 5 Ma. Our results show support for a 9 million year old split between a small-bodied group and large-bodied group of tamarins. These phylogenetic data, in concert with the consistent difference in body size between the two groups and geographical patterns (small-bodied tamarins and large-bodied tamarins have an unusually high degree of geographic overlap for congeners) lend support to our suggestion to split Saguinus into two genera, and we propose the use of distinct generic names; Leontocebus and Saguinus, respectively.

geiger v2.0: an expanded suite of methods for fitting macroevolutionary models to phylogenetic trees

SUMMARY

Phylogenetic comparative methods are essential for addressing evolutionary hypotheses with interspecific data. The scale and scope of such data have increased dramatically in the past few years. Many existing approaches are either computationally infeasible or inappropriate for data of this size. To address both of these problems, we present geiger v2.0, a complete overhaul of the popular R package geiger. We have reimplemented existing methods with more efficient algorithms and have developed several new approaches for accomodating heterogeneous models and data types.

AVAILABILITY AND IMPLEMENTATION

 This R package is available on the CRAN repository http://cran.r-project.org/web/packages/geiger/. All source code is also available on github http://github.com/mwpennell/geiger-v2. geiger v2.0 depends on the ape package.

CONTACT

mwpennell@gmail.com

SUPPLEMENTARY INFORMATION

 Supplementary data are available at Bioinformatics online.

Body fineness ratio as a predictor of maximum prolonged-swimming speed in coral reef fishes

The ability to sustain high swimming speeds is believed to be an important factor affecting resource acquisition in fishes. While we have gained insights into how fin morphology and motion influences swimming performance in coral reef fishes, the role of other traits, such as body shape, remains poorly understood. We explore the ability of two mechanistic models of the causal relationship between body fineness ratio and endurance swimming-performance to predict maximum prolonged-swimming speed (Umax ) among 84 fish species from the Great Barrier Reef, Australia. A drag model, based on semi-empirical data on the drag of rigid, submerged bodies of revolution, was applied to species that employ pectoral-fin propulsion with a rigid body at U max. An alternative model, based on the results of computer simulations of optimal shape in self-propelled undulating bodies, was applied to the species that swim by body-caudal-fin propulsion at Umax . For pectoral-fin swimmers, Umax increased with fineness, and the rate of increase decreased with fineness, as predicted by the drag model. While the mechanistic and statistical models of the relationship between fineness and Umax were very similar, the mechanistic (and statistical) model explained only a small fraction of the variance in Umax . For body-caudal-fin swimmers, we found a non-linear relationship between fineness and Umax , which was largely negative over most of the range of fineness. This pattern fails to support either predictions from the computational models or standard functional interpretations of body shape variation in fishes. Our results suggest that the widespread hypothesis that a more optimal fineness increases endurance-swimming performance via reduced drag should be limited to fishes that swim with rigid bodies.

A multi-locus molecular timescale for the origin and diversification of eels (Order: Anguilliformes)

Anguilliformes are an ecologically diverse group of predominantly marine fishes whose members are easily recognized by their extremely elongate bodies, and universal lack of pelvic fins. Recent studies based on mitochondrial loci, including full mitogenomes, have called into question the monophyly of both the Anguilliformes, which appear to be paraphyletic without the inclusion of the Saccopharyngiformes (gulper eels and allies), as well as other more commonly known eel families (e.g., Congridae, Serrivomeridae). However, no study to date has investigated anguilliform interrelationships using nuclear loci. Here we present a new phylogenetic hypothesis for the Anguilliformes based on five markers (the nuclear loci Early Growth Hormone 3, Myosin Heavy Polypeptide 6 and Recombinase Activating Gene 1, as well as the mitochondrial genes Cytochrome b and Cytochrome Oxidase I). Our sampling spans 148 species and includes 19 of the 20 extant families of anguilliforms and saccopharyngiforms. Maximum likelihood analysis reveals that saccopharyngiform eels are deeply nested within the anguilliforms, and supports the non-monophyly of Congridae and Nettastomatidae, as well as that of Derichthyidae and Chlopsidae. Our analyses suggest that Protanguilla may be the sister group of the Synaphobranchidae, though the recent hypothesis that this species is the sister group to all other anguilliforms cannot be rejected. The molecular phylogeny, time-calibrated using a Bayesian relaxed clock approach and seven fossil calibration points, reveals a Late Cretaceous origin of this expanded anguilliform clade (stem age ~116 Ma, crown age ~99 Ma). Most major (family level) lineages originated between the end of the Cretaceous and Early Eocene, suggesting that anguilliform radiation may have been facilitated by the recovery of marine ecosystems following the KP extinction.