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

Phylogenomic analyses of the pantropical Platycerium Desv. (Platycerioideae) reveal their complex evolution and historical biogeography

Platycerium is a genus of pantropical epiphytic ferns consisting of ca. 18 species and are highly sought after by horticultural enthusiasts. Although the monophyly of this genus has been well supported in previous molecular studies, as an intercontinentally disjunct genus, the origin and distribution pattern of Platycerium were elusive and controversial. This is mainly due to limited taxon sampling, a plastid representing only a single coalescent history, the lack of fossil evidence, and so on. Here, by utilizing genome-skimming sequencing, transcriptome sequencing, and flow cytometry, we integrated chloroplast genomes, data of single-copy nuclear genes, ploidy levels, morphology, and geographic distribution to understand the species phylogeny and the evolutionary and biogeographic history of Platycerium. Our major results include: (1) based on both plastid and nuclear datasets, Platycerium is consistently resolved into three fully supported clades: the Afro-American (AA) clade, the Javan-Australian (JA) clade, and the Malayan-Asian (MA) clade. The AA clade and MA clade are further divided into three and two subclades, respectively; (2) a large amount of gene tree conflict, as well as cytonuclear discordance, was found and can be explained by hybridization and incomplete lineage sorting, and most of the hybridization hypotheses represented ancient hybridization events; (3) through molecular dating, the crown age of Platycerium is determined to be at approximately 32.79 Ma based on the plastid dataset or 29.08 Ma based on the nuclear dataset in the Middle Oligocene; (4) ancestral area reconstruction analysis from different datasets showed that Platycerium most likely originated from Indochina; (5) current distribution patterns are resultant from long-distance dispersals, ancient orogeny, and an ancient climate event; and (6) species diversification was driven by polyploidization, dispersal, and hybridization. This study presented here will help understand the evolution of tropical plant flora and provide a reference for the cultivation and breeding of staghorn ferns.

Wing pattern diversity in Eunica butterflies (Nymphalidae: Biblidinae): phylogenetic analysis implies decoupled adaptive trends in dorsal sexual dimorphism and ventral eyespot evolution

Butterfly eyespots are wing patterns reminiscent of vertebrate eyes, formed by concentric rings of contrastingly coloured scales. Eyespots are usually located close to the wing margin and often regarded as the single most conspicuous pattern element of butterfly wing colour displays. Recent efforts to understand the processes involved in the formation of eyespots have been driven mainly by evo-devo approaches focused on model species. However, patterns of change implied by phylogenetic relationships can also inform hypotheses about the underlying developmental mechanisms associated with the formation or disappearance of eyespots, and the limits of phenotypic diversity occurring in nature. Here we present a combined evidence phylogenetic hypothesis for the genus Eunica, a prominent member of diverse Neotropical butterfly communities, that features notable variation among species in eyespot patterns on the ventral hind wing surface. The data matrix consists of one mitochondrial gene region (COI), four nuclear gene regions (GAPDH, RPS5, EF1a and Wingless) and 68 morphological characters. A combined cladistic analysis with all the characters concatenated produced a single most parsimonious tree that, although fully resolved, includes many nodes with modest branch support. The phylogenetic hypothesis presented corroborates a previously proposed morphological trend leading to the loss of eyespots, together with an increase in the size of the conserved eyespots, relative to outgroup taxa. Furthermore, wing colour pattern dimorphism and the presence of androconia suggest that the most remarkable instances of sexual dimorphism are present in the species of Eunica with the most derived eyespot patterns, and are in most cases accompanied by autapomorphic combinations of scent scales and "hair pencils". We discuss natural and sexual selection as potential adaptive explanations for dorsal and ventral wing patterns.

Ancient hybridization in Curcuma (Zingiberaceae)-Accelerator or brake in lineage diversifications?

Hybridization is a widespread phenomenon in the evolution of plants and exploring its role is crucial to understanding diversification processes of many taxonomic groups. Recently, more attention is focused on the role of ancient hybridization that has repeatedly been shown as triggers of evolutionary radiation, although in some cases, it can prevent further diversification. The causes, frequency, and consequences of ancient hybridization remain to be explored. Here, we present an account of several events of ancient hybridization in turmeric, the economically important plant genus Curcuma (Zingiberaceae), which harbors about 130 known species. We analyzed 1094 targeted low-copy genes and plastomes obtained by next-generation sequencing of 37 species of Curcuma, representing the known genetic diversity and spanning the geographical distribution of the genus. Using phylogenetic network analysis, we show that the entire genus Curcuma as well as its most speciose lineage arose via introgression from the genus Pyrgophyllum and one of the extinct lineages, respectively. We also document a single event of ancient hybridization, with C. vamana as a product, that represents an evolutionary dead end. We further discuss distinct circumstances of those hybridization events that deal mainly with (in)congruence in chromosome counts of the parental lineages.

Response to Tcherkez and Farquhar: Rubisco adaptation is more limited by phylogenetic constraint than by catalytic trade-off

Rubisco is the primary entry point for carbon into the biosphere. It has been widely proposed that rubisco is highly constrained by catalytic trade-offs due to correlations between the enzyme's kinetic traits across species. In previous work, we have shown that the strength of these correlations, and thus the strength of catalytic trade-offs, have been overestimated due to the presence of phylogenetic signal in the kinetic trait data (Bouvier et al., 2021). We demonstrated that only the trade-offs between the Michaelis constant for CO2 and carboxylase turnover, and between the Michaelis constants for CO2 and O2 were robust to phylogenetic effects. We further demonstrated that phylogenetic constraints have limited rubisco adaptation to a greater extent than the combined action of catalytic trade-offs. Recently, however, our claims have been contested by Tcherkez and Farquhar (2021), who have argued that the phylogenetic signal we detect in rubisco kinetic traits is an artefact of species sampling, the use of rbcL-based trees for phylogenetic inference, laboratory-to-laboratory variability in kinetic measurements, and homoplasy of the C4 trait. In the present article, we respond to these criticisms on a point-by-point basis and conclusively show that all are unfounded. As such, we stand by our original conclusions. Namely, although rubisco kinetic evolution has been limited by biochemical trade-offs, these are not absolute and have been previously overestimated due to phylogenetic biases. Instead, rubisco adaptation has in fact been more limited by phylogenetic constraint.

