An ultraconserved element probe set for velvet worms (Onychophora)

Onychophora are cryptic, soil-dwelling invertebrates known for their biogeographic affinities, diversity of reproductive modes, close phylogenetic relationship to arthropods, and peculiar prey capture mechanism. The 216 valid species of Onychophora are grouped into two families - Peripatopsidae and Peripatidae - and apart from a few relationships among major lineages within these two families, a stable phylogenetic backbone for the phylum has yet to be resolved. This has hindered our understanding of onychophoran biogeographic patterns, evolutionary history, and systematics. Neopatida, the Neotropical clade of peripatids, has proved particularly difficult, with recalcitrant nodes and low resolution, potentially due to rapid radiation of the group during the Cretaceous. Previous studies have had to compromise between number of loci and number of taxa due to limitations of Sanger sequencing and phylotranscriptomics, respectively. Additionally, aspects of their genome size and structure have made molecular phylogenetics difficult and data matrices have been affected by missing data. To address these issues, we leveraged recent, published transcriptomes and the first high quality genome for the phylum and designed a high affinity ultraconserved element (UCE) probe set for Onychophora. This new probe set, consisting of ∼ 20,000 probes that target 1,465 loci across both families, has high locus recovery and phylogenetic utility. Phylogenetic analyses recovered the monophyly of major clades of Onychophora and revealed a novel lineage from the Neotropics that challenges our current understanding of onychophoran biogeographic endemicity. This new resource could drastically increase the power of molecular datasets and potentially allow access to genomic scale data from archival museum specimens to further tackle the issues exasperating onychophoran systematics.

Plotting for change: an analytical framework to aid decisions on which lineages are candidate species in phylogenomic species discovery

A recent study argued that coalescent-based models of species delimitation mostly delineate population structure, not species, and called for the validation of candidate species using biological information additional to the genetic information, such as phenotypic or ecological data. Here, we introduce a framework to interrogate genomic datasets and coalescent-based species trees for the presence of candidate species in situations where additional biological data are unavailable, unobtainable or uninformative. For de novo genomic studies of species boundaries, we propose six steps: (1) visualize genetic affinities among individuals to identify both discrete and admixed genetic groups from first principles and to hold aside individuals involved in contemporary admixture for independent consideration; (2) apply phylogenetic techniques to identify lineages; (3) assess diagnosability of those lineages as potential candidate species; (4) interpret the diagnosable lineages in a geographical context (sympatry, parapatry, allopatry); (5) assess significance of difference or trends in the context of sampling intensity; and (6) adopt a holistic approach to available evidence to inform decisions on species status in the difficult cases of allopatry. We adopt this approach to distinguish candidate species from within-species lineages for a widespread species complex of Australian freshwater fishes (Retropinna spp.). Our framework addresses two cornerstone issues in systematics that are often not discussed explicitly in genomic species discovery: diagnosability and how to determine it, and what criteria should be used to decide whether diagnosable lineages are conspecific or represent different species.

Comparing Ultraconserved Elements and Exons for Phylogenomic Analyses of Middle American Cichlids: When Data Agree to Disagree

Choosing among types of genomic markers to be used in a phylogenomic study can have a major influence on the cost, design, and results of a study. Yet few attempts have been made to compare categories of next-generation sequence markers limiting our ability to compare the suitability of these different genomic fragment types. Here, we explore properties of different genomic markers to find if they vary in the accuracy of component phylogenetic trees and to clarify the causes of conflict obtained from different data sets or inference methods. As a test case, we explore the causes of discordance between phylogenetic hypotheses obtained using a novel data set of ultraconserved elements (UCEs) and a recently published exon data set of the cichlid tribe Heroini. Resolving relationships among heroine cichlids has historically been difficult, and the processes of colonization and diversification in Middle America and the Greater Antilles are not yet well understood. Despite differences in informativeness and levels of gene tree discordance between UCEs and exons, the resulting phylogenomic hypotheses generally agree on most relationships. The independent data sets disagreed in areas with low phylogenetic signal that were overwhelmed by incomplete lineage sorting and nonphylogenetic signals. For UCEs, high levels of incomplete lineage sorting were found to be the major cause of gene tree discordance, whereas, for exons, nonphylogenetic signal is most likely caused by a reduced number of highly informative loci. This paucity of informative loci in exons might be due to heterogeneous substitution rates that are problematic to model (i.e., computationally restrictive) resulting in systematic errors that UCEs (being less informative individually but more uniform) are less prone to. These results generally demonstrate the robustness of phylogenomic methods to accommodate genomic markers with different biological and phylogenetic properties. However, we identify common and unique pitfalls of different categories of genomic fragments when inferring enigmatic phylogenetic relationships.

Behavioral, morphological, and ecological trait evolution in two clades of New World Sparrows (Aimophila and Peucaea, Passerellidae)

The New World sparrows (Passerellidae) are a large, diverse group of songbirds that vary in morphology, behavior, and ecology. Thus, they are excellent for studying trait evolution in a phylogenetic framework. We examined lability versus conservatism in morphological and behavioral traits in two related clades of sparrows (Aimophila, Peucaea), and assessed whether habitat has played an important role in trait evolution. We first inferred a multi-locus phylogeny which we used to reconstruct ancestral states, and then quantified phylogenetic signal among morphological and behavioral traits in these clades and in New World sparrows more broadly. Behavioral traits have a stronger phylogenetic signal than morphological traits. Specifically, vocal duets and song structure are the most highly conserved traits, and nesting behavior appears to be maintained within clades. Furthermore, we found a strong correlation between open habitat and unpatterned plumage, complex song, and ground nesting. However, even within lineages that share the same habitat type, species vary in nesting, plumage pattern, song complexity, and duetting. Our findings highlight trade-offs between behavior, morphology, and ecology in sparrow diversification.

