Multiple mechanisms explain loss of anthocyanins from betalain-pigmented Caryophyllales, including repeated wholesale loss of a key anthocyanidin synthesis enzyme.

In this study, we investigate the genetic mechanisms responsible for the loss of anthocyanins in betalain-pigmented Caryophyllales, considering our hypothesis of multiple transitions to betalain pigmentation. Utilizing transcriptomic and genomic datasets across 357 species and 31 families, we scrutinize 18 flavonoid pathway genes and six regulatory genes spanning four transitions to betalain pigmentation. We examined evidence for hypotheses of wholesale gene loss, modified gene function, altered gene expression, and degeneration of the MBW (MYB-bHLH-WD40) trasnscription factor complex, within betalain-pigmented lineages. Our analyses reveal that most flavonoid synthesis genes remain conserved in betalain-pigmented lineages, with the notable exception of TT19 orthologs, essential for the final step in anthocyanidin synthesis, which appear to have been repeatedly and entirely lost. Additional late-stage flavonoid pathway genes upstream of TT19 also manifest strikingly reduced expression in betalain-pigmented species. Additionally, we find repeated loss and alteration in the MBW transcription complex essential for canonical anthocyanin synthesis. Consequently, the loss and exclusion of anthocyanins in betalain-pigmented species appear to be orchestrated through several mechanisms: loss of a key enzyme, downregulation of synthesis genes, and degeneration of regulatory complexes. These changes have occurred iteratively in Caryophyllales, often coinciding with evolutionary transitions to betalain pigmentation.

MycoRed: Betalain pigments enable in vivo real-time visualisation of arbuscular mycorrhizal colonisation.

Arbuscular mycorrhiza (AM) are mutualistic interactions formed between soil fungi and plant roots. AM symbiosis is a fundamental and widespread trait in plants with the potential to sustainably enhance future crop yields. However, improving AM fungal association in crop species requires a fundamental understanding of host colonisation dynamics across varying agronomic and ecological contexts. To this end, we demonstrate the use of betalain pigments as in vivo visual markers for the occurrence and distribution of AM fungal colonisation by Rhizophagus irregularis in Medicago truncatula and Nicotiana benthamiana roots. Using established and novel AM-responsive promoters, we assembled multigene reporter constructs that enable the AM-controlled expression of the core betalain synthesis genes. We show that betalain colouration is specifically induced in root tissues and cells where fungal colonisation has occurred. In a rhizotron setup, we also demonstrate that betalain staining allows for the noninvasive tracing of fungal colonisation along the root system over time. We present MycoRed, a useful innovative method that will expand and complement currently used fungal visualisation techniques to advance knowledge in the field of AM symbiosis.

The link between ancient whole-genome duplications and cold adaptations in the Caryophyllaceae.

PREMISE: The Caryophyllaceae (the carnation family) have undergone multiple transitions into colder climates and convergence on cushion plant adaptation, indicating that they may provide a natural system for cold adaptation research. Previous research has suggested that putative ancient whole-genome duplications (WGDs) are correlated with niche shifts into colder climates across the Caryophyllales. Here, we explored the genomic changes potentially involved in one of these discovered shifts in the Caryophyllaceae. METHODS: We constructed a data set combining 26 newly generated transcriptomes with 45 published transcriptomes, including 11 cushion plant species across seven genera. With this data set, we inferred a dated phylogeny for the Caryophyllaceae and mapped ancient WGDs and gene duplications onto the phylogeny. We also examined functional groups enriched for gene duplications related to the climatic shift. RESULTS: The ASTRAL topology was mostly congruent with the current consensus of relationships within the family. We inferred 15 putative ancient WGDs in the family, including eight that have not been previously published. The oldest ancient WGD (ca. 64.4-56.7 million years ago), WGD1, was found to be associated with a shift into colder climates by previous research. Gene regions associated with ubiquitination were overrepresented in gene duplications retained after WGD1 and those convergently retained by cushion plants in Colobanthus and Eremogone, along with other functional annotations. CONCLUSIONS: Gene family expansions induced by ancient WGDs may have contributed to the shifts to cold climatic niches in the Caryophyllaceae. Transcriptomic data are crucial resources that help unravel heterogeneity in deep-time evolutionary patterns in plants.

A mycorrhiza-associated receptor-like kinase with an ancient origin in the green lineage.