Weighted ASTRID: fast and accurate species trees from weighted internode distances

BACKGROUND

Species tree estimation is a basic step in many biological research projects, but is complicated by the fact that gene trees can differ from the species tree due to processes such as incomplete lineage sorting (ILS), gene duplication and loss (GDL), and horizontal gene transfer (HGT), which can cause different regions within the genome to have different evolutionary histories (i.e., "gene tree heterogeneity"). One approach to estimating species trees in the presence of gene tree heterogeneity resulting from ILS operates by computing trees on each genomic region (i.e., computing "gene trees") and then using these gene trees to define a matrix of average internode distances, where the internode distance in a tree T between two species x and y is the number of nodes in T between the leaves corresponding to x and y. Given such a matrix, a tree can then be computed using methods such as neighbor joining. Methods such as ASTRID and NJst (which use this basic approach) are provably statistically consistent, very fast (low degree polynomial time) and have had high accuracy under many conditions that makes them competitive with other popular species tree estimation methods. In this study, inspired by the very recent work of weighted ASTRAL, we present weighted ASTRID, a variant of ASTRID that takes the branch uncertainty on the gene trees into account in the internode distance.

RESULTS

Our experimental study evaluating weighted ASTRID typically shows improvements in accuracy compared to the original (unweighted) ASTRID, and shows competitive accuracy against weighted ASTRAL, the state of the art. Our re-implementation of ASTRID also improves the runtime, with marked improvements on large datasets.

CONCLUSIONS

Weighted ASTRID is a new and very fast method for species tree estimation that typically improves upon ASTRID and has comparable accuracy to weighted ASTRAL, while remaining much faster. Weighted ASTRID is available at https://github.com/RuneBlaze/internode .

Plastid phylogenomics provides novel insights into the infrafamilial relationship of Polypodiaceae

The polygrammoids (Polypodiaceae) are the most species-rich and diversified epiphytic fern lineages, and hold an important role to understand the deep diverging events and rapid adaptation to changing environments in the plant tree of life. Despite progress in the phylogeny of this group of ferns in previous multilocus phylogenetic studies, uncertainty remains especially in backbone relationships among closely related clades, and the phylogenetic placement of recalcitrant species or lineages. Here, we investigated the deep phylogenetic relationships within Polypodiaceae by sampling all major lineages and using 81 plastid genomes (plastomes), of which 70 plastomes were newly sequenced with high-throughput sequencing technology. Based on parsimony, maximum-likelihood, Bayesian and multispecies coalescent analyses of genome skimming data, we achieved a better resolution of the backbone phylogeny of Polypodiaceae. Using simulated data matrices, we detected that potential phylogenetic artefacts, such as long-branch attraction and insufficient taxonomic sampling, may have a confounding impact on the incongruence of phylogenetic inferences. Furthermore, our phylogenetic analyses offer greater resolution than previous multilocus studies, providing a robust framework for future phylogenetic implications on the subfamilial taxonomy of Polypodiaceae. Our phylogenomic study not only demonstrates the advantage of a character-rich plastome dataset for resolving the recalcitrant lineages that have undergone rapid radiation, but also sheds new light on integrative explorations understanding the evolutionary history of large fern groups in the genomic era.

Sweet or salty? The origin of freshwater gastrotrichs (Gastrotricha, Chaetonotida) revealed by molecular phylogenetic analysis

Chaetonotidae is the most diverse and widely distributed family of the order Chaetonotida (Gastrotricha) and includes both marine and freshwater species. Although the family is regarded as a sister taxon to the exclusively marine Xenotrichulidae, the type of environment, marine or freshwater, where Chaetonotidae originated is still not known. Here, we reconstructed the phylogeny of the family based on molecular sequence data and mapped both morphological and ecological characters to determine the ancestral environment of the first members of the family. Our results revealed that the freshwater genus Bifidochaetus is the earliest branching lineage in the paraphyletic Chaetonotidae (encompassing Dasydytidae and Neogosseidae). Moreover, we reconstructed Lepidochaetus-Cephalionotus clade as a monophyletic sister group to the remaining chaetonotids, which supports Kisielewski's morphological based hypothesis concerning undifferentiated type of body scales as a most primary character in Chaetonotidae. We also found that reversals to marine habitats occurred independently in different Chaetonotidae lineages, thus marine species in the genera Heterolepidoderma, Halichaetonotus, Aspidiophorus and subgenera Chaetonotus (Schizochaetonotus) or Chaetonotus (Marinochaetus) should be assumed as having secondarily invaded the marine environment. Character mapping revealed a series of synapomorphies that define the clade that includes Chaetonotidae (with Dasydytidae and Neogosseidae), the most important of which may be those linked to reproduction.

Miocene Climate and Habitat Change Drove Diversification in Bicyclus, Africa's Largest Radiation of Satyrine Butterflies