Phylogenomic data resolve higher-level relationships within South American Liolaemus lizards

Phylogenomic approaches now generate hundreds of loci representative of the whole genome that can be used for phylogenetic analyses. The South American lizard genus Liolaemus is the most species-rich vertebrate radiation from temperate zones (more than 265 described species), yet most higher-level phylogenetic relationships within Liolaemus remain poorly resolved. In this study, we used 584 nuclear loci collected using targeted sequenced capture to estimate the phylogenetic relationships among 26 species representing the two subgenera within Liolaemus (Eulaemus + Liolaemus), and all major groups within Eulaemus. Previous molecular and morphological-based phylogenetic analyses of Eulaemus based on a limited number of characters resolved few higher-level relationships, although one point of agreement is that the early divergence within Eulaemus corresponds to the lineomaculatus section, followed by the diversification of eight main clades that are strongly supported and recognized. Liolaemus probably experienced relatively rapid divergences during parts of its evolutionary history, and a phylogenomic approach was used to resolve the relationships among the major groups. The new analyses presented here support the division of Liolaemus into two subgenera, and resolve relationships among many of the major clades of Eulaemus with strong support. A Bayesian divergence dating analysis using 44 protein-coding genes provides an estimation of the split of the two Liolaemus subgenera of approximately 19,7 ma (95% HPD = 16,94-23,04), while diversification within Eulaemus started at 15,05 ma (95% HPD = 12,94 - 17,59) among the L. lineomaculatus and the L. montanus series by Mid Miocene. A novel phylogenetic network analyses for SNP data identified two hybridizing edges among different groups of Eulaemus at different points in time. Having a solid phylogenetic hypothesis of the main Eulaemus clades opens new opportunities to test a variety of macroevolutionary questions for this unique radiation.

Probabilistic Species Tree Distances: Implementing the Multispecies Coalescent to Compare Species Trees Within the Same Model-Based Framework Used to Estimate Them

Despite the ubiquitous use of statistical models for phylogenomic and population genomic inferences, this model-based rigor is rarely applied to post hoc comparison of trees. In a recent study, Garba et al. derived new methods for measuring the distance between two gene trees computed as the difference in their site pattern probability distributions. Unlike traditional metrics that compare trees solely in terms of geometry, these measures consider gene trees and associated parameters as probabilistic models that can be compared using standard information theoretic approaches. Consequently, probabilistic measures of phylogenetic tree distance can be far more informative than simply comparisons of topology and/or branch lengths alone. However, in their current form, these distance measures are not suitable for the comparison of species tree models in the presence of gene tree heterogeneity. Here, we demonstrate an approach for how the theory of Garba et al. (2018), which is based on gene tree distances, can be extended naturally to the comparison of species tree models. Multispecies coalescent (MSC) models parameterize the discrete probability distribution of gene trees conditioned upon a species tree with a particular topology and set of divergence times (in coalescent units), and thus provide a framework for measuring distances between species tree models in terms of their corresponding gene tree topology probabilities. We describe the computation of probabilistic species tree distances in the context of standard MSC models, which assume complete genetic isolation postspeciation, as well as recent theoretical extensions to the MSC in the form of network-based MSC models that relax this assumption and permit hybridization among taxa. We demonstrate these metrics using simulations and empirical species tree estimates and discuss both the benefits and limitations of these approaches. We make our species tree distance approach available as an R package called pSTDistanceR, for open use by the community.

Description and phylogeny of a new species of Liolaemus (Iguania: Liolaemidae) endemic to the south of the Plurinational State of Bolivia

The Liolaemus montanus group is a diverse group of lizards that ranges from central Peru to southwestern Mendoza, Argentina, including much of the Plurinational State of Bolivia ("Bolivia") and Chile. The species of this group mainly inhabit high elevation areas with cold temperatures. In the last years, several species of this group have been described, mostly in Argentina and Chile. In Bolivia, there are at least thirteen valid species belonging to the L. montanus group. In this study, we describe a new species of the L. montanus group with a marked endemism in the Cordillera de Sama of the Tarija Department, Bolivia, and a combination of unique character states that allows its formal description as a new species. The phylogenetic relationships based on analysis of 159 morphological characters suggest that it belongs to the L. montanus group, and that it is closest to Liolaemus pulcherrimus, which is found allopatrically in a small area of the Jujuy Province, Argentina. The multivariate analyses of 66 morphological characters support the phylogenetic relationships. Statistical analyses of inter-species comparisons of morphological characters are not considered the only methods due to the non-independence of some characters states among species; thus, a phylogenetic analysis is recommended. The detailed revision of specimens of the L. montanus group held in the collections of Bolivia is filling major geographic gaps and improving our understanding of the phylogenetic and biogeographic relationships of this widely distributed group of South American lizards.

Species delimitation and evolutionary reconstruction within an integrative taxonomic framework: A case study on Rhinolophus macrotis complex (Chiroptera: Rhinolophidae)

Species delimitation and evolutionary reconstruction remain challenging for non-model species that have experienced reticulate evolution and exhibit conflicting patterns of differentiation among multiple lines of evidence, such as mitochondrial and nuclear data and phenotypes. Here, we applied an integrative taxonomic approach to a case study of Rhinolophus macrotis complex, whose taxonomic status remains controversial, to provide insight into the systematics and evolutionary history of these species. By integrating traditional genetic markers with different modes of inheritance, genome-wide SNPs as well as phenotypic characteristics, we clarified the presence of three closely related species, R. episcopus, R. siamensis, and R. osgoodi, within this complex, and proposed a new taxonomic treatment for R. osgoodi. Our results suggested that hybridization and introgression are the main causes of low mtDNA divergence in these species. Combined with the demographic inference, we deduced that glacial-interglacial cycles drove geographic isolation and secondary contacts of these species, then promoted hybridization and lineage fusion among them, finally resulting in a reticulate evolutionary pattern. Overall, our study highlights the importance of combining multiple types of data to delimit species, especially those with conserved morphology, and to reveal the sophisticated processes of speciation.

