Do missing data influence the accuracy of divergence-time estimation with BEAST?

A time-calibrated salamander phylogeny including 765 species and 503 genes

Recent time-calibrated amphibian phylogenies agree on the family-level relationships among extant salamanders but had disparate sampling regimes and inferred very different divergence times. For example, a recent phylogenomic study based on 220 nuclear loci had limited taxon sampling (41 species) and estimated relatively young divergence dates, whereas a more extensive supermatrix study based on 15 genes and 481 species estimated dates that were 22-45 million years older for major clades. Here, we combined phylogenomic and supermatrix approaches to estimate the largest salamander phylogeny to date based on molecular markers. Our matrix contained 765 salamander species and 503 genes (with 92.3% missing data overall). We included 284 more species than the previous largest salamander phylogeny (59% increase) and sampled approximately 93% of all currently described salamander species. Our dating analyses incorporated more than twice as many fossil calibration points within salamanders as previous studies. Maximum-likelihood estimates of tree topology yielded family-level relationships that were consistent with earlier studies. Nearly all species were placed in the expected genera, despite extensive missing data in many species. Bootstrap support was generally high across the tree but was poor in some clades where sampling of genes was limited (e.g., among some bolitoglossine salamanders). The dating analyses yielded age estimates for major clades that were generally intermediate between those from the previous phylogenomic and supermatrix analyses. We also provide a set of 200 time-calibrated trees for use in comparative analyses.

A Phylogenomic Backbone for Acoelomorpha Inferred From Transcriptomic Data

Xenacoelomorpha are mostly microscopic, morphologically simple worms, lacking many structures typical of other bilaterians. Xenacoelomorphs-which include three main groups, namely Acoela, Nemertodermatida, and Xenoturbella-have been proposed to be an early diverging Bilateria, sister to protostomes and deuterostomes, but other phylogenomic analyses have recovered this clade nested within the deuterostomes, as sister to Ambulacraria. The position of Xenacoelomorpha within the metazoan tree has understandably attracted a lot of attention, overshadowing the study of phylogenetic relationships within this group. Given that Xenoturbella includes only six species whose relationships are well understood, we decided to focus on the most speciose Acoelomorpha (Acoela + Nemertodermatida). Here, we have sequenced 29 transcriptomes, doubling the number of sequenced species, to infer a backbone tree for Acoelomorpha based on genomic data. The recovered topology is mostly congruent with previous studies. The most important difference is the recovery of Paratomella as the first off-shoot within Acoela, dramatically changing the reconstruction of the ancestral acoel. Besides, we have detected incongruence between the gene trees and the species tree, likely linked to incomplete lineage sorting, and some signal of introgression between the families Dakuidae and Mecynostomidae, which hampers inferring the correct placement of this family and, particularly, of the genus Notocelis. We have also used this dataset to infer for the first time diversification times within Acoelomorpha, which coincide with known bilaterian diversification and extinction events. Given the importance of morphological data in acoelomorph phylogenetics, we tested several partitions and models. Although morphological data failed to recover a robust phylogeny, phylogenetic placement has proven to be a suitable alternative when a reference phylogeny is available.

The Phylogenetic Relationships of Major Lizard Families Using Mitochondrial Genomes and Selection Pressure Analyses in Anguimorpha

Anguimorpha, within the order Squamata, represents a group with distinct morphological and behavioral characteristics in different ecological niches among lizards. Within Anguimorpha, there is a group characterized by limb loss, occupying lower ecological niches, concentrated within the subfamily Anguinae. Lizards with limbs and those without exhibit distinct locomotor abilities when adapting to their habitats, which in turn necessitate varying degrees of energy expenditure. Mitochondria, known as the metabolic powerhouses of cells, play a crucial role in providing approximately 95% of an organism's energy. Functionally, mitogenomes (mitochondrial genomes) can serve as a valuable tool for investigating potential adaptive evolutionary selection behind limb loss in reptiles. Due to the variation of mitogenome structures among each species, as well as its simple genetic structure, maternal inheritance, and high evolutionary rate, the mitogenome is increasingly utilized to reconstruct phylogenetic relationships of squamate animals. In this study, we sequenced the mitogenomes of two species within Anguimorpha as well as the mitogenomes of two species in Gekkota and four species in Scincoidea. We compared these data with the mitogenome content and evolutionary history of related species. Within Anguimorpha, between the mitogenomes of limbless and limbed lizards, a branch-site model analysis supported the presence of 10 positively selected sites: Cytb protein (at sites 183 and 187), ND2 protein (at sites 90, 155, and 198), ND3 protein (at site 21), ND5 protein (at sites 12 and 267), and ND6 protein (at sites 72 and 119). These findings suggested that positive selection of mitogenome in limbless lizards may be associated with the energy requirements for their locomotion. Additionally, we acquired data from 205 mitogenomes from the NCBI database. Bayesian inference (BI) and Maximum Likelihood (ML) trees were constructed using the 13 mitochondrial protein-coding genes (PCGs) and two rRNAs (12S rRNA and 16S rRNA) from 213 mitogenomes. Our phylogenetic tree and the divergence time estimates for Squamata based on mitogenome data are consistent with results from previous studies. Gekkota was placed at the root of Squamata in both BI and ML trees. However, within the Toxicofera clade, due to long-branch attraction, Anguimorpha and (Pleurodonta + (Serpentes + Acrodonta)) were closely related groupings, which might indicate errors and also demonstrate that mitogenome-based phylogenetic trees may not effectively resolve long-branch attraction issues. Additionally, we reviewed the origin and diversification of Squamata throughout the Mesozoic era, suggesting that Squamata originated in the Late Triassic (206.05 Mya), with the diversification of various superfamilies occurring during the Cretaceous period. Future improvements in constructing squamate phylogenetic relationships using mitogenomes will rely on identifying snake and acrodont species with slower evolutionary rates, ensuring comprehensive taxonomic coverage of squamate diversity, and increasing the number of genes analyzed.

A supermatrix phylogeny of the world's bees (Hymenoptera: Anthophila)

The increasing availability of large molecular phylogenies has provided new opportunities to study the evolution of species traits, their origins and diversification, and biogeography; yet there are limited attempts to synthesise existing phylogenetic information for major insect groups. Bees (Hymenoptera: Anthophila) are a large group of insect pollinators that have a worldwide distribution, and a wide variation in ecology, morphology, and life-history traits, including sociality. For these reasons, as well as their major economic importance as pollinators, numerous molecular phylogenetic studies of family and genus-level relationships have been published, providing an opportunity to assemble a bee 'tree-of-life'. We used publicly available genetic sequence data, including phylogenomic data, reconciled to a taxonomic database, to produce a concatenated supermatrix phylogeny for the Anthophila comprising 4,586 bee species, representing 23% of species and 82% of genera. At family, subfamily, and tribe levels, support for expected relationships was robust, but between and within some genera relationships remain uncertain. Within families, sampling of genera ranged from 67 to 100% but species coverage was lower (17-41%). Our phylogeny mostly reproduces the relationships found in recent phylogenomic studies with a few exceptions. We provide a summary of these differences and the current state of molecular data available and its gaps. We discuss the advantages and limitations of this bee supermatrix phylogeny (available online at beetreeoflife.org), which may enable new insights into long standing questions about evolutionary drivers in bees, and potentially insects more generally.

