Spurious 99% bootstrap and jackknife support for unsupported clades

Phylogenomic data resolved the deep relationships of Gymnogynoideae (Selaginellaceae)

The unresolved phylogenetic framework within the Selaginellaceae subfamily Gymnogynoideae (ca. 130 species) has hindered our comprehension of the diversification and evolution of Selaginellaceae, one of the most important lineages in land plant evolution. Here, based on plastid and nuclear data extracted from genomic sequencing of more than 90% species of all genera except two in Gymnogynoideae, a phylogenomic study focusing on the contentious relationships among the genera in Gymnogynoideae was conducted. Our major results included the following: (1) Only single-copy region (named NR) and only one ribosomal operon was firstly found in Afroselaginella among vascular plants, the plastome structure of Gymnogynoideae is diverse among the six genera, and the direct repeats (DR) type is inferred as the ancestral state in the subfamily; (2) The first strong evidence was found to support Afroselaginella as a sister to Megaloselaginella. Alternative placements of Ericetorum and Gymnogynum were detected, and their relationships were investigated by analyzing the variation of phylogenetic signals; and (3) The most likely genus-level relationships in Gymnogynoideae might be: ((Bryodesma, Lepidoselaginella), (((Megaloselaginella, Afroselaginella), Ericetorum), Gymnogynum)), which was supported by maximum likelihood phylogeny based on plastid datasets, maximum likelihood, and Bayesian inference based on SCG dataset and concatenated nuclear and plastid datasets and the highest proportion of phylogenetic signals of plastid genes.

Plastid phylogenomics and fossil evidence provide new insights into the evolutionary complexity of the 'woody clade' in Saxifragales

BACKGROUND

The "woody clade" in Saxifragales (WCS), encompassing four woody families (Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae), is a phylogenetically recalcitrant node in the angiosperm tree of life, as the interfamilial relationships of the WCS remain contentious. Based on a comprehensive sampling of WCS genera, this study aims to recover a robust maternal backbone phylogeny of the WCS by analyzing plastid genome (plastome) sequence data using Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods, and to explore the possible causes of the phylogenetic recalcitrance with respect to deep relationships within the WCS, in combination with molecular and fossil evidence.

RESULTS

Although the four WCS families were identically resolved as monophyletic, the MP analysis recovered different tree topologies for the relationships among Altingiaceae, Cercidiphyllaceae, and Daphniphyllaceae from the ML and BI phylogenies. The fossil-calibrated plastome phylogeny showed that the WCS underwent a rapid divergence of crown groups in the early Cretaceous (between 104.79 and 100.23 Ma), leading to the origin of the stem lineage ancestors of Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae within a very short time span (∼4.56 Ma). Compared with the tree topology recovered in a previous study based on nuclear genome data, cytonuclear discordance regarding the interfamilial relationships of the WCS was detected.

CONCLUSIONS

Molecular and fossil evidence imply that the early divergence of the WCS might have experienced radiative diversification of crown groups, extensive extinctions at the genus and species levels around the Cretaceous/Paleocene boundary, and ancient hybridization. Such evolutionarily complex events may introduce biases in topological estimations within the WCS due to incomplete lineage sorting, cytonuclear discordance, and long-branch attraction, potentially impacting the accurate reconstruction of deep relationships.

Phylogeny and biogeography of the northern temperate genus Dracocephalum s.l. (Lamiaceae)

