Large Multi-locus Datasets Research Articles

The evolution of multi-component weapons in the superfamily of leaf-footed bugs.

Sexually selected weapons, such as the antlers of deer, claws of crabs, and tusks of beaked whales, are strikingly diverse across taxa and even within groups of closely related species. Phylogenetic comparative studies have typically taken a simplified approach to investigate the evolution of weapon diversity, examining the gains and losses of entire weapons, major shifts in size or type, or changes in location. Less understood is how individual weapon components evolve and assemble into a complete weapon. We addressed this question by examining weapon evolution in the diverse, multi-component hind-leg and body weapons of leaf-footed bugs, superfamily Coreoidea (Hemiptera: Heteroptera). Male leaf-footed bugs use their morphological weapons to fight for access to mating territories. We used a large multilocus dataset comprised of ultraconserved element loci for 248 species and inferred evolutionary transitions among component states using ancestral state estimation. Our results suggest that weapons added components over time with some evidence of a cyclical evolutionary pattern-gains of components followed by losses and then gains again. Furthermore, our best estimate indicated that certain trait combinations evolved repeatedly across the phylogeny, suggesting that they function together in battle or that they are genetically correlated. This work reveals the remarkable and dynamic evolution of weapon form in the leaf-footed bugs and provides insights into weapon assembly and disassembly over evolutionary time.

Quantitatively defining species boundaries with more efficiency and more biological realism

We introduce a widely applicable species delimitation method based on the multispecies coalescent model that is more efficient and more biologically realistic than existing methods. We extend a threshold-based method to allow the ancestral speciation rate to vary through time as a smooth piecewise function. Furthermore, we introduce the cutting-edge proposal kernels of StarBeast3 to this model, thus enabling rapid species delimitation on large molecular datasets and allowing the use of relaxed molecular clock models. We validate these methods with genomic sequence data and SNP data, and show they are more efficient than existing methods at achieving parameter convergence during Bayesian MCMC. Lastly, we apply these methods to two datasets (Hemidactylus and Galagidae) and find inconsistencies with the published literature. Our methods are powerful for rapid quantitative testing of species boundaries in large multilocus datasets and are implemented as an open source BEAST 2 package called SPEEDEMON.

Geography is more important than life history in the recent diversification of the tiger salamander complex

The North American tiger salamander species complex, including its best-known species, the Mexican axolotl, has long been a source of biological fascination. The complex exhibits a wide range of variation in developmental life history strategies, including populations and individuals that undergo metamorphosis; those able to forego metamorphosis and retain a larval, aquatic lifestyle (i.e., paedomorphosis); and those that do both. The evolution of a paedomorphic life history state is thought to lead to increased population genetic differentiation and ultimately reproductive isolation and speciation, but the degree to which it has shaped population- and species-level divergence is poorly understood. Using a large multilocus dataset from hundreds of samples across North America, we identified genetic clusters across the geographic range of the tiger salamander complex. These clusters often contain a mixture of paedomorphic and metamorphic taxa, indicating that geographic isolation has played a larger role in lineage divergence than paedomorphosis in this system. This conclusion is bolstered by geography-informed analyses indicating no effect of life history strategy on population genetic differentiation and by model-based population genetic analyses demonstrating gene flow between adjacent metamorphic and paedomorphic populations. This fine-scale genetic perspective on life history variation establishes a framework for understanding how plasticity, local adaptation, and gene flow contribute to lineage divergence. Many members of the tiger salamander complex are endangered, and the Mexican axolotl is an important model system in regenerative and biomedical research. Our results chart a course for more informed use of these taxa in experimental, ecological, and conservation research.

Phylogenomics and biogeography of the world's thrushes (Aves, Turdus): new evidence for a more parsimonious evolutionary history.

To elucidate the relationships and spatial range evolution across the world of the bird genus Turdus (Aves), we produced a large genomic dataset comprising ca 2 million nucleotides for ca 100 samples representing 53 species, including over 2000 loci. We estimated time-calibrated maximum-likelihood and multispecies coalescent phylogenies and carried out biogeographic analyses. Our results indicate that there have been considerably fewer trans-oceanic dispersals within the genus Turdus than previously suggested, such that the Palaearctic clade did not originate in America and the African clade was not involved in the colonization of the Americas. Instead, our findings suggest that dispersal from the Western Palaearctic via the Antilles to the Neotropics might have occurred in a single event, giving rise to the rich Neotropical diversity of Turdus observed today, with no reverse dispersals to the Palaearctic or Africa. Our large multilocus dataset, combined with dense species-level sampling and analysed under probabilistic methods, brings important insights into historical biogeography and systematics, even in a scenario of fast and spatially complex diversification.