Frog phylogeny: A time-calibrated, species-level tree based on hundreds of loci and 5,242 species

Large-scale, time-calibrated phylogenies from supermatrix studies have become crucial for evolutionary and ecological studies in many groups of organisms. However, in frogs (anuran amphibians), there is a serious problem with existing supermatrix estimates. Specifically, these trees are based on a limited number of loci (15 or fewer), and the higher-level relationships estimated are discordant with recent phylogenomic estimates based on much larger numbers of loci. Here, we attempted to rectify this problem by generating an expanded supermatrix and combining this with data from phylogenomic studies. To assist in aligning ribosomal sequences for this supermatrix, we developed a new program (TaxonomyAlign) to help perform taxonomy-guided alignments. The new combined matrix contained 5,242 anuran species with data from 307 markers, but with 95% missing data overall. This dataset represented a 71% increase in species sampled relative to the previous largest supermatrix analysis of anurans (adding 2,175 species). Maximum-likelihood analyses generated a tree in which higher-level relationships (and estimated clade ages) were generally concordant with those from phylogenomic analyses but were more discordant with the previous largest supermatrix analysis. We found few obvious problems arising from the extensive missing data in most species. We also generated a set of 100 time-calibrated trees for use in comparative analyses. Overall, we provide an improved estimate of anuran phylogeny based on the largest number of combined taxa and markers to date. More broadly, we demonstrate the potential to combine phylogenomic and supermatrix analyses in other groups of organisms.

Parallel evolution, atavism, and extensive introgression explain the radiation of Epimedium sect. Diphyllon (Berberidaceae) in southern East Asia

East Asia is the richest region of plant biodiversity in the northern temperate zone, and its radiation provides key insights for understanding rapid speciation, including evolutionary patterns and processes. However, it is challenging to investigate the recent evolutionary radiation among plants because of the lack of genetic divergence, phenotypic convergence, and interspecific gene flow. Epimedium sect. Diphyllon is a rarely studied plant lineage endemic to East Asia, especially highly diversified in its southern part. In this study, we report a robust phylogenomic analysis based on genotyping-by-sequencing data of this lineage. The results revealed a clear biogeographic pattern for Epimedium sect. Diphyllon with recognition into two major clades corresponding to the Sino-Himalayan and Sino-Japanese subkingdoms of East Asian Flora and rapid diversification of the extant species dated to the Pleistocene. Evolutionary radiation of Epimedium sect. Diphyllon is characterized by recent and predominant parallel evolution and atavism between the two subkingdom regions, with extensive reticulating hybridization within each region during the course of diversification in southern East Asia. A parallel-atavism-introgression hypothesis is referred to in explaining the radiation of plant diversity in southern East Asia, which represents a potential model for the rapid diversification of plants under global climate cooling in the late Tertiary. Our study advances our understanding of the evolutionary processes of plant radiation in East Asia as well as in other biodiversity hotspot regions.

Species assemblages of insular Caribbean Sticta (lichenized Ascomycota: Peltigerales) over ecological and evolutionary time scales

Phylogenetic approaches to macroevolution have provided unique insight into evolutionary relationships, ancestral ranges, and diversification patterns for many taxa. Similar frameworks have also been developed to assess how environmental and/or spatial variables shape species diversity and distribution patterns at different spatial/temporal scales, but studies implementing these are still scarce for many groups, including lichens. Here, we combine phylogeny-based ancestral range reconstruction and diversification analysis with community phylogenetics to reconstruct evolutionary origins and assess patterns of taxonomic and phylogenetic relatedness between island communities of the lichenized fungal genus Sticta in the Caribbean. Sampling was carried out in the Greater Antilles (Cuba, Jamaica, Dominican Republic, and Puerto Rico) and Lesser Antilles (Dominica, Guadeloupe, and Martinique). Data for six molecular loci were obtained for 64 candidate Caribbean species and used to perform both macroevolutionary phylogenetics, which also included worldwide taxa, and phylobetadiversity analyses, which emphasized island-level communities. Our work uncovered high levels of island endemism (∼59%) in Caribbean Sticta. We estimate initial colonization of the region occurred about 19 Mya from a South American ancestor. Reverse migration events by Caribbean lineages to South America were also inferred. We found no evidence for increased diversification rates associated with range expansion into the Caribbean. Taxonomic and phylogenetic turnover between island-level communities was most strongly correlated with environmental variation rather than with geographic distance. We observed less dissimilarity among communities from the Dominican Republic and Jamaica than between these islands and the Lesser Antilles/Puerto Rico. High levels of hidden diversity and endemism in Caribbean Sticta reaffirm that islands are crucial for the maintenance of global biodiversity of lichenized fungi. Altogether, our findings suggest that strong evolutionary links exist between Caribbean and South American biotas but at regional scales, species assemblages exhibit complex taxonomic and phylogenetic relationships that are determined by local environments and shared evolutionary histories.

Revised molecular phylogeny, global biogeography, and diversification of palms subfamily Coryphoideae (Arecaceae) based on low copy nuclear and plastid regions

Coryphoideae are palmate-leaved palms from the family Arecaceae consisting of 46 genera representing 421 species. Although several phylogenetic analyses based on different genomic regions have been carried out on Coryphoideae, a fully resolved molecular phylogenetic tree has not been reported yet. To achieve this, we applied two phylogenetic reconstruction methods: Maximum Likelihood and Bayesian Inference, using amplified sampling by retrieving chloroplast and nuclear DNA sequences from NCBI and adding newly produced sequences from Indian accession into the dataset. The same dataset (chloroplast + nuclear DNA sequences) was used to estimate divergence times and the evolutionary history of Coryphoideae with a Bayesian uncorrelated, lognormal relaxed-clock approach and a Statistical Divergence-Vicariance Analysis method, respectively. The phylogenetic analyses based on a combined chloroplast and nuclear DNA sequence dataset showed well-resolved relationships within the subfamily. Both phylogenetic trees divide Coryphoideae into two main groups: CSPT (Crysophileae, Sabaleae, Phoeniceae, and Trachycarpeae) and the Syncarpous group. These main groups are segregated into eight tribes (Trachycarpeae, Phoeniceae, Sabaleae, Crysophileae, Borasseae, Corypheae, Caryoteae, and Chuniophoeniceae) and four subtribes (Rhapidine, Livistoninae, Hyphaeninae, and Lataniinae) with strong support-values. Most previously unresolved and doubtful relationships within tribes Trachycarpeae and Crysophilieae are now resolved and well-supported. The reconstructed phylogenetic trees support all previous systematic revisions of the subfamily. All Indian sampled species of Arenga, Bentinckia, Hyphaene, and Trachycarpus show close relation with their respective congeneric species. Molecular dating results and integration of biogeography suggest that Coryphoideae originated in Laurasia at ~95.12 Ma and then diverged into the tropical and subtropical regions of the whole world. This study offers the correct combination of nuclear and plastid regions to test the current and future systematic revisions.

