Plant Tree Of Life Research Articles

Mapping of specialized metabolite terms onto a plant phylogeny using text mining and large language models.

Plants produce a staggering array of chemicals that are the basis for organismal function and important human nutrients and medicines. However, it is poorly defined how these compounds evolved and are distributed across the plant kingdom, hindering a systematic view and understanding of plant chemical diversity. Recent advances in plant genome/transcriptome sequencing have provided a well-defined molecular phylogeny of plants, on which the presence of diverse natural products can be mapped to systematically determine their phylogenetic distribution. Here, we built a proof-of-concept workflow where previously reported diverse tyrosine-derived plant natural products were mapped onto the plant tree of life. Plant chemical-species associations were mined from literature, filtered, evaluated through manual inspection of over 2500 scientific articles, and mapped onto the plant phylogeny. The resulting "phylochemical" map confirmed several highly lineage-specific compound class distributions, such as betalain pigments and Amaryllidaceae alkaloids. The map also highlighted several lineages enriched in dopamine-derived compounds, including the orders Caryophyllales, Liliales, and Fabales. Additionally, the application of large language models, using our manually curated data as a ground truth set, showed that post-mining processing can largely be automated with a low false-positive rate, critical for generating a reliable phylochemical map. Although a high false-negative rate remains a challenge, our study demonstrates that combining text mining with language model-based processing can generate broader phylochemical maps, which will serve as a valuable community resource to uncover key evolutionary events that underlie plant chemical diversity and enable system-level views of nature's millions of years of chemical experimentation.

Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

IRplus: An Augmented Tool to Detect Inverted Repeats in Plastid Genomes.

High-throughput sequencing methods have increased the accessibility of plastid genomes, which are crucial for clarifying phylogenetic relationships. Current large sequencing efforts require software tools for routine display of their distinctive quadripartite structure, which is denoted by four junction sites. By concentrating on these junctions and their close vicinity, IRscope has emerged as the standard tool for detection of this structure and creating simplified comparative graphical maps of plastid genomes. Here, we provide an augmented version (IRplus) that encompasses a novel set of functions such as integrated error detection, flexible color schemes, and an upgraded method to detect inverted repeats in genomic sequences. Spanning across the plant tree of life, IRplus allows the quick visualization of various sets of plastid genomes and features, next to smooth interoperability with other widely used annotation file formats and platforms. The IRplus can be accessed at https://irscope.shinyapps.io/IRplus/, and source codes are freely available at https://github.com/AmiryousefiLab/IRplus.

Differential timing of gene expression and recruitment in independent origins of CAM in the Agavoideae (Asparagaceae).

Crassulacean acid metabolism (CAM) photosynthesis has evolved repeatedly across the plant tree of life, however our understanding of the genetic convergence across independent origins remains hampered by the lack of comparative studies. Here, we explore gene expression profiles in eight species from the Agavoideae (Asparagaceae) encompassing three independent origins of CAM. Using comparative physiology and transcriptomics, we examined the variable modes of CAM in this subfamily and the changes in gene expression across time of day and between well watered and drought-stressed treatments. We further assessed gene expression and the molecular evolution of genes encoding phosphoenolpyruvate carboxylase (PPC), an enzyme required for primary carbon fixation in CAM. Most time-of-day expression profiles are largely conserved across all eight species and suggest that large perturbations to the central clock are not required for CAM evolution. By contrast, transcriptional response to drought is highly lineage specific. Yucca and Beschorneria have CAM-like expression of PPC2, a copy of PPC that has never been shown to be recruited for CAM in angiosperms. Together the physiological and transcriptomic comparison of closely related C3 and CAM species reveals similar gene expression profiles, with the notable exception of differential recruitment of carboxylase enzymes for CAM function.

Exploring Zosterophyll Relationships within a More Broadly Sampled Character Space: A Focus on Anatomy

