Gene tree quality affects empirical coalescent branch length estimation

Assessing effects of gene tree error in coalescent analyses have widely ignored coalescent branch lengths (CBLs) despite their potential utility in estimating ancestral population demographics and detecting species tree anomaly zones. However, the ability of coalescent methods to obtain accurate estimates remains largely unexplored. Errors in gene trees should lead to underestimates of the true CBL, and for a given set of comparisons, longer CBLs should be more accurate. Here, we furthered our empirical understanding of how error in gene tree quality (i.e., locus informativeness and gene tree resolution) affect CBLs using four datasets comprised of ultraconserved elements (UCE) or exons for clades that exhibit wide ranges of branch lengths. For each dataset, we compared the impact of locus informativeness (assessed using number of parsimony-informative sites) and gene tree resolution on CBL estimates. Our results, in general, showed that CBLs were drastically shorter when estimates included low informative loci. Gene tree resolution also had an impact on UCE datasets, with polytomous gene trees producing longer branches than randomly resolved gene trees. However, resolution did not appear to affect CBL estimates from the more informative exon datasets. Thus, as expected, gene tree quality affects CBL estimates, though this can generally be minimized by using moderate filtering to select more informative loci and/or by allowing polytomies in gene trees. These approaches, as well as additional contributions to improve CBL estimation, should lead to CBLs that are useful for addressing evolutionary and biological questions.     

The impact of anchored phylogenomics and taxon sampling on phylogenetic inference in narrow‐mouthed frogs (Anura,Microhylidae)

Despite considerable progress in unravelling the phylogenetic relationships of microhylid frogs, relationships among subfamilies remain largely unstable and many genera are not demonstrably monophyletic. Here, we used five alternative combinations of DNA sequence data (ranging from seven loci for 48 taxa to up to 73 loci for as many as 142 taxa) generated using the anchored phylogenomics sequencing method (66 loci, derived from conserved genome regions, for 48 taxa) and Sanger sequencing (seven loci for up to 142 taxa) to tackle this problem. We assess the effects of character sampling, taxon sampling, analytical methods and assumptions in phylogenetic inference of microhylid frogs. The phylogeny of microhylids shows high susceptibility to different analytical methods and datasets used for the analyses. Clades inferred from maximum‐likelihood are generally more stable across datasets than those inferred from parsimony. Parsimony trees inferred within a tree‐alignment framework are generally better resolved and better supported than those inferred within a similarity‐alignment framework, even under the same cost matrix (equally weighted) and same treatment of gaps (as a fifth nucleotide state). We discuss potential causes for these differences in resolution and clade stability among discovery operations. We also highlight the problem that commonly used algorithms for model‐based analyses do not explicitly model insertion and deletion events (i.e. gaps are treated as missing data). Our results corroborate the monophyly of Microhylidae and most currently recognized subfamilies but fail to provide support for relationships among subfamilies. Several taxonomic updates are provided, including naming of two new subfamilies, both monotypic.     

Phylogenomic analysis of the beetle suborder Adephaga with comparison of tailored and generalized ultraconserved element probe performance

Adephaga is the second largest suborder of beetles (Coleoptera) and they serve as important arthropod predators in both aquatic and terrestrial ecosystems. The suborder is divided into Geadephaga comprising terrestrial families and Hydradephaga for aquatic lineages. Despite numerous studies, phylogenetic relationships among the adephagan families and monophyly of the Hydradephaga itself remain in question. Here we conduct a comprehensive phylogenomic analysis of the suborder using ultraconserved elements (UCEs). This study presents the first in vitro test of a newly developed UCE probe set customized for use within Adephaga that includes both probes tailored specifically for the suborder, alongside generalized Coleoptera probes previously found to work in adephagan taxa. We assess the utility of the entire probe set, as well as comparing the tailored and generalized probes alone for reconstructing evolutionary relationships. Our analyses recovered strong support for the paraphyly of Hydradephaga with whirligig beetles (Gyrinidae) placed as sister to all other adephagan families. Geadephaga was strongly supported as monophyletic and placed sister to a clade composed of Haliplidae + Dytiscoidea. Monophyly of Dytiscoidea was strongly supported with relationships among the dytiscoid families resolved and strongly supported. Relationships among the subfamilies of Dytiscidae were strongly supported but largely incongruent with prior phylogenetic estimates for the family. The results of our UCE probe comparison showed that tailored probes alone outperformed generalized probes alone, as well as the full combined probe set (containing both types of probes), under decreased taxon sampling. When taxon sampling was increased, the full combined probe set outperformed both tailored probes and generalized probes alone. This study provides further evidence that UCE probe sets customized for a focal group result in a greater number of recovered loci and substantially improve phylogenomic analysis.     

Resolution of a concatenation/coalescence kerfuffle: partitioned coalescence support and a robust family‐level tree for Mammalia

