Quantification of congruence among gene trees with polytomies using overall success of resolution for phylogenomic coalescent analyses.

Gene-tree-inference error can cause species-tree-inference artefacts in summary phylogenomic coalescent analyses. Here we integrate two ways of accommodating these inference errors: collapsing arbitrarily or dubiously resolved gene-tree branches, and subsampling gene trees based on their pairwise congruence. We tested the effect of collapsing gene-tree branches with 0% approximate-likelihood-ratio-test (SH-like aLRT) support in likelihood analyses and strict consensus trees for parsimony, and then subsampled those partially resolved trees based on congruence measures that do not penalize polytomies. For this purpose we developed a new TNT script for congruence sorting (congsort), and used it to calculate topological incongruence for eight phylogenomic datasets using three distance measures: standard Robinson-Foulds (RF) distances; overall success of resolution (OSR), which is based on counting both matching and contradicting clades; and RF contradictions, which only counts contradictory clades. As expected, we found that gene-tree incongruence was often concentrated in clades that are arbitrarily or dubiously resolved and that there was greater congruence between the partially collapsed gene trees and the coalescent and concatenation topologies inferred from those genes. Coalescent branch lengths typically increased as the most incongruent gene trees were excluded, although branch supports typically did not. We investigated two successful and complementary approaches to prioritizing genes for investigation of alignment or homology errors. Coalescent-tree clades that contradicted concatenation-tree clades were generally less robust to gene-tree subsampling than congruent clades. Our preferred approach to collapsing likelihood gene-tree clades (0% SH-like aLRT support) and subsampling those trees (OSR) generally outperformed competing approaches for a large fungal dataset with respect to branch lengths, support and congruence. We recommend widespread application of this approach (and strict consensus trees for parsimony-based analyses) for improving quantification of gene-tree congruence/conflict, estimating coalescent branch lengths, testing robustness of coalescent analyses to gene-tree-estimation error, and improving topological robustness of summary coalescent analyses. This approach is quick and easy to implement, even for huge datasets.

Benefits of alignment quality-control processing steps and an Angiosperms353 phylogenomics pipeline applied to the Celastrales.

We examined the impact of successive alignment quality-control steps on downstream phylogenomic analyses. We applied a recently published phylogenomics pipeline that was developed for the Angiosperms353 target-sequence-capture probe set to the flowering plant order Celastrales. Our final dataset consists of 158 species, including at least one exemplar from all 109 currently recognized Celastrales genera. We performed nine quality-control steps and compared the inferred resolution, branch support, and topological congruence of the inferred gene and species trees with those generated after each of the first six steps. We describe and justify each of our quality-control steps, including manual masking, in detail so that they may be readily applied to other lineages. We found that highly supported clades could generally be relied upon even if stringent orthology and alignment quality-control measures had not been applied. But separate instances were identified, for both concatenation and coalescence, wherein a clade was highly supported before manual masking but then subsequently contradicted. These results are generally reassuring for broad-scale analyses that use phylogenomics pipelines, but also indicate that we cannot rely exclusively on these analyses to conclude how challenging phylogenetic problems are best resolved.

Gene-tree misrooting drives conflicts in phylogenomic coalescent analyses of palaeognath birds.

Phylogenomic analyses of ancient rapid radiations can produce conflicting results that are driven by differential sampling of taxa and characters as well as the limitations of alternative analytical methods. We re-examine basal relationships of palaeognath birds (ratites and tinamous) using recently published datasets of nucleotide characters from 20,850 loci as well as 4301 retroelement insertions. The original studies attributed conflicting resolutions of rheas in their inferred coalescent and concatenation trees to concatenation failing in the anomaly zone. By contrast, we find that the coalescent-based resolution of rheas is premised upon extensive gene-tree estimation errors. Furthermore, retroelement insertions contain much more conflict than originally reported and multiple insertion loci support the basal position of rheas found in concatenation trees, while none were reported in the original publication. We demonstrate how even remarkable congruence in phylogenomic studies may be driven by long-branch misplacement of a divergent outgroup, highly incongruent gene trees, differential taxon sampling that can result in gene-tree misrooting errors that bias species-tree inference, and gross homology errors. What was previously interpreted as broad, robustly supported corroboration for a single resolution in coalescent analyses may instead indicate a common bias that taints phylogenomic results across multiple genome-scale datasets. The updated retroelement dataset now supports a species tree with branch lengths that suggest an ancient anomaly zone, and both concatenation and coalescent analyses of the huge nucleotide datasets fail to yield coherent, reliable results in this challenging phylogenetic context.