Premise of research. Zosterophylls are a major constituent of Siluro-Devonian tracheophyte floras worldwide and gave rise to the morphologically diverse lycopsid clade. Despite this pivotal position in the plant tree of life, relationships within the group have remained poorly understood. Furthermore, anatomical data have been undersampled in previous phylogenetic analyses including zosterophylls. With this in mind, we addressed the relationships among zosterophylls using a data set that maximized anatomical characters and employing phylogenetic and phenetic methods.Methodology. A matrix of 40 total characters (12 of which are anatomical) and 19 zosterophyll species with known anatomy, representing 17 genera, was compiled and analyzed using parsimony-constrained phylogenetic methods, clustering, and ordination methods.Pivotal results. Phenetic analyses show sensitivity to taxon sampling and support the placement of Renalia among the zosterophylls but do not support taxonomic inferences strongly congruent with those supported by phylogenetic analyses. Phylogenetic analyses consistently recover two major clades: one lacks internal resolution and comprises the bulk of the zosterophyll taxa included in the analyses; the other clade includes the zosterophyll Ventarura and the lycopsid Sengelia, often accompanied by Discalis and Trichopherophyton. Analyses using subsets of characters (only morphology or anatomy) recover trees that differ from those obtained using the total set of characters (morphology and anatomy).Conclusions. The incongruence between the results of the total character analyses and those using only morphology or only anatomy highlights the importance of broadening the sampling of morphological character space. Because both anatomy and morphology are part of the identity and evolutionary history of a species, the relationships recovered by the inclusion of both morphological and anatomical characters are more likely to reflect natural evolutionary relationships. Breadth of character sampling and not the amount of phylogenetic resolution should be the primary criterion for selecting among alternative hypotheses of relationships.

Macroevolutionary consequences of mast seeding.

Masting characterizes large, intermittent and highly synchronous seeding events among individual plants and is found throughout the plant Tree of Life (ToL). Although masting can increase plant fitness, little is known about whether it results in evolutionary changes across entire clades, such as by promoting speciation or enhanced trait selection. Here, we tested if masting has macroevolutionary consequences by combining the largest existing dataset of population-level reproductive time series and time-calibrated phylogenetic tree of vascular plants. We found that the coefficient of variation (CVp) of reproductive output for 307 species covaried with evolutionary history, and more so within clades than expected by random. Speciation rates estimated at the species level were highest at intermediate values of CVp and regional-scale synchrony (Sr) in seed production, that is, there were unimodal correlations. There was no support for monotonic correlations between either CVp or Sr and rates of speciation or seed size evolution. These results were robust to different sampling decisions, and we found little bias in our dataset compared with the wider plant ToL. While masting is often adaptive and encompasses a rich diversity of reproductive behaviours, we suggest it may have few consequences beyond the species level.This article is part of the theme issue ‘The ecology and evolution of synchronized seed production in plants’.

Trajectories of Homoeolog-Specific Expression in Allotetraploid Tragopogon castellanus Populations of Independent Origins.

Polyploidization can have a significant ecological and evolutionary impact by providing substantially more genetic material that may result in novel phenotypes upon which selection may act. While the effects of polyploidization are broadly reviewed across the plant tree of life, the reproducibility of these effects within naturally occurring, independently formed polyploids is poorly characterized. The flowering plant genus Tragopogon (Asteraceae) offers a rare glimpse into the intricacies of repeated allopolyploid formation with both nascent (< 90 years old) and more ancient (mesopolyploids) formations. Neo- and mesopolyploids in Tragopogon have formed repeatedly and have extant diploid progenitors that facilitate the comparison of genome evolution after polyploidization across a broad span of evolutionary time. Here, we examine four independently formed lineages of the mesopolyploid Tragopogon castellanus for homoeolog expression changes and fractionation after polyploidization. We show that expression changes are remarkably similar among these independently formed polyploid populations with large convergence among expressed loci, moderate convergence among loci lost, and stochastic silencing. We further compare and contrast these results for T. castellanus with two nascent Tragopogon allopolyploids. While homoeolog expression bias was balanced in both nascent polyploids and T. castellanus, the degree of additive expression was significantly different, with the mesopolyploid populations demonstrating more non-additive expression. We suggest that gene dosage and expression noise minimization may play a prominent role in regulating gene expression patterns immediately after allopolyploidization as well as deeper into time, and these patterns are conserved across independent polyploid lineages.

