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

Repeated shifts out of tropical climates preceded by whole genome duplication

While flowering plants have diversified in virtually every terrestrial clime, climate constrains the distribution of individual lineages. Overcoming climatic constraints may be associated with diverse evolutionary phenomena including whole genome duplication (WGD), gene-tree conflict, and life-history changes. Climatic shifts may also have facilitated increases in flowering plant diversification rates. We investigate climatic shifts in the flowering plant order Ericales, which consists of c. 14 000 species with diverse climatic tolerances. We estimate phylogenetic trees from transcriptomic data, 64 chloroplast loci, and Angiosperms353 nuclear loci that, respectively, incorporate 147, 4508, and 2870 Ericales species. We use these phylogenetic trees to analyse how climatic shifts are associated with WGD, gene-tree conflict, life-history, and diversification rates. Early branches in the phylogenetic trees are extremely short, and have high levels of gene-tree conflict and at least one WGD. On lineages descended from these early branches, there is a significant association between climatic shifts (primarily out of tropical climates), further WGDs, and life-history. Extremely short early branches, and their associated gene-tree conflict and WGDs, appear to underpin the explosive origin of numerous species rich Ericales clades. The evolution of diverse climatic tolerances in these species rich clades is tightly associated with WGD and life-history.

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.

An updated phylogeny and adaptive evolution within Amaranthaceae s.l. inferred from multiple phylogenomic datasets

Amaranthaceae s.l. is a widely distributed family consisting of over 170 genera and 2000 species. Previous molecular phylogenetic studies have shown that Amaranthaceae s.s. and traditional Chenopodiaceae form a monophyletic group (Amaranthaceae s.l.), however, the relationships within this evolutionary branch have yet to be fully resolved. In this study, we assembled the complete plastomes and full-length ITS of 21 Amaranthaceae s.l. individuals and compared them with 38 species of Amaranthaceae s.l. Through plastome structure and sequence alignment analysis, we identified a reverse complementary region approximately 5200 bp long in the genera Atriplex and Chenopodium. Adaptive evolution analysis revealed significant positive selection in eight genes, which likely played a driving role in the evolution of Amaranthaceae s.l., as demonstrated by partitioned evolutionary analysis. Furthermore, we found that about two-thirds of the examined species lack the ycf15 gene, potentially associated with natural selection pressures from their adapted habitats. The phylogenetic tree indicated that some genera (Chenopodium, Halogeton, and Subtr. Salsolinae) are paraphyletic lineages. Our results strongly support the clustering of Amaranthaceae s.l. with monophyletic traditional Chenopodiaceae (Clades I and II) and Amaranthaceae s.s. After a comprehensive analysis, we determined that cytonuclear conflict, gene selection by adapted habitats, and incomplete lineage sorting (ILS) events were the primary reasons for the inconsistent phylogeny of Amaranthaceae s.l. During the last glacial period, certain species within Amaranthaceae s.l. underwent adaptations to different environments and began to differentiate rapidly. Since then, these species may have experienced morphological and genetic changes distinct from those of other genera due to intense selection pressure.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Medullary bundles in Caryophyllales: form, function, and evolution

The occurrence of conducting vascular tissue in the pith (CVTP) of tracheophytes is noteworthy. Medullary bundles, one of the remarkable examples of CVTP, evolved multiple times across angiosperms, notably in the Caryophyllales. Yet, information on the occurrence of medullary bundles is fragmented, hampering our understanding of their structure-function relationships, and evolutionary implications. Using three plastid molecular markers (matK, rbcL, and rps16 intron), a phylogeny is constructed for 561 species of Caryophyllales, and anatomical data are assembled for 856 species across 40 families to investigate the diversity of medullary bundles, their function, evolution, and diversification dynamics. Additionally, correlated evolution between medullary bundles and successive cambia was tested. Medullary bundles are ancestrally absent in Caryophyllales and evolved in core and noncore families. They are structurally diverse (e.g. number, arrangement, and types of bundles) and functionally active throughout the plant's lifespan, providing increased hydraulic conductivity, especially in herbaceous plants. Acquisition of medullary bundles does not explain diversification rate heterogeneity but is correlated to a higher diversification rate. Disparate developmental pathways were found leading to rampant convergent evolution of CVTP in Caryophyllales. These findings indicate the diversification of medullary bundles and vascular tissues as another central theme for functional and comparative molecular studies in Caryophyllales.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Whole Genome Duplication Events Likely Contributed to the Aquatic Adaptive Evolution of Parkerioideae

As the only aquatic lineage of Pteridaceae, Parkerioideae is distinct from many xeric-adapted species of the family and consists of the freshwater Ceratopteris species and the only mangrove ferns from the genus Acrostichum. Previous studies have shown that whole genome duplication (WGD) has occurred in Parkerioideae at least once and may have played a role in their adaptive evolution; however, more in-depth research regarding this is still required. In this study, comparative and evolutionary transcriptomics analyses were carried out to identify WGDs and explore their roles in the environmental adaptation of Parkerioideae. Three putative WGD events were identified within Parkerioideae, two of which were specific to Ceratopteris and Acrostichum, respectively. The functional enrichment analysis indicated that the lineage-specific WGD events have played a role in the adaptation of Parkerioideae to the low oxygen concentrations of aquatic habitats, as well as different aquatic environments of Ceratopteris and Acrostichum, such as the adaptation of Ceratopteris to reduced light levels and the adaptation of Acrostichum to high salinity. Positive selection analysis further provided evidence that the putative WGD events may have facilitated the adaptation of Parkerioideae to changes in habitat. Moreover, the gene family analysis indicated that the plasma membrane H+-ATPase (AHA), vacuolar H+-ATPase (VHA), and suppressor of K+ transport growth defect 1 (SKD1) may have been involved in the high salinity adaptation of Acrostichum. Our study provides new insights into the evolution and adaptations of Parkerioideae in different aquatic environments.

Subgenome dominance shapes novel gene evolution in the decaploid pitcher plant Nepenthes gracilis