Compared to other regions, the drivers of diversification in Africa are poorly understood. We studied a radiation of insects with over 100 species occurring in a wide range of habitats across the Afrotropics to investigate the fundamental evolutionary processes and geological events that generate and maintain patterns of species richness on the continent. By investigating the evolutionary history of Bicyclus butterflies within a phylogenetic framework, we inferred the group's origin at the Oligo-Miocene boundary from ancestors in the Congolian rainforests of central Africa. Abrupt climatic fluctuations during the Miocene (ca. 19-17 Ma) likely fragmented ancestral populations, resulting in at least eight early-divergent lineages. Only one of these lineages appears to have diversified during the drastic climate and biome changes of the early Miocene, radiating into the largest group of extant species. The other seven lineages diversified in forest ecosystems during the late Miocene and Pleistocene when climatic conditions were more favourable-warmer and wetter. Our results suggest changing Neogene climate, uplift of eastern African orogens, and biotic interactions might have had different effects on the various subclades of Bicyclus, producing one of the most spectacular butterfly radiations in Africa.

Testing the utility of alternative metrics of branch support to address the ancient evolutionary radiation of tunas, stromateoids, and allies (Teleostei: Pelagiaria)

The use of high-throughput sequencing technologies to produce genome-scale datasets was expected to settle some long-standing controversies across the Tree of Life, particularly in areas where short branches occur at deep timescales. Instead, these datasets have often yielded many well-supported but conflicting topologies, and highly variable gene-tree distributions. A variety of branch-support metrics beyond the nonparametric bootstrap are now available to assess how robust a phylogenetic hypothesis may be, as well as new methods to quantify gene-tree discordance. We applied multiple branch support metrics to an ancient group of marine fishes (Teleostei: Pelagiaria) whose interfamilial relationships have proven difficult to resolve due to a rapid accumulation of lineages very early in its history. We analyzed hundreds of loci including published UCE data and newly generated exonic data along with their flanking regions to represent all 16 extant families for more than 150 out of 284 valid species in the group. Branch support was lower for interfamilial relationships (except the SH-like aLRT and aBayes methods) regardless of the type of marker used. Several nodes that were highly supported with bootstrap had very low site and gene-tree concordance, revealing underlying conflict. Despite this conflict, we were able to identify four consistent interfamilial clades, each comprised of two or three families. Combining exons with their flanking regions also produced increased branch lengths in the deep branches of the pelagiarian tree. Our results demonstrate the limitations of employing current metrics of branch support and species-tree estimation when assessing the confidence of ancient evolutionary radiations and emphasize the necessity to embrace alternative measurements to explore phylogenetic uncertainty and discordance in phylogenomic datasets.

Molecular phylogenetics of sub-Saharan African natricine snakes, and the biogeographic origins of the Seychelles endemic Lycognathophis seychellensis

Phylogenetic relationships of sub-Saharan African natricine snakes are understudied and poorly understood, which in turn has precluded analyses of the historical biogeography of the Seychelles endemic Lycognathophis seychellensis. We inferred the phylogenetic relationships of Seychelles and mainland sub-Saharan natricines by analysing a multilocus DNA sequence dataset for three mitochondrial (mt) and four nuclear (nu) genes. The mainland sub-Saharan natricines and L. seychellensis comprise a well-supported clade. Two maximally supported sets of relationships within this clade are (Limnophis,Natriciteres) and (Afronatrix,(Hydraethiops,Helophis)). The relationships of L. seychellensis with respect to these two lineages are not clearly resolved by analysing concatenated mt and nu data. Analysed separately, nu data best support a sister relationship of L. seychellensis with (Afronatrix,(Hydraethiops,Helophis)) and mt data best support a sister relationship with all mainland sub-Saharan natricines. Methods designed to cope with incomplete lineage sorting strongly favour the former hypothesis. Genetic variation among up to 33 L. seychellensis from five Seychelles islands is low. Fossil calibrated divergence time estimates support an overseas dispersal of the L. seychellensis lineage to the Seychelles from mainland Africa ca. 43-25 million years before present (Ma), rather than this taxon being a Gondwanan relic.

Swimming through the sands of the Sahara and Arabian deserts: Phylogeny of sandfish skinks (Scincidae, Scincus) reveals a recent and rapid diversification

Large parts of the Sahara Desert and Arabia are covered by sand seas and sand dunes, which are inhabited by specialized animal communities. For example, many lizards have developed adaptations to life in loose sand, including sand-swimming behavior. The best-known sand swimmers of the Saharo-Arabia are the sandfish skinks (genus Scincus). Although there are currently only four Scincus species recognized, their phylogenetic relationships have not yet been addressed in detail. We use eight genetic markers (three mitochondrial, five nuclear) and a complete sampling of species to infer the relationships within the genus. We employ multiple phylogenetic approaches to reconstruct the evolutionary history of these skinks and to assess the level of reticulation at the onset of their radiation. Our results indicate the presence of five strongly supported species-level lineages, four represented by the currently recognized species and the fifth by S. scincus conirostris, which does not form a clade with S. scincus. Based on these results we elevate the Iranian and northern Arabian S. conirostris to the species level. The two Saharan species, S. albifasciatus and S. scincus, are sister in all analyses. Deeper relationships within the genus, however, remained largely unresolved despite the extensive genetic data set. This basal polytomy, together with the fact that we detected no sign of hybridization in the history of the genus, indicates that the diversification of the five Scincus species was rapid, burst-like, and not followed by secondary hybridization events. Divergence time estimations show a Middle Pliocene crown radiation of the genus (3.3 Mya). We hypothesize that the aridification of the Saharo-Arabia that began in the Late Miocene triggered the initial diversification of Scincus, and that the subsequent expansion of sand deserts enabled their dispersal over the large Saharan and Arabian range. We discuss the evolution of body form in sand swimming lizards and ponder how Scincus retained their fully limbed morphology despite being sand swimmers that are typically limbless.

Molecular phylogenetic inference of the howler monkey radiation (Primates: Alouatta)