Collapsing dubiously resolved gene-tree branches in phylogenomic coalescent analyses.

In two-step coalescent analyses of phylogenomic data, gene-tree topologies are treated as fixed prior to species-tree inference. Although all gene-tree conflict is assumed to be caused by lineage sorting when applying these methods, in empirical datasets much of the conflict can be caused by estimation error. Weakly supported and even arbitrarily resolved clades are important sources of this estimation error for gene trees inferred from few informative characters relative to the number of sampled terminals, and the resulting extraneous conflict among gene trees can negatively impact species-tree inference. In this study, we quantified the relative severity of alternative methods for collapsing gene-tree branches for seven empirical datasets and quantified their effects on species-tree inference. The branch-collapsing methods that we employed were based on the strict consensus of optimal topologies, various bootstrap thresholds, and 0% approximate likelihood ratio test (SH-like aLRT) support. Up to 86% of internal gene-tree branches are dubiously or arbitrarily resolved in reanalyses of these published phylogenomic datasets, and collapsing these branches increased inferred species-tree coalescent branch lengths by up to 455%. For two datasets, the longer inferred branch lengths sometimes impacted inference of anomaly-zone conditions. Although branch-collapsing methods did not consistently affect the species-tree topology, they often increased branch support. The more severe and clearly justified gene-tree branch-collapsing methods, which we recommend be broadly applied for two-step coalescent analyses, are use of the strict consensus in parsimony analyses and the collapse clades with 0% SH-like aLRT support in likelihood analyses. Collapsing dubiously or arbitrarily resolved branches in gene trees sometimes improved congruence between coalescent-based results and concatenation trees. In such cases, we contend that the resolution provided by concatenation should be preferred and that incomplete lineage sorting is a poor explanation for the initial conflict between phylogenetic approaches.

Treatment for Severe Coronavirus Disease 2019 With the Seraph-100 Microbind Affinity Blood Filter.

To determine whether Seraph-100 (Exthera Medical Corporation, Martinez, CA) treatment provides clinical benefit for severe coronavirus disease 2019 cases that require mechanical ventilation and vasopressor support. DATA SOURCES: The first two patients in the United States treated with the novel Seraph-100 device. These cases were reviewed by the Food and Drug Administration prior to granting an emergency use authorization for treatment of coronavirus disease 2019. STUDY SELECTION: Case series. DATA EXTRACTION: Vasopressor dose, mean arterial pressure, temperature, interleukin-6, C-reactive protein, and other biomarker levels were documented both before and after Seraph-100 treatments. DATA SYNTHESIS: Vasopressor dose, temperature, interleukin-6, and C-reactive protein levels declined after Seraph-100 treatments. Severe acute respiratory syndrome coronavirus 2 viremia was confirmed in the one patient tested and cleared by the completion of treatments. CONCLUSIONS: Seraph-100 use may improve hemodynamic stability in coronavirus disease 2019 cases requiring mechanical ventilation and vasopressor support. These findings warrant future study of a larger cohort with the addition of mortality and total hospital day outcomes.

Mitochondrial fatty acid ß-oxidation is required for storage-lipid catabolism in a marine diatom.