Plastid phylogenomics provides novel insights into the infrafamilial relationship of Polypodiaceae.

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

R2R3-MYBGene Evolution in Plants, Incorporating Ferns into the Story

Premise of research. MYB transcription factors are one of the largest families of genes in land plants. They are involved in many regulatory processes, such as primary and secondary metabolism, cell fate and identity, developmental processes, and responses to biotic and abiotic stresses. Phylogenetic studies have mainly focused on R2R3-MYB evolution in seed plants; however, a comprehensive sampling across land plants is lacking, as none have included ferns. Methodology. To better understand the evolution of R2R3-MYB genes in land plants, we surveyed the R2R3-MYB gene sequences from six genomes belonging to the main lineages of land plants (bryophytes, lycophytes, ferns, angiosperms, and gymnosperms). In addition, we searched transcriptome sequences of selected fern and lycophyte transcriptomes. We assembled a nucleotide matrix of the MYB domain and conducted maximum likelihood analyses to infer phylogenetic relationships and gene-tree reconciliation analyses to infer gene duplications. We labeled the main clades with the known MYB functions onto the resulting tree. To detect reported as well as new conserved motifs, we ran protein motif analyses through MEME. Pivotal results. Our results support previous studies indicating that R2R3-MYBs in seed plants are more diverse than in any other lineage of vascular plants. Most of the well-supported clades recovered are inferred to be already present in the ancestor of land plants. We found that ferns have numerous copies of R2R3-MYBs, although not as many as seed plants. Protein motif analyses revealed that R2R3-MYB motifs are highly conserved across land plants, suggesting that most R2R3-MYB orthologues, including those of ferns, might have DNA-binding capabilities and could be involved in regulatory processes that are similar to those of their angiosperm homologues. Our analyses showed that there are no fern orthologues of the angiosperm MYB leaf developmental genes ARP, raising interesting questions about the evolution and origin of leaves in ferns and seed plants. Conclusions. Our results provide the phylogenetic context for research on the genetic and functional evolution of an important gene family of developmental and metabolic regulators across the plant tree of life.

Polyploids increase overall diversity despite higher turnover than diploids in the Brassicaceae.

Although polyploidy is widespread across the plant Tree of Life, its long-term evolutionary significance is still poorly understood. Here, we examine the effects of polyploidy in explaining the large-scale evolutionary patterns within angiosperms by focusing on a single family exhibiting extensive interspecific variation in chromosome numbers. We inferred ploidy from haploid chromosome numbers for 80% of species in the most comprehensive species-level chronogram for the Brassicaceae. After evaluating a total of 94 phylogenetic models of diversification, we found that ploidy influences diversification rates across the Brassicaceae. We also found that despite diversifying at a similar rate to diploids, polyploids have played a significant role in driving present-day differences in species richness among clades. Overall, in addition to highlighting the complexity in the evolutionary consequences of polyploidy, our results suggest that rare successful polyploids persist while significantly contributing to the long-term evolution of clades. Our findings further indicate that polyploidy has played a major role in driving the long-term evolution of the Brassicaceae and highlight the potential of polyploidy in shaping present-day diversity patterns across the plant Tree of Life.

Biogeography and phylogeny of masting: do global patterns fit functional hypotheses?

Interannual variability of seed crops (CVp) has profound consequences for plant populations and food webs, where high CVp is termed 'masting'. Here we ask: is global variation in CVp better predicted by plant or habitat differences consistent with adaptive economies of scale, in which flower and seed benefits increase disproportionately during mast years; or by passive mechanisms, in which seed production responds to variation in resource availability associated with climate variability? To address this question, we compiled a dataset for phylogenetic comparative analysis of long-term fruit/seed production for plants comprising 920 time series spanning 311 plant species. Factors associated with both adaptive benefits of CVp (wind pollination and seed dispersal) and climatic variability (variability of summer precipitation) were among the best predictors of global variation in CVp. We observed a hump-shaped relationship between CVp and latitude and intermediate phylogenetic and geographic signals in CVp. CVp is patterned nonrandomly across the globe and over the plant tree of life, where high CVp is associated with species benefiting from economies of scale of seed or flower production and with species that experience variable rainfall over summer months when seeds usually mature.