The northern temperate genus Dracocephalum consists of approximately 70 species mainly distributed in the steppe-desert biomes of Central and West Asia and the alpine region of the Qinghai-Tibetan Plateau (QTP). Previous work has shown that Dracocephalum is not monophyletic and might include Hyssopus and Lallemantia. This study attempts to clarify the phylogenetic relationships, diversification patterns, and the biogeographical history of the three genera (defined as Dracocephalum s.l.). Based on a sampling of 66 taxa comprising more than 80% from extant species of Dracocephalum s.l., morphological, phylogenetic (maximum parsimony, likelihood, and Bayesian inference based on nuclear ITS and ETS, plastid rpl32-trnL, trnL-trnF, ycf1, and ycf1-rps15, and two low-copy nuclear markers AT3G09060 and AT1G09680), molecular dating, diversification, and ancestral range estimation analyses were carried out. Our results demonstrate that both Hyssopus and Lallemantia are embedded within Dracocephalum and nine well-supported clades can be recognized within Dracocephalum s.l. Analyses of divergence times suggest that the genus experienced an early rapid radiation during the middle to late Miocene with major lineages diversifying within a relatively narrow timescale. Ancestral area reconstruction analyses indicate that Dracocephalum s.l. originated in Central and West Asia and southern Siberia, and dispersed from Central and West Asia into the QTP and adjacent areas twice independently during the Pliocene. The aridification of the Asian interior possibly promoted the rapid radiation of Dracocephalum within this region, and the uplift of the QTP appears to have triggered the dispersal and recent rapid diversification of the genus in the QTP and adjacent regions. Combining molecular phylogenetic and morphological evidence, a revised infrageneric classification of Dracocephalum s.l. is proposed, which recognizes nine sections within the genus.

An empirical test of the relationship between the bootstrap and likelihood ratio support in maximum likelihood phylogenetic analysis

In maximum likelihood (ML), the support for a clade can be calculated directly as the likelihood ratio (LR) or log-likelihood difference (S, LLD) of the best trees with and without the clade of interest. However, bootstrap (BS) clade frequencies are more pervasive in ML phylogenetics and are almost universally interpreted as measuring support. In addition to theoretical arguments against that interpretation, BS has several undesirable attributes for a support measure. For example, it does not vary in proportion to optimality or identify clades that are rejected by the evidence and can be overestimated due to missing data. Nevertheless, if BS is a reliable predictor of S, then it might be an efficient indirect method of measuring support-an attractive possibility, given the speed of many BS implementations. To assess the relationship between S and BS, we analyzed 106 empirical datasets retrieved from TreeBASE. Also, to evaluate the degree to which S and BS are affected by the number of replicates during suboptimal tree searches for S and pseudoreplicates during BS estimation, we randomly selected 5 of the 106 datasets and analyzed them using variable numbers of replicates and pseudoreplicates, respectively. The correlation between S and BS was extremely weak in the datasets we analyzed. Increasing the number of replicates during tree search decreased the estimated values of S for most clades, but the magnitude of change was small. In contrast, although increasing pseudoreplicates affected BS values for only approximately 40% of clades, values both increased and decreased, and they did so at much greater magnitudes. Increasing replicates/pseudoreplicates affected the rank order of clades in each tree for both S and BS. Our findings show decisively that BS is not an efficient indirect method of measuring support and suggest that even quite superficial searches to calculate S provide better estimates of support.

The FUSED LEAVES1-ADHERENT1 regulatory module is required for maize cuticle development and organ separation

All aerial epidermal cells in land plants are covered by the cuticle, an extracellular hydrophobic layer that provides protection against abiotic and biotic stresses and prevents organ fusion during development. Genetic and morphological analysis of the classic maize adherent1 (ad1) mutant was combined with genome-wide binding analysis of the maize MYB transcription factor FUSED LEAVES1 (FDL1), coupled with transcriptional profiling of fdl1 mutants. We show that AD1 encodes an epidermally-expressed 3-KETOACYL-CoA SYNTHASE (KCS) belonging to a functionally uncharacterized clade of KCS enzymes involved in cuticular wax biosynthesis. Wax analysis in ad1 mutants indicates that AD1 functions in the formation of very-long-chain wax components. We demonstrate that FDL1 directly binds to CCAACC core motifs present in AD1 regulatory regions to activate its expression. Over 2000 additional target genes of FDL1, including many involved in cuticle formation, drought response and cell wall organization, were also identified. Our results identify a regulatory module of cuticle biosynthesis in maize that is conserved across monocots and eudicots, and highlight previously undescribed factors in lipid metabolism, transport and signaling that coordinate organ development and cuticle formation.