PhyloMAd: efficient assessment of phylogenomic model adequacy.

Statistical phylogenetic inference plays an important role in evolutionary biology. The accuracy of phylogenetic methods relies on having suitable models of the evolutionary process. Various tools allow comparisons of candidate phylogenetic models, but assessing the absolute performance of models remains a considerable challenge. We introduce PhyloMAd, a user-friendly application for assessing the adequacy of commonly used models of nucleotide substitution and among-lineage rate variation. Our software implements a fast, likelihood-based method of model assessment that is tractable for analyses of large multi-locus datasets. PhyloMAd provides a means of informing model improvement, or selecting data to enhance the evolutionary signal in phylogenomic analyses. PhyloMAd, together with a manual, a tutorial and the source code, are freely available from the GitHub repository github.com/duchene/phylomad.

Filtering of target sequence capture individuals facilitates species tree construction in the plant subtribe Iochrominae (Solanaceae)

Advances in sequencing technology have made it possible to produce large multi-locus datasets required for species tree analyses. One challenge with constructing high throughput sequencing datasets, however, is that missing information is propagated at different steps in the sequence preparation process. To date, species tree studies have focused on filtering and removing errors that occur at particular loci. Given the way that high throughput sequencing datasets are constructed, however, large amounts of error or ambiguity may also manifest across individuals. Here we use a novel tree-based multivariate clustering method to identify and remove individuals with low phylogenetic signal in a nuclear sequence capture dataset for the Iochrominae clade (Solanaceae). Our results suggest that the low quality tips are the result of the library preparation process (e.g. unequal pooling) rather than poor capture due to phylogenetic distance from the reference species. After implementing the clustering approach and removing low quality tips, we construct an Iochrominae species tree that resolves a number of unknown relationships. We propose this pipeline as a valuable tool for species tree reconstruction with phylogenomic datasets containing variable levels of missing data.

Population genomic data delineate conservation units in mottled ducks (Anas fulvigula)

Large multi-locus data sets provide an unprecedented opportunity for detecting fine-scale population structure and identifying populations that warrant independent conservation and management. In this study, we used 3161 polymorphic loci obtained from double-digest restriction-site associated DNA sequencing (ddRAD-seq) to examine genetic differentiation between populations of mottled ducks (Anas fulvigula) and their close relative, the mallard (A. platyrhynchos). Florida (FL) and western gulf coast (WGC) mottled ducks group into discrete populations in ordination analyses and can be assigned unambiguously to their population of origin based on multi-locus genotypes. Moreover, we find limited evidence of recent admixture between populations, and the range gap separating them is a significant barrier to gene flow. We found no evidence of diversifying selection in the sample of loci we analyzed; instead, the genome-wide pattern of differentiation between mottled duck populations is consistent with historical isolation and divergence by neutral genetic drift. Thus, our results indicate that FL and WGC mottled ducks are discrete populations on separate evolutionary pathways and should be managed as independent conservation units, and perhaps recognized as distinct subspecies given observed differences in plumage and morphology. In contrast, the recently established South Atlantic coastal zone population is indistinguishable from the WGC population but shows evidence of a small amount of admixture from FL mottled ducks. Our study demonstrates a quantitative, operational approach for delineating discrete populations for conservation using ddRAD-seq data; multi-locus genotypes allow for a robust test of historical isolation and evolutionary independence on much shorter evolutionary timescales than the overly conservative standard of reciprocal monophyly in mtDNA.

Resolving the higher-order phylogenetic relationships of the circumtropical Mabuya group (Squamata: Scincidae): An out-of-Asia diversification