Photoautotrophic growth in nature requires the accumulation of energy-containing molecules via photosynthesis during daylight to fuel nighttime catabolism. Many diatoms store photosynthate as the neutral lipid triacylglycerol (TAG). While the pathways of diatom fatty acid and TAG synthesis appear to be well conserved with plants, the pathways of TAG catabolism and downstream fatty acid ß-oxidation have not been characterised in diatoms. We identified a putative mitochondria-targeted, bacterial-type acyl-CoA dehydrogenase (PtMACAD1) that is present in Stramenopile and Hacrobian eukaryotes, but not found in plants, animals or fungi. Gene knockout, protein-YFP tags and physiological assays were used to determine PtMACAD1's role in the diatom Phaeodactylum tricornutum. PtMACAD1 is located in the mitochondria. Absence of PtMACAD1 led to no consumption of TAG at night and slower growth in light : dark cycles compared with wild-type. Accumulation of transcripts encoding peroxisomal-based ß-oxidation did not change in response to day : night cycles or to PtMACAD1 knockout. Mutants also hyperaccumulated TAG after the amelioration of N limitation. We conclude that diatoms utilise mitochondrial ß-oxidation; this is in stark contrast to the peroxisomal-based pathways observed in plants and green algae. We infer that this pattern is caused by retention of catabolic pathways from the host during plastid secondary endosymbiosis.

Divergence and support among slightly suboptimal likelihood gene trees.

Contemporary phylogenomic studies frequently incorporate two-step coalescent analyses wherein the first step is to infer individual-gene trees, generally using maximum-likelihood implemented in the popular programs PhyML or RAxML. Four concerns with this approach are that these programs only present a single fully resolved gene tree to the user despite potential for ambiguous support, insufficient phylogenetic signal to fully resolve each gene tree, inexact computer arithmetic affecting the reported likelihood of gene trees, and an exclusive focus on the most likely tree while ignoring trees that are only slightly suboptimal or within the error tolerance. Taken together, these four concerns are sufficient for RAxML and PhyML users to be suspicious of the resulting (perhaps over-resolved) gene-tree topologies and (perhaps unjustifiably high) bootstrap support for individual clades. In this study, we sought to determine how frequently these concerns apply in practice to contemporary phylogenomic studies that use RAxML for gene-tree inference. We did so by re-analyzing 100 genes from each of ten studies that, taken together, are representative of many empirical phylogenomic studies. Our seven findings are as follows. First, the few search replicates that are frequently applied in phylogenomic studies are generally insufficient to find the optimal gene-tree topology. Second, there is often more topological variation among slightly suboptimal gene trees relative to the best-reported tree than can be safely ignored. Third, the Shimodaira-Hasegawa-like approximate likelihood ratio test is highly effective at identifying dubiously supported clades and outperforms the alternative approaches of relying on bootstrap support or collapsing minimum-length branches. Fourth, the bootstrap can, but rarely does, indicate high support for clades that are not supported amongst slightly suboptimal trees. Fifth, increasing the accuracy by which RAxML optimizes model-parameter values generally has a nominal effect on selection of optimal trees. Sixth, tree searches using the GTRCAT model were generally less effective at finding optimal known trees than those using the GTRGAMMA model. Seventh, choice of gene-tree sampling strategy can affect inferred coalescent branch lengths, species-tree topology and branch support.

ILS-Aware Analysis of Low-Homoplasy Retroelement Insertions: Inference of Species Trees and Introgression Using Quartets.

DNA sequence alignments have provided the majority of data for inferring phylogenetic relationships with both concatenation and coalescent methods. However, DNA sequences are susceptible to extensive homoplasy, especially for deep divergences in the Tree of Life. Retroelement insertions have emerged as a powerful alternative to sequences for deciphering evolutionary relationships because these data are nearly homoplasy-free. In addition, retroelement insertions satisfy the "no intralocus-recombination" assumption of summary coalescent methods because they are singular events and better approximate neutrality relative to DNA loci commonly sampled in phylogenomic studies. Retroelements have traditionally been analyzed with parsimony, distance, and network methods. Here, we analyze retroelement data sets for vertebrate clades (Placentalia, Laurasiatheria, Balaenopteroidea, Palaeognathae) with 2 ILS-aware methods that operate by extracting, weighting, and then assembling unrooted quartets into a species tree. The first approach constructs a species tree from retroelement bipartitions with ASTRAL, and the second method is based on split-decomposition with parsimony. We also develop a Quartet-Asymmetry test to detect hybridization using retroelements. Both ILS-aware methods recovered the same species-tree topology for each data set. The ASTRAL species trees for Laurasiatheria have consecutive short branch lengths in the anomaly zone whereas Palaeognathae is outside of this zone. For the Balaenopteroidea data set, which includes rorquals (Balaenopteridae) and gray whale (Eschrichtiidae), both ILS-aware methods resolved balaeonopterids as paraphyletic. Application of the Quartet-Asymmetry test to this data set detected 19 different quartets of species for which historical introgression may be inferred. Evidence for introgression was not detected in the other data sets.