Receptor-like kinases (RLKs) are key cell signaling components. The rice ARBUSCULAR RECEPTOR-LIKE KINASE 1 (OsARK1) regulates the arbuscular mycorrhizal (AM) association postarbuscule development and belongs to an undefined subfamily of RLKs. Our phylogenetic analysis revealed that ARK1 has an ancient paralogue in spermatophytes, ARK2 Single ark2 and ark1/ark2 double mutants in rice showed a nonredundant AM symbiotic function for OsARK2 Global transcriptomics identified a set of genes coregulated by the two RLKs, suggesting that OsARK1 and OsARK2 orchestrate symbiosis in a common pathway. ARK lineage proteins harbor a newly identified SPARK domain in their extracellular regions, which underwent parallel losses in ARK1 and ARK2 in monocots. This protein domain has ancient origins in streptophyte algae and defines additional overlooked groups of putative cell surface receptors.

Two independently evolved natural mutations additively deregulate TyrA enzymes and boost tyrosine production in planta.

l-Tyrosine is an essential amino acid for protein synthesis and is also used in plants to synthesize diverse natural products. Plants primarily synthesize tyrosine via TyrA arogenate dehydrogenase (TyrAa or ADH), which are typically strongly feedback inhibited by tyrosine. However, two plant lineages, Fabaceae (legumes) and Caryophyllales, have TyrA enzymes that exhibit relaxed sensitivity to tyrosine inhibition and are associated with elevated production of tyrosine-derived compounds, such as betalain pigments uniquely produced in core Caryophyllales. Although we previously showed that a single D222N substitution is primarily responsible for the deregulation of legume TyrAs, it is unknown when and how the deregulated Caryophyllales TyrA emerged. Here, through phylogeny-guided TyrA structure-function analysis, we found that functionally deregulated TyrAs evolved early in the core Caryophyllales before the origin of betalains, where the E208D amino acid substitution in the active site, which is at a different and opposite location from D222N found in legume TyrAs, played a key role in the TyrA functionalization. Unlike legumes, however, additional substitutions on non-active site residues further contributed to the deregulation of TyrAs in Caryophyllales. The introduction of a mutation analogous to E208D partially deregulated tyrosine-sensitive TyrAs, such as Arabidopsis TyrA2 (AtTyrA2). Moreover, the combined introduction of D222N and E208D additively deregulated AtTyrA2, for which the expression in Nicotiana benthamiana led to highly elevated accumulation of tyrosine in planta. The present study demonstrates that phylogeny-guided characterization of key residues underlying primary metabolic innovations can provide powerful tools to boost the production of essential plant natural products.

Are seven amino acid substitutions sufficient to explain the evolution of high l-DOPA 4,5-dioxygenase activity leading to betalain pigmentation? Revisiting the gain-of-function mutants of Bean et al. (2018).

This work revisits a publication by Bean et al. (2018) that reports seven amino acid substitutions are essential for the evolution of l-DOPA 4,5-dioxygenase (DODA) activity in Caryophyllales. In this study, we explore several concerns which led us to replicate the analyses of Bean et al. (2018). Our comparative analyses, with structural modelling, implicate numerous residues additional to those identified by Bean et al. (2018), with many of these additional residues occurring around the active site of BvDODAα1. We therefore replicated the analyses of Bean et al. (2018) to re-observe the effect of their original seven residue substitutions in a BvDODAα2 background, that is the BvDODAα2-mut3 variant. Multiple in vivo assays, in both Saccharomyces cerevisiae and Nicotiana benthamiana, did not result in visible DODA activity in BvDODAα2-mut3, with betalain production always 10-fold below BvDODAα1. In vitro assays also revealed substantial differences in both catalytic activity and pH optima between BvDODAα1, BvDODAα2 and BvDODAα2-mut3 proteins, explaining their differing performance in vivo. In summary, we were unable to replicate the in vivo analyses of Bean et al. (2018), and our quantitative in vivo and in vitro analyses suggest a minimal effect of these seven residues in altering catalytic activity of BvDODAα2. We conclude that the evolutionary pathway to high DODA activity is substantially more complex than implied by Bean et al. (2018).

The evidence for anthocyanins in the betalain-pigmented genus Hylocereus is weak.