Historical demography of four gecko species specializing in boulder cave habitat: Implications in the evolutionary dead end hypothesis and conservation

Specialization in narrow ecological niches may not only help species to survive in competitive or unique environments but also contribute to their extermination over evolutionary time. Although the "evolutionary dead end" hypothesis has long been debated, empirical evidence from species with detailed information on niche specialization and evolutionary history remains rare. Here we use a group of four closely related Cnemaspis gecko species that depend highly on granite boulder caves in the Mekong Delta to investigate the potential impact of ecological specialization on their evolution and population dynamics. Isolated by unsuitable floodplain habitats, these boulder-dwelling geckos are among the most narrowly distributed Squamata in the world. We applied several coalescence-based approaches combined with the RAD-seq technique to estimate their divergence times, gene flow and demographic fluctuations during the speciation and population differentiation processes. Our results reveal long-term population shrinkage in the four geckos and limited gene flow during their divergence. The results suggest that the erosion and fragmentation of the granite boulder hills have greatly impacted population divergence and declines. The habitat specialization of these geckos has led to fine-scaled speciation in these granite rocky hills; in contrast, specialization might also have pushed these species toward the edge of extinction. Our study also emphasizes the conservation urgency of these vulnerable, cave-dependent geckos.

Fish diversification at the pace of geomorphological changes: evolutionary history of western Iberian Leuciscinae (Teleostei: Leuciscidae) inferred from multilocus sequence data

The evolutionary history of western Iberian Leuciscinae, obligatory freshwater fish, is directly linked to the evolution of the hydrographic network of the Iberian Peninsula after its isolation from the rest of Europe, which involved dramatic rearrangements such as the transition from endorheic lakes to open basins draining to the Atlantic. Previous phylogenetic research on western Iberian leuciscines, using mainly mitochondrial DNA and more recently one or two nuclear genes, has found contradictory results and there remain many unresolved issues regarding species relationships, taxonomy, and evolutionary history. Moreover, there is a lack of integration between phylogenetic and divergence time estimates and information on the timing of geomorphological changes and paleobasin rearrangements in the Iberian Peninsula. This study presents the first comprehensive fossil-calibrated multilocus coalescent species tree of western Iberian Leuciscinae (including 14 species of Achondrostoma, Iberochondrostoma, Pseudochondrostoma and Squalius endemic to the Iberian Peninsula, seven of which endemic to Portugal) based on seven nuclear genes, and from which we infer their biogeographic history by comparing divergence time estimates to known dated geological events. The phylogenetic pattern suggests slow-paced evolution of leuciscines during the Early-Middle Miocene endorheic phase of the main Iberian river basins, with the shift to exorheism in the late Neogene-Quaternary allowing westward dispersals that resulted in many cladogenetic events and a high rate of endemism in western Iberia. The results of this study also: (i) confirm the paraphyly of S. pyrenaicus with respect to S. carolitertii, and thus the possible presence of a new taxon in the Portuguese Tagus currently assigned to S. pyrenaicus; (ii) support the taxonomic separation of the Guadiana and Sado populations of S. pyrenaicus; (iii) show the need for further population sampling and taxonomic research to clarify the phylogenetic status of A. arcasii from the Minho basin and of the I. lusitanicum populations in the Sado and Tagus basins; and (iv) indicate that A. occidentale, I. olisiponensis and P. duriensis are the most ancient lineages within their respective genera.

Uncovering hidden specific diversity of Andean glassfrogs of the Centrolene buckleyi species complex (Anura: Centrolenidae)

The glassfrog Centrolene buckleyi has been recognized as a species complex. Herein, using coalescence-based species delimitation methods, we evaluate the specific diversity within this taxon. Four coalescence approaches (generalized mixed Yule coalescents, Bayesian general mixed Yule-coalescent, Poisson tree processes, and Bayesian Poisson tree processes) were consistent with the delimitation results, identifying four lineages within what is currently recognized as C. buckleyi. We propose three new candidate species that should be tested with nuclear markers, morphological, and behavioral data. In the meantime, for conservation purposes, candidate species should be considered evolutionary significant units, in light of observed population crashes in the C. buckleyi species complex. Finally, our results support the validity of C. venezuelense, formerly considered as a subspecies of C. buckleyi.

Hidden diversity within the depauperate genera of the snake tribe Lampropeltini (Serpentes, Colubridae)

Accurate representation of lineage diversity through complete taxon sampling is crucial to understanding the evolution of biodiversity, particularly when using molecular phylogenetics to estimate evolutionary relationships. In this interest, taxonomic diversity is often used as a proxy for lineage diversity even though the two concepts are not synonymous. We explore this within the snake tribe Lampropeltini which includes some of the most conspicuous and heavily studied snakes in North America. Both the taxonomy and hypothesized relationships within this tribe have been in flux. The number of species has increased from 23 to 51 over the last thirty years, predominately within three of the nine genera (Lampropeltis, Pantherophis, Pituophis). The remaining six depauperate genera (Arizona, Bogertophis, Cemophora, Pseudelaphe, Rhinocheilus, and Senticolis) have been poorly represented in phylogenetic studies. To estimate evolutionary relationships and determine if the dichotomy in depauperate and speciose genera within Lampropeltini is a function of taxon sampling or truly represents the lineage diversity, we estimated the phylogeny of this group using nuclear and mitochondrial loci in a concatenated and coalescent framework with the largest sampling of the six depauperate genera to date. In addition, we estimated the divergence dates among the genera to assess whether the instability of Lampropeltini phylogenetic relationships is due to an adaptive radiation. While some nodes still remain unresolved, the generic-level relationships we recovered agree with those of a recent next-generation study that used a much larger set of loci for fewer individuals. We also tested two putative species, Arizona pacata and Pseudelaphe phaescens, for the first time phylogenetically and find evidence that they are distinct lineages. Overall, we find that the taxonomic and genetic diversity are not correlated in Lampropeltini and that representing putative diversity in phylogenies will lead to a better estimate of evolutionary histories, especially in groups with complex radiations.