The FUSED LEAVES1-ADHERENT1 regulatory module is required for maize cuticle development and organ separation

All aerial epidermal cells in land plants are covered by the cuticle, an extracellular hydrophobic layer that provides protection against abiotic and biotic stresses and prevents organ fusion during development. Genetic and morphological analysis of the classic maize adherent1 (ad1) mutant was combined with genome-wide binding analysis of the maize MYB transcription factor FUSED LEAVES1 (FDL1), coupled with transcriptional profiling of fdl1 mutants. We show that AD1 encodes an epidermally-expressed 3-KETOACYL-CoA SYNTHASE (KCS) belonging to a functionally uncharacterized clade of KCS enzymes involved in cuticular wax biosynthesis. Wax analysis in ad1 mutants indicates that AD1 functions in the formation of very-long-chain wax components. We demonstrate that FDL1 directly binds to CCAACC core motifs present in AD1 regulatory regions to activate its expression. Over 2000 additional target genes of FDL1, including many involved in cuticle formation, drought response and cell wall organization, were also identified. Our results identify a regulatory module of cuticle biosynthesis in maize that is conserved across monocots and eudicots, and highlight previously undescribed factors in lipid metabolism, transport and signaling that coordinate organ development and cuticle formation.

Phylogeny and diversification history of the large Neotropical genus Philodendron (Araceae): Accelerated speciation in a lineage dominated by epiphytes

PREMISE OF THE STUDY

Philodendron is a large genus of ~560 species and among the most conspicuous epiphytic components of Neotropical forests, yet its phylogenetic relationships, timing of divergence, and diversification history have remained unclear. We present a comprehensive phylogenetic study for Philodendron and investigate its diversification, including divergence-time estimates and diversification rate shift analyses.

METHODS

We performed the largest phylogenetic reconstruction for Philodendron to date, including 125 taxa with a combined dataset of three plastid regions (petD, rpl16, and trnK/matK). We estimated divergence times using Bayesian evolutionary analysis sampling trees and inferred shifts in diversification rates using Bayesian analysis of macroevolutionary mixtures.

KEY RESULTS

We found that Philodendron, its three subgenera, and the closely related genus Adelonema are monophyletic. Within Philodendron subgenus Philodendron, 12 statistically well-supported clades are recognized. The genus Philodendron originated ~25 mya and a diversification rate upshift was detected at the origin of subgenus Philodendron ~12 mya.

CONCLUSIONS

Philodendron is a species-rich Neotropical lineage that diverged from Adelonema during the late Oligocene. Within Philodendron, the three subgenera currently accepted are recovered in two lineages: one contains the subgenera Meconostigma and Pteromischum and the other contains subgenus Philodendron. The lineage containing subgenera Meconostigma and Pteromischum underwent a consistent diversification rate. By contrast, a diversification rate upshift occurred within subgenus Philodendron ~12 mya. This diversification rate upshift is associated with the species radiation of the most speciose subgenus within Philodendron. The sections accepted within subgenus Philodendron are not congruent with the clades recovered. Instead, the clades are geographically defined.

The prevalence of terraced treescapes in analyses of phylogenetic data sets

BACKGROUND

The pattern of data availability in a phylogenetic data set may lead to the formation of terraces, collections of equally optimal trees. Terraces can arise in tree space if trees are scored with parsimony or with partitioned, edge-unlinked maximum likelihood. Theory predicts that terraces can be large, but their prevalence in contemporary data sets has never been surveyed. We selected 26 data sets and phylogenetic trees reported in recent literature and investigated the terraces to which the trees would belong, under a common set of inference assumptions. We examined terrace size as a function of the sampling properties of the data sets, including taxon coverage density (the proportion of taxon-by-gene positions with any data present) and a measure of gene sampling "sufficiency". We evaluated each data set in relation to the theoretical minimum gene sampling depth needed to reduce terrace size to a single tree, and explored the impact of the terraces found in replicate trees in bootstrap methods.