Despite an abundance of phylogenetic studies focused on intrageneric relationships of members of the Mabuya group, the intergeneric relationships have remained difficult to resolve. The most-persistent unresolved regions of the phylogeny of the group include: (1) the placement of the Middle-Eastern Trachylepis with respect to the Afro-Malagasy Trachylepis and its taxonomic status; (2) the phylogenetic position of the Cape Verdean Chioninia within the larger Mabuya group; (3) support for the placement of Dasia with respect to the entire group; and (4) the phylogenetic placement of Eutropis novemcarinata with respect to other Eutropis and Dasia. In this study, we include representatives of all these taxa as well as African Eumecia and Neotropical Mabuya. We seek to address these phylogenetic and systematic issues by generating a well-resolved and supported phylogeny for the Mabuya group as a whole that can be used to develop a stable taxonomy and reconstruct the geographic patterns of diversification within the group. To meet these goals, we built a large multi-locus dataset of 11 markers (nine nuclear and two mitochondrial), and performed concatenated and species tree analyses to generate a well-supported phylogeny for the group. Statistical topology tests reject the monophyly of Middle-Eastern Trachylepis with Afro-Malagasy Trachylepis, and to reflect monophyly we place the Middle-Eastern species into a previously described genus, Heremites. Cape-Verdean Chioninia are resolved as the strongly supported sister-group to Afro-Malagasy Trachylepis. Monophyly of the Southeast-Asian genera, Eutropis and Dasia, is not supported, with a clade composed of Dasia+Eutropis novemcarinata more closely related to the rest of the Mabuya group than to the remaining Eutropis. The phylogenetic position of E. novemcarinata renders Eutropis polyphyletic, and we therefore describe and place E. novemcarinata into a new monotypic genus, Toenayar, to preserve monophyly among the genera. In light of these novel findings, we review and discuss the historical biogeography of the entire Mabuya group.

The art and science of species delimitation

Abstract DNA-based approaches to systematics have changed dramatically during the last two decades with the rise of DNA barcoding methods and newer multi-locus methods for species delimitation. During the last half-decade, partly driven by the new sequencing technologies, the focus has shifted to multi-locus sequence data and the identification of species within the framework of the multi-species coalescent (MSC). In this paper, I discuss model-based Bayesian methods for species delimitation that have been developed in recent years using the MSC. Several approximate methods for species delimitation (and their limitations) are also discussed. Explicit species delimitation models have the advantage of clarifying more precisely what is being delimited and what assumptions we are making in doing so. Moreover, the methods can be very powerful when applied to large multi-locus datasets and thus take full advantage of data generated using today’s technologies.

Forked Tongues Revisited: Molecular Apomorphies Support Morphological Hypotheses of Squamate Evolution

Incongruence between morphological and molecular-based phylogenetic hypotheses has been reported across a wide range of taxa. Specifically, morphological and molecular hypotheses of squamate phylogeny have been consistently incongruent and have been notoriously difficult to reconcile. With the ever-rising popular view of the superiority of molecular data over morphological data, studies have been biased toward seeking explanations for homoplasy in the morphological data. We propose considering approaches that do not make a priori assumptions about the superiority of one type of data over another. We present two case studies using the proto-oncogene c-mos, as well as a large multi-locus dataset, to examine apomorphy distributions across hypotheses. This approach reveals novel insights into data incongruence currently plaguing squamate phylogeny, with direct implications for other scenarios of phylogenetic incongruence.

Phylogeographic inference using Bayesian model comparison across a fragmented chorus frog species complex.

Fragmented species complexes provide an interesting system for investigating biogeographic history and the present distribution of genetic variation. Recent advances in sequencing technology and statistical phylogeography enable the collection and rigorous analysis of large multilocus data sets, but designing studies that produce meaningful phylogeographic inferences remains challenging. We implemented a Bayesian model comparison approach to investigate previous biogeographic hypotheses while simultaneously inferring the presence of genetic structure in a chorus frog species complex. The Illinois chorus frog (Pseudacris illinoensis), originally described as a subspecies of the broadly distributed Strecker's chorus frog (Pseudacris streckeri), occurs in small, disjunct regions associated with scarce sand prairie habitats that have been impacted by human development. We used high-throughput sequencing to develop and collect a multitiered genetic data set comprised of three different marker types (23 anonymous nuclear sequence loci, four mitochondrial genes and 14 microsatellite loci) designed to address questions across different evolutionary timescales. Phylogenetic analyses uncovered a deep divergence between populations in the Edwards Plateau of central Texas and all other P.streckeri/P. illinoensis populations, but suggest the disjunct distribution of P.illinoensis occurred more recently. Our best-supported migration model is consistent with the hypothesis that central Texas represented a refugium from which populations expanded via multiple routes. This model also indicates that disjunct northern and southern regions of P.illinoensis should be considered genetically distinct management units. Our study provides an evolutionary context for future studies and conservation efforts in P.illinoensis and demonstrates the utility of model-based approaches for phylogeographic inference.