Adenine·cytosine substitutions are an alternative pathway of compensatory mutation in angiosperm ITS2.

Compensatory mutations are crucial for functional RNA because they maintain RNA configuration and thus function. Compensatory mutation has traditionally been considered to be a two-step substitution through the GU-base-pair intermediate. We tested for an alternative AC-mediated compensatory mutation (ACCM). We investigated ACCMs by using a comprehensive sampling of ribosomal internal transcribed spacer 2 (ITS2) from 3934 angiosperm species in 80 genera and 55 families. We predicted ITS2 consensus secondary structures by using LocARNA for structure-based alignment and partitioning paired and unpaired regions. We examined and compared the substitution rates and frequencies among base pairs by using RNA-specific models. Base-pair states of ACCMs were mapped onto the inferred phylogenetic trees to infer their evolution. All types of compensatory mutations involving the AC intermediate were observed, but the most frequent substitutions were with AU or GC pairs, which are part of the AU-AC-GC pathway. Compared with the GU intermediate, AC had a lower frequency and higher mutability. Within the AU-AC-GC pathway, the AU-AC substitution rate was much slower than the AC-GC substitution rate. No consistently higher overall rate was identified for either pathway among all 80 sampled lineages, though compensatory mutations through the AC intermediate averaged about half that through the GU intermediate. These results demonstrate an alternative compensatory mutation between AU and GC that helps address the controversial inference of inferred simultaneous double substitutions.

Alternative analyses of compensatory base changes in an ITS2 phylogeny of Corydalis (Papaveraceae).

BACKGROUND AND AIMS: Compensatory base changes (CBCs) that occur in stems of ribosomal internal transcribed spacer 2 (ITS2) can have important phylogenetic implications because they are not expected to occur within a single species and also affect selection of appropriate DNA substitution models. These effects have been demonstrated when studying ancient lineages. Here we examine these effects to quantify their importance within a more recent lineage by using both DNA- and RNA-specific models. METHODS: We examined the phylogenetic implications of the CBC process by using a comprehensive sampling of ITS2 from ten closely related species of Corydalis. We predicted ITS2 secondary structures by using homology modelling, which was then used for a structure-based alignment. Paired and unpaired regions were analysed separately and in combination by using both RNA-specific substitution models and conventional DNA models. We mapped all base-pair states of CBCs on the phylogenetic tree to infer their evolution and relative timing. KEY RESULTS: Our results indicate that selection acted to increase the thermodynamic stability of the secondary structure. Thus, the unpaired and paired regions did not evolve under a common substitution model. Only two CBCs occurred within the lineage sampled and no striking differences in topology or support for the shared clades were found between trees constructed using DNA- or RNA-specific substitution models. CONCLUSIONS: Although application of RNA-specific substitution models remains preferred over more conventional DNA models, we infer that application of conventional DNA models is unlikely to be problematic when conducting phylogenetic analyses of ITS2 within closely related lineages wherein few CBCs are observed. Each of the two CBCs was found within the same lineages but was not observed within a given species, which supports application of the CBC species concept.

Partitioned coalescence support reveals biases in species-tree methods and detects gene trees that determine phylogenomic conflicts.