Recent phylogenetic analyses of a large dataset for mammalian families (169 taxa, 26 loci) portray contrasting results. Supermatrix (concatenation) methods support a generally robust tree with only a few inconsistently resolved polytomies, whereas MP‐EST coalescence analysis of the same dataset yields a weakly supported tree that conflicts with many traditionally recognized clades. Here, we evaluate this discrepancy via improved coalescence analyses with reference to the rich history of phylogenetic studies on mammals. This integration clearly demonstrates that both supermatrix and coalescence analyses of just 26 loci yield a congruent, well‐supported phylogenetic hypothesis for Mammalia. Discrepancies between published studies are explained by implementation of overly simple DNA substitution models, inadequate tree‐search routines and limitations of the MP‐EST method. We develop a simple measure, partitioned coalescence support (PCS), which summarizes the distribution of support and conflict among gene trees for a given clade. Extremely high PCS scores for outlier gene trees at two nodes in the mammalian tree indicate a troubling bias in the MP‐EST method. We conclude that in this age of phylogenomics, a solid understanding of systematics fundamentals, choice of valid methodology and a broad knowledge of a clade's taxonomic history are still required to yield coherent phylogenetic inferences.     

Molecular and chemical characters to evaluate species status of two cuckoo bumblebees: Bombus barbutellus and Bombus maxillosus (Hymenoptera,Apidae, Bombini)

Many methods, based on morphological, molecular or chemical characters, have been used to address the question of species taxonomic status. Integrative taxonomy aims to define stronger supported taxonomic hypotheses by considering complementary datasets from different characters. By following an integrative approach, the present study includes molecular, chemical and morphological criteria to establish the taxonomic status of two rare and doubtful cuckoo bumblebee taxa: Bombus (Psithyrus) barbutellus and Bombus (Psithyrus) maxillosus. These two sympatric taxa are discriminated by few morphological criteria (mainly wing darkness and hair length). We used these morphological character diagnoses to establish an a priori status of our samples (23 specimens). We developed a combined molecular dataset from one nuclear gene, elongation factor 1α (EF‐1α), and one mitochondrial gene, cytochrome c oxidase subunit I (COI), spanning 1623 bp, and a chemical dataset of sexual marking pheromones (73 compounds). The molecular data were subjected to maximum‐likelihood and Bayesian phylogenetic inference under partitioned model and maximum parsimony. The chemical data were analysed by clustering and the two‐group k‐means method to test divergences between the two species. The resulting phylogenetic trees show no consistent divergence between the two taxa. Moreover, we found no divergence in the sexual marking pheromones in the clustering and two‐group k‐means analyses. These converging results support the conspecificity of both taxa. Nonetheless, our determinations using the traditional morphological criteria separated our samples into two taxa. We conclude that the morphological criteria seem to relate to intraspecific variations: B. maxillosus is regarded as a syn.n. of B. barbutellus.     

A phylogeny of megachiropteran bats (Mammalia: Chiroptera: Pteropodidae) based on direct optimization analysis of one nuclear and four mitochondrial genes

The phylogeny of megachiropteran bats (Mammalia: Chiroptera: Pteropodidae) has been investigated using several different molecular datasets. These studies differed widely in taxonomic and locus sampling, and their results tended to lack resolution of internal nodes and were themselves largely incongruent. To address this, we assembled a data set of 5 loci (up to 3.5 kbp from 12S rDNA, 16S rDNA, tDNA‐valine, cytochrome b, and the nuclear gene c‐mos) for 43 species of megachiropterans and 6 microchiropteran outgroups. We analyzed these data with direct optimization under equal costs for substitutions and indels. We used POY in a parallel setting, and searches consisted of replicated swapping + refinements (ratcheting, tree fusing, and iterative pass optimization). Our results indicate that Megachiroptera and all recognized genera (including Pteropus) are monophyletic, and that Melonycteris is the sister group of the clade containing all the other genera. Clades previously proposed using molecular data, as well as many new and traditional groups, were well‐supported, and various sources suggest that the degree of conflict with morphological data may be considerably less marked than previously supposed. Analysis of individual loci suffer 70% loss in the number of compatible groups recovered across all analyses with respect to combined analyses. Our results indicate that, within Megachiroptera, nectarivory and cave‐dwelling originated several times, but echolocation (used for obstacle detection) evolved only once. Megachiropterans likely originated in SE Asia‐Melanesia, and colonized Africa at least four times.     

Phylogenomic analysis suggests Coreidae and Alydidae (Hemiptera: Heteroptera) are not monophyletic