Does Evolutionary History Correlate with Contemporary Extinction Risk by Influencing Range Size Dynamics?

Extinction threatens many species yet is predicted by few factors across the plant tree of life (ToL). Taxon age is one factor that may associate with extinction if occupancy of geographic and adaptive zones varies with time, but evidence for such an association has been equivocal. Age-dependent occupancy can also influence diversification rates and thus extinction risk where new taxa have small range and population sizes. To test how age, diversification, and range size were correlated with extinction, we analyzed 639 well-sampled genera representing 8,937 species from across the plant ToL. We found a greater proportion of species were threatened by contemporary extinction in younger and faster-diversifying genera. When we directly tested how range size mediated this pattern in two large, well-sampled groups, our results varied. In conifers, potential range size was smaller in older species and was correlated with higher extinction risk. Age on its own had no direct effect on extinction when accounting for its influence on range size. In palm species, age was neither directly nor indirectly correlated with extinction risk. Our results suggest that range size dynamics may explain differing patterns of extinction risk across the ToL, with consequences for biodiversity conservation.

Correlates of hybridization in plants.

Hybridization is a biological phenomenon increasingly recognized as an important evolutionary process in both plants and animals, as it is linked to speciation, radiation, extinction, range expansion and invasion, and allows for increased trait diversity in agricultural and horticultural systems. Estimates of hybridization frequency vary across taxonomic groups, but causes of this variation are unknown. Here, we ask on a global scale whether hybridization is linked to any of 11 traits related to plant life history, reproduction, genetic predisposition, and environment or opportunity. Given that hybridization is not evenly distributed across the plant tree of life, we use phylogenetic generalized least squares regression models and phylogenetic path analysis to detect statistical associations between hybridization and plant traits at both the family and genus levels. We find that perenniality and woodiness are each weakly associated with an increased frequency of hybridization in univariate analyses, but path analysis suggests that the direct linkage is between perenniality and increased hybridization (with woodiness having only an indirect relationship with hybridization via perenniality). Weak associations between higher rates of hybridization and higher outcrossing rates, abiotic pollination syndromes, vegetative reproductive modes, larger genomes, and less variable genome sizes are detectable in some cases but not others. We argue that correlational evidence at the global scale, such as that presented here, provides a robust framework for forming hypotheses to examine and test drivers of hybridization at a more mechanistic level.

Geophytism in monocots leads to higher rates of diversification.

Geophytes, plants with buds on underground structures, are found throughout the plant tree of life. These below ground structures allow plants to inhabit highly seasonal and disturbance-prone environments across ecosystems. Past researchers have hypothesised that the bulbous, cormous and tuberous habits promote diversification, but this had yet to be tested. Using a comprehensive monocot data set of almost 13000 taxa, we investigated the effects of the geophytic habit on diversification using both state-dependent and state-independent models. We found that geophytes exhibit increased rates of diversification relative to nongeophytes. State-dependent analyses recovered higher yet similar rates of diversification for bulbous, cormous and tuberous taxa compared with rhizomatous and nongeophytic taxa. However, the state-independent model returned no difference in rates among the different traits. Geophytism shows higher rates of diversification relative to nongeophytes but we found little support for the hypothesis that the evolution of the bulb, corm or tuber appears to provide a diversification increase relative to rhizomatous and nongeophytic taxa. Our broad-scale analysis highlights the overall evolutionary importance of the geophytic habit (i.e. belowground bud placement). However, our results also suggest that belowground morphological diversity alone cannot explain this rate increase. In order to further test the evolutionary significance of these underground structures, future studies should consider these in combination with other biotic and abiotic factors.

Methods for exploring the plant tree of life

Methods for exploring the plant tree of life

Chromploid: An R package for chromosome number evolution across the plant tree of life.