Phylogenomic analysis on the exceptionally diverse fish clade Gobioidei (Actinopterygii: Gobiiformes) and data-filtering based on molecular clocklikeness

The use of genome-scale data to infer phylogenetic relationships has gained in popularity in recent years due to the progress made in target-gene capture and sequencing techniques. Data filtering, the approach of excluding data inconsistent with the model from analyses, presumably could alleviate problems caused by systematic errors in phylogenetic inference. Different data filtering criteria, such as those based on evolutionary rate and molecular clocklikeness as well as others have been proposed for selecting useful phylogenetic markers, yet few studies have tested these criteria using phylogenomic data. We developed a novel set of single-copy nuclear coding markers to capture thousands of target genes in gobioid fishes, a species-rich lineages of vertebrates, and tested the effects of data-filtering methods based on substitution rate and molecular clocklikeness while attempting to control for the compounding effects of missing data and variation in locus length. We found that molecular clocklikeness was a better predictor than overall substitution rate for phylogenetic usefulness of molecular markers in our study. In addition, when the 100 best ranked loci for our predictors were concatenated and analyzed using maximum likelihood, or combined in a coalescent-based species-tree analysis, the resulting trees showed a well-resolved topology of Gobioidei that mostly agrees with previous studies. However, trees generated from the 100 least clocklike frequently recovered conflicting, and in some cases clearly erroneous topologies with strong support, thus indicating strong systematic biases in those datasets. Collectively these results suggest that data filtering has the potential improve the performance of phylogenetic inference when using both a concatenation approach as well as methods that rely on input from individual gene trees (i.e. coalescent species-tree approaches), which may be preferred in scenarios where incomplete lineage sorting is likely to be an issue.

Phylogenomic analysis of the Chilean clade of Liolaemus lizards (Squamata: Liolaemidae) based on sequence capture data

The genus Liolaemus is one of the most ecologically diverse and species-rich genera of lizards worldwide. It currently includes more than 250 recognized species, which have been subject to many ecological and evolutionary studies. Nevertheless, Liolaemus lizards have a complex taxonomic history, mainly due to the incongruence between morphological and genetic data, incomplete taxon sampling, incomplete lineage sorting and hybridization. In addition, as many species have restricted and remote distributions, this has hampered their examination and inclusion in molecular systematic studies. The aims of this study are to infer a robust phylogeny for a subsample of lizards representing the Chilean clade (subgenus Liolaemus sensu stricto), and to test the monophyly of several of the major species groups. We use a phylogenomic approach, targeting 541 ultra-conserved elements (UCEs) and 44 protein-coding genes for 16 taxa. We conduct a comparison of phylogenetic analyses using maximum-likelihood and several species tree inference methods. The UCEs provide stronger support for phylogenetic relationships compared to the protein-coding genes; however, the UCEs outnumber the protein-coding genes by 10-fold. On average, the protein-coding genes contain over twice the number of informative sites. Based on our phylogenomic analyses, all the groups sampled are polyphyletic. Liolaemus tenuis tenuis is difficult to place in the phylogeny, because only a few loci (nine) were recovered for this species. Topologies or support values did not change dramatically upon exclusion of L. t. tenuis from analyses, suggesting that missing data did not had a significant impact on phylogenetic inference in this data set. The phylogenomic analyses provide strong support for sister group relationships between L. fuscus, L. monticola, L. nigroviridis and L. nitidus, and L. platei and L. velosoi. Despite our limited taxon sampling, we have provided a reliable starting hypothesis for the relationships among many major groups of the Chilean clade of Liolaemus that will help future work aimed at resolving the Liolaemus phylogeny.

Gene flow during glacial habitat shifts facilitates character displacement in a Neotropical flycatcher radiation

BACKGROUND

Pleistocene climatic fluctuations are known to be an engine of biotic diversification at higher latitudes, but their impact on highly diverse tropical areas such as the Andes remains less well-documented. Specifically, while periods of global cooling may have led to fragmentation and differentiation at colder latitudes, they may - at the same time - have led to connectivity among insular patches of montane tropical habitat with unknown consequences on diversification. In the present study we utilized ~5.5 kb of DNA sequence data from eight nuclear loci and one mitochondrial gene alongside diagnostic morphological and bioacoustic markers to test the effects of Pleistocene climatic fluctuations on diversification in a complex of Andean tyrant-flycatchers of the genus Elaenia.

RESULTS

Population genetic and phylogenetic approaches coupled with coalescent simulations demonstrated disparate levels of gene flow between the taxon chilensis and two parapatric Elaenia taxa predominantly during the last glacial period but not thereafter, possibly on account of downward shifts of montane forest habitat linking the populations of adjacent ridges. Additionally, morphological and bioacoustic analyses revealed a distinct pattern of character displacement in coloration and vocal traits between the two sympatric taxa albiceps and pallatangae, which were characterized by a lack of gene flow.

CONCLUSION

Our study demonstrates that global periods of cooling are likely to have facilitated gene flow among Andean montane Elaenia flycatchers that are more isolated from one another during warm interglacial periods such as the present era. We also identify a hitherto overlooked case of plumage and vocal character displacement, underpinning the complexities of gene flow patterns caused by Pleistocene climate change across the Andes.

Integrative species delimitation of the widespread North American jumping mice (Zapodinae)