A protocol for targeted enrichment of intron-containing sequence markers for recent radiations: A phylogenomic example from Heuchera (Saxifragaceae).

Premise of the study:Phylogenetic inference is moving to large multilocus data sets, yet there remains uncertainty in the choice of marker and sequencing method at low taxonomic levels. To address this gap, we present a method for enriching long loci spanning intron-exon boundaries in the genus Heuchera.Methods:Two hundred seventy-eight loci were designed using a splice-site prediction method combining transcriptomic and genomic data. Biotinylated probes were designed for enrichment of these loci. Reference-based assembly was performed using genomic references; additionally, chloroplast and mitochondrial genomes were used as references for off-target reads. The data were aligned and subjected to coalescent and concatenated phylogenetic analyses to demonstrate support for major relationships.Results:Complete or nearly complete (>99%) sequences were assembled from essentially all loci from all taxa. Aligned introns showed a fourfold increase in divergence as opposed to exons. Concatenated analysis gave decisive support to all nodes, and support was also high and relationships mostly similar in the coalescent analysis. Organellar phylogenies were also well-supported and conflicted with the nuclear signal.Discussion:Our approach shows promise for resolving a recent radiation. Enrichment for introns is highly successful with little or no sequencing dropout at low taxonomic levels despite higher substitution and indel frequencies, and should be exploited in studies of species complexes.

Multilocus assessment of phylogenetic relationships in Alytes (Anura, Alytidae)

With the advent of large multilocus datasets, molecular systematics is experiencing very rapid progress, but important challenges remain regarding data analysis and interpretation. Midwife toads (genus Alytes) exemplify two of the most widespread problems for accurate phylogenetic reconstruction: discerning the causes of discordance between gene trees, and resolving short internodes produced during rapid, successive splitting events. The three species in subgenus Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), the sister group to A. obstetricans, have disjunct and highly restricted geographical ranges, which are thought to result from old vicariant events affecting their common ancestor, but their phylogenetic relationships are still unresolved. In this study we re-address the phylogeny of Alytes with a special focus on the relationships in Baleaphryne with a multilocus dataset including >9000 base pairs of mitochondrial DNA and four nuclear markers (3142bp) in all recognized taxa, including all subspecies of A. obstetricans. Both concatenation and species tree analyses suggest that A. muletensis, endemic to the Balearic island of Mallorca, is the sister taxon to a clade comprising the southeastern Iberian endemic A. dickhilleni and the North African A. maurus. This scenario is consistent with palaeogeological evidence associated with the fragmentation of the Betic-Rifean Massif, followed by the opening of the Strait of Gibraltar. On the other hand, analyses of intraspecific variation in A. obstetricans are inconclusive regarding relationships between major clades and conflict with current subspecific taxonomy.

A comparison of supermatrix and supertree methods for multilocus phylogenetics using organismal datasets.

It has been proposed that supertree approaches should be applied to large multilocus datasets to achieve computational tractability. Large datasets such as those derived from phylogenomics studies can be broken into many locus-specific tree searches and the resulting trees can be stitched together via a supertree method. Using simulated data, workers have reported that they can rapidly construct a supertree that is comparable to the results of heuristic tree search on the entire dataset. To test this assertion with organismal data, we compare tree length under the parsimony criterion and computational time for 20 multilocus datasets using supertree (SuperFine and SuperTriplets) and supermatrix (heuristic search in TNT) approaches. Tree length and computational times were compared among methods using the Wilcoxon matched-pairs signed rank test. Supermatrix searches produced significantly shorter trees than either supertree approach (SuperFine or SuperTriplets; P<0.0002 in both cases). Moreover, the processing time of supermatrix search was significantly lower than SuperFine+locus-specific search (P<0.01) but roughly equivalent to that of SuperTriplets+locus-specific search (P>0.4, not significant). In conclusion, we show by using real rather than simulated data that there is no basis, either in time tractability or in tree length, for use of supertrees over heuristic tree search using a supermatrix for phylogenomics.

Concatenation and Concordance in the Reconstruction of Mouse Lemur Phylogeny: An Empirical Demonstration of the Effect of Allele Sampling in Phylogenetics