Subgenome dominance after whole-genome duplication generates distinction in gene number and expression at the level of chromosome sets, but it remains unclear how this process may be involved in evolutionary novelty. Here we generated a chromosome-scale genome assembly of the Asian pitcher plant Nepenthes gracilis to analyse how its novel traits (dioecy and carnivorous pitcher leaves) are linked to genomic evolution. We found a decaploid karyotype and a clear indication of subgenome dominance. A male-linked and pericentromerically located region on the putative sex chromosome was identified in a recessive subgenome and was found to harbour three transcription factors involved in flower and pollen development, including a likely neofunctionalized LEAFY duplicate. Transcriptomic and syntenic analyses of carnivory-related genes suggested that the paleopolyploidization events seeded genes that subsequently formed tandem clusters in recessive subgenomes with specific expression in the digestive zone of the pitcher, where specialized cells digest prey and absorb derived nutrients. A genome-scale analysis suggested that subgenome dominance likely contributed to evolutionary innovation by permitting recessive subgenomes to diversify functions of novel tissue-specific duplicates. Our results provide insight into how polyploidy can give rise to novel traits in divergent and successful high-ploidy lineages.

Gene duplications facilitate C4-CAM compatibility in common purslane

Common purslane (Portulaca oleracea) integrates both C4 and crassulacean acid metabolism (CAM) photosynthesis pathways and is a promising model plant to explore C4-CAM plasticity. Here, we report a high-quality chromosome-level genome of nicotinamide adenine dinucleotide (NAD)-malic enzyme (ME) subtype common purslane that provides evidence for 2 rounds of whole-genome duplication (WGD) with an ancient WGD (P-β) in the common ancestor to Portulacaceae and Cactaceae around 66.30 million years ago (Mya) and another (Po-α) specific to common purslane lineage around 7.74 Mya. A larger number of gene copies encoding key enzymes/transporters involved in C4 and CAM pathways were detected in common purslane than in related species. Phylogeny, conserved functional site, and collinearity analyses revealed that the Po-α WGD produced the phosphoenolpyruvate carboxylase-encoded gene copies used for photosynthesis in common purslane, while the P-β WGD event produced 2 ancestral genes of functionally differentiated (C4- and CAM-specific) beta carbonic anhydrases involved in the C4 + CAM pathways. Additionally, cis-element enrichment analysis in the promoters showed that CAM-specific genes have recruited both evening and midnight circadian elements as well as the Abscisic acid (ABA)-independent regulatory module mediated by ethylene-response factor cis-elements. Overall, this study provides insights into the origin and evolutionary process of C4 and CAM pathways in common purslane, as well as potential targets for engineering crops by integrating C4 or CAM metabolism.

Phylogenomic Analyses Reveal an Allopolyploid Origin of Core Didymocarpinae (Gesneriaceae) Followed by Rapid Radiation

Allopolyploid plants have long been regarded as possessing genetic advantages under certain circumstances due to the combined effects of their hybrid origins and duplicated genomes. However, the evolutionary consequences of allopolyploidy in lineage diversification remain to be fully understood. Here, we investigate the evolutionary consequences of allopolyploidy using 138 transcriptomic sequences of Gesneriaceae, including 124 newly sequenced, focusing particularly on the largest subtribe Didymocarpinae. We estimated the phylogeny of Gesneriaceae using concatenated and coalescent-based methods based on five different nuclear matrices and 27 plastid genes, focusing on relationships among major clades. To better understand the evolutionary affinities in this family, we applied a range of approaches to characterize the extent and cause of phylogenetic incongruence. We found that extensive conflicts between nuclear and chloroplast genomes and among nuclear genes were caused by both incomplete lineage sorting (ILS) and reticulation, and we found evidence of widespread ancient hybridization and introgression. Using the most highly supported phylogenomic framework, we revealed multiple bursts of gene duplication throughout the evolutionary history of Gesneriaceae. By incorporating molecular dating and analyses of diversification dynamics, our study shows that an ancient allopolyploidization event occurred around the Oligocene-Miocene boundary, which may have driven the rapid radiation of core Didymocarpinae.

Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat

Common buckwheat, Fagopyrum esculentum, is an orphan crop domesticated in southwest China that exhibits heterostylous self-incompatibility. Here we present chromosome-scale assemblies of a self-compatible F. esculentum accession and a self-compatible wild relative, Fagopyrum homotropicum, together with the resequencing of 104 wild and cultivated F. esculentum accessions. Using these genomic data, we report the roles of transposable elements and whole-genome duplications in the evolution of Fagopyrum. In addition, we show that (1) the breakdown of heterostyly occurs through the disruption of a hemizygous gene jointly regulating the style length and female compatibility and (2) southeast Tibet was involved in common buckwheat domestication. Moreover, we obtained mutants conferring the waxy phenotype for the first time in buckwheat. These findings demonstrate the utility of our F. esculentum assembly as a reference genome and promise to accelerate buckwheat research and breeding.

Phylotranscriptomics supports numerous polyploidization events and phylogenetic relationships in Nicotiana

Introduction

Nicotiana L. (Solanaceae) is of great scientific and economic importance, and polyploidization has been pivotal in shaping this genus. Despite many previous studies on the Nicotiana phylogenetic relationship and hybridization, evidence from whole genome data is still lacking.

Methods

In this study, we obtained 995 low-copy genes and plastid transcript fragments from the transcriptome datasets of 26 Nicotiana species, including all sections. We reconstructed the phylogenetic relationship and phylogenetic network of diploid species.

Results

The incongruence among gene trees showed that the formation of N. sylvestris involved incomplete lineage sorting. The nuclear-plastid discordance and nuclear introgression absence indicated that organelle capture from section Trigonophyllae was involved in forming section Petunioides. Furthermore, we analyzed the evolutionary origin of polyploid species and dated the time of hybridization events based on the analysis of PhyloNet, sequence similarity search, and phylogeny of subgenome approaches. Our results highly evidenced the hybrid origins of five polyploid sections, including sections Nicotiana, Repandae, Rusticae, Polydicliae, and Suaveolentes. Notably, we provide novel insights into the hybridization event of section Polydicliae and Suaveolentes. The section Polydicliae formed from a single hybridization event between maternal progenitor N. attenuata and paternal progenitor N. undulata; the N. sylvestris (paternal progenitor) and the N. glauca (maternal progenitor) were involved in the formation of section Suaveolentes.

Discussion

This study represents the first exploration of Nicotiana polyploidization events and phylogenetic relationships using the high-throughput RNA-seq approach. It will provide guidance for further studies in molecular systematics, population genetics, and ecological adaption studies in Nicotiana and other related species.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

Inference of Ancient Polyploidy from Genomic Data