Howler monkeys (Alouatta), comprising between nine and 14 species and ranging from southern Mexico to northern Argentina, are the most widely distributed platyrrhines. Previous phylogenetic studies of howlers have used chromosomal and morphological characters and a limited number of molecular markers; however, branching patterns conflict between studies or remain unresolved. We performed a new phylogenetic analysis of Alouatta using both concatenated and coalescent-based species tree approaches based on 14 unlinked non-coding intergenic nuclear regions. Our taxon sampling included five of the seven South American species (Alouatta caraya, Alouatta belzebul, Alouatta guariba, Alouatta seniculus, Alouatta sara) and the two recognized species from Mesoamerica (Alouatta pigra, Alouatta palliata). Similarly to previous studies, our phylogenies supported a Mesoamerican clade and a South American clade. For the South American howlers, both methods recovered the Atlantic Forest endemic A. guariba as sister to all remaining South American species, albeit with moderate support. Moreover, we found no support for the previously proposed sister relationship between A. guariba and A. belzebul. For the first time, a clade composed of A. sara and A. caraya was identified. The relationships among the other South American howlers, however, were not fully supported. Our estimates for divergence times within Alouatta are generally older compared to estimates in earlier studies. However, they conform to recent studies proposing a Miocene age for the Isthmus of Panama and for the uplift of the northern Andes. Our results also point to an early genetic isolation of A. guariba in the Atlantic Forest, in agreement with the hypothesis of biotic exchange across South American rain forests in the Miocene. Collectively, these findings contribute to a better understanding of the diversification processes among howler monkey species; however, they also suggest that further comprehension of the evolutionary history of the Alouatta radiation will rely on broadened taxonomic, geographic, and genomic sampling.

Do Alignment and Trimming Methods Matter for Phylogenomic (UCE) Analyses?

Alignment is a crucial issue in molecular phylogenetics because different alignment methods can potentially yield very different topologies for individual genes. But it is unclear if the choice of alignment methods remains important in phylogenomic analyses, which incorporate data from hundreds or thousands of genes. For example, problematic biases in alignment might be multiplied across many loci, whereas alignment errors in individual genes might become irrelevant. The issue of alignment trimming (i.e., removing poorly aligned regions or missing data from individual genes) is also poorly explored. Here, we test the impact of 12 different combinations of alignment and trimming methods on phylogenomic analyses. We compare these methods using published phylogenomic data from ultraconserved elements (UCEs) from squamate reptiles (lizards and snakes), birds, and tetrapods. We compare the properties of alignments generated by different alignment and trimming methods (e.g., length, informative sites, missing data). We also test whether these data sets can recover well-established clades when analyzed with concatenated (RAxML) and species-tree methods (ASTRAL-III), using the full data ($\sim $5000 loci) and subsampled data sets (10% and 1% of loci). We show that different alignment and trimming methods can significantly impact various aspects of phylogenomic data sets (e.g., length, informative sites). However, these different methods generally had little impact on the recovery and support values for well-established clades, even across very different numbers of loci. Nevertheless, our results suggest several "best practices" for alignment and trimming. Intriguingly, the choice of phylogenetic methods impacted the phylogenetic results most strongly, with concatenated analyses recovering significantly more well-established clades (with stronger support) than the species-tree analyses. [Alignment; concatenated analysis; phylogenomics; sequence length heterogeneity; species-tree analysis; trimming].

Complete mitochondrial genomes provide current refined phylogenomic hypotheses for relationships among ten Hirundo species

Hirundo is the most species-rich genus of the passerine swallow family (Hirundinidae) and has a cosmopolitan distribution. Here we report the complete, annotated mitochondrial genomes for 25 individuals from 10 of the 14 extant Hirundo species; these include representatives from four subspecies of the barn swallow, H. rustica. Mitogenomes were conserved in size, ranging from 18,500 to 18,700 base pairs. They all contained 13 protein-coding regions, 22 tRNAs, a control region, and large and small ribosomal subunits. Phylogenetic analysis resolved most of the relationships between the studied species and subspecies which were largely consistent with previously published trees. Several new relationships were observed within the phylogeny that could have only been discovered with the increased amount of genetic material. This study represents the largest Hirundo mitochondrial phylogeny to date, and could serve as a vital tool for other studies focusing on the evolution of the Hirundo genus.

Multilocus phylogeny of Paratelmatobiinae (Anura: Leptodactylidae) reveals strong spatial structure and previously unknown diversity in the Atlantic Forest hotspot

The Brazilian Atlantic Forest harbors high levels of anuran diversity and endemism, including several taxa restricted to small geographic ranges. Here, we provide a multilocus phylogeny for Paratelmatobiinae, a leptodactylid subfamily composed of small-ranged species distributed in the Brazilian Atlantic Forest and in the campo rupestre ecosystem. We performed Bayesian inference and maximum likelihood analyses using three mitochondrial and five nuclear markers, and a matrix comprising a broad taxonomic sampling. We then delimitated independently evolving lineages within the group. We recovered Paratelmatobiinae and each of its four genera as monophyletic and robustly supported. Five putatively new species included in our analyses were unambiguously supported in the phylogenetic trees and delimitation analyses. We also recovered other deeply divergent and geographically structured lineages within the four genera of Paratelmatobiinae. Our estimation of divergence times indicates that diversification in the subfamily began in the Eocene and continued until the Pleistocene. We discuss possible scenarios of diversification for the four genera of Paratelmatobiinae, and outline the implications of our findings for taxonomy and conservation.

Phylogeny of Diplazium (Athyriaceae) revisited: Resolving the backbone relationships based on plastid genomes and phylogenetic tree space analysis

