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.

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).

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.

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.

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.

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.

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.

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.

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.

RNA-dependent RNA polymerase 1 in potato (Solanum tuberosum) and its relationship to other plant RNA-dependent RNA polymerases.

Cellular RNA-dependent RNA polymerases (RDRs) catalyze synthesis of double-stranded RNAs that can serve to initiate or amplify RNA silencing. Arabidopsis thaliana has six RDR genes; RDRs 1, 2 and 6 have roles in anti-viral RNA silencing. RDR6 is constitutively expressed but RDR1 expression is elevated following plant treatment with defensive phytohormones. RDR1 also contributes to basal virus resistance. RDR1 has been studied in several species including A. thaliana, tobacco (Nicotiana tabacum), N. benthamiana, N. attenuata and tomato (Solanum lycopersicum) but not to our knowledge in potato (S. tuberosum). StRDR1 was identified and shown to be salicylic acid-responsive. StRDR1 transcript accumulation decreased in transgenic potato plants constitutively expressing a hairpin construct and these plants were challenged with three viruses: potato virus Y, potato virus X, and tobacco mosaic virus. Suppression of StRDR1 gene expression did not increase the susceptibility of potato to these viruses. Phylogenetic analysis of RDR genes present in potato and in a range of other plant species identified a new RDR gene family, not present in potato and found only in Rosids (but apparently lost in the Rosid A. thaliana) for which we propose the name RDR7.

Rosid radiation and the rapid rise of angiosperm-dominated forests.

The rosid clade (70,000 species) contains more than one-fourth of all angiosperm species and includes most lineages of extant temperate and tropical forest trees. Despite progress in elucidating relationships within the angiosperms, rosids remain the largest poorly resolved major clade; deep relationships within the rosids are particularly enigmatic. Based on parsimony and maximum likelihood (ML) analyses of separate and combined 12-gene (10 plastid genes, 2 nuclear; >18,000 bp) and plastid inverted repeat (IR; 24 genes and intervening spacers; >25,000 bp) datasets for >100 rosid species, we provide a greatly improved understanding of rosid phylogeny. Vitaceae are sister to all other rosids, which in turn form 2 large clades, each with a ML bootstrap value of 100%: (i) eurosids I (Fabidae) include the nitrogen-fixing clade, Celastrales, Huaceae, Zygophyllales, Malpighiales, and Oxalidales; and (ii) eurosids II (Malvidae) include Tapisciaceae, Brassicales, Malvales, Sapindales, Geraniales, Myrtales, Crossosomatales, and Picramniaceae. The rosid clade diversified rapidly into these major lineages, possibly over a period of <15 million years, and perhaps in as little as 4 to 5 million years. The timing of the inferred rapid radiation of rosids [108 to 91 million years ago (Mya) and 107-83 Mya for Fabidae and Malvidae, respectively] corresponds with the rapid rise of angiosperm-dominated forests and the concomitant diversification of other clades that inhabit these forests, including amphibians, ants, placental mammals, and ferns.