Next‐generation sequencing technologies (NGS) allow systematists to amass a wealth of genomic data from non‐model species for phylogenetic resolution at various temporal scales. However, phylogenetic inference for many lineages dominated by non‐model species has not yet benefited from NGS, which can complement Sanger sequencing studies. One such lineage, whose phylogenetic relationships remain uncertain, is the diverse, agriculturally important and charismatic Coreoidea (Hemiptera: Heteroptera). Given the lack of consensus on higher‐level relationships and the importance of a robust phylogeny for evolutionary hypothesis testing, we use a large data set comprised of hundreds of ultraconserved element (UCE) loci to infer the phylogeny of Coreoidea (excluding Stenocephalidae and Hyocephalidae), with emphasis on the families Coreidae and Alydidae. We generated three data sets by including alignments that contained loci sampled for at least 50%, 60%, or 70% of the total taxa, and inferred phylogeny using maximum likelihood and summary coalescent methods. Twenty‐six external morphological features used in relatively comprehensive phylogenetic analyses of coreoids were also re‐evaluated within our molecular phylogenetic framework. We recovered 439–970 loci per species (16%–36% of loci targeted) and combined this with previously generated UCE data for 12 taxa. All data sets, regardless of analytical approach, yielded topologically similar and strongly supported trees, with the exception of outgroup relationships and the position of Hydarinae. We recovered a monophyletic Coreoidea, with Rhopalidae highly supported as the sister group to Alydidae + Coreidae. Neither Alydidae nor Coreidae were monophyletic; the coreid subfamilies Hydarinae and Pseudophloeinae were recovered as more closely related to Alydidae than to other coreid subfamilies. Coreinae were paraphyletic with respect to Meropachyinae. Most morphological traits were homoplastic with several clades defined by few, if any, synapomorphies. Our results demonstrate the utility of phylogenomic approaches in generating robust hypotheses for taxa with long‐standing phylogenetic problems and highlight that novel insights may come from such approaches.     

Phylogenetic incongruence between nuclear and mitochondrial markers in the Asian colobines and the evolution of the langurs and leaf monkeys

Evidence of incongruence between mitochondrial and nuclear gene trees is now becoming documented with increasing frequency. Among the Old World monkeys, this discordance has been well demonstrated in the Cercopithecinae, but has not yet been investigated in the Colobinae. The mitochondrial relationships between the colobine genera have recently been clarified and cluster Presbytis and Trachypithecus as sister taxa to the exclusion of Semnopithecus. This is incongruent with previous morphological hypotheses that suggest the latter two are sister taxa, and perhaps even congeneric. In addition to analyzing a previously published 10,896 bp mitochondrial dataset, we sequenced and analyzed a 4297 bp fragment of the X-chromosome in order to test the competing mitochondrial and morphological phylogenetic hypotheses. The results from the mitochondrial dataset again support a Presbytis+Trachypithecus group while the X-chromosomal dataset supported a Semnopithecus+Trachypithecus group. A Shimodaira-Hasegawa test performed on both datasets indicates that the mitochondrial and X-chromosomal trees are significantly better at explaining their respective datasets than alternative topologies (p<0.05). We suggest that differential lineage sorting or ancient hybridization may be the cause of this strong discordance between the mitochondrial and X-chromosomal markers in these taxa.     

Ultraconserved elements show utility in phylogenetic inference of Adephaga (Coleoptera) and suggest paraphyly of ‘Hydradephaga’

The beetle suborder Adephaga has been the subject of many phylogenetic reconstructions utilizing a variety of data sources and inference methods. However, no strong consensus has yet emerged on the relationships among major adephagan lineages. Ultraconserved elements (UCEs) have proved useful for inferring difficult or unresolved phylogenies at varying timescales in vertebrates, arachnids and Hymenoptera. Recently, a UCE bait set was developed for Coleoptera using polyphagan genomes and a member of the order Strepsiptera as an outgroup. Here, we examine the utility of UCEs for reconstructing the phylogeny of adephagan families, in the first in vitro application a UCE bait set in Coleoptera. Our final dataset included 305 UCE loci for 18 representatives of all adephagan families except Aspidytidae, and two polyphagan outgroups, with a total concatenated length of 83 547 bp. We inferred trees using maximum likelihood analyses of the concatenated UCE alignment and coalescent species tree methods (astral ii , ASTRID, svdquartets ). Although the coalescent species tree methods had poor resolution and weak support, concatenated analyses produced well‐resolved, highly supported trees. Hydradephaga was recovered as paraphyletic, with Gyrinidae sister to Geadephaga and all other adephagans. Haliplidae was recovered as sister to Dytiscoidea, with Hygrobiidae and Amphizoidae successive sisters to Dytiscidae. Finally, Noteridae was recovered as monophyletic and sister to Meruidae. Given the success of UCE data for resolving phylogenetic relationships within Adephaga, we suggest the potential for further resolution of relationships within Adephaga using UCEs with improved taxon sampling, and by developing Adephaga‐specific probes.     