Whole-genome sequence data have revealed that numerous eukaryotic organisms derive from distant polyploid ancestors, even when these same organisms are genetically and karyotypically diploid. Such ancient whole-genome duplications (WGDs) have been important for long-term genome evolution and are often speculatively associated with important evolutionary events such as key innovations, adaptive radiations, or survival after mass extinctions. Clearly, reliable methods for unveiling ancient WGDs are key toward furthering understanding of the long-term evolutionary significance of polyploidy. In this chapter, we describe a set of basic established comparative genomics approaches for the inference of ancient WGDs from genomic data based on empirical age distributions and collinearity analyses, explain the principles on which they are based, and illustrate a basic workflow using the software "wgd," geared toward a typical exploratory analysis of a newly obtained genome sequence.

The nearly complete assembly of the Cercis chinensis genome and Fabaceae phylogenomic studies provide insights into new gene evolution

Fabaceae is a large family of angiosperms with high biodiversity that contains a variety of economically important crops and model plants for the study of biological nitrogen fixation. Polyploidization events have been extensively studied in some Fabaceae plants, but the occurrence of new genes is still concealed, owing to a lack of genomic information on certain species of the basal clade of Fabaceae. Cercis chinensis (Cercidoideae) is one such species; it diverged earliest from Fabaceae and is essential for phylogenomic studies and new gene predictions in Fabaceae. To facilitate genomic studies on Fabaceae, we performed genome sequencing of C. chinensis and obtained a 352.84 Mb genome, which was further assembled into seven pseudochromosomes with 30 612 predicted protein-coding genes. Compared with other legume genomes, that of C. chinensis exhibits no lineage-specific polyploidization event. Further phylogenomic analyses of 22 legumes and 11 other angiosperms revealed that many gene families are lineage specific before and after the diversification of Fabaceae. Among them, dozens of genes are candidates for new genes that have evolved from intergenic regions and are thus regarded as de novo-originated genes. They differ significantly from established genes in coding sequence length, exon number, guanine-cytosine content, and expression patterns among tissues. Functional analysis revealed that many new genes are related to asparagine metabolism. This study represents an important advance in understanding the evolutionary pattern of new genes in legumes and provides a valuable resource for plant phylogenomic studies.

Phylotranscriptomics of liverworts: revisiting the backbone phylogeny and ancestral gene duplications

BACKGROUND AND AIMS

With some 7300 extant species, liverworts (Marchantiophyta) represent one of the major land plant lineages. The backbone relationships, such as the phylogenetic position of Ptilidiales, and the occurrence and timing of whole-genome duplications, are still contentious.

METHODS

Based on analyses of the newly generated transcriptome data for 38 liverworts and complemented with those publicly available, we reconstructed the evolutionary history of liverworts and inferred gene duplication events along the 55 taxon liverwort species tree.

KEY RESULTS

Our phylogenomic study provided an ordinal-level liverwort nuclear phylogeny and identified extensive gene tree conflicts and cyto-nuclear incongruences. Gene duplication analyses based on integrated phylogenomics and Ks distributions indicated no evidence of whole-genome duplication events along the backbone phylogeny of liverworts.

CONCLUSIONS

With a broadened sampling of liverwort transcriptomes, we re-evaluated the backbone phylogeny of liverworts, and provided evidence for ancient hybridizations followed by incomplete lineage sorting that shaped the deep evolutionary history of liverworts. The lack of whole-genome duplication during the deep evolution of liverworts indicates that liverworts might represent one of the few major embryophyte lineages whose evolution was not driven by whole-genome duplications.

Genomic analyses provide insights into the evolution and salinity adaptation of halophyte Tamarix chinensis

BACKGROUND

The woody halophyte Tamarix chinensis is a pioneer tree species in the coastal wetland ecosystem of northern China, exhibiting high resistance to salt stress. However, the genetic information underlying salt tolerance in T. chinensis remains to be seen. Here we present a genomic investigation of T. chinensis to elucidate the underlying mechanism of its high resistance to salinity.

RESULTS

Using a combination of PacBio and high-throughput chromosome conformation capture data, a chromosome-level T. chinensis genome was assembled with a size of 1.32 Gb and scaffold N50 of 110.03 Mb. Genome evolution analyses revealed that T. chinensis significantly expanded families of HAT and LIMYB genes. Whole-genome and tandem duplications contributed to the expansion of genes associated with the salinity adaptation of T. chinensis. Transcriptome analyses were performed on root and shoot tissues during salt stress and recovery, and several hub genes responding to salt stress were identified. WRKY33/40, MPK3/4, and XBAT31 were critical in responding to salt stress during early exposure, while WRKY40, ZAT10, AHK4, IRX9, and CESA4/8 were involved in responding to salt stress during late stress and recovery. In addition, PER7/27/57/73 encoding class III peroxidase and MCM3/4/5/7 encoding DNA replication licensing factor maintained up/downregulation during salt stress and recovery stages.

CONCLUSIONS

The results presented here reveal the genetic mechanisms underlying salt adaptation in T. chinensis, thus providing important genomic resources for evolutionary studies on tamarisk and plant salt tolerance genetic improvement.

Evolution of complex genome architecture in gymnosperms

Gymnosperms represent an ancient lineage that diverged from early spermatophytes during the Devonian. The long fossil records and low diversity in living species prove their complex evolutionary history, which included ancient radiations and massive extinctions. Due to their ultra-large genome size, the whole-genome assembly of gymnosperms has only generated in the past 10 years and is now being further expanded into more taxonomic representations. Here, we provide an overview of the publicly available gymnosperm genome resources and discuss their assembly quality and recent findings in large genome architectures. In particular, we describe the genomic features most related to changes affecting the whole genome. We also highlight new realizations relative to repetitive sequence dynamics, paleopolyploidy, and long introns. Based on the results of relevant genomic studies of gymnosperms, we suggest additional efforts should be made toward exploring the genomes of medium-sized (5–15 gigabases) species. Lastly, more comparative analyses among high-quality assemblies are needed to understand the genomic shifts and the early species diversification of seed plants.

Species Tree Estimation and the Impact of Gene Loss Following Whole-Genome Duplication