Reference-Based RADseq Unravels the Evolutionary History of Polar Species in 'the Crux Lichenologorum' Genus Usnea (Parmeliaceae, Ascomycota)

Nearly 90% of fungal diversity, one of the most speciose branches in the tree of life, remains undescribed. Lichenized fungi as symbiotic associations are still a challenge for species delimitation, and current species diversity is vastly underestimated. The ongoing democratization of Next-Generation Sequencing is turning the tables. Particularly, reference-based RADseq allows for metagenomic filtering of the symbiont sequence and yields robust phylogenomic trees of closely related species. We implemented reference-based RADseq to disentangle the evolution of neuropogonoid lichens, which inhabit harsh environments and belong to Usnea (Parmeliaceae, Ascomycota), one of the most taxonomically intriguing genera within lichenized fungi. Full taxon coverage of neuropogonoid lichens was sampled for the first time, coupled with phenotype characterizations. More than 20,000 loci of 126 specimens were analyzed through concatenated and coalescent-based methods, including time calibrations. Our analysis addressed the major taxonomic discussions over recent decades. Subsequently, two species are newly described, namely U. aymondiana and U. fibriloides, and three species names are resurrected. The late Miocene and Pliocene-Pleistocene boundary is inferred as the timeframe for neuropogonoid lichen diversification. Ultimately, this study helped fill the gap of fungal diversity by setting a solid backbone phylogeny which raises new questions about which factors may trigger complex evolutionary scenarios.

Illusion of flight? Absence, evidence and the age of winged insects

The earliest fossils of winged insects (Pterygota) are mid-Carboniferous (latest Mississippian, 328–324 Mya), but estimates of their age based on fossil-calibrated molecular phylogenetic studies place their origin at 440–370 Mya during the Silurian or Devonian. This discrepancy would require that winged insects evaded fossilization for at least the first ~50 Myr of their history. Here, we examine the plausibility of such a gap in the fossil record, and possible explanations for it, based on comparisons with the fossil records of other arthropod groups, the distribution of first occurrence dates of pterygote families, phylogenetically informed simulations of the fossilization of Palaeozoic insects, and re-analysis of data presented by Misof and colleagues using updated fossil calibrations under a variety of prior probability settings. We do not find support for the mechanisms previously suggested to account for such an extended gap in the pterygote fossil record, including sampling bias, preservation bias, and body size. We suggest that inference of an early origin of Pterygota long prior to their first appearance in the fossil record is probably an analytical artefact of taxon sampling and choice of fossil calibration points, possibly compounded by heterogeneity in rates of sequence evolution or speciation, including radiations or ‘bursts’ during their early history.

Molecular phylogenetics reveals the evolutionary history of marine fishes (Actinopterygii) endemic to the subtropical islands of the Southwest Pacific

Remote oceanic islands of the Pacific host elevated levels of actinopterygian (ray-finned fishes) endemism. Characterizing the evolutionary histories of these endemics has provided insight into the generation and maintenance of marine biodiversity in many regions. The subtropical islands of Lord Howe, Norfolk, and Rangitāhua (Kermadec) in the Southwest Pacific are yet to be comprehensively studied. Here, we characterize the spatio-temporal diversification of marine fishes endemic to these Southwest Pacific islands by combining molecular phylogenies and the geographic distribution of species. We built Bayesian ultrametric trees based on open-access and newly generated sequences for five mitochondrial and ten nuclear loci, and using fossil data for time calibration. We present the most comprehensive phylogenies to date for marine ray-finned fish genera, comprising 34 species endemic to the islands, including the first phylogenetic placements for 11 endemics. Overall, our topologies confirm the species status of all endemics, including three undescribed taxa. Our phylogenies highlight the predominant affinity of these endemics with the Australian fish fauna (53%), followed by the East Pacific (15%), and individual cases where the closest sister taxon of our endemic is found in the Northwest Pacific and wider Indo-Pacific. Nonetheless, for a quarter of our focal endemics, their geographic affinity remains unresolved due to sampling gaps within their genera. Our divergence time estimates reveal that the majority of endemic lineages (67.6%) diverged after the emergence of Lord Howe (6.92 Ma), the oldest subtropical island in the Southwest Pacific, suggesting that these islands have promoted diversification. However, divergence ages of some endemics pre-date the emergence of the islands, suggesting they may have originated outside of these islands, or, in some cases, ages may be overestimated due to unsampled taxa. To fully understand the role of the Southwest Pacific subtropical islands as a 'cradle' for diversification, our study advocates for further regional surveys focused on tissue collection for DNA analysis.

Variation in the substitution rates among the human mitochondrial haplogroup U sublineages

Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from ¹⁴C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.

A snapshot of progenitor-derivative speciation in Iberodes (Boraginaceae)

Traditional classification of speciation modes has focused on physical barriers to gene flow. Allopatric speciation with complete reproductive isolation is viewed as the most common mechanism of speciation. Parapatry and sympatry, by contrast, entail speciation in the face of ongoing gene flow, making them more difficult to detect. The genus Iberodes (Boraginaceae, NW Europe) comprises five species with contrasting morphological traits, habitats, and species distributions. Based on the predominance of narrow and geographically distant endemic species, we hypothesized that geographic barriers were responsible for most speciation events in Iberodes. We undertook an integrative study including: (i) phylogenomics through restriction-site associated DNA sequencing, (ii) genetic structure analyses, (iii) demographic modeling, (iv) morphometrics, and (v) climatic niche modeling and niche overlap analysis. Results revealed a history of recurrent progenitor-derivative speciation manifested by a paraphyletic pattern of nested species differentiation. Budding speciation mediated by ecological differentiation is suggested for the coastal lineage, deriving from the inland widespread I. linifolia during Late Pliocene. Meanwhile, geographic isolation followed by niche shifts are suggested for the more recent differentiation of the coastland taxa. Our work provides a model for distinguishing speciation via ecological differentiation of peripheral, narrowly endemic I. kuzinskyanae and I. littoralis from a widespread extant ancestor, I. linifolia. Ultimately, our results illustrate a case of Pliocene speciation in the probable absence of geographic barriers and get away from the traditional cladistic perspective of speciation as producing two species from an extinct ancestor, thus reminding us that phylogenetic trees tell only part of the story.

Phylogenetic value of jaw elements of lacertid lizards (Squamata: Lacertoidea): a case study with Oligocene material from France