Delimiting species can be challenging, but is a key step for the critical examination of evolutionary history and for prioritizing conservation efforts. Because systematic relationships are often determined iteratively using tests based on taxonomy, such methods can fail to detect cryptic variation and result in biased conclusions. Conversely, discovery-based approaches provide a powerful way to define operational taxonomic units and test species boundaries. We compare both approaches (taxonomy-based delimitation - TBD and discovery-based delimitation - DBD) within North American jumping mice (Zapodinae) using broad sampling, multilocus analyses, and ecological tests. This group diversified through the dynamic glacial-interglacial periods of the Quaternary and phylogeographic tests reveal 28 lineages that correspond poorly with current taxonomy (4 species, 32 nominal subspecies). However, neither the 4-species or 28-lineage hypotheses are optimal for species-level classification. Rather, information theoretic approaches (Bayes Factors) indicate a 15-species hypothesis is best for characterizing genetic variation in this group, with subsequent iterative pairwise ecological tests failing to confirm four species pairs. Taken together, evolutionary and ecological tests capture divergence among 11 putative species that, if upheld by additional tests, will lead to taxonomic revision and reevaluation of conservation plans.

The birth of aposematism: High phenotypic divergence and low genetic diversity in a young clade of poison frogs

Rapid radiation coupled with low genetic divergence often hinders species delimitation and phylogeny estimation even if putative species are phenotypically distinct. Some aposematic species, such as poison frogs (Dendrobatidae), have high levels of intraspecific color polymorphism, which can lead to overestimation of species when phenotypic divergence primarily guides species delimitation. We explored this possibility in the youngest origin of aposematism (3-7 MYA) in poison frogs, Epipedobates, by comparing genetic divergence with color and acoustic divergence. We found low genetic divergence (2.6% in the 16S gene) despite substantial differences in color and acoustic signals. While chemical defense is inferred to have evolved in the ancestor of Epipedobates, aposematic coloration evolved at least twice or was lost once in Epipedobates, suggesting that it is evolutionarily labile. We inferred at least one event of introgression between two cryptically colored species with adjacent ranges (E. boulengeri and E. machalilla). We also find evidence for peripheral isolation resulting in phenotypic divergence and potential speciation of the aposematic E. tricolor from the non-aposematic E. machalilla. However, we were unable to estimate a well-supported species tree or delimit species using multispecies coalescent models. We attribute this failure to factors associated with recent speciation including mitochondrial introgression, incomplete lineage sorting, and too few informative molecular characters. We suggest that species delimitation within young aposematic lineages such as Epipedobates will require genome-level molecular studies. We caution against relying solely on DNA barcoding for species delimitation or identification and highlight the value of phenotypic divergence and natural history in delimiting species.

Phylogeography of the heavily poached African common pangolin (Pholidota, Manis tricuspis) reveals six cryptic lineages as traceable signatures of Pleistocene diversification

Knowledge on faunal diversification in African rainforests remains scarce. We used phylogeography to assess (i) the role of Pleistocene climatic oscillations in the diversification of the African common pangolin (Manis tricuspis) and (ii) the utility of our multilocus approach for taxonomic delineation and trade tracing of this heavily poached species. We sequenced 101 individuals for two mitochondrial DNA (mtDNA), two nuclear DNA and one Y-borne gene fragments (totalizing 2602 bp). We used a time-calibrated, Bayesian inference phylogenetic framework and conducted character-based, genetic and phylogenetic delineation of species hypotheses within African common pangolins. We identified six geographic lineages partitioned into western Africa, Ghana, the Dahomey Gap, western central Africa, Gabon and central Africa, all diverging during the Middle to Late Pleistocene. MtDNA (cytochrome b + control region) was the sole locus to provide diagnostic characters for each of the six lineages. Tree-based Bayesian delimitation methods using single- and multilocus approaches gave high support for 'species' level recognition of the six African common pangolin lineages. Although the diversification of African common pangolins occurred during Pleistocene cyclical glaciations, causative correlation with traditional rainforest refugia and riverine barriers in Africa was not straightforward. We conclude on the existence of six cryptic lineages within African common pangolins, which might be of major relevance for future conservation strategies. The high discriminative power of the mtDNA markers used in this study should allow an efficient molecular tracing of the regional origin of African common pangolin seizures.

Resolving evolutionary relationships in lichen-forming fungi using diverse phylogenomic datasets and analytical approaches

Evolutionary histories are now being inferred from unprecedented, genome-scale datasets for a broad range of organismal groups. While phylogenomic data has helped in resolving a number of difficult, long-standing questions, constructing appropriate datasets from genomes is not straightforward, particularly in non-model groups. Here we explore the utility of phylogenomic data to infer robust phylogenies for a lineage of closely related lichen-forming fungal species. We assembled multiple, distinct nuclear phylogenomic datasets, ranging from ca. 25 Kb to 16.8 Mb and inferred topologies using both concatenated gene tree approaches and species tree methods based on the multispecies coalescent model. In spite of evidence for rampant incongruence among individual loci, these genome-scale datasets provide a consistent, well-supported phylogenetic hypothesis using both concatenation and multispecies coalescent approaches (ASTRAL-II and SVDquartets). However, the popular full hierarchical coalescent approach implemented in *BEAST provided inconsistent inferences, both in terms of nodal support and topology, with smaller subsets of the phylogenomic data. While comparable, well-supported topologies can be accurately inferred with only a small fraction of the overall genome, consistent results across a variety of datasets and methodological approaches provide reassurance that phylogenomic data can effectively be used to provide robust phylogenies for closely related lichen-forming fungal lineages.

Computational Performance and Statistical Accuracy of *BEAST and Comparisons with Other Methods

