Comparing Likelihood Ratios to Understand Genome-Wide Variation in Phylogenetic Support.

Genomic data have only sometimes brought resolution to the tree of life. Large phylogenomic studies can reach conflicting conclusions about important relationships, with mutually exclusive hypotheses receiving strong support. Reconciling such differences requires a detailed understanding of how phylogenetic signal varies among data sets. Two complementary strategies for better understanding phylogenomic conflicts are to examine support on a locus-by-locus basis and use support values that capture a larger range of variation in phylogenetic information, such as likelihood ratios. Likelihood ratios can be calculated using either maximum or marginal likelihoods. Despite being conceptually similar, differences in how these ratios are calculated and interpreted have not been closely examined in phylogenomics. Here, we compare the behavior of maximum and marginal likelihood ratios when evaluating alternate resolutions of recalcitrant relationships among major squamate lineages. We find that these ratios are broadly correlated between loci, but the correlation is driven by extreme values. As a consequence, the proportion of loci that support a hypothesis can change depending on which ratio is used and whether smaller values are discarded. In addition, maximum likelihood ratios frequently exhibit identical support for alternate hypotheses, making conflict resolution a challenge. We find surprising support for a sister relationship between snakes and iguanians across four different phylogenomic data sets in contrast to previous empirical studies. [Bayes factors; likelihood ratios; marginal likelihood; maximum likelihood; phylogenomics; squamates.].

Locomotory mode transitions alter phenotypic evolution and lineage diversification in an ecologically rich clade of mammals.

The relationship between organismal function and form is a cornerstone of biology because functional diversity is key to generating and maintaining ecological diversity. Morphological changes often occur in unison with behavioral or ecological transitions, and this process may foster diversification, but alternately could trap a species on an adaptive peak. We estimated the most comprehensive phylogenetic hypothesis of Murinae, a young (∼15 million years) and diverse (∼700 species) clade of mammals. We then tested for correlated evolution among four morphological traits with potential links to locomotor modes (Arboreal, General, Terrestrial, and Amphibious), then investigated the effects of locomotion on morphological and lineage diversification. We found unique combinations of trait values for each locomotor mode, including strong covariance between the tail and hindfoot lengths of specialized Arboreal and ecologically flexible General species. Low diversification rates and long branch lengths suggest that specialized lineages represent stable evolutionary "cul-de-sacs." General species, characterized by the classic "rat-like" body plan and broad locomotor abilities, have narrow optimal trait values and slow phenotypic evolution, but high lineage diversification rates. Our findings suggest that versatile, generalist forms act as seeds of species diversity and morphological specialization, which together build ecologically diverse radiations.

Experimental removal of introduced slider turtles offers new insight into competition with a native, threatened turtle.

The red-eared slider turtle (Trachemys scripta elegans; RES) is often considered one of the world's most invasive species. Results from laboratory and mesocosm experiments suggest that introduced RES outcompete native turtles for key ecological resources, but such experiments can overestimate the strength of competition. We report on the first field experiment with a wild turtle community, involving introduced RES and a declining native species of conservation concern, the western pond turtle (Emys marmorata; WPT). Using a before/after experimental design, we show that after removing most of an introduced RES population, the remaining RES dramatically shifted their spatial basking distribution in a manner consistent with strong intraspecific competition. WPT also altered their spatial basking distribution after the RES removal, but in ways inconsistent with strong interspecific competition. However, we documented reduced levels of WPT basking post-removal, which may reflect a behavioral shift attributable to the lower density of the turtle community. WPT body condition also increased after we removed RES, consistent with either indirect or direct competition between WPT and RES and providing the first evidence that RES can compete with a native turtle in the wild. We conclude that the negative impacts on WPT basking by RES in natural contexts are more limited than suggested by experiments with captive turtles, although wild WPT do appear to compete for food with introduced RES. Our results highlight the importance of manipulative field experiments when studying biological invasions, and the potential value of RES removal as a management strategy for WPT.

P3: Phylogenetic Posterior Prediction in RevBayes.

Tests of absolute model fit are crucial in model-based inference because poorly structured models can lead to biased parameter estimates. In Bayesian inference, posterior predictive simulations can be used to test absolute model fit. However, such tests have not been commonly practiced in phylogenetic inference due to a lack of convenient and flexible software. Here, we describe our newly implemented tests of model fit using posterior predictive testing, based on both data- and inference-based test statistics, in the phylogenetics software RevBayes. This new implementation makes a large spectrum of models available for use through a user-friendly and flexible interface.

Phylogenomic analyses of 539 highly informative loci dates a fully resolved time tree for the major clades of living turtles (Testudines).

Accurate time-calibrated phylogenies are the centerpiece of many macroevolutionary studies, and the relationship between the size and scale of molecular data sets and the density and accuracy of fossil calibrations is a key element of time tree studies. Here, we develop a target capture array specifically for living turtles, compare its efficiency to an ultraconserved element (UCE) dataset, and present a time-calibrated molecular phylogeny based on 539 nuclear loci sequenced from 26 species representing the breadth of living turtle diversity plus outgroups. Our gene array, based on three fully sequenced turtle genomes, is 2.4 times more variable across turtles than a recently published UCE data set for an identical subset of 13 species, confirming that taxon-specific arrays return more informative data per sequencing effort than UCEs. We used our genomic data to estimate the ages of living turtle clades including a mid-late Triassic origin for crown turtles and a mid-Carboniferous split of turtles from their sister group, Archosauria. By specifically excluding several of the earliest potential crown turtle fossils and limiting the age of fossil calibration points to the unambiguous crown lineage Caribemys oxfordiensis from the Late Jurassic (Oxfordian, 163.5-157.3Ma) we corroborate a relatively ancient age for living turtles. We also provide novel age estimates for five of the ten testudine families containing more than a single species, as well as several intrafamilial clades. Most of the diversity of crown turtles appears to date to the Paleogene, well after the Cretaceous-Paleogene mass extinction 66mya.

Exon capture optimization in amphibians with large genomes.

Gathering genomic-scale data efficiently is challenging for nonmodel species with large, complex genomes. Transcriptome sequencing is accessible for organisms with large genomes, and sequence capture probes can be designed from such mRNA sequences to enrich and sequence exonic regions. Maximizing enrichment efficiency is important to reduce sequencing costs, but relatively few data exist for exon capture experiments in nonmodel organisms with large genomes. Here, we conducted a replicated factorial experiment to explore the effects of several modifications to standard protocols that might increase sequence capture efficiency for amphibians and other taxa with large, complex genomes. Increasing the amounts of c0 t-1 repetitive sequence blocker and individual input DNA used in target enrichment reactions reduced the rates of PCR duplication. This reduction led to an increase in the percentage of unique reads mapping to target sequences, essentially doubling overall efficiency of the target capture from 10.4% to nearly 19.9% and rendering target capture experiments more efficient and affordable. Our results indicate that target capture protocols can be modified to efficiently screen vertebrates with large genomes, including amphibians.