Whole-genome duplication (WGD) occurs broadly and repeatedly across the history of eukaryotes and is recognized as a prominent evolutionary force, especially in plants. Immediately following WGD, most genes are present in two copies as paralogs. Due to this redundancy, one copy of a paralog pair commonly undergoes pseudogenization and is eventually lost. When speciation occurs shortly after WGD; however, differential loss of paralogs may lead to spurious phylogenetic inference resulting from the inclusion of pseudoorthologs–paralogous genes mistakenly identified as orthologs because they are present in single copies within each sampled species. The influence and impact of including pseudoorthologs versus true orthologs as a result of gene extinction (or incomplete laboratory sampling) are only recently gaining empirical attention in the phylogenomics community. Moreover, few studies have yet to investigate this phenomenon in an explicit coalescent framework. Here, using mathematical models, numerous simulated data sets, and two newly assembled empirical data sets, we assess the effect of pseudoorthologs on species tree estimation under varying degrees of incomplete lineage sorting (ILS) and differential gene loss scenarios following WGD. When gene loss occurs along the terminal branches of the species tree, alignment-based (BPP) and gene-tree-based (ASTRAL, MP-EST, and STAR) coalescent methods are adversely affected as the degree of ILS increases. This can be greatly improved by sampling a sufficiently large number of genes. Under the same circumstances, however, concatenation methods consistently estimate incorrect species trees as the number of genes increases. Additionally, pseudoorthologs can greatly mislead species tree inference when gene loss occurs along the internal branches of the species tree. Here, both coalescent and concatenation methods yield inconsistent results. These results underscore the importance of understanding the influence of pseudoorthologs in the phylogenomics era. [Coalescent method; concatenation method; incomplete lineage sorting; pseudoorthologs; single-copy gene; whole-genome duplication.]

The genome of the recretohalophyte Limonium bicolor provides insights into salt gland development and salinity adaptation during terrestrial evolution

Halophytes have evolved specialized strategies to cope with high salinity. The extreme halophyte sea lavender (Limonium bicolor) lacks trichomes but possesses salt glands on its epidermis that can excrete harmful ions, such as sodium, to avoid salt damage. Here, we report a high-quality, 2.92-Gb, chromosome-scale L. bicolor genome assembly based on a combination of Illumina short reads, single-molecule, real-time long reads, chromosome conformation capture (Hi-C) data, and Bionano genome maps, greatly enriching the genomic information on recretohalophytes with multicellular salt glands. Although the L. bicolor genome contains genes that show similarity to trichome fate genes from Arabidopsis thaliana, it lacks homologs of the decision fate genes GLABRA3, ENHANCER OF GLABRA3, GLABRA2, TRANSPARENT TESTA GLABRA2, and SIAMESE, providing a molecular explanation for the absence of trichomes in this species. We identified key genes (LbHLH and LbTTG1) controlling salt gland development among classical trichome homologous genes and confirmed their roles by showing that their mutations markedly disrupted salt gland initiation, salt secretion, and salt tolerance, thus offering genetic support for the long-standing hypothesis that salt glands and trichomes may share a common origin. In addition, a whole-genome duplication event occurred in the L. bicolor genome after its divergence from Tartary buckwheat and may have contributed to its adaptation to high salinity. The L. bicolor genome resource and genetic evidence reported in this study provide profound insights into plant salt tolerance mechanisms that may facilitate the engineering of salt-tolerant crops.

Phylogenomic analyses in Phrymaceae reveal extensive gene tree discordance in relationships among major clades

PREMISE

Phylogenomic datasets using genomes and transcriptomes provide rich opportunities beyond resolving bifurcating phylogenetic relationships. Monkeyflower (Phrymaceae) is a model system for evolutionary ecology. However, it lacks a well-supported phylogeny as a basis for a stable taxonomy and for macroevolutionary comparisons.

METHODS

We sampled 24 genomes and transcriptomes in Phrymaceae and closely related families, including eight newly sequenced transcriptomes. We reconstructed the phylogeny using IQ-TREE and ASTRAL, evaluated gene tree discordance using PhyParts, Quartet Sampling, and a cloudogram, and carried out reticulation analyses using PhyloNet and HyDe. We searched for whole genome duplication (WGD) events using chromosome numbers, synonymous distances, and gene duplication events as evidence.

RESULTS

Most gene trees support the monophyly of Phrymaceae and each of its tribes. Most gene trees also support tribe Mimuleae being sister to Phrymeae + Diplaceae + Leucocarpeae, with extensive gene tree discordance among the latter three. Despite the discordance, the monophyly of Mimulus s.l. is rejected, and no individual reticulation event among the Phrymaceae tribes is well-supported. Reticulation likely occurred among Erythranthe bicolor and closely related species. No ancient WGD was detected in Phrymaceae. Instead, small-scale duplications are among potential drivers of macroevolutionary diversification of Phrymaceae.

CONCLUSIONS

We show that analysis of reticulate evolution is sensitive to taxon sampling and methods used. We also demonstrate that phylogenomic datasets using genomes and transcriptomes present rich opportunities to investigate gene family evolution and genome duplication events involved in lineage diversification and adaptation.

Phylogenomic Analyses of Alismatales Shed Light into Adaptations to Aquatic Environments

Land plants first evolved from freshwater algae, and flowering plants returned to water as early as the Cretaceous and multiple times subsequently. Alismatales is the largest clade of aquatic angiosperms including all marine angiosperms, as well as terrestrial plants. We used Alismatales to explore plant adaptations to aquatic environments by analyzing a data set that included 95 samples (89 Alismatales species) covering four genomes and 91 transcriptomes (59 generated in this study). To provide a basis for investigating adaptations, we assessed phylogenetic conflict and whole-genome duplication (WGD) events in Alismatales. We recovered a relationship for the three main clades in Alismatales as (Tofieldiaceae, Araceae) + core Alismatids. We also found phylogenetic conflict among the three main clades that was best explained by incomplete lineage sorting and introgression. Overall, we identified 18 putative WGD events across Alismatales. One of them occurred at the most recent common ancestor of core Alismatids, and three occurred at seagrass lineages. We also found that lineage and life-form were both important for different evolutionary patterns for the genes related to freshwater and marine adaptation. For example, several light- or ethylene-related genes were lost in the seagrass Zosteraceae, but are present in other seagrasses and freshwater species. Stomata-related genes were lost in both submersed freshwater species and seagrasses. Nicotianamine synthase genes, which are important in iron intake, expanded in both submersed freshwater species and seagrasses. Our results advance the understanding of the adaptation to aquatic environments and WGDs using phylogenomics.

Global Patterns of Subgenome Evolution in Organelle-Targeted Genes of Six Allotetraploid Angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary "mismatches" between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

A well-supported nuclear phylogeny of Poaceae and implications for the evolution of C4 photosynthesis