Under the multispecies coalescent model of molecular evolution, gene trees have independent evolutionary histories within a shared species tree. In comparison, supermatrix concatenation methods assume that gene trees share a single common genealogical history, thereby equating gene coalescence with species divergence. The multispecies coalescent is supported by previous studies which found that its predicted distributions fit empirical data, and that concatenation is not a consistent estimator of the species tree. *BEAST, a fully Bayesian implementation of the multispecies coalescent, is popular but computationally intensive, so the increasing size of phylogenetic data sets is both a computational challenge and an opportunity for better systematics. Using simulation studies, we characterize the scaling behavior of *BEAST, and enable quantitative prediction of the impact increasing the number of loci has on both computational performance and statistical accuracy. Follow-up simulations over a wide range of parameters show that the statistical performance of *BEAST relative to concatenation improves both as branch length is reduced and as the number of loci is increased. Finally, using simulations based on estimated parameters from two phylogenomic data sets, we compare the performance of a range of species tree and concatenation methods to show that using *BEAST with tens of loci can be preferable to using concatenation with thousands of loci. Our results provide insight into the practicalities of Bayesian species tree estimation, the number of loci required to obtain a given level of accuracy and the situations in which supermatrix or summary methods will be outperformed by the fully Bayesian multispecies coalescent.

Comparing species tree estimation with large anchored phylogenomic and small Sanger-sequenced molecular datasets: an empirical study on Malagasy pseudoxyrhophiine snakes

BACKGROUND

Using molecular data generated by high throughput next generation sequencing (NGS) platforms to infer phylogeny is becoming common as costs go down and the ability to capture loci from across the genome goes up. While there is a general consensus that greater numbers of independent loci should result in more robust phylogenetic estimates, few studies have compared phylogenies resulting from smaller datasets for commonly used genetic markers with the large datasets captured using NGS. Here, we determine how a 5-locus Sanger dataset compares with a 377-locus anchored genomics dataset for understanding the evolutionary history of the pseudoxyrhophiine snake radiation centered in Madagascar. The Pseudoxyrhophiinae comprise ~86 % of Madagascar's serpent diversity, yet they are poorly known with respect to ecology, behavior, and systematics. Using the 377-locus NGS dataset and the summary statistics species-tree methods STAR and MP-EST, we estimated a well-supported species tree that provides new insights concerning intergeneric relationships for the pseudoxyrhophiines. We also compared how these and other methods performed with respect to estimating tree topology using datasets with varying numbers of loci.

METHODS

Using Sanger sequencing and an anchored phylogenomics approach, we sequenced datasets comprised of 5 and 377 loci, respectively, for 23 pseudoxyrhophiine taxa. For each dataset, we estimated phylogenies using both gene-tree (concatenation) and species-tree (STAR, MP-EST) approaches. We determined the similarity of resulting tree topologies from the different datasets using Robinson-Foulds distances. In addition, we examined how subsets of these data performed compared to the complete Sanger and anchored datasets for phylogenetic accuracy using the same tree inference methodologies, as well as the program *BEAST to determine if a full coalescent model for species tree estimation could generate robust results with fewer loci compared to the summary statistics species tree approaches. We also examined the individual gene trees in comparison to the 377-locus species tree using the program MetaTree.

RESULTS

Using the full anchored dataset under a variety of methods gave us the same, well-supported phylogeny for pseudoxyrhophiines. The African pseudoxyrhophiine Duberria is the sister taxon to the Malagasy pseudoxyrhophiines genera, providing evidence for a monophyletic radiation in Madagascar. In addition, within Madagascar, the two major clades inferred correspond largely to the aglyphous and opisthoglyphous genera, suggesting that feeding specializations associated with tooth venom delivery may have played a major role in the early diversification of this radiation. The comparison of tree topologies from the concatenated and species-tree methods using different datasets indicated the 5-locus dataset cannot beused to infer a correct phylogeny for the pseudoxyrhophiines under any method tested here and that summary statistics methods require 50 or more loci to consistently recover the species-tree inferred using the complete anchored dataset. However, as few as 15 loci may infer the correct topology when using the full coalescent species tree method *BEAST. MetaTree analyses of each gene tree from the Sanger and anchored datasets found that none of the individual gene trees matched the 377-locus species tree, and that no gene trees were identical with respect to topology.

CONCLUSIONS

Our results suggest that ≥50 loci may be necessary to confidently infer phylogenies when using summaryspecies-tree methods, but that the coalescent-based method *BEAST consistently recovers the same topology using only 15 loci. These results reinforce that datasets with small numbers of markers may result in misleading topologies, and further, that the method of inference used to generate a phylogeny also has a major influence on the number of loci necessary to infer robust species trees.

Extreme food-plant specialisation in Megabombus bumblebees as a product of long tongues combined with short nesting seasons

Megabombus bumblebees have unusually long tongues and are generally more specialised than other bumblebees in their choice of food plants. The phylogeny of Megabombus bumblebees shows that speciation was concentrated in two periods. Speciation in the first period (ca 4.25–1.5 Ma) is associated with the late rise of the Hengduan Mountains at the eastern end of the Qinghai-Tibetan plateau. Speciation in the second period (1.2–0.3 Ma) is associated with climatic cooling in the northern forests. The most extreme food-specialist species belong to the second period, which may point to climate as a factor in specialisation. These extreme specialist species occur either in the far north (Bombus consobrinus), or at high elevations (Bombus gerstaeckeri), in situations where long tongues coincide with the shortest nesting seasons. Species with the longest tongues but occurring further south (even at high elevations) use a broader range of food plants.

Assessing species boundaries using multilocus species delimitation in a morphologically conserved group of neotropical freshwater fishes, the Poecilia sphenops species complex (Poeciliidae)