Despite progress in resolving the phylogeny of twinsorus ferns (Diplazium) based on multilocus phylogenetic studies, uncertainty remains especially for deep, or backbone relationships among closely related clades, suggesting a classic case of rapid evolutionary radiation. Here, we investigated the deep phylogenetic relationships within Diplazium by sampling all major clades and using 51 plastid genomes (plastomes), of which 38 were newly sequenced with high-throughput sequencing technology, resulting more than 127,000 informative sites. Using parsimony, maximum likelihood and Bayesian analyses of plastome sequences, we largely resolved the backbone of the phylogeny of Diplazium with strong support. However, we also detected phylogenetic incongruence among different datasets and moderately to poorly supported relationships, particularly at several extremely short internal branches. By using phylogenetic tree space and topology-clustering analyses, we provide evidence that conflicting phylogenetic signals can be found across the trees estimated from individual chloroplast protein-coding genes, which may underlie the difficulty of systematics of Diplazium. Furthermore, our phylogenetic estimate offers more resolution over previous multilocus analyses, providing a framework for future taxonomic revisions of sectional classification of Diplazium. Our study demonstrates the advantage of a character-rich plastome dataset, combining the comparison of different phylogenetic methods, for resolving the recalcitrant lineages that have undergone rapid radiation and dramatic changes in evolutionary rates.

Multiple historical processes obscure phylogenetic relationships in a taxonomically difficult group (Lobariaceae, Ascomycota)

In the age of next-generation sequencing, the number of loci available for phylogenetic analyses has increased by orders of magnitude. But despite this dramatic increase in the amount of data, some phylogenomic studies have revealed rampant gene-tree discordance that can be caused by many historical processes, such as rapid diversification, gene duplication, or reticulate evolution. We used a target enrichment approach to sample 400 single-copy nuclear genes and estimate the phylogenetic relationships of 13 genera in the lichen-forming family Lobariaceae to address the effect of data type (nucleotides and amino acids) and phylogenetic reconstruction method (concatenation and species tree approaches). Furthermore, we examined datasets for evidence of historical processes, such as rapid diversification and reticulate evolution. We found incongruence associated with sequence data types (nucleotide vs. amino acid sequences) and with different methods of phylogenetic reconstruction (species tree vs. concatenation). The resulting phylogenetic trees provided evidence for rapid and reticulate evolution based on extremely short branches in the backbone of the phylogenies. The observed rapid and reticulate diversifications may explain conflicts among gene trees and the challenges to resolving evolutionary relationships. Based on divergence times, the diversification at the backbone occurred near the Cretaceous-Paleogene (K-Pg) boundary (65 Mya) which is consistent with other rapid diversifications in the tree of life. Although some phylogenetic relationships within the Lobariaceae family remain with low support, even with our powerful phylogenomic dataset of up to 376 genes, our use of target-capturing data allowed for the novel exploration of the mechanisms underlying phylogenetic and systematic incongruence.

The bizarre skull of Xenotyphlops sheds light on synapomorphies of Typhlopoidea

The emerging picture of non-monophyly of scolecophidian snakes is increasingly indicative that fossorial lifestyle, myrmecophagous diet, and miniaturisation are powerful drivers of morphological evolution in squamate skulls. We provide a detailed description of the skull of Xenotyphlops grandidieri, with reference to the skulls of other scolecophidian snakes. The skull, which shows dramatic ventral inflection of the snout complex, is remarkably bizarre, and the mouth opening is more ventrally oriented than in other typhlopoids. The eyes are strongly reduced, and the enlarged and rather flat anterior head shield is covered in numerous sensillae. We put forward several potential explanations for the evolution of these unusual modifications. On the other hand, Xenotyphlops shares numerous synapomorphies with other typhlopoid snakes, including the highly specialized jaw mechanism. We argue that the key differences between the jaw mechanisms of Leptotyphlopidae, Anomalepididae, and Typhlopoidea provide compelling evidence for a strong role of convergence in the evolution of the scolecophidian bauplan, and these clades therefore cannot be interpreted as representative of ancestral anatomy or ecology among snakes.

Evolution of nuchal glands, unusual defensive organs of Asian natricine snakes (Serpentes: Colubridae), inferred from a molecular phylogeny

A large body of evidence indicates that evolutionary innovations of novel organs have facilitated the subsequent diversification of species. Investigation of the evolutionary history of such organs should provide important clues for understanding the basis for species diversification. An Asian natricine snake, Rhabdophis tigrinus, possesses a series of unusual organs, called nuchal glands, which contain cardiotonic steroid toxins known as bufadienolides. Rhabdophis tigrinus sequesters bufadienolides from its toad prey and stores them in the nuchal glands as a defensive mechanism. Among more than 3,500 species of snakes, only 17 Asian natricine species are known to possess nuchal glands or their homologues. These 17 species belong to three nominal genera, Balanophis, Macropisthodon, and Rhabdophis. In Macropisthodon and Rhabdophis, however, species without nuchal glands also exist. To infer the evolutionary history of the nuchal glands, we investigated the molecular phylogenetic relationships among Asian natricine species with and without nuchal glands, based on variations in partial sequences of Mt-CYB, Cmos, and RAG1 (total 2,767 bp). Results show that all species with nuchal glands belong to a single clade (NGC). Therefore, we infer that the common ancestor of this clade possessed nuchal glands with no independent origins of the glands within the members. Our results also imply that some species have secondarily lost the glands. Given the estimated divergence time of related species, the ancestor of the nuchal gland clade emerged 19.18 mya. Our study shows that nuchal glands are fruitful subjects for exploring the evolution of novel organs. In addition, our analysis indicates that reevaluation of the taxonomic status of the genera Balanophis and Macropisthodon is required. We propose to assign all species belonging to the NGC to the genus Rhabdophis, pending further study.

Plastome sequences and exploration of tree-space help to resolve the phylogeny of riceflowers (Thymelaeaceae: Pimelea)