Several extinct species are known from the family Lacertidae, but due to poor preservation, many of them are based on single bones. Here, we compare phylogenetic signals of disarticulated premaxillae, maxillae and dentaries of lacertids from four French Oligocene localities (Coderet, La Colombière, Roqueprune 2, Mas de Got B). We identified five morphotypes among the premaxillae, six among the maxillae, and ten among the dentaries. These morphotypes were scored as individual taxa per locality into three separate character matrices with the same 246 characters, one matrix for each jaw element. Subsequently, the phylogenetic position of the morphotypes was tested using maximum parsimony. The consensus trees with the dentaries and the maxillae found a large polytomy including all taxa except the outgroup taxon Gekko gecko. The consensus tree with the premaxillae showed a considerably more resolved topology but found all morphotype taxa outside Lacertidae. In a second step, we compared the constitution of our three datasets and the morphotype taxa. Our results suggest that a combination of convergent characters and missing data led to the outgroup position of the premaxilla morphotype taxa. The poor resolution of the maxillae strict consensus is likely a consequence of their fragmentary preservation. For the dentaries, a high amount of missing data due to the high number of morphotype taxa most likely caused the poor tree resolution. Indeed, tests with fewer morphotypes found tree resolutions comparable to the premaxilla data. When linking the morphotypes, five possible lacertid "species" were found. Comparison with already known French Oligocene lacertid species points to a slightly higher species richness of Lacertidae at that time than known before. Reliable species classification based on phylogeny only seems possible when combining the jaw elements or in association with other cranial and postcranial material, putting some doubt on species identifications based on single bones.

Divergence time estimation of Galliformes based on the best gene shopping scheme of ultraconserved elements

BACKGROUND

Divergence time estimation is fundamental to understanding many aspects of the evolution of organisms, such as character evolution, diversification, and biogeography. With the development of sequence technology, improved analytical methods, and knowledge of fossils for calibration, it is possible to obtain robust molecular dating results. However, while phylogenomic datasets show great promise in phylogenetic estimation, the best ways to leverage the large amounts of data for divergence time estimation has not been well explored. A potential solution is to focus on a subset of data for divergence time estimation, which can significantly reduce the computational burdens and avoid problems with data heterogeneity that may bias results.

RESULTS

In this study, we obtained thousands of ultraconserved elements (UCEs) from 130 extant galliform taxa, including representatives of all genera, to determine the divergence times throughout galliform history. We tested the effects of different "gene shopping" schemes on divergence time estimation using a carefully, and previously validated, set of fossils. Our results found commonly used clock-like schemes may not be suitable for UCE dating (or other data types) where some loci have little information. We suggest use of partitioning (e.g., PartitionFinder) and selection of tree-like partitions may be good strategies to select a subset of data for divergence time estimation from UCEs. Our galliform time tree is largely consistent with other molecular clock studies of mitochondrial and nuclear loci. With our increased taxon sampling, a well-resolved topology, carefully vetted fossil calibrations, and suitable molecular dating methods, we obtained a high quality galliform time tree.

CONCLUSIONS

We provide a robust galliform backbone time tree that can be combined with more fossil records to further facilitate our understanding of the evolution of Galliformes and can be used as a resource for comparative and biogeographic studies in this group.

Redefining the Distributional Boundaries and Phylogenetic Relationships for Ctenomids From Central Argentina

With about 68 recognized living species, subterranean rodents of the genus Ctenomys are found in a multiplicity of habitats, from the dunes of the Atlantic coast to the Andes Mountains, including environments ranging from humid steppes of Pampas to the dry deserts of Chaco region. However, this genus needs an exhaustive reevaluation of its systematic and phylogenetic relationships regarding the different groups that compose it. This knowledge is essential to propose biodiversity conservation strategies both at species level and at higher hierarchical levels. In order to clarify the taxonomy and the recent evolutionary history from populations of Ctenomys in the Pampas region, Argentina, phylogenetic relationships among them were evaluated using mitochondrial DNA sequences: gene encoding cytochrome b protein (1,140 bp) and the non-coding D-loop region (434 bp). To infer the divergence times inside the Ctenomys clade, a Bayesian calibrate tree using fossil remains data from different families within Caviomorpha was performed at first. Secondly, that calibration data was used as priors in a new Bayesian phylogenetic inference within the genus Ctenomys. This phylogenetic tree emphasized on species currently distributed on the Pampas region, more precisely considering both the talarum and mendocinus groups. Bayesian inferences (BI) were integrated with the results of a Maximum Likelihood approach (ML). Based on these results, the distributional limits of the mendocinus and talarum groups appear to be related to the physiognomy of the Pampas region soils. On the other hand, the validity of C. pundti complex as a differentiated species of C. talarum is debated. According to previous evidence from morphological and chromosomal studies, these results show a very low divergence between those species that originally were classified within the talarum group. Mitochondrial DNA sequences from populations associated with these putative species have not recovered as reciprocal monophyletic groups in the phylogenetic analyses. In conclusion, C. talarum and C. pundti complex might be considered as the same biological species, or lineages going through a recent or incipient differentiation process. The results obtained in this study have important implications for conservation policies and practices, since both species are currently categorized as Vulnerable and Endangered, respectively.

UCE Phylogenomics, detection of a putative hybrid population, and one older mitogenomic node age of Batrachuperus salamanders

The prevalence of incomplete lineage sorting complicates the examination of hybridization and species-level paraphyly with gene trees of a small number of loci. In Asian mountain salamanders of the genus Batrachuperus, possible hybridization and species paraphyly had been identified by utilizing mitochondrial genealogy and fixed allozyme differences. Here we sampled 2909 UCEs in 44 local populations from all six Batrachuperus species, inferred gene and species trees, compared them with mitochondrial and allozyme results, and examined the potential hybridization and species paraphyly. The clustering pattern of single-locus trees, increased proportion of heterozygous SNPs, allele frequency-based migration edge estimation, and intrapopulation long branches (as expected from an increase of genetic lineage and nucleotide diversity) support that an eastern B. karlschmidti population has experienced admixture with B. tibetanus. On the 2909-UCE concatenated and species trees, lower nodal supports were observed when similar proportions of loci agreed with alternative topologies, i.e., a reciprocal monophyly between a Pengxian lineage and the remainder of B. pinchonii (0.379) or a paraphyly of the latter with respect to Pengxian (0.362). The UCE phylogenomics agreed with the relatively recent groupings in the allozyme dendrogram. Despite incomplete lineage sorting, the mitochondrial trees were similar to the UCE trees for deeper relationships of the genus. However, one significant branch-length level discordance was identified. The branch between the common ancestor of B. daochengensis and B. yenyuanensis and common ancestor of the genus was approximately three times shorter on the mitochondrial tree than on the UCE tree, suggesting that the split of the mitochondrial lineages was likely a few million years earlier than the split of species. This finding supports considering possible ancestral polymorphism when interpreting different divergence dates estimated from mitochondrial and genome-wide data.

DNA barcodes combined with multi-locus data of representative taxa can generate reliable higher-level phylogenies