Poaceae (the grasses) includes rice, maize, wheat, and other crops, and is the most economically important angiosperm family. Poaceae is also one of the largest plant families, consisting of over 11 000 species with a global distribution that contributes to diverse ecosystems. Poaceae species are classified into 12 subfamilies, with generally strong phylogenetic support for their monophyly. However, many relationships within subfamilies, among tribes and/or subtribes, remain uncertain. To better resolve the Poaceae phylogeny, we generated 342 transcriptomic and seven genomic datasets; these were combined with other genomic and transcriptomic datasets to provide sequences for 357 Poaceae species in 231 genera, representing 45 tribes and all 12 subfamilies. Over 1200 low-copy nuclear genes were retrieved from these datasets, with several subsets obtained using additional criteria, and used for coalescent analyses to reconstruct a Poaceae phylogeny. Our results strongly support the monophyly of 11 subfamilies; however, the subfamily Puelioideae was separated into two non-sister clades, one for each of the two previously defined tribes, supporting a hypothesis that places each tribe in a separate subfamily. Molecular clock analyses estimated the crown age of Poaceae to be ∼101 million years old. Ancestral character reconstruction of C₃/C₄ photosynthesis supports the hypothesis of multiple independent origins of C₄ photosynthesis. These origins are further supported by phylogenetic analysis of the ppc gene family that encodes the phosphoenolpyruvate carboxylase, which suggests that members of three paralogous subclades (ppc-aL1a, ppc-aL1b, and ppc-B2) were recruited as functional C₄ppc genes. This study provides valuable resources and a robust phylogenetic framework for evolutionary analyses of the grass family.

Evolutionary Genetics of Cacti: Research Biases, Advances and Prospects

Here, we present a review of the studies of evolutionary genetics (phylogenetics, population genetics, and phylogeography) using genetic data as well as genome scale assemblies in Cactaceae (Caryophyllales, Angiosperms), a major lineage of succulent plants with astonishing diversity on the American continent. To this end, we performed a literature survey (1992-2021) to obtain detailed information regarding key aspects of studies investigating cactus evolution. Specifically, we summarize the advances in the following aspects: molecular markers, species delimitation, phylogenetics, hybridization, biogeography, and genome assemblies. In brief, we observed substantial growth in the studies conducted with molecular markers in the past two decades. However, we found biases in taxonomic/geographic sampling and the use of traditional markers and statistical approaches. We discuss some methodological and social challenges for engaging the cactus community in genomic research. We also stressed the importance of integrative approaches, coalescent methods, and international collaboration to advance the understanding of cactus evolution.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

Analysis of Paralogs in Target Enrichment Data Pinpoints Multiple Ancient Polyploidy Events in Alchemilla s.l. (Rosaceae)

Target enrichment is becoming increasingly popular for phylogenomic studies. Although baits for enrichment are typically designed to target single-copy genes, paralogs are often recovered with increased sequencing depth, sometimes from a significant proportion of loci, especially in groups experiencing whole-genome duplication (WGD) events. Common approaches for processing paralogs in target enrichment data sets include random selection, manual pruning, and mainly, the removal of entire genes that show any evidence of paralogy. These approaches are prone to errors in orthology inference or removing large numbers of genes. By removing entire genes, valuable information that could be used to detect and place WGD events is discarded. Here, we used an automated approach for orthology inference in a target enrichment data set of 68 species of Alchemilla s.l. (Rosaceae), a widely distributed clade of plants primarily from temperate climate regions. Previous molecular phylogenetic studies and chromosome numbers both suggested ancient WGDs in the group. However, both the phylogenetic location and putative parental lineages of these WGD events remain unknown. By taking paralogs into consideration and inferring orthologs from target enrichment data, we identified four nodes in the backbone of Alchemilla s.l. with an elevated proportion of gene duplication. Furthermore, using a gene-tree reconciliation approach, we established the autopolyploid origin of the entire Alchemilla s.l. and the nested allopolyploid origin of four major clades within the group. Here, we showed the utility of automated tree-based orthology inference methods, previously designed for genomic or transcriptomic data sets, to study complex scenarios of polyploidy and reticulate evolution from target enrichment data sets.[Alchemilla; allopolyploidy; autopolyploidy; gene tree discordance; orthology inference; paralogs; Rosaceae; target enrichment; whole genome duplication.].

Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.

On the potential of Angiosperms353 for population genomic studies

PREMISE

The successful application of universal targeted sequencing markers, such as those developed for the Angiosperms353 probe set, within populations could reduce or eliminate the need for specific marker development, while retaining the benefits of full-gene sequences in population-level analyses. However, whether the Angiosperms353 markers provide sufficient variation within species to calculate demographic parameters is untested.

METHODS

Using herbarium specimens from a 50-year-old floristic survey in Texas, we sequenced 95 samples from 24 species using the Angiosperms353 probe set. Our data workflow calls variants within species and prepares data for population genetic analysis using standard metrics. In our case study, gene recovery was affected by genomic library concentration only at low concentrations and displayed limited phylogenetic bias.

RESULTS

We identified over 1000 segregating variants with zero missing data for 92% of species and demonstrate that Angiosperms353 markers contain sufficient variation to estimate pairwise nucleotide diversity (π)-typically between 0.002 and 0.010, with most variation found in flanking non-coding regions. In a subset of variants that were filtered to reduce linkage, we uncovered high heterozygosity in many species, suggesting that denser sampling within species should permit estimation of gene flow and population dynamics.

DISCUSSION

Angiosperms353 should benefit conservation genetic studies by providing universal repeatable markers, low missing data, and haplotype information, while permitting inclusion of decades-old herbarium specimens.

On the Origin of Carnivory: Molecular Physiology and Evolution of Plants on an Animal Diet

Charles Darwin recognized that carnivorous plants thrive in nutrient-poor soil by capturing animals. Although the concept of botanical carnivory has been known for nearly 150 years, its molecular mechanisms and evolutionary origins have not been well understood until recently. In the last decade, technical advances have fueled the genome and transcriptome sequencings of active and passive hunters, leading to a better understanding of the traits associated with the carnivorous syndrome, from trap leaf development and prey digestion to nutrient absorption, exemplified, for example, by the Venus flytrap (Dionaea muscipula), pitcher plant (Cephalotus follicularis), and bladderwort (Utricularia gibba). The repurposing of defense-related genes is an important trend in the evolution of plant carnivory. In this review, using the Venus flytrap as a representative of the carnivorous plants, we summarize the molecular mechanisms underlying their ability to attract, trap, and digest prey and discuss the origins of plant carnivory in relation to their genomic evolution.