Accurately delimiting species is fundamentally important for understanding species diversity and distributions and devising effective strategies to conserve biodiversity. However, species delimitation is problematic in many taxa, including 'non-adaptive radiations' containing morphologically cryptic lineages. Fortunately, coalescent-based species delimitation methods hold promise for objectively estimating species limits in such radiations, using multilocus genetic data. Using coalescent-based approaches, we delimit species and infer evolutionary relationships in a morphologically conserved group of Central American freshwater fishes, the Poecilia sphenops species complex. Phylogenetic analyses of multiple genetic markers (sequences of two mitochondrial DNA genes and five nuclear loci) from 10/15 species and genetic lineages recognized in the group support the P. sphenops species complex as monophyletic with respect to outgroups, with eight mitochondrial 'major-lineages' diverged by ≥2% pairwise genetic distances. From general mixed Yule-coalescent models, we discovered (conservatively) 10 species within our concatenated mitochondrial DNA dataset, 9 of which were strongly supported by subsequent multilocus Bayesian species delimitation and species tree analyses. Results suggested species-level diversity is underestimated or overestimated by at least ~15% in different lineages in the complex. Nonparametric statistics and coalescent simulations indicate genealogical discordance among our gene tree results has mainly derived from interspecific hybridization in the nuclear genome. However, mitochondrial DNA show little evidence for introgression, and our species delimitation results appear robust to effects of this process. Overall, our findings support the utility of combining multiple lines of genetic evidence and broad phylogeographical sampling to discover and validate species using coalescent-based methods. Our study also highlights the importance of testing for hybridization versus incomplete lineage sorting, which aids inference of not only species limits but also evolutionary processes influencing genetic diversity.

A tribal level phylogeny of Lake Tanganyika cichlid fishes based on a genomic multi-marker approach

The species-flocks of cichlid fishes in the East African Great Lakes Victoria, Malawi and Tanganyika constitute the most diverse extant adaptive radiations in vertebrates. Lake Tanganyika, the oldest of the lakes, harbors the morphologically and genetically most diverse assemblage of cichlids and contains the highest number of endemic cichlid genera of all African lakes. Based on morphological grounds, the Tanganyikan cichlid species have been grouped into 12-16 distinct lineages, so-called tribes. While the monophyly of most of the tribes is well established, the phylogenetic relationships among the tribes remain largely elusive. Here, we present a new tribal level phylogenetic hypothesis for the cichlid fishes of Lake Tanganyika that is based on the so far largest set of nuclear markers and a total alignment length of close to 18kb. Using next-generation amplicon sequencing with the 454 pyrosequencing technology, we compiled a dataset consisting of 42 nuclear loci in 45 East African cichlid species, which we subjected to maximum likelihood and Bayesian inference phylogenetic analyses. We analyzed the entire concatenated dataset and each marker individually, and performed a Bayesian concordance analysis and gene tree discordance tests. Overall, we find strong support for a position of the Oreochromini, Boulengerochromini, Bathybatini and Trematocarini outside of a clade combining the substrate spawning Lamprologini and the mouthbrooding tribes of the 'H-lineage', which are both strongly supported to be monophyletic. The Eretmodini are firmly placed within the 'H-lineage', as sister-group to the most species-rich tribe of cichlids, the Haplochromini. The phylogenetic relationships at the base of the 'H-lineage' received less support, which is likely due to high speciation rates in the early phase of the radiation. Discordance among gene trees and marker sets further suggests the occurrence of past hybridization and/or incomplete lineage sorting in the cichlid fishes of Lake Tanganyika.

When do species-tree and concatenated estimates disagree? An empirical analysis with higher-level scincid lizard phylogeny

Simulation studies suggest that coalescent-based species-tree methods are generally more accurate than concatenated analyses. However, these species-tree methods remain impractical for many large datasets. Thus, a critical but unresolved issue is when and why concatenated and coalescent species-tree estimates will differ. We predict such differences for branches in concatenated trees that are short, weakly supported, and have conflicting gene trees. We test these predictions in Scincidae, the largest lizard family, with data from 10 nuclear genes for 17 ingroup taxa and 44 genes for 12 taxa. We support our initial predictions, andsuggest that simply considering uncertainty in concatenated trees may sometimes encompass the differences between these methods. We also found that relaxed-clock concatenated trees can be surprisingly similar to the species-tree estimate. Remarkably, the coalescent species-tree estimates had slightly lower support values when based on many more genes (44 vs. 10) and a small (∼30%) reduction in taxon sampling. Thus, taxon sampling may be more important than gene sampling when applying species-tree methods to deep phylogenetic questions. Finally, our coalescent species-tree estimates tentatively support division of Scincidae into three monophyletic subfamilies, a result otherwise found only in concatenated analyses with extensive species sampling.

A time-calibrated, multi-locus phylogeny of piranhas and pacus (Characiformes: Serrasalmidae) and a comparison of species tree methods

The phylogeny of piranhas, pacus, and relatives (family Serrasalmidae) was inferred on the basis of DNA sequences from eleven gene fragments that include the mitochondrial control region plus 10 nuclear genes (two exons and eight introns). The new data were obtained for a representative sampling of 53 specimens, collected from all major South American rivers, accounting for over 40% of the valid species and all genera excluding Utiaritichthys. Two fossil calibration points and relaxed-clock Bayesian analyses were used to estimate the timing of diversification. The new multilocus dataset also is used to compare several species-tree approaches against the results obtained using the concatenated alignment analyzed under maximum likelihood and Bayesian inference. Individual gene trees showed substantial topological discordance, but analyses based on concatenation and Bayesian and maximum likelihood-based species trees approaches converged onto a single phylogeny. The resulting phylogenetic hypothesis is robust and supports a division of the family into three major clades, consistent with previous results based on mitochondrial DNA alone. The earliest branching event separated a "pacu" clade (Colossoma, Mylossoma and Piaractus) from the rest of the family in the Late Cretaceous (over 68 Ma). The other two clades, that contain most of the diversity, are formed by the "true piranhas" (Metynnis, Pygopristis, Pygocentrus, Pristobrycon, Catoprion, and Serrasalmus) and the Myleus-like pacus (the Myleus clade). The "true" piranha clade originated during the Eocene (∼53 Ma) but the most recent diversification of flesh-eating piranhas within the genera Serrasalmus and Pygocentrus did not start until the Miocene (∼17 Ma). A comparison of species tree approaches indicates that most methods tested are consistent with results obtained by concatenation, suggesting that the gene-tree incongruence observed is mild and will not produce misleading results under simple concatenation analysis. Non-monophyly of several genera (Pristobrycon, Tometes, Myloplus, Mylesinus) and putative species (Serrasalmus rhombeus) was obtained, suggesting that further study of this family is necessary.