Genomic datasets sometimes support conflicting phylogenetic relationships when different tree-building methods are applied. Coherent interpretations of such results are enabled by partitioning support for controversial relationships among the constituent genes of a phylogenomic dataset. For the supermatrix (=concatenation) approach, several methods that measure the distribution of support and conflict among loci were introduced over 15 years ago. More recently, partitioned coalescence support (PCS) was developed for phylogenetic coalescence methods that account for incomplete lineage sorting and use the summed fits of gene trees to estimate the species tree. Here, we automate computation of PCS to permit application of this index to genome-scale matrices that include hundreds of loci. Reanalyses of four phylogenomic datasets for amniotes, land plants, skinks, and angiosperms demonstrate how PCS scores can be used to: (1) compare conflicting results favored by alternative coalescence methods, (2) identify outlier gene trees that have a disproportionate influence on the resolution of contentious relationships, (3) assess the effects of missing data in species-tree analysis, and (4) clarify biases in commonly-implemented coalescence methods and support indices. We show that key phylogenomic conclusions from these analyses often hinge on just a few gene trees and that results can be driven by specific biases of a particular coalescence method and/or the differential weight placed on gene trees with high versus low taxon sampling. The attribution of exceptionally high weight to some gene trees and very low weight to other gene trees counters the basic logic of phylogenomic coalescence analysis; even clades in species trees with high support according to commonly used indices (likelihood-ratio test, bootstrap, Bayesian local posterior probability) can be unstable to the removal of only one or two gene trees with high PCS. Computer simulations cannot adequately describe all of the contingencies and complexities of empirical genetic data. PCS scores complement simulation work by providing specific insights into a particular dataset given the assumptions of the phylogenetic coalescence method that is applied. In combination with standard measures of nodal support, PCS provides a more complete understanding of the overall genomic evidence for contested evolutionary relationships in species trees.

Gene-wise resampling outperforms site-wise resampling in phylogenetic coalescence analyses.