Here we respond to Zhou (BMC Genomics 21:734, 2020) "Combined Transcriptome and Metabolome analysis of Pitaya fruit unveiled the mechanisms underlying peel and pulp color formation" published in BMC Genomics. Given the evolutionary conserved anthocyanin biosynthesis pathway in betalain-pigmented species, we are open to the idea that species with both anthocyanins and betalains might exist. However, in absence of LC-MS/MS spectra, apparent lack of biological replicates, and no comparison to authentic standards, the findings of Zhou (BMC Genomics 21:734, 2020) are not a strong basis to propose the presence of anthocyanins in betalain-pigmented pitaya. In addition, our re-analysis of the datasets indicates the misidentification of important genes and the omission of key flavonoid and anthocyanin synthesis genes ANS and DFR. Finally, our re-analysis of the RNA-Seq dataset reveals no correlation between anthocyanin biosynthesis gene expression and pigment status.

Conical petal epidermal cells, regulated by the MYB transcription factor MIXTA, have an ancient origin within the angiosperms.

Conical epidermal cells occur on the tepals (perianth organs, typically petals and/or sepals) of the majority of animal-pollinated angiosperms, where they play both visual and tactile roles in pollinator attraction, providing grip to foraging insects, and enhancing colour, temperature, and hydrophobicity. To explore the evolutionary history of conical epidermal cells in angiosperms, we surveyed the tepal epidermis in representative species of the ANA-grade families, the early-diverging successive sister lineages to all other extant angiosperms, and analysed the function of a candidate regulator of cell outgrowth from Cabomba caroliniana (Nymphaeales). We identified conical cells in at least two genera from different families (Austrobaileya and Cabomba). A single SBG9 MYB gene was isolated from C. caroliniana and found to induce strong differentiation of cellular outgrowth, including conical cells, when ectopically expressed in Nicotiana tabacum. Ontogenetic analysis and quantitative reverse transcription-PCR established that CcSBG9A1 is spatially and temporally expressed in a profile which correlates with a role in conical cell development. We conclude that conical or subconical cells on perianth organs are ancient within the angiosperms and most probably develop using a common genetic programme initiated by a SBG9 MYB transcription factor.

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l.

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.].

The report of anthocyanins in the betalain-pigmented genus Hylocereus is not well evidenced and is not a strong basis to refute the mutual exclusion paradigm.

Here we respond to the paper entitled "Contribution of anthocyanin pathways to fruit flesh coloration in pitayas" (Fan et al., BMC Plant Biol 20:361, 2020). In this paper Fan et al. 2020 propose that the anthocyanins can be detected in the betalain-pigmented genus Hylocereus, and suggest they are responsible for the colouration of the fruit flesh. We are open to the idea that, given the evolutionary maintenance of fully functional anthocyanin synthesis genes in betalain-pigmented species, anthocyanin pigmentation might co-occur with betalain pigments, as yet undetected, in some species. However, in absence of the LC-MS/MS spectra and co-elution/fragmentation of the authentic standard comparison, the findings of Fan et al. 2020 are not credible. Furthermore, our close examination of the paper, and re-analysis of datasets that have been made available, indicate numerous additional problems. Namely, the failure to detect betalains in an untargeted metabolite analysis, accumulation of reported anthocyanins that does not correlate with the colour of the fruit, absence of key anthocyanin synthesis genes from qPCR data, likely mis-identification of key anthocyanin genes, unreproducible patterns of correlated RNAseq data, lack of gene expression correlation with pigmentation accumulation, and putative transcription factors that are weak candidates for transcriptional up-regulation of the anthocyanin pathway.

The land plant-specific MIXTA-MYB lineage is implicated in the early evolution of the plant cuticle and the colonization of land.

The evolution of a lipid-based cuticle on aerial plant surfaces that protects against dehydration is considered a fundamental innovation in the colonization of the land by the green plants. However, key evolutionary steps in the early regulation of cuticle synthesis are still poorly understood, owing to limited studies in early-diverging land plant lineages. Here, we characterize a land plant specific subgroup 9 R2R3 MYB transcription factor MpSBG9, in the early-diverging land plant model Marchantia polymorpha, that is homologous to MIXTA proteins in vascular plants. The MpSBG9 functions as a key regulator of cuticle biosynthesis by preferentially regulating expression of orthologous genes for cutin formation, but not wax biosynthesis genes. The MpSBG9 also promotes the formation of papillate cells on the adaxial surface of M. polymorpha, which is consisitent with its canonical role in vascular plants. Our observations imply conserved MYB transcriptional regulation in the control of the cutin biosynthesis pathway as a core genetic network in the common ancestor of all land plants, implicating the land plant-specific MIXTA MYB lineage in the early origin and evolution of the cuticle.