Data sets comprising small numbers of genetic markers are not always able to resolve phylogenetic relationships. This has frequently been the case in molecular systematic studies of plants, with many analyses being based on sequence data from only two or three chloroplast genes. An example of this comes from the riceflowers Pimelea Banks & Sol. ex Gaertn. (Thymelaeaceae), a large genus of flowering plants predominantly distributed in Australia. Despite the considerable morphological variation in the genus, low sequence divergence in chloroplast markers has led to the phylogeny of Pimelea remaining largely uncertain. In this study, we resolve the backbone of the phylogeny of Pimelea in comprehensive Bayesian and maximum-likelihood analyses of plastome sequences from 41 taxa. However, some relationships received only moderate to poor support, and the Pimelea clade contained extremely short internal branches. By using topology-clustering analyses, we demonstrate that conflicting phylogenetic signals can be found across the trees estimated from individual chloroplast protein-coding genes. A relaxed-clock dating analysis reveals that Pimelea arose in the mid-Miocene, with most divergences within the genus occurring during a subsequent rapid diversification. Our new phylogenetic estimate offers better resolution and is more strongly supported than previous estimates, providing a platform for future taxonomic revisions of both Pimelea and the broader subfamily. Our study has demonstrated the substantial improvements in phylogenetic resolution that can be achieved using plastome-scale data sets in plant molecular systematics.

Conflicting Evolutionary Histories of the Mitochondrial and Nuclear Genomes in New World Myotis Bats

The rapid diversification of Myotis bats into more than 100 species is one of the most extensive mammalian radiations available for study. Efforts to understand relationships within Myotis have primarily utilized mitochondrial markers and trees inferred from nuclear markers lacked resolution. Our current understanding of relationships within Myotis is therefore biased towards a set of phylogenetic markers that may not reflect the history of the nuclear genome. To resolve this, we sequenced the full mitochondrial genomes of 37 representative Myotis, primarily from the New World, in conjunction with targeted sequencing of 3648 ultraconserved elements (UCEs). We inferred the phylogeny and explored the effects of concatenation and summary phylogenetic methods, as well as combinations of markers based on informativeness or levels of missing data, on our results. Of the 294 phylogenies generated from the nuclear UCE data, all are significantly different from phylogenies inferred using mitochondrial genomes. Even within the nuclear data, quartet frequencies indicate that around half of all UCE loci conflict with the estimated species tree. Several factors can drive such conflict, including incomplete lineage sorting, introgressive hybridization, or even phylogenetic error. Despite the degree of discordance between nuclear UCE loci and the mitochondrial genome and among UCE loci themselves, the most common nuclear topology is recovered in one quarter of all analyses with strong nodal support. Based on these results, we re-examine the evolutionary history of Myotis to better understand the phenomena driving their unique nuclear, mitochondrial, and biogeographic histories.

Novel Approaches for Phylogenetic Inference from Morphological Data and Total-Evidence Dating in Squamate Reptiles (Lizards, Snakes, and Amphisbaenians)

Here, I combine previously underutilized models and priors to perform more biologically realistic phylogenetic inference from morphological data, with an example from squamate reptiles. When coding morphological characters, it is often possible to denote ordered states with explicit reference to observed or hypothetical ancestral conditions. Using this logic, we can integrate across character-state labels and estimate meaningful rates of forward and backward transitions from plesiomorphy to apomorphy. I refer to this approach as MkA, for “asymmetric.” The MkA model incorporates the biological reality of limited reversal for many phylogenetically informative characters, and significantly increases likelihoods in the empirical data sets. Despite this, the phylogeny of Squamata remains contentious. Total-evidence analyses using combined morphological and molecular data and the MkA approach tend toward recent consensus estimates supporting a nested Iguania. However, support for this topology is not unambiguous across data sets or analyses, and no mechanism has been proposed to explain the widespread incongruence between partitions, or the hidden support for various topologies in those partitions. Furthermore, different morphological data sets produced by different authors contain both different characters and different states for the same or similar characters, resulting in drastically different placements for many important fossil lineages. Effort is needed to standardize ontology for morphology, resolve incongruence, and estimate a robust phylogeny. The MkA approach provides a preliminary avenue for investigating morphological evolution while accounting for temporal evidence and asymmetry in character-state changes.

Evolutionary Dynamics of the W Chromosome in Caenophidian Snakes

The caenophidian (assigned also as “advanced”) snakes are traditionally viewed as a group of reptiles with a limited karyotypic variation and stable ZZ/ZW sex chromosomes. The W chromosomes of the caenophidian snakes are heterochromatic, and pioneering studies demonstrated that they are rich in repetitive elements. However, a comparative study of the evolutionary dynamics of the repetitive content of the W chromosome across the whole lineage is missing. Using molecular-cytogenetic techniques, we explored the distribution of four repetitive motifs (microsatellites GATA, GACA, AG and telomeric-like sequences), which are frequently accumulated in differentiated sex chromosomes in vertebrates, in the genomes of 13 species of the caenophidian snakes covering a wide phylogenetic spectrum of the lineage. The results demonstrate a striking variability in the morphology and the repetitive content of the W chromosomes even between closely-related species, which is in contrast to the homology and long-term stability of the gene content of the caenophidian Z chromosome. We uncovered that the tested microsatellite motifs are accumulated on the degenerated, heterochromatic W chromosomes in all tested species of the caenophidian snakes with the exception of the Javan file snake representing a basal clade. On the other hand, the presence of the accumulation of the telomeric-like sequences on the caenophidian W chromosome is evolutionary much less stable. Moreover, we demonstrated that large accumulations of telomeric-like motifs on the W chromosome contribute to sexual differences in the number of copies of the telomeric and telomeric-like repeats estimated by quantitative PCR, which might be confusing and incorrectly interpreted as sexual differences in telomere length.