RESULTS

Terraces were identified in nearly all data sets with taxon coverage densities < 0.90. They were not found, however, in high-coverage-density (i.e., ≥ 0.94) transcriptomic and genomic data sets. The terraces could be very large, and size varied inversely with taxon coverage density and with gene sampling sufficiency. Few data sets achieved a theoretical minimum gene sampling depth needed to reduce terrace size to a single tree. Terraces found during bootstrap resampling reduced overall support.

CONCLUSIONS

If certain inference assumptions apply, trees estimated from empirical data sets often belong to large terraces of equally optimal trees. Terrace size correlates to data set sampling properties. Data sets seldom include enough genes to reduce terrace size to one tree. When bootstrap replicate trees lie on a terrace, statistical support for phylogenetic hypotheses may be reduced. Although some of the published analyses surveyed were conducted with edge-linked inference models (which do not induce terraces), unlinked models have been used and advocated. The present study describes the potential impact of that inference assumption on phylogenetic inference in the context of the kinds of multigene data sets now widely assembled for large-scale tree construction.

Why are there so many sedges? Sumatroscirpeae, a missing piece in the evolutionary puzzle of the giant genus Carex (Cyperaceae)

For over a century, the origins and mechanisms underlying the diversification of the enormous temperate genus Carex (>2100 species; Cariceae, Cyperaceae) have remained largely speculative. Characteristics such as its diverse ecology, varied biogeography, and intriguing cytology have made Carex a powerful model for studying plant evolution, but its uncertain sister-group relationships hinder its use in studies that depend on accurate ancestral state estimates and biogeographic inferences. To identify the sister to Carex, we estimated the phylogeny of all genera in the Cariceae-Dulichieae-Scirpeae clade (CDS) using three plastid and two nuclear ribosomal markers. Ancestral state reconstructions of key characters were made, and a time-calibrated tree estimated. Carex is strongly supported as sister to the rare East Asian Sumatroscirpus, sole genus of a new tribe, Sumatroscirpeae trib. nov. Believed to be unique to Carex, the perigynium (prophyllar bract enclosing a flower) is in fact a synapomorphy shared with this small tribe (∼4 species) that appeared 36 Mya. We thus suggest the initial key innovation in the remarkable diversification of Carex is not the perigynium, but could be the release of mechanical constraints on perigynia through spikelet truncation, resulting in novel adaptive morphologies. Monoecy, chromosomal change, and rapid inflorescence development enabling phenological isolation may also be involved. The tiny tribe Sumatroscirpeae will provide unprecedented insights into the inflorescence homology, evolution, diversification, and biogeographic history of its sister-group Carex, one of the world's most diverse plant lineages.

Can quartet analyses combining maximum likelihood estimation and Hennigian logic overcome long branch attraction in phylogenomic sequence data?

Systematic biases such as long branch attraction can mislead commonly relied upon model-based (i.e. maximum likelihood and Bayesian) phylogenetic methods when, as is usually the case with empirical data, there is model misspecification. We present PhyQuart, a new method for evaluating the three possible binary trees for any quartet of taxa. PhyQuart was developed through a process of reciprocal illumination between a priori considerations and the results of extensive simulations. It is based on identification of site-patterns that can be considered to support a particular quartet tree taking into account the Hennigian distinction between apomorphic and plesiomorphic similarity, and employing corrections to the raw observed frequencies of site-patterns that exploit expectations from maximum likelihood estimation. We demonstrate through extensive simulation experiments that, whereas maximum likeilihood estimation performs well in many cases, it can be outperformed by PhyQuart in cases where it fails due to extreme branch length asymmetries producing long-branch attraction artefacts where there is only very minor model misspecification.