Taxa are frequently labeled incertae sedis when their placement is debated at ranks above the species level, such as their subgeneric, generic, or subtribal placement. This is a pervasive problem in groups with complex systematics due to difficulties in identifying suitable synapomorphies. In this study, we propose combining DNA barcodes with a multilocus backbone phylogeny in order to assign taxa to genus or other higher-level categories. This sampling strategy generates molecular matrices containing large amounts of missing data that are not distributed randomly: barcodes are sampled for all representatives, and additional markers are sampled only for a small percentage. We investigate the effects of the degree and randomness of missing data on phylogenetic accuracy using simulations for up to 100 markers in 1000-tips trees, as well as a real case: the subtribe Polyommatina (Lepidoptera: Lycaenidae), a large group including numerous species with unresolved taxonomy. Our simulation tests show that when a strategic and representative selection of species for higher-level categories has been made for multigene sequencing (approximately one per simulated genus), the addition of this multigene backbone DNA data for as few as 5–10% of the specimens in the total data set can produce high-quality phylogenies, comparable to those resulting from 100% multigene sampling. In contrast, trees based exclusively on barcodes performed poorly. This approach was applied to a 1365-specimen data set of Polyommatina (including ca. 80% of described species), with nearly 8% of representative species included in the multigene backbone and the remaining 92% included only by mitochondrial COI barcodes, a phylogeny was generated that highlighted potential misplacements, unrecognized major clades, and placement for incertae sedis taxa. We use this information to make systematic rearrangements within Polyommatina, and to describe two new genera. Finally, we propose a systematic workflow to assess higher-level taxonomy in hyperdiverse groups. This research identifies an additional, enhanced value of DNA barcodes for improvements in higher-level systematics using large data sets. [Birabiro; DNA barcoding; incertae sedis; Kipepeo; Lycaenidae; missing data; phylogenomic; phylogeny; Polyommatina; supermatrix; systematics; taxonomy]

Phylogenetics of mud snakes (Squamata: Serpentes: Homalopsidae): A paradox of both undescribed diversity and taxonomic inflation

Mud snakes (Serpentes: Homalopsidae) are a family of 55 described, mainly aquatic, species primarily distributed throughout mainland Southeast Asia and the Indo-Australian Archipelago. Although they have been the focus of prior research, the basic relationships amongst genera and species remain poorly known. We used a combined mitochondrial and nuclear gene dataset to infer their phylogenetic relationships, using the highest levels of taxon and geographic sampling for any homalopsid phylogeny to date (62% generic and 62% species coverage; 140 individuals). Our results recover two reciprocally monophyletic groups: the fangless Brachyorrhos and its sister clade comprised of all rear-fanged homalopsids. Most genera and interspecific relationships were monophyletic and strongly supported, but intergeneric relationships and intraspecific population structure lack support. We find evidence of both undescribed diversity as well as cases of taxonomic inflation within several species. Tree-based species delimitation approaches (mPTP) support potential new candidate species as distinct from their conspecifics and also suggest that many named taxa may not be distinct species. Divergence date estimation and lineage-through-time analyses indicate lower levels of speciation in the Eocene, with a subsequent burst in diversification in the Miocene. Homalopsids may have diversified most rapidly during the Pliocene and Pleistocene, possibly in relation to tectonic shifts and sea-level fluctuations that took place in Sundaland and the Sahul Shelf. Our analyses provide new insights on homalopsid taxonomy, a baseline phylogeny for the family, and further biogeographic implications demonstrating how dynamic tectonics and Quaternary sea level changes may have shaped a widespread, diverse family of snakes.

Rapid radiation of angraecoids (Orchidaceae, Angraecinae) in tropical Africa characterised by multiple karyotypic shifts under major environmental instability

Angraecoid orchids present a remarkable diversity of chromosome numbers, which makes them a highly suitable system for exploring the impact of karyotypic changes on cladogenesis, diversification and morphological differentiation. We compiled an annotated cytotaxonomic checklist for 126 species of Angraecinae, which was utilised to reconstruct chromosomal evolution using a newly-produced, near-comprehensive phylogenetic tree that includes 245 angraecoid taxa. In tandem with this improved phylogenetic framework, using combined Bayesian, maximum likelihood and parsimony approaches on ITS-1 and five plastid markers, we propose a new cladistic nomenclature for the angraecoids, and we estimate a new timeframe for angraecoid radiation based on a secondary calibration, and calculate diversification rates using a Bayesian approach. Coincident divergence dates between clades with identical geographical distributions in the angraecoids and the pantropical orchid genus Bulbophyllum suggest that the same events may have intervened in the dispersal of these two epiphytic groups between Asia, continental Africa, Madagascar and the Neotropics. The major angraecoid lineages probably began to differentiate in the Middle Miocene, and most genera and species emerged respectively around the Late Miocene-Pliocene boundary and the Pleistocene. Ancestral state reconstruction using maximum likelihood estimation revealed an eventful karyotypic history dominated by descending dysploidy. Karyotypic shifts seem to have paralleled cladogenesis in continental tropical Africa, where approximately 90% of the species have descended from at least one inferred transition from n = 17-18 to n = 25 during the Middle Miocene Climatic Transition, followed by some clade-specific descending and ascending dysploidy from the Late Miocene to the Pleistocene. Conversely, detected polyploidy is restricted to a few species lineages mostly originating during the Pleistocene. No increases in net diversification could be related to chromosome number changes, and the apparent net diversification was found to be highest in Madagascar, where karyotypic stasis predominates. Finally, shifts in chromosome number appear to have paralleled the evolution of rostellum structure, leaflessness, and conspicuous changes in floral colour.

Evolution in the Model Genus Antirrhinum Based on Phylogenomics of Topotypic Material

Researchers in phylogenetic systematics typically choose a few individual representatives of every species for sequencing based on convenience (neighboring populations, herbarium specimens, samples provided by experts, garden plants). However, few studies are based on original material, type material or topotypic material (living specimens from the locality where the type material was collected). The use of type or topotypic material in phylogenetic studies is paramount particularly when taxonomy is complex, such as that of Antirrhinum (Plantaginaceae). In this paper, we used topotypic materials of Antirrhinum at the species level (34 species proposed by previous authors), 87 specimens representing the species distributions and >50,000 informative nucleotide characters (from ∼4,000 loci) generated by the genotyping-by-sequencing (GBS) technique: (i) to test two explicit taxonomic hypotheses widely followed by local taxonomic treatments; (ii) to robustly estimate phylogenetic relationships; (iii) to investigate the evolution of key morphological characters and biogeographic centers of differentiation. Two GBS phylogenies based on two datasets (87 localities and 34 topotypic specimens) revealed that: (1) Sutton's (1988) taxonomic account is the most congruent with phylogenetic results, whereas division of Antirrhinum into three major clades disagrees with Rothmaler's (1956) infrageneric classification; (2) monophyly of populations currently included in the same species is primarily supported; (3) the historically recognized Antirrhinum majus group is not monophyletic; (4) sister-group relationships are robust for eight species pairs; (5) the evolutionary radiation of 26 species since the Pliocene is underpinned given a high rate of diversification (0.54 spp. Myr-1); (6) a geographic pattern of speciation is reconstructed, with northern Iberia as the center of early diversification followed by more recent speciation in southeastern Iberia; and (7) multiple acquisitions of key taxonomic characters in the course of Antirrhinum diversification are strongly supported, with no evidence of hybridization between major clades. Our results also suggest incipient speciation in some geographic areas and point to future avenues of research in evolution and systematics of Antirrhinum.

Mitochondrial genomics reveals the evolutionary history of the porpoises (Phocoenidae) across the speciation continuum