Evolution of Portulacineae Marked by Gene Tree Conflict and Gene Family Expansion Associated with Adaptation to Harsh Environments.

Several plant lineages have evolved adaptations that allow survival in extreme and harsh environments including many families within the plant clade Portulacineae (Caryophyllales) such as the Cactaceae, Didiereaceae, and Montiaceae. Here, using newly generated transcriptomic data, we reconstructed the phylogeny of Portulacineae and examined potential correlates between molecular evolution and adaptation to harsh environments. Our phylogenetic results were largely congruent with previous analyses, but we identified several early diverging nodes characterized by extensive gene tree conflict. For particularly contentious nodes, we present detailed information about the phylogenetic signal for alternative relationships. We also analyzed the frequency of gene duplications, confirmed previously identified whole genome duplications (WGD), and proposed a previously unidentified WGD event within the Didiereaceae. We found that the WGD events were typically associated with shifts in climatic niche but did not find a direct association with WGDs and diversification rate shifts. Diversification shifts occurred within the Portulacaceae, Cactaceae, and Anacampserotaceae, and whereas these did not experience WGDs, the Cactaceae experienced extensive gene duplications. We examined gene family expansion and molecular evolutionary patterns with a focus on genes associated with environmental stress responses and found evidence for significant gene family expansion in genes with stress adaptation and clades found in extreme environments. These results provide important directions for further and deeper examination of the potential links between molecular evolutionary patterns and adaptation to harsh environments.

Evolution of l-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales.

The evolution of l-DOPA 4,5-dioxygenase activity, encoded by the gene DODA, was a key step in the origin of betalain biosynthesis in Caryophyllales. We previously proposed that l-DOPA 4,5-dioxygenase activity evolved via a single Caryophyllales-specific neofunctionalisation event within the DODA gene lineage. However, this neofunctionalisation event has not been confirmed and the DODA gene lineage exhibits numerous gene duplication events, whose evolutionary significance is unclear. To address this, we functionally characterised 23 distinct DODA proteins for l-DOPA 4,5-dioxygenase activity, from four betalain-pigmented and five anthocyanin-pigmented species, representing key evolutionary transitions across Caryophyllales. By mapping these functional data to an updated DODA phylogeny, we then explored the evolution of l-DOPA 4,5-dioxygenase activity. We find that low l-DOPA 4,5-dioxygenase activity is distributed across the DODA gene lineage. In this context, repeated gene duplication events within the DODA gene lineage give rise to polyphyletic occurrences of elevated l-DOPA 4,5-dioxygenase activity, accompanied by convergent shifts in key functional residues and distinct genomic patterns of micro-synteny. In the context of an updated organismal phylogeny and newly inferred pigment reconstructions, we argue that repeated convergent acquisition of elevated l-DOPA 4,5-dioxygenase activity is consistent with recurrent specialisation to betalain synthesis in Caryophyllales.

Plastid phylogenomic insights into the evolution of Caryophyllales.

The Caryophyllales includes 40 families and 12,500 species, representing a large and diverse clade of angiosperms. Collectively, members of the clade grow on all continents and in all terrestrial biomes and often occupy extreme habitats (e.g., xeric, salty). The order is characterized by many taxa with unusual adaptations including carnivory, halophytism, and multiple origins of C4 photosynthesis. However, deep phylogenetic relationships within the order have long been problematic due to putative rapid divergence. To resolve the deep-level relationships of Caryophyllales, we performed phylogenomic analyses of all 40 families of Caryophyllales. We time-calibrated the molecular phylogeny of this clade, and evaluated putative correlations among plastid structural changes and rates of molecular substitution. We recovered a well-resolved and well-supported phylogeny of the Caryophyllales that was largely congruent with previous estimates of this order. Our results provide improved support for the phylogenetic position of several key families within this clade. The crown age of Caryophyllales was estimated at ca. 114.4 million years ago (Ma), with periods of rapid divergence in the mid-Cretaceous. A strong, positive correlation between nucleotide substitution rate and plastid structural changes was detected. Our study highlights the importance of broad taxon sampling in phylogenomic inference and provides a firm basis for future investigations of molecular, morphological, and ecophysiological evolution in Caryophyllales.