Premise of the StudyPolyploidy has profound evolutionary consequences for land plants. Despite the availability of large phylogenetic and chromosomal data sets, estimating the rates of polyploidy and chromosomal evolution across the tree of life remains a challenging, computationally complex problem. We introduce the R package chromploid, which allows scientists to perform inference of chromosomal evolution rates across large phylogenetic trees.Methods and Resultschromploid is an open‐source package in the R environment that calculates the likelihood function of models of chromosome evolution. Models of discrete character evolution can be customized using chromploid. We demonstrate the performance of the BiChroM model, testing for associations between rates of chromosome doubling (as a proxy for polyploidy) and a binary phenotypic character, within chromploid using simulations and empirical data from Solanum. In simulations, estimated chromosome‐doubling rates were unbiased and the variance decreased with larger trees, but distinguishing small differences in rates of chromosome doubling, even from large data sets, remains challenging. In the Solanum data set, a custom model of chromosome number evolution demonstrated higher rates of chromosome doubling in herbaceous species compared to woody.Conclusionschromploid enables researchers to perform robust likelihood‐based inferences using complex models of chromosome number evolution across large phylogenies.

A novel experimental system using the liverwort Marchantia polymorpha and its fungal endophytes reveals diverse and context-dependent effects.

Fungal symbioses are ubiquitous in plants, but their effects have mostly been studied in seed plants. This study aimed to assess the diversity of fungal endophyte effects in a bryophyte and identify factors contributing to the variability of outcomes in these interactions. Fungal endophyte cultures and axenic liverwort clones were isolated from wild populations of the liverwort, Marchantia polymorpha. These collections were combined in a gnotobiotic system to test the effects of fungal isolates on the growth rates of hosts under laboratory conditions. Under the experimental conditions, fungi isolated from M.polymorpha ranged from aggressively pathogenic to strongly growth-promoting, but the majority of isolates caused no detectable change in host growth. Growth promotion by selected fungi depended on nutrient concentrations and was inhibited by coinoculation with multiple fungi. The M.polymorpha endophyte system expands the resources for this model liverwort. The experiments presented here demonstrate a wealth of diversity in fungal interactions even in a host reported to lack standard mycorrhizal symbiosis. In addition, they show that some known pathogens of vascular plants live in M.polymorpha and can confer benefits to this nonvascular host. This highlights the importance of studying endophyte effects across the plant tree of life.

JSE IBC Special Issue on Frontiers in Plant Systematics and Evolution

JSE IBC Special Issue on Frontiers in Plant Systematics and Evolution

A new resource for the development of SSR markers: Millions of loci from a thousand plant transcriptomes.

Premise of the study:The One Thousand Plant Transcriptomes Project (1KP, 1000+ assembled plant transcriptomes) provides an enormous resource for developing microsatellite loci across the plant tree of life. We developed loci from these transcriptomes and tested their utility.Methods and Results:Using software packages and custom scripts, we identified microsatellite loci in 1KP transcriptomes. We assessed the potential for cross-amplification and whether loci were biased toward exons, as compared to markers derived from genomic DNA. We characterized over 5.7 million simple sequence repeat (SSR) loci from 1334 plant transcriptomes. Eighteen percent of loci substantially overlapped with open reading frames (ORFs), and electronic PCR revealed that over half the loci would amplify successfully in conspecific taxa. Transcriptomic SSRs were approximately three times more likely to map to translated regions than genomic SSRs.Conclusions:We believe microsatellites still have a place in the genomic age—they remain effective and cost-efficient markers. The loci presented here are a valuable resource for researchers.

Genome-wide identification, evolutionary, and expression analyses of histone H3 variants in plants.

Histone variants alter the nucleosome structure and play important roles in chromosome segregation, transcription, DNA repair, and sperm compaction. Histone H3 is encoded by many genes in most eukaryotic species and is the histone that contains the largest variety of posttranslational modifications. Compared with the metazoan H3 variants, little is known about the complex evolutionary history of H3 variants proteins in plants. Here, we study the identification, evolutionary, and expression analyses of histone H3 variants from genomes in major branches in the plant tree of life. Firstly we identified all the histone three related (HTR) genes from the examined genomes, then we classified the four groups variants: centromeric H3, H3.1, H3.3 and H3-like, by phylogenetic analysis, intron information, and alignment. We further demonstrated that the H3 variants have evolved under strong purifying selection, indicating the conservation of HTR proteins. Expression analysis revealed that the HTR has a wide expression profile in maize and rice development and plays important roles in development.