Historical variation in food resources is expected to be a major driver of cetacean evolution, especially for the smallest species like porpoises. Despite major conservation issues among porpoise species (e.g., vaquita and finless), their evolutionary history remains understudied. Here, we reconstructed their evolutionary history across the speciation continuum. Phylogenetic analyses of 63 mitochondrial genomes suggest that porpoises radiated during the deep environmental changes of the Pliocene. However, all intra-specific subdivisions were shaped during the Quaternary glaciations. We observed analogous evolutionary patterns in both hemispheres associated with convergent evolution to coastal versus oceanic environments. This suggests that similar mechanisms are driving species diversification in northern (harbor and Dall's) and southern species (spectacled and Burmeister's). In contrast to previous studies, spectacled and Burmeister's porpoises shared a more recent common ancestor than with the vaquita that diverged from southern species during the Pliocene. The low genetic diversity observed in the vaquita carried signatures of a very low population size since the last 5,000 years. Cryptic lineages within Dall's, spectacled and Pacific harbor porpoises suggest a richer evolutionary history than previously suspected. These results provide a new perspective on the mechanisms driving diversification in porpoises and an evolutionary framework for their conservation.

Integrative taxonomy of extant maras supports the recognition of the genera Pediolagus and Dolichotis within the Dolichotinae (Rodentia, Caviidae)

The subfamily Dolichotinae (Rodentia, Caviidae) includes two living species (Dolichotis patagonum and D. salinicola) of cursorial, long-legged rodents that inhabit semiarid thorn scrubs and shrublands in southern South America. Some authors consider that the morphological differences between D. patagonum and D. salinicola warrant circumscription of the latter in its own genus, Pediolagus. Based on a phylogeny and divergence times grounded on molecular data, as well as a qualitative and quantitative assessment of morphological differences, the distinction between Dolichotis and Pediolagus is equal to or greater than differences seen for other intergeneric comparisons within Cavioidea. Based on these results, we argue that the taxonomic designations of Dolichotis and Pediolagus should be retained.

The evolution of autotomy in leaf-footed bugs

Sacrificing body parts is one of many behaviors that animals use to escape predation. This trait, termed autotomy, is classically associated with lizards. However, several other taxa also autotomize, and this trait has independently evolved multiple times throughout Animalia. Despite having multiple origins and being an iconic antipredatory trait, much remains unknown about the evolution of autotomy. Here, we combine morphological, behavioral, and genomic data to investigate the evolution of autotomy within leaf-footed bugs and allies (Insecta: Hemiptera: Coreidae + Alydidae). We found that the ancestor of leaf-footed bugs autotomized and did so slowly; rapid autotomy (<2 min) then arose multiple times. The ancestor likely used slow autotomy to reduce the cost of injury or to escape nonpredatory entrapment but could not use autotomy to escape predation. This result suggests that autotomy to escape predation is a co-opted benefit (i.e., exaptation), revealing one way that sacrificing a limb to escape predation may arise. In addition to identifying the origins of rapid autotomy, we also show that across species variation in the rates of autotomy can be explained by body size, distance from the equator, and enlargement of the autotomizable appendage.

Analyses of Plastome Sequences Improve Phylogenetic Resolution and Provide New Insight Into the Evolutionary History of Asian Sonerileae/Dissochaeteae

Sonerileae/Dissochaeteae (Melastomataceae) comprises ca. 50 genera, two thirds of which occur in Southeast Asia. Phylogenetic relationships within this clade remain largely unclear, which hampers our understanding of its origin, evolution, and biogeography. Here, we explored the use of chloroplast genomes in phylogenetic reconstruction of Sonerileae/Dissochaeteae, by sampling 138 species and 23 genera in this clade. A total of 151 complete plastid genomes were assembled for this study. Plastid genomic data provided better support for the backbone of the Sonerileae/Dissochaeteae phylogeny, and also for relationships among most closely related species, but failed to resolve the short internodes likely resulted from rapid radiation. Trees inferred from plastid genome and nrITS sequences were largely congruent regarding the major lineages of Sonerileae/Dissochaeteae. The present analyses recovered 15 major lineages well recognized in both nrITS and plastid phylogeny. Molecular dating and biogeographical analyses indicated a South American origin for Sonerileae/Dissochaeteae during late Eocene (stem age: 34.78 Mya). Two dispersal events from South America to the Old World were detected in late Eocene (33.96 Mya) and Mid Oligocene (28.33 Mya) respectively. The core Asian clade began to diversify around early Miocene in Indo-Burma and dispersed subsequently to Malesia and Sino-Japanese regions, possibly promoted by global temperature changes and East Asian monsoon activity. Our analyses supported previous hypothesis that Medinilla reached Madagascar by transoceanic dispersal in Miocene. In addition, generic limits of some genera concerned were discussed.

Post K-Pg diversification of the mammalian order Eulipotyphla as suggested by phylogenomic analyses of ultra-conserved elements

The origin of the mammalian order Eulipotyphla has been debated intensively with arguments around whether they began diversifying before or after the Cretaceous-Palaeogene (K-Pg) boundary at 66 Ma. Here, we used an in-solution nucleotide capture method and next generation DNA sequencing to determine the sequence of hundreds of ultra-conserved elements (UCEs), and conducted phylogenomic and molecular dating analyses for the four extant eulipotyphlan lineages-Erinaceidae, Solenodontidae, Soricidae, and Talpidae. Concatenated maximum-likelihood analyses with single or partitioned models and a coalescent species-tree analysis showed that divergences among the four major eulipotyphlan lineages occurred within a short period of evolutionary time, but did not resolve the interrelationships among them. Alternative suboptimal phylogenetic hypotheses received consistently the same amount of support from different UCE loci, and were not significantly different from the maximum likelihood tree topology, suggesting the prevalence of stochastic lineage sorting. Molecular dating analyses that incorporated among-lineage evolutionary rate differences supported a scenario where the four eulipotyphlan families diversified between 57.8 and 63.2 Ma. Given short branch lengths with low support values, traces of rampant genome-wide stochastic lineage sorting, and post K-Pg diversification, we concluded that the crown eulipotyphlan lineages arose through a rapid diversification after the K-Pg boundary when novel niches were created by the mass extinction of species.

Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

Molecular phylogeny reveals the gradual evolutionary transition to soft-bodiedness in click-beetles and identifies sub-Saharan Africa as a cradle of diversity for Drilini (Coleoptera: Elateridae)

Drilini are soft-bodied predatory click-beetles (Elateridae: Agrypninae) with incompletely metamorphosed females. Due to divergent morphology, their classification has been contentious. We present the first densely sampled molecular phylogeny of Drilini based on nuclear and mitochondrial markers. Altogether, 44 species, representing all genera, were analysed using maximum likelihood and the Bayesian approach. Molecular analyses recovered five major clades that were also well supported by morphology. Afrotropical lineages mark deep splits. Most Palearctic species belong to a terminal clade. A few species of the predominantly Afrotropical Selasia are distributed from Arabia up to the Himalayas and Thailand. The origin of Drilini is dated to the Late Eocene (~35.5 Mya) and rapid radiation is identified from the Eocene/Oligocene to the Middle Miocene. We describe the gradual transformation of male morphological traits, e.g. the level of sclerotization, structure of mouthparts, loss of the thoracic interlocking mechanism, shortened elytra and expanded larviform abdomen. Five new genera, Austroselasia, Habeshaselasia, Illubaboria, Malacodrilus and Mashaselasia are proposed. Latoselasia, previously a subgenus of Wittmerselasia, is elevated to the genus rank. Five new species, Habeshaselasia nekemtensis, H. illubaborensis, Illubaboria bicolor, Malacodrilus hajeki and Mashaselasia aethiopica are described.