Linear Plasmids and the Rate of Sequence Evolution in Plant Mitochondrial Genomes

The mitochondrial genomes of flowering plants experience frequent insertions of foreign sequences, including linear plasmids that also exist in standalone forms within mitochondria, but the history and phylogenetic distribution of plasmid insertions is not well known. Taking advantage of the increased availability of plant mitochondrial genome sequences, we performed phylogenetic analyses to reconstruct the evolutionary history of these plasmids and plasmid-derived insertions. Mitochondrial genomes from multiple land plant lineages (including liverworts, lycophytes, ferns, and gymnosperms) include fragmented remnants from ancient plasmid insertions. Such insertions are much more recent and widespread in angiosperms, in which approximately 75% of sequenced mitochondrial genomes contain identifiable plasmid insertions. Although conflicts between plasmid and angiosperm phylogenies provide clear evidence of repeated horizontal transfers, we were still able to detect significant phylogenetic concordance, indicating that mitochondrial plasmids have also experienced sustained periods of (effectively) vertical transmission in angiosperms. The observed levels of sequence divergence in plasmid-derived genes suggest that nucleotide substitution rates in these plasmids, which often encode their own viral-like DNA polymerases, are orders of magnitude higher than in mitochondrial chromosomes. Based on these results, we hypothesize that the periodic incorporation of mitochondrial genes into plasmids contributes to the remarkable heterogeneity in substitution rates among genes that has recently been discovered in some angiosperm mitochondrial genomes. In support of this hypothesis, we show that the recently acquired ψtrnP-trnW gene region in a maize linear plasmid is evolving significantly faster than homologous sequences that have been retained in the mitochondrial chromosome in closely related grasses.

Biogeography of the Malagasy Celastraceae: Multiple independent origins followed by widespread dispersal of genera from Madagascar

Of the 97 currently recognized genera of Celastraceae, 19 are native to Madagascar, including six endemics. In this study we conducted the most thorough phylogenetic analysis of Celastraceae yet completed with respect to both character and taxon sampling, and include representatives of five new endemic genera. Fifty-one new accessions, together with 328 previously used accessions of Celastrales, were sampled for morphological characters, two rDNA gene regions, and two plastid gene regions. The endemic Malagasy genera are resolved in two separate lineages-Xenodrys by itself and all other endemic genera in a clade that also includes four lineages inferred to have dispersed from Madagascar: Brexia madagascariensis (Mascarene Islands, coastal Africa), Elaeodendron (West Indies, Africa to New Caledonia), and Pleurostylia (Africa to New Caledonia). Of the 12 extant Malagasy Celastraceae lineages identified, eight are clearly of African origin. The origins of the remaining four lineages are less clear, but reasonable possibilities include America, Eurasia, Africa, southern India, Malesia, and Australia. Based on 95% credible age intervals from fossil-calibrated molecular dating, all 12 extant Malagasy Celastraceae lineages appear to have arisen following dispersal after the separation of Madagascar from other landmasses within the last 70 million years.

Testing reticulate versus coalescent origins of Erica lusitanica using a species phylogeny of the northern heathers (Ericeae, Ericaceae)