The evolutionary history of Darwin's finches: speciation, gene flow, and introgression in a fragmented landscape

Many classic examples of adaptive radiations take place within fragmented systems such as islands or mountains, but the roles of mosaic landscapes and variable gene flow in facilitating species diversification is poorly understood. Here we combine phylogenetic and landscape genetic approaches to understand diversification in Darwin's finches, a model adaptive radiation. We combined sequence data from 14 nuclear introns, mitochondrial markers, and microsatellite variation from 51 populations of all 15 recognized species. Phylogenetic species-trees recovered seven major finch clades: ground, tree, vegetarian, Cocos Island, grey and green warbler finches, and a distinct clade of sharp-beaked ground finches (Geospiza cf. difficilis) basal to all ground and tree finches. The ground and tree finch clades lack species-level phylogenetic structure. Interisland gene flow and interspecies introgression vary geographically in predictable ways. First, several species exhibit concordant patterns of population divergence across the channel separating the Galápagos platform islands from the separate volcanic province of northern islands. Second, peripheral islands have more admixed populations while central islands maintain more distinct species boundaries. This landscape perspective highlights a likely role for isolation of peripheral populations in initial divergence, and demonstrates that peripheral populations may maintain genetic diversity through outbreeding during the initial stages of speciation.

How do you solve a problem like Letharia? A new look at cryptic species in lichen-forming fungi using Bayesian clustering and SNPs from multilocus sequence data

The inclusion of molecular data is increasingly an integral part of studies assessing species boundaries. Analyses based on predefined groups may obscure patterns of differentiation, and population assignment tests provide an alternative for identifying population structure and barriers to gene flow. In this study, we apply population assignment tests implemented in the programs STRUCTURE and BAPS to single nucleotide polymorphisms from DNA sequence data generated for three previous studies of the lichenized fungal genus Letharia. Previous molecular work employing a gene genealogical approach circumscribed six species-level lineages within the genus, four putative lineages within the nominal taxon L. columbiana (Nutt.) J.W. Thomson and two sorediate lineages. We show that Bayesian clustering implemented in the program STRUCTURE was generally able to recover the same six putative Letharia lineages. Population assignments were largely consistent across a range of scenarios, including: extensive amounts of missing data, the exclusion of SNPs from variable markers, and inferences based on SNPs from as few as three gene regions. While our study provided additional evidence corroborating the six candidate Letharia species, the equivalence of these genetic clusters with species-level lineages is uncertain due, in part, to limited phylogenetic signal. Furthermore, both the BAPS analysis and the ad hoc ΔK statistic from results of the STRUCTURE analysis suggest that population structure can possibly be captured with fewer genetic groups. Our findings also suggest that uneven sampling across taxa may be responsible for the contrasting inferences of population substructure. Our results consistently supported two distinct sorediate groups, ‘L. lupina’ and L. vulpina, and subtle morphological differences support this distinction. Similarly, the putative apotheciate species ‘L. lucida’ was also consistently supported as a distinct genetic cluster. However, additional studies will be required to elucidate the relationships of other L. columbiana s.l. populations with the two sorediate genetic clusters.

An empirical evaluation of two-stage species tree inference strategies using a multilocus dataset from North American pines

Background

As it becomes increasingly possible to obtain DNA sequences of orthologous genes from diverse sets of taxa, species trees are frequently being inferred from multilocus data. However, the behavior of many methods for performing this inference has remained largely unexplored. Some methods have been proven to be consistent given certain evolutionary models, whereas others rely on criteria that, although appropriate for many parameter values, have peculiar zones of the parameter space in which they fail to converge on the correct estimate as data sets increase in size.

Results

Here, using North American pines, we empirically evaluate the behavior of 24 strategies for species tree inference using three alternative outgroups (72 strategies total). The data consist of 120 individuals sampled in eight ingroup species from subsection Strobus and three outgroup species from subsection Gerardianae, spanning ∼47 kilobases of sequence at 121 loci. Each “strategy” for inferring species trees consists of three features: a species tree construction method, a gene tree inference method, and a choice of outgroup. We use multivariate analysis techniques such as principal components analysis and hierarchical clustering to identify tree characteristics that are robustly observed across strategies, as well as to identify groups of strategies that produce trees with similar features. We find that strategies that construct species trees using only topological information cluster together and that strategies that use additional non-topological information (e.g., branch lengths) also cluster together. Strategies that utilize more than one individual within a species to infer gene trees tend to produce estimates of species trees that contain clades present in trees estimated by other strategies. Strategies that use the minimize-deep-coalescences criterion to construct species trees tend to produce species tree estimates that contain clades that are not present in trees estimated by the Concatenation, RTC, SMRT, STAR, and STEAC methods, and that in general are more balanced than those inferred by these other strategies.

Conclusions

When constructing a species tree from a multilocus set of sequences, our observations provide a basis for interpreting differences in species tree estimates obtained via different approaches that have a two-stage structure in common, one step for gene tree estimation and a second step for species tree estimation. The methods explored here employ a number of distinct features of the data, and our analysis suggests that recovery of the same results from multiple methods that tend to differ in their patterns of inference can be a valuable tool for obtaining reliable estimates.