Diversification rates and species richness across the Tree of Life

Species richness varies dramatically among clades across the Tree of Life, by over a million-fold in some cases (e.g. placozoans versus arthropods). Two major explanations for differences in richness among clades are the clade-age hypothesis (i.e. species-rich clades are older) and the diversification-rate hypothesis (i.e. species-rich clades diversify more rapidly, where diversification rate is the net balance of speciation and extinction over time). Here, we examine patterns of variation in diversification rates across the Tree of Life. We address how rates vary across higher taxa, whether rates within higher taxa are related to the subclades within them, and how diversification rates of clades are related to their species richness. We find substantial variation in diversification rates, with rates in plants nearly twice as high as in animals, and rates in some eukaryotes approximately 10-fold faster than prokaryotes. Rates for each kingdom-level clade are then significantly related to the subclades within them. Although caution is needed when interpreting relationships between diversification rates and richness, a positive relationship between the two is not inevitable. We find that variation in diversification rates seems to explain most variation in richness among clades across the Tree of Life, in contrast to the conclusions of previous studies.

Paleoclimate determines diversification patterns in the fossorial snake family Uropeltidae Cuvier, 1829

Understanding how and why diversification rates vary across evolutionary time is central to understanding how biodiversity is generated and maintained. Recent mathematical models that allow estimation of diversification rates across time from reconstructed phylogenies have enabled us to make inferences on how biodiversity copes with environmental change. Here, we explore patterns of temporal diversification in Uropeltidae, a diverse fossorial snake family. We generate a time-calibrated phylogenetic hypothesis for Uropeltidae and show a significant correlation between diversification rate and paleotemperature during the Cenozoic. We show that the temporal diversification pattern of this group is punctuated by one rate shift event with a decrease in diversification and turnover rate between ca. 11Ma to present, but there is no strong support for mass extinction events. The analysis indicates higher turnover during periods of drastic climatic fluctuations and reduced diversification rates associated with contraction and fragmentation of forest habitats during the late Miocene. Our study highlights the influence of environmental fluctuations on diversification rates in fossorial taxa such as uropeltids, and raises conservation concerns related to present rate of climate change.

Scent gland constituents of the Middle American burrowing python, Loxocemus bicolor (Serpentes: Loxocemidae)

Analysis by gas chromatography/mass spectrometry of the scent gland secretions of male and female Middle American burrowing pythons (Loxocemus bicolor) revealed the presence of over 300 components including cholesterol, fatty acids, glyceryl monoalkyl ethers, and alcohols. The fatty acids, over 100 of which were identified, constitute most of the compounds in the secretions and show the greatest structural diversity. They include saturated and unsaturated, unbranched and mono-, di-, and trimethyl-branched compounds ranging in carbon-chain length from 13 to 24. The glyceryl monoethers possess saturated or unsaturated, straight or methyl-branched alkyl chains ranging in carbon-chain length from 13 to 24. Alcohols, which have not previously been reported from the scent glands, possess straight, chiefly saturated carbon chains ranging in length from 13 to 24. Sex or individual differences in secretion composition were not observed. Compounds in the scent gland secretions of L. bicolor may deter offending arthropods, such as ants.

Variable wing venation in Agathiphaga (Lepidoptera: Agathiphagidae) is key to understanding the evolution of basal moths

Details of the ancestral groundplan of wing venation in moths remain uncertain, despite approximately a century of study. Here, we describe a 3-branched subcostal vein, a 5-branched medial vein and a 2-branched cubitus posterior vein on the forewing of Agathiphaga vitiensis Dumbleton 1952 from Vanuatu. Such veins had not previously been described in any Lepidoptera. Because wing veins are typically lost during lepidopteran evolutionary history, rarely-if ever-to be regained, the venation of A. vitiensis probably represents the ancestral character state for moths. Wing venation is often used to identify fossil insects as moths, because wing scales are not always preserved; the presence of a supposedly trichopteran 3-branched subcostal vein in crown Lepidoptera may decrease the certainty with which certain amphiesmenopteran fossils from the Mesozoic can be classified. And because plesiomorphic veins can influence the development of lepidopteran wing patterns even if not expressed in the adult wing, the veins described here may determine the location of wing pattern elements in many lepidopteran taxa.

A Species-Level Phylogeny of Extant Snakes with Description of a New Colubrid Subfamily and Genus

BACKGROUND

With over 3,500 species encompassing a diverse range of morphologies and ecologies, snakes make up 36% of squamate diversity. Despite several attempts at estimating higher-level snake relationships and numerous assessments of generic- or species-level phylogenies, a large-scale species-level phylogeny solely focusing on snakes has not been completed. Here, we provide the largest-yet estimate of the snake tree of life using maximum likelihood on a supermatrix of 1745 taxa (1652 snake species + 7 outgroup taxa) and 9,523 base pairs from 10 loci (5 nuclear, 5 mitochondrial), including previously unsequenced genera (2) and species (61).

RESULTS

Increased taxon sampling resulted in a phylogeny with a new higher-level topology and corroborate many lower-level relationships, strengthened by high nodal support values (> 85%) down to the species level (73.69% of nodes). Although the majority of families and subfamilies were strongly supported as monophyletic with > 88% support values, some families and numerous genera were paraphyletic, primarily due to limited taxon and loci sampling leading to a sparse supermatrix and minimal sequence overlap between some closely-related taxa. With all rogue taxa and incertae sedis species eliminated, higher-level relationships and support values remained relatively unchanged, except in five problematic clades.

CONCLUSION