Whilst most of the immense species richness of heathers (Calluna, Daboecia and Erica: Ericeae; Ericaceae) is endemic to Africa, particularly the Cape Floristic Region, the oldest lineages are found in the Northern Hemisphere. We present phylogenetic hypotheses for the major clades of Ericeae represented by multiple accessions of all northern Erica species and placeholder taxa for the large nested African/Madagascan clade. We identified consistent, strongly supported conflict between gene trees inferred from ITS and chloroplast DNA sequences with regard to the position of Erica lusitanica. We used coalescent simulations to test whether this conflict could be explained by coalescent stochasticity, as opposed to reticulation (e.g. hybridisation), given estimates of clade ages, generation time and effective population sizes (Ne). A standard approach, comparing overall differences between real and simulated trees, could not clearly reject coalescence. However, additional simulations showed that at the (higher) Ne necessary to explain conflict in E. lusitanica, further topological conflict would also be expected. Ancient hybridisation between ancestors of northern species is therefore a plausible scenario to explain the origin of E. lusitanica, and its morphological similarities to E. arborea. Assuming either process influences the results of species tree and further evolutionary inference. The coalescence scenario is equivocal with regard the standing hypothesis of stepping stone dispersal of Erica from Europe into Africa; whereas reticulate evolution in E. lusitanica would imply that the colonisation of Tropical East Africa by E. arborea instead occurred independently of dispersals within the rest of the African/Madagascan clade.

Phylogenetic relationships in Epidendroideae (Orchidaceae), one of the great flowering plant radiations: progressive specialization and diversification

BACKGROUND AND AIMS

The largest subfamily of orchids, Epidendroideae, represents one of the most significant diversifications among flowering plants in terms of pollination strategy, vegetative adaptation and number of species. Although many groups in the subfamily have been resolved, significant relationships in the tree remain unclear, limiting conclusions about diversification and creating uncertainty in the classification. This study brings together DNA sequences from nuclear, plastid and mitochrondrial genomes in order to clarify relationships, to test associations of key characters with diversification and to improve the classification.

METHODS

Sequences from seven loci were concatenated in a supermatrix analysis for 312 genera representing most of epidendroid diversity. Maximum-likelihood and parsimony analyses were performed on this matrix and on subsets of the data to generate trees and to investigate the effect of missing values. Statistical character-associated diversification analyses were performed.

KEY RESULTS

Likelihood and parsimony analyses yielded highly resolved trees that are in strong agreement and show significant support for many key clades. Many previously proposed relationships among tribes and subtribes are supported, and some new relationships are revealed. Analyses of subsets of the data suggest that the relatively high number of missing data for the full analysis is not problematic. Diversification analyses show that epiphytism is most strongly associated with diversification among epidendroids, followed by expansion into the New World and anther characters that are involved with pollinator specificity, namely early anther inflexion, cellular pollinium stalks and the superposed pollinium arrangement.

CONCLUSIONS

All tested characters show significant association with speciation in Epidendroideae, suggesting that no single character accounts for the success of this group. Rather, it appears that a succession of key features appeared that have contributed to diversification, sometimes in parallel.

Platyzoan paraphyly based on phylogenomic data supports a noncoelomate ancestry of spiralia

Based on molecular data three major clades have been recognized within Bilateria: Deuterostomia, Ecdysozoa, and Spiralia. Within Spiralia, small-sized and simply organized animals such as flatworms, gastrotrichs, and gnathostomulids have recently been grouped together as Platyzoa. However, the representation of putative platyzoans was low in the respective molecular phylogenetic studies, in terms of both, taxon number and sequence data. Furthermore, increased substitution rates in platyzoan taxa raised the possibility that monophyletic Platyzoa represents an artifact due to long-branch attraction. In order to overcome such problems, we employed a phylogenomic approach, thereby substantially increasing 1) the number of sampled species within Platyzoa and 2) species-specific sequence coverage in data sets of up to 82,162 amino acid positions. Using established and new measures (long-branch score), we disentangled phylogenetic signal from misleading effects such as long-branch attraction. In doing so, our phylogenomic analyses did not recover a monophyletic origin of platyzoan taxa that, instead, appeared paraphyletic with respect to the other spiralians. Platyhelminthes and Gastrotricha formed a monophylum, which we name Rouphozoa. To the exclusion of Gnathifera, Rouphozoa and all other spiralians represent a monophyletic group, which we name Platytrochozoa. Platyzoan paraphyly suggests that the last common ancestor of Spiralia was a simple-bodied organism lacking coelomic cavities, segmentation, and complex brain structures, and that more complex animals such as annelids evolved from such a simply organized ancestor. This conclusion contradicts alternative evolutionary scenarios proposing an annelid-like ancestor of Bilateria and Spiralia and several independent events of secondary reduction.