An integrative phylogenomic approach illuminates the evolutionary history of cockroaches and termites (Blattodea)

Phylogenetic relationships among subgroups of cockroaches and termites are still matters of debate. Their divergence times and major phenotypic transitions during evolution are also not yet settled. We addressed these points by combining the first nuclear phylogenomic study of termites and cockroaches with a thorough approach to divergence time analysis, identification of endosymbionts, and reconstruction of ancestral morphological traits and behaviour. Analyses of the phylogenetic relationships within Blattodea robustly confirm previously uncertain hypotheses such as the sister-group relationship between Blaberoidea and remaining Blattodea, and Lamproblatta being the closest relative to the social and wood-feeding Cryptocercus and termites. Consequently, we propose new names for various clades in Blattodea: Cryptocercus + termites = Tutricablattae; Lamproblattidae + Tutricablattae = Kittrickea; and Blattoidea + Corydioidea = Solumblattodea. Our inferred divergence times contradict previous studies by showing that most subgroups of Blattodea evolved in the Cretaceous, reducing the gap between molecular estimates of divergence times and the fossil record. On a phenotypic level, the blattodean ground-plan is for egg packages to be laid directly in a hole while other forms of oviposition, including ovovivipary and vivipary, arose later. Finally, other changes in egg care strategy may have allowed for the adaptation of nest building and other novelties.

An In Silico Comparison of Protocols for Dated Phylogenomics

In the age of genome-scale DNA sequencing, choice of molecular marker arguably remains an important decision in planning a phylogenetic study. Using published genomes from 23 primate species, we make a standardized comparison of four of the most frequently used protocols in phylogenomics, viz., targeted sequence-enrichment using ultraconserved element and exon-capture probes, and restriction-site-associated DNA sequencing (RADseq and ddRADseq). Here, we present a procedure to perform in silico extractions from genomes and create directly comparable data sets for each class of marker. We then compare these data sets in terms of both phylogenetic resolution and ability to consistently and precisely estimate clade ages using fossil-calibrated molecular-clock models. Furthermore, we were also able to directly compare these results to previously published data sets from Sanger-sequenced nuclear exons and mitochondrial genomes under the same analytical conditions. Our results show-although with the exception of the mitochondrial genome data set and the smallest ddRADseq data set-that for uncontroversial nodes all data classes performed equally well, that is they recovered the same well supported topology. However, for one difficult-to-resolve node comprising a rapid diversification, we report well supported but conflicting topologies among the marker classes consistent with the mismodeling of gene tree heterogeneity as demonstrated by species tree analyses of single nucleotide polymorphisms. Likewise, clade age estimates showed consistent discrepancies between data sets under strict and relaxed clock models; for recent nodes, clade ages estimated by nuclear exon data sets were younger than those of the UCE, RADseq and mitochondrial data, but vice versa for the deepest nodes in the primate phylogeny. This observation is explained by temporal differences in phylogenetic informativeness (PI), with the data sets with strong PI peaks toward the present underestimating the deepest node ages. Finally, we conclude by emphasizing that while huge numbers of loci are probably not required for uncontroversial phylogenetic questions-for which practical considerations such as ease of data generation, sharing, and aggregating, therefore become increasingly important-accurately modeling heterogeneous data remains as relevant as ever for the more recalcitrant problems.

What explains high plant richness in East Asia? Time and diversification in the tribe Lysimachieae (Primulaceae)

What causes the disparity in biodiversity among regions is a fundamental question in biogeography, ecology, and evolutionary biology. Evolutionary and biogeographic processes (speciation, extinction, dispersal) directly determine species richness patterns, and can be studied using integrative phylogenetic approaches. However, the strikingly high richness of East Asia relative to other Northern Hemisphere regions remains poorly understood from this perspective. Here, for the first time, we test two general hypotheses (older colonization time, faster diversification rate) to explain this pattern, using the plant tribe Lysimachieae (Primulaceae) as a model system. We generated a new time-calibrated phylogeny for Lysimachieae (13 genes, 126 species), to estimate colonization times and diversification rates for each region and to test the relative importance of these two factors for explaining regional richness patterns. We find that neither time nor diversification rates alone explain richness patterns among regions in Lysimachieae. Instead, a new index that combines both factors explains global richness patterns in the group and their high East Asian biodiversity. Based on our results from Lysimachieae, we suggest that the high richness of plants in East Asia may be explained by a combination of older colonization times and faster diversification rates in this region.

Taxon sampling effects on the quantification and comparison of community phylogenetic diversity

Ecologists are increasingly making use of molecular phylogenies, especially in the fields of community ecology and conservation. However, these phylogenies are often used without full appreciation of their underlying assumptions and uncertainties. A frequent practice in ecological studies is inferring a phylogeny with molecular data from taxa only within the community of interest. These "inferred community phylogenies" are inherently biased in their taxon sampling. Despite the importance of comprehensive sampling in constructing phylogenies, the implications of using inferred community phylogenies in ecological studies have not been examined. Here, we evaluate how taxon sampling affects the quantification and comparison of community phylogenetic diversity using both simulated and empirical data sets. We demonstrate that inferred community trees greatly underestimate phylogenetic diversity and that the probability of incorrectly ranking community diversity can reach up to 25%, depending on the dating methods employed. We argue that to reach reliable conclusions, ecological studies must improve their taxon sampling and generate the best phylogeny possible.

Frequent and parallel habitat transitions as driver of unbounded radiations in the Cape flora

The enormous species richness in the Cape Floristic Region (CFR) of Southern Africa is the result of numerous radiations, but the temporal progression and possible mechanisms of these radiations are still poorly understood. Here, we explore the macroevolutionary dynamics of the Restionaceae, which include 340 species that are found in all vegetation types in the Cape flora and are ecologically dominant in fynbos. Using an almost complete (i.e., 98%) species-level time calibrated phylogeny and models of diversification dynamics, we show that species diversification is constant through the Cenozoic, with no evidence of an acceleration with the onset of the modern winter-wet climate, or a recent density-dependent slowdown. Contrary to expectation, species inhabiting the oldest (montane) and most extensive (drylands) habitats did not undergo higher diversification rates than species in the younger (lowlands) and more restricted (wetland) habitats. We show that the rate of habitat transitions is more closely related to the speciation rate than to time, and that more than a quarter of all speciation events are associated with habitat transitions. This suggests that the unbounded Restionaceae diversification resulted from numerous, parallel, habitat shifts, rather than persistence in a habitat stimulating speciation. We speculate that this could be one of the mechanisms resulting in the hyperdiverse Cape flora.