Nuclear phylotranscriptomics and phylogenomics support numerous polyploidization events and hypotheses for the evolution of rhizobial nitrogen-fixing symbiosis in Fabaceae

Fabaceae are the third largest angiosperm family, with 765 genera and ∼19 500 species. They are important both economically and ecologically, and global Fabaceae crops are intensively studied in part for their nitrogen-fixing ability. However, resolution of the intrasubfamilial Fabaceae phylogeny and divergence times has remained elusive, precluding a reconstruction of the evolutionary history of symbiotic nitrogen fixation in Fabaceae. Here, we report a highly resolved phylogeny using >1500 nuclear genes from newly sequenced transcriptomes and genomes of 391 species, along with other datasets, for a total of 463 legumes spanning all 6 subfamilies and 333 of 765 genera. The subfamilies are maximally supported as monophyletic. The clade comprising subfamilies Cercidoideae and Detarioideae is sister to the remaining legumes, and Duparquetioideae and Dialioideae are successive sisters to the clade of Papilionoideae and Caesalpinioideae. Molecular clock estimation revealed an early radiation of subfamilies near the K/Pg boundary, marked by mass extinction, and subsequent divergence of most tribe-level clades within ∼15 million years. Phylogenomic analyses of thousands of gene families support 28 proposed putative whole-genome duplication/whole-genome triplication events across Fabaceae, including those at the ancestors of Fabaceae and five of the subfamilies, and further analyses supported the Fabaceae ancestral polyploidy. The evolution of rhizobial nitrogen-fixing nodulation in Fabaceae was probed by ancestral character reconstruction and phylogenetic analyses of related gene families and the results support the hypotheses of one or two switch(es) to rhizobial nodulation followed by multiple losses. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Fabaceae.

Phylogenomic conflict coincides with rapid morphological innovation

Evolutionary biologists have long been fascinated with the episodes of rapid phenotypic innovation that underlie the emergence of major lineages. Although our understanding of the environmental and ecological contexts of such episodes has steadily increased, it has remained unclear how population processes contribute to emergent macroevolutionary patterns. One insight gleaned from phylogenomics is that gene-tree conflict, frequently caused by population-level processes, is often rampant during the origin of major lineages. With the understanding that phylogenomic conflict is often driven by complex population processes, we hypothesized that there may be a direct correspondence between instances of high conflict and elevated rates of phenotypic innovation if both patterns result from the same processes. We evaluated this hypothesis in six clades spanning vertebrates and plants. We found that the most conflict-rich regions of these six clades also tended to experience the highest rates of phenotypic innovation, suggesting that population processes shaping both phenotypic and genomic evolution may leave signatures at deep timescales. Closer examination of the biological significance of phylogenomic conflict may yield improved connections between micro- and macroevolution and increase our understanding of the processes that shape the origin of major lineages across the Tree of Life.

Phylotranscriptomic analysis of Dillenia indica L. (Dilleniales, Dilleniaceae) and its systematics implication

The recent massive development in the next-generation sequencing platforms and bioinformatics tools including cloud based computing have proven extremely useful in understanding the deeper-level phylogenetic relationships of angiosperms. The present phylotranscriptomic analyses address the poorly known evolutionary relationships of the order Dilleniales to order of the other angiosperms using the minimum evolution method. The analyses revealed the nesting of the representative taxon of Dilleniales in the MPT but distinct from the representative of the order Santalales, Caryophyllales, Asterales, Cornales, Ericales, Lamiales, Saxifragales, Fabales, Malvales, Vitales and Berberidopsidales.

Secretory structures in plants: Lessons from the Plumbaginaceae on their origin, evolution and roles in stress tolerance

The Plumbaginaceae (non-core Caryophyllales) is a family well known for species adapted to a wide range of arid and saline habitats. Of its salt-tolerant species, at least 45 are in the genus Limonium; two in each of Aegialitis, Limoniastrum and Myriolimon, and one each in Psylliostachys, Armeria, Ceratostigma, Goniolimon and Plumbago. All the halophytic members of the family have salt glands and salt glands are also common in the closely related Tamaricaceae and Frankeniaceae. The halophytic species of the three families can secrete a range of ions (Na⁺ , K⁺ , Ca²⁺ , Mg²⁺ , Cl^- , HCO₃^- , SO₄^2- ) and other elements (As, Cd, Cr, Cu, Fe, Mn, Ni, Pb and Zn). Salt glands are, however, absent in salt-tolerant members of the sister family Polygonaceae. We describe the structure of the salt glands in the three families and consider whether glands might have arisen as a means to avoid the toxicity of Na⁺ and/or Cl^- or to regulate Ca²⁺ concentrations with the leaves. We conclude that the establishment of lineages with salt glands took place after the split between the Polygonaceae and its sister group the Plumbaginaceae.

The Origin of the Legumes is a Complex Paleopolyploid Phylogenomic Tangle Closely Associated with the Cretaceous-Paleogene (K-Pg) Mass Extinction Event

The consequences of the Cretaceous–Paleogene (K–Pg) boundary (KPB) mass extinction for the evolution of plant diversity remain poorly understood, even though evolutionary turnover of plant lineages at the KPB is central to understanding assembly of the Cenozoic biota. The apparent concentration of whole genome duplication (WGD) events around the KPB may have played a role in survival and subsequent diversification of plant lineages. To gain new insights into the origins of Cenozoic biodiversity, we examine the origin and early evolution of the globally diverse legume family (Leguminosae or Fabaceae). Legumes are ecologically (co-)dominant across many vegetation types, and the fossil record suggests that they rose to such prominence after the KPB in parallel with several well-studied animal clades including Placentalia and Neoaves. Furthermore, multiple WGD events are hypothesized to have occurred early in legume evolution. Using a recently inferred phylogenomic framework, we investigate the placement of WGDs during early legume evolution using gene tree reconciliation methods, gene count data and phylogenetic supernetwork reconstruction. Using 20 fossil calibrations we estimate a revised timeline of legume evolution based on 36 nuclear genes selected as informative and evolving in an approximately clock-like fashion. To establish the timing of WGDs we also date duplication nodes in gene trees. Results suggest either a pan-legume WGD event on the stem lineage of the family, or an allopolyploid event involving (some of) the earliest lineages within the crown group, with additional nested WGDs subtending subfamilies Papilionoideae and Detarioideae. Gene tree reconciliation methods that do not account for allopolyploidy may be misleading in inferring an earlier WGD event at the time of divergence of the two parental lineages of the polyploid, suggesting that the allopolyploid scenario is more likely. We show that the crown age of the legumes dates to the Maastrichtian or early Paleocene and that, apart from the Detarioideae WGD, paleopolyploidy occurred close to the KPB. We conclude that the early evolution of the legumes followed a complex history, in which multiple auto- and/or allopolyploidy events coincided with rapid diversification and in association with the mass extinction event at the KPB, ultimately underpinning the evolutionary success of the Leguminosae in the Cenozoic. [Allopolyploidy; Cretaceous–Paleogene (K–Pg) boundary; Fabaceae, Leguminosae; paleopolyploidy; phylogenomics; whole genome duplication events]