Resolving ancient radiations: can complete plastid gene sets elucidate deep relationships among the tropical gingers (Zingiberales)?

BACKGROUND AND AIMS

Zingiberales comprise a clade of eight tropical monocot families including approx. 2500 species and are hypothesized to have undergone an ancient, rapid radiation during the Cretaceous. Zingiberales display substantial variation in floral morphology, and several members are ecologically and economically important. Deep phylogenetic relationships among primary lineages of Zingiberales have proved difficult to resolve in previous studies, representing a key region of uncertainty in the monocot tree of life.

METHODS

Next-generation sequencing was used to construct complete plastid gene sets for nine taxa of Zingiberales, which were added to five previously sequenced sets in an attempt to resolve deep relationships among families in the order. Variation in taxon sampling, process partition inclusion and partition model parameters were examined to assess their effects on topology and support.

KEY RESULTS

Codon-based likelihood analysis identified a strongly supported clade of ((Cannaceae, Marantaceae), (Costaceae, Zingiberaceae)), sister to (Musaceae, (Lowiaceae, Strelitziaceae)), collectively sister to Heliconiaceae. However, the deepest divergences in this phylogenetic analysis comprised short branches with weak support. Additionally, manipulation of matrices resulted in differing deep topologies in an unpredictable fashion. Alternative topology testing allowed statistical rejection of some of the topologies. Saturation fails to explain observed topological uncertainty and low support at the base of Zingiberales. Evidence for conflict among the plastid data was based on a support metric that accounts for conflicting resampled topologies.

CONCLUSIONS

Many relationships were resolved with robust support, but the paucity of character information supporting the deepest nodes and the existence of conflict suggest that plastid coding regions are insufficient to resolve and support the earliest divergences among families of Zingiberales. Whole plastomes will continue to be highly useful in plant phylogenetics, but the current study adds to a growing body of literature suggesting that they may not provide enough character information for resolving ancient, rapid radiations.

Decisive data sets in phylogenomics: lessons from studies on the phylogenetic relationships of primarily wingless insects

Phylogenetic relationships of the primarily wingless insects are still considered unresolved. Even the most comprehensive phylogenomic studies that addressed this question did not yield congruent results. To get a grip on these problems, we here analyzed the sources of incongruence in these phylogenomic studies by using an extended transcriptome data set. Our analyses showed that unevenly distributed missing data can be severely misleading by inflating node support despite the absence of phylogenetic signal. In consequence, only decisive data sets should be used which exclusively comprise data blocks containing all taxa whose relationships are addressed. Additionally, we used Four-cluster Likelihood Mapping (FcLM) to measure the degree of congruence among genes of a data set, as a measure of support alternative to bootstrap. FcLM showed incongruent signal among genes, which in our case is correlated neither with functional class assignment of these genes nor with model misspecification due to unpartitioned analyses. The herein analyzed data set is the currently largest data set covering primarily wingless insects, but failed to elucidate their interordinal phylogenetic relationships. Although this is unsatisfying from a phylogenetic perspective, we try to show that the analyses of structure and signal within phylogenomic data can protect us from biased phylogenetic inferences due to analytical artifacts.

Phylogeny of the Asparagales based on three plastid and two mitochondrial genes

PREMISE OF THE STUDY

The Asparagales, with ca. 40% of all monocotyledons, include a host of commercially important ornamentals in families such as Orchidaceae, Alliaceae, and Iridaceae, and several important crop species in genera such as Allium, Aloe, Asparagus, Crocus, and Vanilla. Though the order is well defined, the number of recognized families, their circumscription, and relationships are somewhat controversial.

METHODS