The systematics and speciation literature is rich with discussion relating to the potential for gene tree/species tree discordance. Numerous mechanisms have been proposed to generate discordance, including differential selection, long-branch attraction, gene duplication, genetic introgression, and/or incomplete lineage sorting. For speciose clades in which divergence has occurred recently and rapidly, recovering the true species tree can be particularly problematic due to incomplete lineage sorting. Unfortunately, the availability of multilocus or "phylogenomic" data sets does not simply solve the problem, particularly when the data are analyzed with standard concatenation techniques. In our study, we conduct a phylogenetic study for a nearly complete species sample of the dwarf and mouse lemur clade, Cheirogaleidae. Mouse lemurs (genus, Microcebus) have been intensively studied over the past decade for reasons relating to their high level of cryptic species diversity, and although there has been emerging consensus regarding the evolutionary diversity contained within the genus, there is no agreement as to the inter-specific relationships within the group. We attempt to resolve cheirogaleid phylogeny, focusing especially on the mouse lemurs, by employing a large multilocus data set. We compare the results of Bayesian concordance methods with those of standard gene concatenation, finding that though concatenation yields the strongest results as measured by statistical support, these results are found to be highly misleading. By employing an approach where individual alleles are treated as operational taxonomic units, we show that phylogenetic results are substantially influenced by the selection of alleles in the concatenation process.

IS A NEW AND GENERAL THEORY OF MOLECULAR SYSTEMATICS EMERGING?

The advent and maturation of algorithms for estimating species trees-phylogenetic trees that allow gene tree heterogeneity and whose tips represent lineages, populations and species, as opposed to genes-represent an exciting confluence of phylogenetics, phylogeography, and population genetics, and ushers in a new generation of concepts and challenges for the molecular systematist. In this essay I argue that to better deal with the large multilocus datasets brought on by phylogenomics, and to better align the fields of phylogeography and phylogenetics, we should embrace the primacy of species trees, not only as a new and useful practical tool for systematics, but also as a long-standing conceptual goal of systematics that, largely due to the lack of appropriate computational tools, has been eclipsed in the past few decades. I suggest that phylogenies as gene trees are a "local optimum" for systematics, and review recent advances that will bring us to the broader optimum inherent in species trees. In addition to adopting new methods of phylogenetic analysis (and ideally reserving the term "phylogeny" for species trees rather than gene trees), the new paradigm suggests shifts in a number of practices, such as sampling data to maximize not only the number of accumulated sites but also the number of independently segregating genes; routinely using coalescent or other models in computer simulations to allow gene tree heterogeneity; and understanding better the role of concatenation in influencing topologies and confidence in phylogenies. By building on the foundation laid by concepts of gene trees and coalescent theory, and by taking cues from recent trends in multilocus phylogeography, molecular systematics stands to be enriched. Many of the challenges and lessons learned for estimating gene trees will carry over to the challenge of estimating species trees, although adopting the species tree paradigm will clarify many issues (such as the nature of polytomies and the star tree paradox), raise conceptually new challenges, or provide new answers to old questions.

Branch Lengths, Support, and Congruence: Testing the Phylogenomic Approach with 20 Nuclear Loci in Snakes

Many authors have claimed that short branches in the Tree of Life will be very difficult to resolve with strong support, even with the large multilocus data sets now made possible by genomic resources. Short branches may be especially problematic because the underlying gene trees are expected to have discordant phylogenetic histories when the time between branching events is very short. Although there are many examples of short branches that are difficult to resolve, surprisingly, no empirical studies have systematically examined the relationships between branch lengths, branch support, and congruence among genes. Here, we examine these fundamental relationships quantitatively using a data set of 20 nuclear loci for 50 species of snakes (representing most traditionally recognized families). A combined maximum likelihood analysis of the 20 loci gives strong support for 69% of the nodes, but many remain weakly supported, with bootstrap values for 20% ranging from 21% to 66%. For the combined-data tree, we find significant correlations between the length of a branch, levels of bootstrap support, and the proportion of genes that are congruent with that branch in the separate analyses of each gene. We also find that strongly supported conflicts between gene trees over the resolution of individual branches are common (roughly 35% of clades), especially for shorter branches. Overall, our results support the hypothesis that short branches may be very difficult to confidently resolve, even with large, multilocus data sets. Nevertheless, our study provides strong support for many clades, including several that were controversial or poorly resolved in previous studies of snake phylogeny.