Our analyses resulted in new topologies at higher- and lower-levels; resolved several previous topological issues; established novel paraphyletic affiliations; designated a new subfamily, Ahaetuliinae, for the genera Ahaetulla, Chrysopelea, Dendrelaphis, and Dryophiops; and appointed Hemerophis (Coluber) zebrinus to a new genus, Mopanveldophis. Although we provide insight into some distinguished problematic nodes, at the deeper phylogenetic scale, resolution of these nodes may require sampling of more slowly-evolving nuclear genes.

Colubrid Venom Composition: An -Omics Perspective

Snake venoms have been subjected to increasingly sensitive analyses for well over 100 years, but most research has been restricted to front-fanged snakes, which actually represent a relatively small proportion of extant species of advanced snakes. Because rear-fanged snakes are a diverse and distinct radiation of the advanced snakes, understanding venom composition among "colubrids" is critical to understanding the evolution of venom among snakes. Here we review the state of knowledge concerning rear-fanged snake venom composition, emphasizing those toxins for which protein or transcript sequences are available. We have also added new transcriptome-based data on venoms of three species of rear-fanged snakes. Based on this compilation, it is apparent that several components, including cysteine-rich secretory proteins (CRiSPs), C-type lectins (CTLs), CTLs-like proteins and snake venom metalloproteinases (SVMPs), are broadly distributed among "colubrid" venoms, while others, notably three-finger toxins (3FTxs), appear nearly restricted to the Colubridae (sensu stricto). Some putative new toxins, such as snake venom matrix metalloproteinases, are in fact present in several colubrid venoms, while others are only transcribed, at lower levels. This work provides insights into the evolution of these toxin classes, but because only a small number of species have been explored, generalizations are still rather limited. It is likely that new venom protein families await discovery, particularly among those species with highly specialized diets.

Phylogenomic analyses reveal novel relationships among snake families

Snakes are a diverse and important group of vertebrates. However, relationships among the major groups of snakes have remained highly uncertain, with recent studies hypothesizing very different (and typically weakly supported) relationships. Here, we address family-level snake relationships with new phylogenomic data from 3776 nuclear loci from ultraconserved elements (1.40million aligned base pairs, 52% missing data overall) sampled from 29 snake species that together represent almost all families, a dataset ∼100 times larger than used in previous studies. We found relatively strong support from species-tree analyses (NJst) for most relationships, including three largely novel clades: (1) a clade uniting the boas, pythons and their relatives, (2) a clade placing cylindrophiids and uropeltids with this clade, and (3) a clade uniting bolyeriids (Round Island boas) with pythonids and their relatives (xenopeltids and loxocemids). Relationships among families of advanced snakes (caenophidians) were also strongly supported. The results show the potential for phylogenomic analyses to resolve difficult groups, but also show a surprising sensitivity of the analyses to the inclusion or exclusion of outgroups.

Phylogenomics of a rapid radiation: is chromosomal evolution linked to increased diversification in north american spiny lizards (Genus Sceloporus)?

Background

Resolving the short phylogenetic branches that result from rapid evolutionary diversification often requires large numbers of loci. We collected targeted sequence capture data from 585 nuclear loci (541 ultraconserved elements and 44 protein-coding genes) to estimate the phylogenetic relationships among iguanian lizards in the North American genus Sceloporus. We tested for diversification rate shifts to determine if rapid radiation in the genus is correlated with chromosomal evolution.

Results

The phylogenomic trees that we obtained for Sceloporus using concatenation and coalescent-based species tree inference provide strong support for the monophyly and interrelationships among nearly all major groups. The diversification analysis supported one rate shift on the Sceloporus phylogeny approximately 20–25 million years ago that is associated with the doubling of the speciation rate from 0.06 species/million years (Ma) to 0.15 species/Ma. The posterior probability for this rate shift occurring on the branch leading to the Sceloporus species groups exhibiting increased chromosomal diversity is high (posterior probability = 0.997).

Conclusions

Despite high levels of gene tree discordance, we were able to estimate a phylogenomic tree for Sceloporus that solves some of the taxonomic problems caused by previous analyses of fewer loci. The taxonomic changes that we propose using this new phylogenomic tree help clarify the number and composition of the major species groups in the genus. Our study provides new evidence for a putative link between chromosomal evolution and the rapid divergence and radiation of Sceloporus across North America.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0628-x) contains supplementary material, which is available to authorized users.

The effects of subsampling gene trees on coalescent methods applied to ancient divergences

Gene-tree-estimation error is a major concern for coalescent methods of phylogenetic inference. We sampled eight empirical studies of ancient lineages with diverse numbers of taxa and genes for which the original authors applied one or more coalescent methods. We found that the average pairwise congruence among gene trees varied greatly both between studies and also often within a study. We recommend that presenting plots of pairwise congruence among gene trees in a dataset be treated as a standard practice for empirical coalescent studies so that readers can readily assess the extent and distribution of incongruence among gene trees. ASTRAL-based coalescent analyses generally outperformed MP-EST and STAR with respect to both internal consistency (congruence between analyses of subsamples of genes with the complete dataset of all genes) and congruence with the concatenation-based topology. We evaluated the approach of subsampling gene trees that are, on average, more congruent with other gene trees as a method to reduce artifacts caused by gene-tree-estimation errors on coalescent analyses. We suggest that this method is well suited to testing whether gene-tree-estimation error is a primary cause of incongruence between concatenation- and coalescent-based results, to reconciling conflicting phylogenetic results based on different coalescent methods, and to identifying genes affected by artifacts that may then be targeted for reciprocal illumination. We provide scripts that automate the process of calculating pairwise gene-tree incongruence and subsampling trees while accounting for differential taxon sampling among genes. Finally, we assert that multiple tree-search replicates should be implemented as a standard practice for empirical coalescent studies that apply MP-EST.