Squamate Conserved Loci (SqCL): A unified set of conserved loci for phylogenomics and population genetics of squamate reptiles

The identification of conserved loci across genomes, along with advances in target capture methods and high‐throughput sequencing, has helped spur a phylogenomics revolution by enabling researchers to gather large numbers of homologous loci across clades of interest with minimal upfront investment in locus design. Target capture for vertebrate animals is currently dominated by two approaches—anchored hybrid enrichment (AHE) and ultraconserved elements (UCE)—and both approaches have proven useful for addressing questions in phylogenomics, phylogeography and population genomics. However, these two sets of loci have minimal overlap with each other; moreover, they do not include many traditional loci that that have been used for phylogenetics. Here, we combine across UCE, AHE and traditional phylogenetic gene locus sets to generate the Squamate Conserved Loci set, a single integrated probe set that can generate high‐quality and highly complete data across all three loci types. We use these probes to generate data for 44 phylogenetically disparate taxa that collectively span approximately 33% of terrestrial vertebrate diversity. Our results generated an average of 4.29 Mb across 4709 loci per individual, of which an average of 2.99 Mb was sequenced to high enough coverage (≥10×) to use for population genetic analyses. We validate the utility of these loci for both phylogenomic and population genomic questions, provide a comparison among these locus sets of their relative usefulness and suggest areas for future improvement.     

Digest: Resolving phylogenomic conflicts in characiform fishes†

How can taxonomists best resolve the challenge of curating and analyzing large phylogenomic datasets that produce incongruent but highly supported topologies? Betancur‐R et al. used a recently established hypothesis‐testing procedure on a large dataset of genes and species to study the evolutionary relationships of characiform fishes, finding that past conclusions of non‐monophyly may have been problematic and establishing monophyly with high confidence. The new findings highlight the importance of using dense taxon sampling to resolve conflicting relationships with phylogenomic data.     

Influences of species, latitudes and methodologies on estimates of phenological response to global warming

New analyses are presented addressing the global impacts of recent climate change on phenology of plant and animal species. A meta‐analysis spanning 203 species was conducted on published datasets from the northern hemisphere. Phenological response was examined with respect to two factors: distribution of species across latitudes and taxonomic affiliation or functional grouping of target species. Amphibians had a significantly stronger shift toward earlier breeding than all other taxonomic/functional groups, advancing more than twice as fast as trees, birds and butterflies. In turn, butterfly emergence or migratory arrival showed three times stronger advancement than the first flowering of herbs, perhaps portending increasing asynchrony in insect–plant interactions. Response was significantly stronger at higher latitudes where warming has been stronger, but latitude explained < 4% of the variation. Despite expectation, latitude was not yet an important predictor of climate change impacts on phenology. The only two previously published estimates of the magnitude of global response are quite different: 2.3 and 5.1 days decade⁻¹ advancement. The scientific community has assumed this difference to be real and has attempted to explain it in terms of biologically relevant phenomena: specifically, differences in distribution of data across latitudes, taxa or time periods. Here, these and other possibilities are explored. All analyses indicate that the difference in estimated response is primarily due to differences between the studies in criteria for incorporating data. It is a clear and automatic consequence of the exclusion by one study of data on ‘stable’ (nonresponsive) species. Once this is accounted for, the two studies support each other, generating similar conclusions despite analyzing substantially nonoverlapping datasets. Analyses here on a new expanded dataset estimate an overall spring advancement across the northern hemisphere of 2.8 days decade⁻¹. This is the first quantitative analysis showing that data‐sampling methodologies significantly impact global (synthetic) estimates of magnitude of global warming response.     

Ultraconserved element (UCE) probe set design: Base genome and initial design parameters critical for optimization

Targeted capture and enrichment approaches have proven effective for phylogenetic study. Ultraconserved elements (UCEs) in particular have exhibited great utility for phylogenomic analyses, with the software package phyluce being among the most utilized pipelines for UCE phylogenomics, including probe design. Despite the success of UCEs, it is becoming increasing apparent that diverse lineages require probe sets tailored to focal taxa in order to improve locus recovery. However, factors affecting probe design and methods for optimizing probe sets to focal taxa remain underexplored. Here, we use newly available beetle (Coleoptera) genomic resources to investigate factors affecting UCE probe set design using phyluce . In particular, we explore the effects of stringency during initial design steps, as well as base genome choice on resulting probe sets and locus recovery. We found that both base genome choice and initial bait design stringency parameters greatly alter the number of resultant probes included in final probe sets and strongly affect the number of loci detected and recovered during in silico testing of these probe sets. In addition, we identify attributes of base genomes that correlated with high performance in probe design. Ultimately, we provide a recommended workflow for using Phyluce to design an optimized UCE probe set that will work across a targeted lineage, and use our findings to develop a new, open‐source UCE probe set for beetles of the suborder Adephaga.