Phylotranscriptomics in Cucurbitaceae Reveal Multiple Whole-Genome Duplications and Key Morphological and Molecular Innovations

The ability of climbing plants to grow upward along others to reach the canopy for photosynthesis is hypothesized as a key innovation in flowering plants. Most members of the Cucurbitaceae, a family containing ∼1000 species and many important crops, are climbers and have characteristic tendrils and pepo fruits. Here, we present 127 newly sequenced transcriptomes and genomes along with other datasets for a total of 136 cucurbits representing all tribes to establish a robust Cucurbitaceae phylogeny containing eight highly resolved major clades. We analyzed whole-genome duplication, diversification dynamics, and ancestral morphologies, and found that after early genome duplication event(s), a burst of diversification and morphological innovations in flower, fruit, and root characters occurred under the climate optimum in the Early Eocene. Species radiation during the Mid-Eocene Climatic Optimum also coincided with several morphological changes shared by 80% of cucurbits. We found that the cucurbit-specific tendril identity gene TEN originated from a paleo-polyploidization event at the origin of the family. Our results support the hypothesis that cucurbit diversifications were probably driven by increased genetic diversity following polyploidizations and by trait morphological innovations under paleo-climate upheavals. Our study provides a phylogenetic framework and new insights into morphological and genomic changes underlying the adaptive evolution of Cucurbitaceae.

Genomes of the Venus Flytrap and Close Relatives Unveil the Roots of Plant Carnivory

Most plants grow and develop by taking up nutrients from the soil while continuously under threat from foraging animals. Carnivorous plants have turned the tables by capturing and consuming nutrient-rich animal prey, enabling them to thrive in nutrient-poor soil. To better understand the evolution of botanical carnivory, we compared the draft genome of the Venus flytrap (Dionaea muscipula) with that of its aquatic sister, the waterwheel plant Aldrovanda vesiculosa, and the sundew Drosera spatulata. We identified an early whole-genome duplication in the family as source for carnivory-associated genes. Recruitment of genes to the trap from the root especially was a major mechanism in the evolution of carnivory, supported by family-specific duplications. Still, these genomes belong to the gene poorest land plants sequenced thus far, suggesting reduction of selective pressure on different processes, including non-carnivorous nutrient acquisition. Our results show how non-carnivorous plants evolved into the most skillful green hunters on the planet.

Gene and genome duplications in the evolution of chemodiversity: perspectives from studies of Lamiaceae

Plants are reservoirs of extreme chemical diversity, yet biosynthetic pathways remain underexplored in the majority of taxa. Access to improved, inexpensive genomic and computational technologies has recently enhanced our understanding of plant specialized metabolism at the biochemical and evolutionary levels including the elucidation of pathways leading to key metabolites. Furthermore, these approaches have provided insights into the mechanisms of chemical evolution, including neofunctionalization and subfunctionalization, structural variation, and modulation of gene expression. The broader utilization of genomic tools across the plant tree of life, and an expansion of genomic resources from multiple accessions within species or populations, will improve our overall understanding of chemodiversity. These data and knowledge will also lead to greater insight into the selective pressures contributing to and maintaining this diversity, which in turn will enable the development of more accurate predictive models of specialized metabolism in plants.

Phylogenomic synteny network analyses reveal ancestral transpositions of auxin response factor genes in plants

BACKGROUND

Auxin response factors (ARFs) have long been a research focus and represent a class of key regulators of plant growth and development. Integrated phylogenomic synteny network analyses were able to provide novel insights into the evolution of the ARF gene family.

RESULTS

Here, more than 3500 ARFs collected from plant genomes and transcriptomes covering major streptophyte lineages were used to reconstruct the broad-scale family phylogeny, where the early origin and diversification of ARF in charophytes was delineated. Based on the family phylogeny, we proposed a unified six-group classification system for angiosperm ARFs. Phylogenomic synteny network analyses revealed the deeply conserved genomic syntenies within each of the six ARF groups and the interlocking syntenic relationships connecting distinct groups. Recurrent duplication events, such as those that occurred in seed plants, angiosperms, core eudicots and grasses contributed to the expansion of ARF genes which facilitated functional diversification. Ancestral transposition activities in important plant families, including crucifers, legumes and grasses, were unveiled by synteny network analyses. Ancestral gene duplications along with transpositions have profound evolutionary significance which may have accelerated the functional diversification process of paralogues.

CONCLUSIONS

The broad-scale family phylogeny in combination with the state-of-art phylogenomic synteny network analyses not only allowed us to infer the evolutionary trajectory of ARF genes across distinct plant lineages, but also facilitated to generate a more robust classification regime for this transcription factor family. Our study provides insights into the evolution of ARFs which will enhance our current understanding of this important transcription factor family.

Nuclear phylogenomic analyses of asterids conflict with plastome trees and support novel relationships among major lineages

PREMISE

Discordance between nuclear and organellar phylogenies (cytonuclear discordance) is a well-documented phenomenon at shallow evolutionary levels but has been poorly investigated at deep levels of plant phylogeny. Determining the extent of cytonuclear discordance across major plant lineages is essential not only for elucidating evolutionary processes, but also for evaluating the currently used framework of plant phylogeny, which is largely based on the plastid genome.

METHODS

We present a phylogenomic examination of a major angiosperm clade (Asteridae) based on sequence data from the nuclear, plastid, and mitochondrial genomes as a means of evaluating currently accepted relationships inferred from the plastome and exploring potential sources of genomic conflict in this group.

RESULTS

We recovered at least five instances of well-supported cytonuclear discordance concerning the placements of major asterid lineages (i.e., Ericales, Oncothecaceae, Aquifoliales, Cassinopsis, and Icacinaceae). We attribute this conflict to a combination of incomplete lineage sorting and hybridization, the latter supported in part by previously inferred whole-genome duplications.

CONCLUSIONS