Genome-wide analyses supported by RNA-Seq reveal non-canonical splice sites in plant genomes.

BACKGROUND: Most eukaryotic genes comprise exons and introns thus requiring the precise removal of introns from pre-mRNAs to enable protein biosynthesis. U2 and U12 spliceosomes catalyze this step by recognizing motifs on the transcript in order to remove the introns. A process which is dependent on precise definition of exon-intron borders by splice sites, which are consequently highly conserved across species. Only very few combinations of terminal dinucleotides are frequently observed at intron ends, dominated by the canonical GT-AG splice sites on the DNA level. RESULTS: Here we investigate the occurrence of diverse combinations of dinucleotides at predicted splice sites. Analyzing 121 plant genome sequences based on their annotation revealed strong splice site conservation across species, annotation errors, and true biological divergence from canonical splice sites. The frequency of non-canonical splice sites clearly correlates with their divergence from canonical ones indicating either an accumulation of probably neutral mutations, or evolution towards canonical splice sites. Strong conservation across multiple species and non-random accumulation of substitutions in splice sites indicate a functional relevance of non-canonical splice sites. The average composition of splice sites across all investigated species is 98.7% for GT-AG, 1.2% for GC-AG, 0.06% for AT-AC, and 0.09% for minor non-canonical splice sites. RNA-Seq data sets of 35 species were incorporated to validate non-canonical splice site predictions through gaps in sequencing reads alignments and to demonstrate the expression of affected genes. CONCLUSION: We conclude that bona fide non-canonical splice sites are present and appear to be functionally relevant in most plant genomes, although at low abundance.

Improved transcriptome sampling pinpoints 26 ancient and more recent polyploidy events in Caryophyllales, including two allopolyploidy events.

Studies of the macroevolutionary legacy of polyploidy are limited by an incomplete sampling of these events across the tree of life. To better locate and understand these events, we need comprehensive taxonomic sampling as well as homology inference methods that accurately reconstruct the frequency and location of gene duplications. We assembled a data set of transcriptomes and genomes from 168 species in Caryophyllales, of which 43 transcriptomes were newly generated for this study, representing one of the most densely sampled genomic-scale data sets available. We carried out phylogenomic analyses using a modified phylome strategy to reconstruct the species tree. We mapped the phylogenetic distribution of polyploidy events by both tree-based and distance-based methods, and explicitly tested scenarios for allopolyploidy. We identified 26 ancient and more recent polyploidy events distributed throughout Caryophyllales. Two of these events were inferred to be allopolyploidy. Through dense phylogenomic sampling, we show the propensity of polyploidy throughout the evolutionary history of Caryophyllales. We also provide a framework for utilizing transcriptome data to detect allopolyploidy, which is important as it may have different macroevolutionary implications compared with autopolyploidy.

Relaxation of tyrosine pathway regulation underlies the evolution of betalain pigmentation in Caryophyllales.

Diverse natural products are synthesized in plants by specialized metabolic enzymes, which are often lineage-specific and derived from gene duplication followed by functional divergence. However, little is known about the contribution of primary metabolism to the evolution of specialized metabolic pathways. Betalain pigments, uniquely found in the plant order Caryophyllales, are synthesized from the aromatic amino acid l-tyrosine (Tyr) and replaced the otherwise ubiquitous phenylalanine-derived anthocyanins. This study combined biochemical, molecular and phylogenetic analyses, and uncovered coordinated evolution of Tyr and betalain biosynthetic pathways in Caryophyllales. We found that Beta vulgaris, which produces high concentrations of betalains, synthesizes Tyr via plastidic arogenate dehydrogenases (TyrAa /ADH) encoded by two ADH genes (BvADHα and BvADHß). Unlike BvADHß and other plant ADHs that are strongly inhibited by Tyr, BvADHα exhibited relaxed sensitivity to Tyr. Also, Tyr-insensitive BvADHα orthologs arose during the evolution of betalain pigmentation in the core Caryophyllales and later experienced relaxed selection and gene loss in lineages that reverted from betalain to anthocyanin pigmentation, such as Caryophyllaceae. These results suggest that relaxation of Tyr pathway regulation increased Tyr production and contributed to the evolution of betalain pigmentation, highlighting the significance of upstream primary metabolic regulation for the diversification of specialized plant metabolism.