Where Fossils Dare and Males Matter: combined morphological and molecular analysis untangles the evolutionary history of the spider ant genus Leptomyrmex Mayr (Hymenoptera : Dolichoderinae)

The distinctive ant genus Leptomyrmex Mayr, 1862 had been thought to be endemic to Australasia for over 150 years, but enigmatic Neotropical fossils have challenged this view for decades. The present study responds to a recent and surprising discovery of extant Leptomyrmex species in Brazil with a thorough evaluation of the Dominican Republic fossil material, which dates to the Miocene. In the first case study of direct fossil inclusion within Formicidae Latreille, 1809, we incorporated both living and the extinct Leptomyrmex species. Through simultaneous analysis of molecular and morphological characters in both Bayesian and parsimony frameworks, we recovered the fossil taxon as sister-group to extant Leptomyrmex in Brazil while considering the influence of taxonomic and character sampling on inferred hypotheses relating to tree topology, biogeography and morphological evolution. We also identified potential loss of signal in the binning of morphological characters and tested the impact of parameterisation on divergence date estimation. Our results highlight the importance of securing sufficient taxon sampling for extant lineages when incorporating fossils and underscore the utility of diverse character sources in accurate placement of fossil terminals. Specifically, we find that fossil placement in this group is influenced by the inclusion of male-based characters and the newly discovered Neotropical ‘Lazarus taxon’.

Genotyping-by-sequencing provides the discriminating power to investigate the subspecies of Daucus carota (Apiaceae)

BACKGROUND

The majority of the subspecies of Daucus carota have not yet been discriminated clearly by various molecular or morphological methods and hence their phylogeny and classification remains unresolved. Recent studies using 94 nuclear orthologs and morphological characters, and studies employing other molecular approaches were unable to distinguish clearly many of the subspecies. Fertile intercrosses among traditionally recognized subspecies are well documented. We here explore the utility of single nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing (GBS) to serve as an effective molecular method to discriminate the subspecies of the D. carota complex.

RESULTS

We used GBS to obtain SNPs covering all nine Daucus carota chromosomes from 162 accessions of Daucus and two related genera. To study Daucus phylogeny, we scored a total of 10,814 or 38,920 SNPs with a maximum of 10 or 30 % missing data, respectively. To investigate the subspecies of D. carota, we employed two data sets including 150 accessions: (i) rate of missing data 10 % with a total of 18,565 SNPs, and (ii) rate of missing data 30 %, totaling 43,713 SNPs. Consistent with prior results, the topology of both data sets separated species with 2n = 18 chromosome from all other species. Our results place all cultivated carrots (D. carota subsp. sativus) in a single clade. The wild members of D. carota from central Asia were on a clade with eastern members of subsp. sativus. The other subspecies of D. carota were in four clades associated with geographic groups: (1) the Balkan Peninsula and the Middle East, (2) North America and Europe, (3) North Africa exclusive of Morocco, and (4) the Iberian Peninsula and Morocco. Daucus carota subsp. maximus was discriminated, but neither it, nor subsp. gummifer (defined in a broad sense) are monophyletic.

CONCLUSIONS

Our study suggests that (1) the morphotypes identified as D. carota subspecies gummifer (as currently broadly circumscribed), all confined to areas near the Atlantic Ocean and the western Mediterranean Sea, have separate origins from sympatric members of other subspecies of D. carota, (2) D. carota subsp. maximus, on two clades with some accessions of subsp. carota, can be distinguished from each other but only with poor morphological support, (3) D. carota subsp. capillifolius, well distinguished morphologically, is an apospecies relative to North African populations of D. carota subsp. carota, (4) the eastern cultivated carrots have origins closer to wild carrots from central Asia than to western cultivated carrots, and (5) large SNP data sets are suitable for species-level phylogenetic studies in Daucus.

A phylogenetic test of sympatric speciation in the Hydrobatinae (Aves: Procellariiformes)

Phylogenetic relationships among species can provide insight into how new species arise. For example, careful consideration of both the phylogenetic and geographic distributions of species in a group can reveal the geographic models of speciation within the group. One such model, sympatric speciation, may be more common than previously thought. The Hydrobatinae (Aves: Procellariformes) is a diverse subfamily of Northern Hemisphere storm-petrels for which the taxonomy is unclear. Previous studies showed that Hydrobates (formally Oceanodroma) castro breeding in the Azores during the cool season is sister species to H. monteiroi, a hot season breeder at the same locations, which suggests sympatric speciation by allochrony. To test whether other species within the subfamily arose via sympatric speciation by allochrony, we sequenced the cytochrome b gene and five nuclear introns to estimate a phylogenetic tree using multispecies coalescent methods, and to test whether species breeding in the same geographic area are monophyletic. We found that speciation within the Hydrobatinae appears to have followed several geographic modes of divergence. Sympatric seasonal species in Japan likely did not arise through sympatric speciation, but allochrony may have played a role in the divergence of H. matsudairae, a cool season breeder, and H. monorhis, a hot season breeder. No other potential cases of sympatric speciation were discovered within the subfamily. Despite breeding in the same geographic area, hydrobatine storm-petrels breeding in Baja California (H. microsoma and H. melania) are each sister to a species breeding off the coast of Peru (H. tethys and H. markhami, respectively). In fact, antitropical sister species appear to have diverged at multiple times, suggesting allochronic divergence might be common. In addition, allopatry has likely played a role in divergence of H. furcata, a north Pacific breeder, and H. pelagius, a north Atlantic breeder. This study demonstrates that a variety of mechanisms of divergence have played a role in generating the diversity of the Hydrobatinae and supports the current taxonomy of the subfamily.

Advances in Time Estimation Methods for Molecular Data

Molecular dating has become central to placing a temporal dimension on the tree of life. Methods for estimating divergence times have been developed for over 50 years, beginning with the proposal of molecular clock in 1962. We categorize the chronological development of these methods into four generations based on the timing of their origin. In the first generation approaches (1960s-1980s), a strict molecular clock was assumed to date divergences. In the second generation approaches (1990s), the equality of evolutionary rates between species was first tested and then a strict molecular clock applied to estimate divergence times. The third generation approaches (since ∼2000) account for differences in evolutionary rates across the tree by using a statistical model, obviating the need to assume a clock or to test the equality of evolutionary rates among species. Bayesian methods in the third generation require a specific or uniform prior on the speciation-process and enable the inclusion of uncertainty in clock calibrations. The fourth generation approaches (since 2012) allow rates to vary from branch to branch, but do not need prior selection of a statistical model to describe the rate variation or the specification of speciation model. With high accuracy, comparable to Bayesian approaches, and speeds that are orders of magnitude faster, fourth generation methods are able to produce reliable timetrees of thousands of species using genome scale data. We found that early time estimates from second generation studies are similar to those of third and fourth generation studies, indicating that methodological advances have not fundamentally altered the timetree of life, but rather have facilitated time estimation by enabling the inclusion of more species. Nonetheless, we feel an urgent need for testing the accuracy and precision of third and fourth generation methods, including their robustness to misspecification of priors in the analysis of large phylogenies and data sets.