Phylogenetic analyses of Asparagales were based on parsimony and maximum likelihood using nucleotide sequence variation in three plastid genes (matK, ndhF, and rbcL) and two mitochondrial genes (atp1 and cob). Branch support was assessed using both jackknife analysis implementing strict-consensus (SC) and bootstrap analysis implementing frequency-within-replicates (FWR). The contribution of edited sites in the mitochondrial genes to topology and branch support was investigated.

KEY RESULTS

The topologies recovered largely agree with previous results, though some clades remain poorly resolved (e.g., Ruscaceae). When the edited sites were included in the analysis, the plastid and mitochondrial genes were highly incongruent. However, when the edited sites were removed, the two partitions became congruent.

CONCLUSIONS

Some deeper nodes in the Asparagales tree remain poorly resolved or unresolved as do the relationships of certain monogeneric families (e.g., Aphyllanthaceae, Ixioliriaceae, Doryanthaceae), whereas support for many families increases. However, the increased support is dominated by plastid data, and the potential influence of mitochondrial and biparentially inherited single or low-copy nuclear genes should be investigated.

Radical instability and spurious branch support by likelihood when applied to matrices with non-random distributions of missing data

Non-random distributions of missing data are a general problem for likelihood-based statistical analyses, including those in a phylogenetic context. Extensive non-randomly distributed missing data are particularly problematic in supermatrix analyses that include many terminals and/or loci. It has been widely reported that missing data can lead to loss of resolution, but only very rarely create misleading or otherwise unsupported results in a parsimony context. Yet this does not hold for all parametric-based analyses because of their assumption of homogeneity across characters and lineages, which can lead to both long-branch attraction and long-branch repulsion. Contrived examples were used to demonstrate that non-random distributions of missing data, even without rate heterogeneity among characters and a well fitting model, can provide misleading likelihood-based topologies and branch-support values that are radically unstable based on slight modifications to character sampling. The same can occur despite complete absence of parsimony-informative characters. Otherwise unsupported resolution and high branch support for these clades were found to occur frequently in 22 empirical examples derived from a published supermatrix. Partitioning characters based on the distribution of missing data helped to decrease, but did not eliminate, these artifacts. These artifacts were exacerbated by low quality tree searches, particularly when holding only a single optimal tree that must be fully resolved.

Phylogeny of Celastraceae tribe Euonymeae inferred from morphological characters and nuclear and plastid genes

The phylogeny of Celastraceae tribe Euonymeae (≈ 230 species in eight genera in both the Old and New Worlds) was inferred using morphological characters together with plastid (matK, trnL-F) and nuclear (ITS and 26S rDNA) genes. Tribe Euonymeae has been defined as those genera of Celastraceae with generally opposite leaves, isomerous carpels, loculicidally dehiscent capsules, and arillate seeds (except Microtropis). Euonymus is the most diverse (129 species) and widely cultivated genus in the tribe. We infer that tribe Euonymeae consists of at least six separate lineages within Celastraceae and that a revised natural classification of the family is needed. Microtropis and Quetzalia are inferred to be distinct sister groups that together are sister to Zinowiewia. The endangered Monimopetalum chinense is an isolated and early derived lineage of Celastraceae that represents an important component of phylogenetic diversity within the family. Hedraianthera is sister to Brassiantha, and we describe a second species (Brassiantha hedraiantheroides A.J. Ford) that represents the first reported occurrence of this genus in Australia. Euonymus globularis, from eastern Australia, is sister to Menepetalum, which is endemic to New Caledonia, and we erect a new genus (Dinghoua R.H. Archer) for it. The Madagascan species of Euonymus are sister to Pleurostylia and recognized as a distinct genus (Astrocassine ined.). Glyptopetalum, Torralbasia, and Xylonymus are all closely related to Euonymus sensu stricto and are questionably distinct from it. Current intrageneric classifications of Euonymus are not completely natural and require revision.