Disparity, diversity, and duplications in the Caryophyllales.

The role played by whole genome duplication (WGD) in plant evolution is actively debated. WGDs have been associated with advantages such as superior colonization, various adaptations, and increased effective population size. However, the lack of a comprehensive mapping of WGDs within a major plant clade has led to uncertainty regarding the potential association of WGDs and higher diversification rates. Using seven chloroplast and nuclear ribosomal genes, we constructed a phylogeny of 5036 species of Caryophyllales, representing nearly half of the extant species. We phylogenetically mapped putative WGDs as identified from analyses on transcriptomic and genomic data and analyzed these in conjunction with shifts in climatic occupancy and lineage diversification rate. Thirteen putative WGDs and 27 diversification shifts could be mapped onto the phylogeny. Of these, four WGDs were concurrent with diversification shifts, with other diversification shifts occurring at more recent nodes than WGDs. Five WGDs were associated with shifts to colder climatic occupancy. While we find that many diversification shifts occur after WGDs, it is difficult to consider diversification and duplication to be tightly correlated. Our findings suggest that duplications may often occur along with shifts in either diversification rate, climatic occupancy, or rate of evolution.

From cacti to carnivores: Improved phylotranscriptomic sampling and hierarchical homology inference provide further insight into the evolution of Caryophyllales.

PREMISE OF THE STUDY: The Caryophyllales contain ~12,500 species and are known for their cosmopolitan distribution, convergence of trait evolution, and extreme adaptations. Some relationships within the Caryophyllales, like those of many large plant clades, remain unclear, and phylogenetic studies often recover alternative hypotheses. We explore the utility of broad and dense transcriptome sampling across the order for resolving evolutionary relationships in Caryophyllales. METHODS: We generated 84 transcriptomes and combined these with 224 publicly available transcriptomes to perform a phylogenomic analysis of Caryophyllales. To overcome the computational challenge of ortholog detection in such a large data set, we developed an approach for clustering gene families that allowed us to analyze >300 transcriptomes and genomes. We then inferred the species relationships using multiple methods and performed gene-tree conflict analyses. KEY RESULTS: Our phylogenetic analyses resolved many clades with strong support, but also showed significant gene-tree discordance. This discordance is not only a common feature of phylogenomic studies, but also represents an opportunity to understand processes that have structured phylogenies. We also found taxon sampling influences species-tree inference, highlighting the importance of more focused studies with additional taxon sampling. CONCLUSIONS: Transcriptomes are useful both for species-tree inference and for uncovering evolutionary complexity within lineages. Through analyses of gene-tree conflict and multiple methods of species-tree inference, we demonstrate that phylogenomic data can provide unparalleled insight into the evolutionary history of Caryophyllales. We also discuss a method for overcoming computational challenges associated with homolog clustering in large data sets.

Widespread paleopolyploidy, gene tree conflict, and recalcitrant relationships among the carnivorous Caryophyllales.

PREMISE OF STUDY: The carnivorous members of the large, hyperdiverse Caryophyllales (e.g., Venus flytrap, sundews, and Nepenthes pitcher plants) represent perhaps the oldest and most diverse lineage of carnivorous plants. However, despite numerous studies seeking to elucidate their evolutionary relationships, the early-diverging relationships remain unresolved. METHODS: To explore the utility of phylogenomic data sets for resolving relationships among the carnivorous Caryophyllales, we sequenced 10 transcriptomes, including all the carnivorous genera except those in the rare West African liana family Dioncophyllaceae. We used a variety of methods to infer the species tree, examine gene tree conflict, and infer paleopolyploidy events. KEY RESULTS: Phylogenomic analyses supported the monophyly of the carnivorous Caryophyllales, with a crown age of 68-83 million years. In contrast to previous analyses, we recovered the remaining noncore Caryophyllales as nonmonophyletic, although the node supporting this relationship contained a significant amount of gene tree discordance. We present evidence that the clade contains at least seven independent paleopolyploidy events, previously unresolved nodes from the literature have high levels of gene tree conflict, and taxon sampling influences topology even in a phylogenomic data set, regardless of the use of coalescent or supermatrix methods. CONCLUSIONS: Our data demonstrate the importance of carefully considering gene tree conflict and taxon sampling in phylogenomic analyses. Moreover, they provide a remarkable example of the propensity for paleopolyploidy in angiosperms, with at least seven such events in a clade of less than 2500 species.