Our results challenge several long-standing hypotheses of asterid relationships and have implications for morphological character evolution and for the importance of ancient whole-genome duplications in early asterid evolution. These findings also highlight the value of reevaluating broad-scale angiosperm and green-plant phylogeny with nuclear genomic data.

A consensus phylogenomic approach highlights paleopolyploid and rapid radiation in the history of Ericales

PREMISE

Large genomic data sets offer the promise of resolving historically recalcitrant species relationships. However, different methodologies can yield conflicting results, especially when clades have experienced ancient, rapid diversification. Here, we analyzed the ancient radiation of Ericales and explored sources of uncertainty related to species tree inference, conflicting gene tree signal, and the inferred placement of gene and genome duplications.

METHODS

We used a hierarchical clustering approach, with tree-based homology and orthology detection, to generate six filtered phylogenomic matrices consisting of data from 97 transcriptomes and genomes. Support for species relationships was inferred from multiple lines of evidence including shared gene duplications, gene tree conflict, gene-wise edge-based analyses, concatenation, and coalescent-based methods, and is summarized in a consensus framework.

RESULTS

Our consensus approach supported a topology largely concordant with previous studies, but suggests that the data are not capable of resolving several ancient relationships because of lack of informative characters, sensitivity to methodology, and extensive gene tree conflict correlated with paleopolyploidy. We found evidence of a whole-genome duplication before the radiation of all or most ericalean families, and demonstrate that tree topology and heterogeneous evolutionary rates affect the inferred placement of genome duplications.

CONCLUSIONS

We provide several hypotheses regarding the history of Ericales, and confidently resolve most nodes, but demonstrate that a series of ancient divergences are unresolvable with these data. Whether paleopolyploidy is a major source of the observed phylogenetic conflict warrants further investigation.

Recurrent genome duplication events likely contributed to both the ancient and recent rise of ferns

Ferns, the second largest group of vascular plants, originated ~400 million years ago (Mya). They became dominant in the ancient Earth landscape before the angiosperms and are still important in current ecosystems. Many ferns have exceptionally high chromosome numbers, possibly resulting from whole-genome duplications (WGDs). However, WGDs have not been investigated molecularly across fern diversity. Here we detected and dated fern WGDs using a phylogenomic approach and by calculating synonymous substitution rates (Ks). We also investigated a possible correlation between proposed WGDs and shifts in species diversification rates. We identified 19 WGDs: three ancient events along the fern phylogenetic backbone that are shared by 66%-97% of extant ferns, with additional lineage-specific WGDs for eight orders, providing strong evidence for recurring genome duplications across fern evolutionary history. We also observed similar Ks peak values for more than half of these WGDs, with multiple WGDs occurring close to the Cretaceous (~145-66 Mya). Despite the repeated WGD events, the biodiversity of ferns declined during the Cretaceous, implying that other factors probably contributed to the floristic turnover from ferns to angiosperms. This study provides molecular evidence for recurring WGDs in ferns and offers important clues to the genomic evolutionary history of ferns.

Has the Polyploid Wave Ebbed?

There was a wave of whole genome duplications (WGD) during and subsequent to the K-Pg interface, which was followed by an increase in the proportion of species that were polyploid. I consider why this wave of polyploid speciation has continued to rise through the divergent evolution of polyploid lineages, and through rounds of homoploid and heteroploid chromosomal change. I also consider why the polyploid speciation wave is likely to rise in the next millennium. I propose that the speed of polyploid genesis through ploidal increase and through diversification among polyploids likely will be greater than the speed of diploid speciation. The increase in polyploid diversity is expected to lag well behind episodes of WGD, owing to the very long period required for species diversification either by lineage splitting or additional rounds of polyploidy, in addition to the long period of genomic adjustment to higher ploidal levels in neopolyploids.

Quartet-Based Computations of Internode Certainty Provide Robust Measures of Phylogenetic Incongruence

Incongruence, or topological conflict, is prevalent in genome-scale data sets. Internode certainty (IC) and related measures were recently introduced to explicitly quantify the level of incongruence of a given internal branch among a set of phylogenetic trees and complement regular branch support measures (e.g., bootstrap, posterior probability) that instead assess the statistical confidence of inference. Since most phylogenomic studies contain data partitions (e.g., genes) with missing taxa and IC scores stem from the frequencies of bipartitions (or splits) on a set of trees, IC score calculation typically requires adjusting the frequencies of bipartitions from these partial gene trees. However, when the proportion of missing taxa is high, the scores yielded by current approaches that adjust bipartition frequencies in partial gene trees differ substantially from each other and tend to be overestimates. To overcome these issues, we developed three new IC measures based on the frequencies of quartets, which naturally apply to both complete and partial trees. Comparison of our new quartet-based measures to previous bipartition-based measures on simulated data shows that: (1) on complete data sets, both quartet-based and bipartition-based measures yield very similar IC scores; (2) IC scores of quartet-based measures on a given data set with and without missing taxa are more similar than the scores of bipartition-based measures; and (3) quartet-based measures are more robust to the absence of phylogenetic signal and errors in phylogenetic inference than bipartition-based measures. Additionally, the analysis of an empirical mammalian phylogenomic data set using our quartet-based measures reveals the presence of substantial levels of incongruence for numerous internal branches. An efficient open-source implementation of these quartet-based measures is freely available in the program QuartetScores (https://github.com/lutteropp/QuartetScores).

Inference of Ancient Whole-Genome Duplications and the Evolution of Gene Duplication and Loss Rates

Gene tree-species tree reconciliation methods have been employed for studying ancient whole-genome duplication (WGD) events across the eukaryotic tree of life. Most approaches have relied on using maximum likelihood trees and the maximum parsimony reconciliation thereof to count duplication events on specific branches of interest in a reference species tree. Such approaches do not account for uncertainty in the gene tree and reconciliation, or do so only heuristically. The effects of these simplifications on the inference of ancient WGDs are unclear. In particular, the effects of variation in gene duplication and loss rates across the species tree have not been considered. Here, we developed a full probabilistic approach for phylogenomic reconciliation-based WGD inference, accounting for both gene tree and reconciliation uncertainty using a method based on the principle of amalgamated likelihood estimation. The model and methods are implemented in a maximum likelihood and Bayesian setting and account for variation of duplication and loss rates across the species tree, using methods inspired by phylogenetic divergence time estimation. We applied our newly developed framework to ancient WGDs in land plants and investigated the effects of duplication and loss rate variation on reconciliation and gene count based assessment of these earlier proposed WGDs.

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.