In summary ("two-step") coalescent analyses of empirical data, researchers typically apply the bootstrap to quantify branch support for clades inferred on the optimal species tree. We tested whether site-wise bootstrap analyses provide consistently more conservative support than gene-wise bootstrap analyses. We did so using data from three empirical phylogenomic studies and employed four coalescent methods (ASTRAL, MP-EST, NJst, and STAR). We demonstrate that application of site-wise bootstrapping generally resulted in gene-trees with substantial additional conflicts relative to the original data and this approach therefore cannot be relied upon to provide conservative support. Instead the site-wise bootstrap can provide high support for apparently incorrect clades. We provide a script (https://github.com/dbsloan/msc_tree_resampling) that implements gene-wise resampling, using either the bootstrap or the jackknife, for use with ASTRAL, MP-EST, NJst, and STAR. We demonstrate that the gene-wise bootstrap outperformed the site-wise bootstrap for the primary focal clades for all four coalescent methods that were applied to all three empirical studies. For summary coalescent analyses we suggest that gene-wise resampling support should be favored over gene + site or site-wise resampling when numerous genes are sampled because site-wise resampling causes substantially greater gene-tree-estimation error.

Narrow thermal tolerance and low dispersal drive higher speciation in tropical mountains.

Species richness is greatest in the tropics, and much of this diversity is concentrated in mountains. Janzen proposed that reduced seasonal temperature variation selects for narrower thermal tolerances and limited dispersal along tropical elevation gradients [Janzen DH (1967) Am Nat 101:233-249]. These locally adapted traits should, in turn, promote reproductive isolation and higher speciation rates in tropical mountains compared with temperate ones. Here, we show that tropical and temperate montane stream insects have diverged in thermal tolerance and dispersal capacity, two key traits that are drivers of isolation in montane populations. Tropical species in each of three insect clades have markedly narrower thermal tolerances and lower dispersal than temperate species, resulting in significantly greater population divergence, higher cryptic diversity, higher tropical speciation rates, and greater accumulation of species over time. Our study also indicates that tropical montane species, with narrower thermal tolerance and reduced dispersal ability, will be especially vulnerable to rapid climate change.

Transcriptome-wide comparison of selenium hyperaccumulator and nonaccumulator Stanleya species provides new insight into key processes mediating the hyperaccumulation syndrome.

To obtain better insight into the mechanisms of selenium hyperaccumulation in Stanleya pinnata, transcriptome-wide differences in root and shoot gene expression levels were investigated in S. pinnata and related nonaccumulator Stanleya elata grown with or without 20 µm selenate. Genes predicted to be involved in sulphate/selenate transport and assimilation or in oxidative stress resistance (glutathione-related genes and peroxidases) were among the most differentially expressed between species; many showed constitutively elevated expression in S. pinnata. A number of defence-related genes predicted to mediate synthesis and signalling of defence hormones jasmonic acid (JA, reported to induce sulphur assimilatory and glutathione biosynthesis genes), salicylic acid (SA) and ethylene were also more expressed in S. pinnata than S. elata. Several upstream signalling genes that up-regulate defence hormone synthesis showed higher expression in S. pinnata than S. elata and might trigger these selenium-mediated defence responses. Thus, selenium hyperaccumulation and hypertolerance in S. pinnata may be mediated by constitutive, up-regulated JA, SA and ethylene-mediated defence systems, associated with elevated expression of genes involved in sulphate/selenate uptake and assimilation or in antioxidant activity. Genes pinpointed in this study may be targets of genetic engineering of plants that may be employed in biofortification or phytoremediation.

Phylogeography of the wild and cultivated stimulant plant qat (Catha edulis, Celastraceae) in areas of historical cultivation.

PREMISE OF THE STUDY: Qat (Catha edulis, Celastraceae) is a woody plant species cultivated for its stimulant alkaloids. Qat is important to the economy and culture in large regions of Ethiopia, Kenya, and Yemen. Despite the importance of this species, the wild origins and dispersal of cultivars have only been described in often contradictory historical documents. We examined the wild origins, human-mediated dispersal, and genetic divergence of cultivated qat compared to wild qat. METHODS: We sampled 17 SSR markers and 1561 wild and cultivated individuals across the historical areas of qat cultivation. KEY RESULTS: On the basis of genetic structure inferred using Bayesian and nonparametric methods, two centers of origin in Kenya and one in Ethiopia were found for cultivated qat. The centers of origin in Ethiopia and northeast of Mt. Kenya are the primary sources of cultivated qat genotypes. Qat cultivated in Yemen is derived from Ethiopian genotypes rather than Yemeni wild populations. Cultivated qat with a wild Kenyan origin has not spread to Ethiopia or Yemen, whereas a small minority of qat cultivated in Kenya originated in Ethiopia. Hybrid genotypes with both Ethiopian and Kenyan parentage are present in northern Kenya. CONCLUSIONS: Ethiopian cultivars have diverged from their wild relatives, whereas Kenyan qat has diverged less. This pattern of divergence could be caused by the extinction of the wild-source qat populations in Ethiopia due to deforestation, undersampling, and/or artificial selection for agronomically important traits.

Employing Two-stage Derivatisation and GC-MS to Assay for Cathine and Related Stimulant Alkaloids across the Celastraceae.

INTRODUCTION: Catha edulis (qat, khat, mirra) is a woody plant species that is grown and consumed in East Africa and Yemen for its stimulant alkaloids cathinone, cathine and norephedrine. Two Celastraceae species, in addition to qat, have been noted for their stimulant properties in ethnobotanical literature. Recent phylogenetic reconstructions place four genera in a clade sister to Catha edulis, and these genera are primary candidates to search for cathine and related alkaloids. OBJECTIVE: Determine if cathine or related alkaloids are present in species of Celastraceae other than Catha edulis. METHODS: Leaf samples from 43 Celastraceae species were extracted in water followed by basification of the aqueous extract and partitioning with methyl-t-butyl ether to provide an alkaloid-enriched fraction. The extract was derivatised in a two-stage process and analysed using GC-MS for the presence of cathine. Related alkaloids and other metabolites in this alkaloid-enriched fraction were tentatively identified. RESULTS: Cathinone, cathine and norephedrine were not detected in any of the 43 Celastraceae species assayed other than Catha edulis. However, the phenylalanine- or tyrosine-derived alkaloid phenylethylamine was identified in five species. Nine species were found to be enriched for numerous sterol- and terpene-like compounds. CONCLUSION: These results indicate that cathine is unique to Catha edulis, and not the compound responsible for the stimulant properties reported in related Celastraceae species. Copyright © 2017 John Wiley & Sons, Ltd.

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

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

Mutually exclusive phylogenomic inferences at the root of the angiosperms: Amborella is supported as sister and Observed Variability is biased.

Identifying the extant sister group to the remaining angiosperms has been a subject of long debate, for which the primary currently competing hypotheses are that Amborella alone is sister or that the clade (Amborella, Nymphaeales) is sister. Both Xi et al. (Syst. Biol., 2014, 63, 919) and Goremykin et al. (Syst. Biol., 2015, 64, 879) identified Amborella as sister in concatenation-based phylogenetic analyses of their 310 nuclear genes and 78 plastid genes, respectively. But after application of Observed Variability-based character subsampling, both papers reported the clade (Amborella, Nymphaeales) as sister. Hence alternative character-sampling strategies may produce highly supported yet mutually exclusive phylogenetic inferences when applied to nuclear and plastid genomic data sets. Edwards et al. (Mol. Phylogenet. Evol., 2016, 94, 447) defended Observed Variability and the (Amborella, Nymphaeales) hypothesis. In this study I respond to Edwards et al.'s (Mol. Phylogenet. Evol., 2016, 94, 447) criticisms of Simmons and Gatesy (Mol. Phylogenet. Evol., 2015, 91, 98) and use Edwards et al.'s (Mol. Phylogenet. Evol., 2016, 94, 447) and Goremykin et al.'s (Syst. Biol., 2015, 64, 879) own data to demonstrate that the best-supported phylogenetic hypothesis is that Amborella alone is sister and that the competing evidence in favour of the (Amborella, Nymphaeales) hypothesis is caused primarily by methodological artifacts (biased character deletion by Observed Variability, MP-EST and STAR generally not being robust to the highly divergent and mis-rooted gene trees that were used).

Biases of tree-independent-character-subsampling methods.

Observed Variability (OV) and Tree Independent Generation of Evolutionary Rates (TIGER) are quick and easy-to-apply tree-independent methods that have been proposed to provide unbiased estimates of each character's rate of evolution and serve as the basis for excluding rapidly evolving characters. Both methods have been applied to multiple phylogenomic datasets, and in many cases the authors considered their trees inferred from the OV- and TIGER-delimited sub-matrices to be better estimates of the phylogeny than their trees based on all characters. In this study we use four sets of simulations and an empirical phylogenomic example to demonstrate that both methods share a systematic bias against characters with more symmetric distributions of character states, against characters with greater observed character-state space, and against large clades in the context of character conflict. As a result these methods can favor convergences and reversals over synapomorphy, exacerbate long-branch attraction, and produce mutually exclusive phylogenetic inferences that are dependent upon differential taxon sampling. We assert that neither OV nor TIGER should be relied upon to increase the ratio of phylogenetic to non-phylogenetic signal in a data matrix. We also assert that skepticism is warranted for empirical phylogenetic results that are based on OV- and/or TIGER-based character deletion wherein a small clade is supported after deletion of characters, yet is contradicted by a larger clade when the entire data matrix was analyzed.

The Complete Plastid Genome of Lagerstroemia fauriei and Loss of rpl2 Intron from Lagerstroemia (Lythraceae).

Lagerstroemia (crape myrtle) is an important plant genus used in ornamental horticulture in temperate regions worldwide. As such, numerous hybrids have been developed. However, DNA sequence resources and genome information for Lagerstroemia are limited, hindering evolutionary inferences regarding interspecific relationships. We report the complete plastid genome of Lagerstroemia fauriei. To our knowledge, this is the first reported whole plastid genome within Lythraceae. This genome is 152,440 bp in length with 38% GC content and consists of two single-copy regions separated by a pair of 25,793 bp inverted repeats. The large single copy and the small single copy regions span 83,921 bp and 16,933 bp, respectively. The genome contains 129 genes, including 17 located in each inverted repeat. Phylogenetic analysis of genera sampled from Geraniaceae, Myrtaceae, and Onagraceae corroborated the sister relationship between Lythraceae and Onagraceae. The plastid genomes of L. fauriei and several other Lythraceae species lack the rpl2 intron, which indicating an early loss of this intron within the Lythraceae lineage. The plastid genome of L. fauriei provides a much needed genetic resource for further phylogenetic research in Lagerstroemia and Lythraceae. Highly variable markers were identified for application in phylogenetic, barcoding and conservation genetic applications.