A reappraisal of the phylogenetic placement of the Aquilegia whole-genome duplication

Lineage-Specific CYP80 Expansion and Benzylisoquinoline Alkaloid Diversity in Early-Diverging Eudicots

Menispermaceae species, as early-diverging eudicots, can synthesize valuable benzylisoquinoline alkaloids (BIAs) like bisbenzylisoquinoline alkaloids (bisBIAs) and sinomenines with a wide range of structural diversity. However, the evolutionary mechanisms responsible for their chemo-diversity are not well understood. Here, a chromosome-level genome assembly of Menispermum dauricum is presented and demonstrated the occurrence of two whole genome duplication (WGD) events that are shared by Ranunculales and specific to Menispermum, providing a model for understanding chromosomal evolution in early-diverging eudicots. The biosynthetic pathway for diverse BIAs in M. dauricum is reconstructed by analyzing the transcriptome and metabolome. Additionally, five catalytic enzymes - one norcoclaurine synthase (NCS) and four cytochrome P450 monooxygenases (CYP450s) - from M. dauricum are responsible for the formation of the skeleton, hydroxylated modification, and C-O/C-C phenol coupling of BIAs. Notably, a novel leaf-specific MdCYP80G10 enzyme that catalyzes C2'-C4a phenol coupling of (S)-reticuline into sinoacutine, the enantiomer of morphinan compounds, with predictable stereospecificity is discovered. Moreover, it is found that Menispermum-specific CYP80 gene expansion, as well as tissue-specific expression, has driven BIA diversity in Menispermaceae as compared to other Ranunculales species. This study sheds light on WGD occurrences in early-diverging eudicots and the evolution of diverse BIA biosynthesis.

Genomic data and ecological niche modeling reveal an unusually slow rate of molecular evolution in the Cretaceous Eupteleaceae

Living fossils are evidence of long-term sustained ecological success. However, whether living fossils have little molecular changes remains poorly known, particularly in plants. Here, we have introduced a novel method that integrates phylogenomic, comparative genomic, and ecological niche modeling analyses to investigate the rate of molecular evolution of Eupteleaceae, a Cretaceous relict angiosperm family endemic to East Asia. We assembled a high-quality chromosome-level nuclear genome, and the chloroplast and mitochondrial genomes of a member of Eupteleaceae (Euptelea pleiosperma). Our results show that Eupteleaceae is most basal in Ranunculales, the earliest-diverging order in eudicots, and shares an ancient whole-genome duplication event with the other Ranunculales. We document that Eupteleaceae has the slowest rate of molecular changes in the observed angiosperms. The unusually low rate of molecular evolution of Eupteleaceae across all three independent inherited genomes and genes within each of the three genomes is in association with its conserved genome architecture, ancestral woody habit, and conserved niche requirements. Our findings reveal the evolution and adaptation of living fossil plants through large-scale environmental change and also provide new insights into early eudicot diversification.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Genome-Wide Identification and Expression Analysis of DWARF53 Gene in Response to GA and SL Related to Plant Height in Banana

Dwarfing is one of the common phenotypic variations in asexually reproduced progeny of banana, and dwarfed banana is not only windproof and anti-fallout but also effective in increasing acreage yield. As a key gene in the strigolactone signalling pathway, DWARF53 (D53) plays an important role in the regulation of the height of plants. In order to gain insight into the function of the banana D53 gene, this study conducted genome-wide identification of banana D53 gene based on the M. acuminata, M. balbisiana and M. itinerans genome database. Analysis of MaD53 gene expression under high temperature, low temperature and osmotic stress based on transcriptome data and RT-qPCR was used to analyse MaD53 gene expression in different tissues as well as in different concentrations of GA and SL treatments. In this study, we identified three MaD53, three MbD53 and two MiD53 genes in banana. Phylogenetic tree analysis showed that D53 Musa are equally related to D53 Asparagales and Poales. Both high and low-temperature stresses substantially reduced the expression of the MaD53 gene, but osmotic stress treatments had less effect on the expression of the MaD53 gene. GR24 treatment did not significantly promote the height of the banana, but the expression of the MaD53 gene was significantly reduced in roots and leaves. GA treatment at 100 mg/L significantly promoted the expression of the MaD53 gene in roots, but the expression of this gene was significantly reduced in leaves. In this study, we concluded that MaD53 responds to GA and SL treatments, but "Yinniaijiao" dwarf banana may not be sensitive to GA and SL.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

Analysis of Carica papaya Informs Lineage-Specific Evolution of the Aquaporin (AQP) Family in Brassicales

Aquaporins (AQPs), a type of intrinsic membrane proteins that transport water and small solutes across biological membranes, play crucial roles in plant growth and development. This study presents a first genome-wide identification and comparative analysis of the AQP gene family in papaya (Carica papaya L.), an economically and nutritionally important fruit tree of tropical and subtropical regions. A total of 29 CpAQP genes were identified, which represent five subfamilies, i.e., nine plasma intrinsic membrane proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). Although the family is smaller than the 35 members reported in Arabidopsis, it is highly diverse, and the presence of CpXIP genes as well as orthologs in Moringa oleifera and Bretschneidera sinensis implies that the complete loss of the XIP subfamily in Arabidopsis is lineage-specific, sometime after its split with papaya but before Brassicaceae-Cleomaceae divergence. Reciprocal best hit-based sequence comparison of 530 AQPs and synteny analyses revealed that CpAQP genes belong to 29 out of 61 identified orthogroups, and lineage-specific evolution was frequently observed in Brassicales. Significantly, the well-characterized NIP3 group was completely lost; lineage-specific loss of the NIP8 group in Brassicaceae occurred sometime before the divergence with Cleomaceae, and lineage-specific loss of NIP2 and SIP3 groups in Brassicaceae occurred sometime after the split with Cleomaceae. In contrast to a predominant role of recent whole-genome duplications (WGDs) on the family expansion in B. sinensis, Tarenaya hassleriana, and Brassicaceae plants, no recent AQP repeats were identified in papaya, and ancient WGD repeats are mainly confined to the PIP subfamily. Subfamily even group-specific evolution was uncovered via comparing exon-intron structures, conserved motifs, the aromatic/arginine selectivity filter, and gene expression profiles. Moreover, down-regulation during fruit ripening and expression divergence of duplicated CpAQP genes were frequently observed in papaya. These findings will not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further studies of AQP genes in papaya and species beyond.

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.

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

Starch and Sucrose Metabolism and Plant Hormone Signaling Pathways Play Crucial Roles in Aquilegia Salt Stress Adaption

Salt stress is one of the main abiotic stresses that strongly affects plant growth. Clarifying the molecular regulatory mechanism in ornamental plants under salt stress is of great significance for the ecological development of saline soil areas. Aquilegia vulgaris is a perennial with a high ornamental and commercial value. To narrow down the key responsive pathways and regulatory genes, we analyzed the transcriptome of A. vulgaris under a 200 mM NaCl treatment. A total of 5600 differentially expressed genes were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis pointed out that starch and sucrose metabolism and plant hormone signal transduction were significantly improved. The above pathways played crucial roles when A. vulgaris was coping with salt stress, and their protein-protein interactions (PPIs) were predicted. This research provides new insights into the molecular regulatory mechanism, which could be the theoretical basis for screening candidate genes in Aquilegia.

Evolutionary analyses and expression patterns of TCP genes in Ranunculales

TCP transcription factors play a role in a large number of developmental processes and are at the crossroads of numerous hormonal biosynthetic and signaling pathways. The complete repertoire of TCP genes has already been characterized in several plant species, but not in any species of early diverging eudicots. We focused on the order Ranunculales because of its phylogenetic position as sister group to all other eudicots and its important morphological diversity. Results show that all the TCP genes expressed in the floral transcriptome of Nigella damascena (Ranunculaceae) are the orthologs of the TCP genes previously identified from the fully sequenced genome of Aquilegia coerulea. Phylogenetic analyses combined with the identification of conserved amino acid motifs suggest that six paralogous genes of class I TCP transcription factors were present in the common ancestor of angiosperms. We highlight independent duplications in core eudicots and Ranunculales within the class I and class II subfamilies, resulting in different numbers of paralogs within the main subclasses of TCP genes. This has most probably major consequences on the functional diversification of these genes in different plant clades. The expression patterns of TCP genes in Nigella damascena were consistent with the general suggestion that CIN and class I TCP genes may have redundant roles or take part in same pathways, while CYC/TB1 genes have more specific actions. Our findings open the way for future studies at the tissue level, and for investigating redundancy and subfunctionalisation in TCP genes and their role in the evolution of morphological novelties.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

Background

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

Results

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

Conclusions

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01420-1.

The slow-evolving Acorus tatarinowii genome sheds light on ancestral monocot evolution

Monocots are one of the most diverse groups of flowering plants, and tracing the evolution of their ancestral genome into modern species is essential for understanding their evolutionary success. Here, we report a high-quality assembly of the Acorus tatarinowii genome, a species that diverged early from all the other monocots. Genome-wide comparisons with a range of representative monocots characterized Acorus as a slowly evolved genome with one whole-genome duplication. Our inference of the ancestral monocot karyotypes provides new insights into the chromosomal evolutionary history assigned to modern species and reveals the probable molecular functions and processes related to the early adaptation of monocots to wetland or aquatic habitats (that is, low levels of inorganic phosphate, parallel leaf venation and ephemeral primary roots). The evolution of ancestral gene order in monocots is constrained by gene structural and functional features. The newly obtained Acorus genome offers crucial evidence for delineating the origin and diversification of monocots, including grasses.

The genome of Corydalis reveals the evolution of benzylisoquinoline alkaloid biosynthesis in Ranunculales

Species belonging to the order Ranunculales have attracted much attention because of their phylogenetic position as a sister group to all other eudicot lineages and their ability to produce unique yet diverse benzylisoquinoline alkaloids (BIAs). The Papaveraceae family in Ranunculales is often used as a model system for studying BIA biosynthesis. Here, we report the chromosome-level genome assembly of Corydalis tomentella, a species of Fumarioideae, one of the two subfamilies of Papaveraceae. Based on comparisons of sequenced Ranunculalean species, we present clear evidence of a shared whole-genome duplication (WGD) event that has occurred before the divergence of Ranunculales but after its divergence from other eudicot lineages. The C. tomentella genome enabled us to integrate isotopic labeling and comparative genomics to reconstruct the BIA biosynthetic pathway for both sanguinarine biosynthesis shared by papaveraceous species and the cavidine biosynthesis that is specific to Corydalis. Also, our comparative analysis revealed that gene duplications, especially tandem gene duplications, underlie the diversification of BIA biosynthetic pathways in Ranunculales. In particular, tandemly duplicated berberine bridge enzyme-like genes appear to be involved in cavidine biosynthesis. In conclusion, our study of the C. tomentella genome provides important insights into the occurrence of WGDs during the early evolution of eudicots, as well as into the evolution of BIA biosynthesis in Ranunculales.

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

Comprehensive Analysis for GRF Transcription Factors in Sacred Lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera) is an aquatic perennial plant with essential food, ornamental, and pharmacological value. Growth-regulating factor (GRF) is a transcription factor (TF) family that plays an important role in regulating the growth and development of plants. In this study, a comprehensive analysis of the GRF family in N. nucifera was performed, and its role in N. nucifera development was studied. A total of eight GRF genes were identified in the N. nucifera genome. Phylogenetic analysis divided the 38 GRF genes into six clades, while the NuGRFs only contained five clades. The analyses of gene structures, motifs, and cis-acting regulatory elements of the GRF gene family were performed. In addition, the chromosome location and collinearity were analyzed. The expression pattern based on transcriptomic data and real-time reverse transcription-quantitative PCR (qRT-PCR) revealed that the GRF genes were expressed in multiple organs and were abundant in actively growing tissues, and the expression levels decreased as the age of N. nucifera increased. Then, 3D structures of the NuGRF proteins were predicted by homology modeling. Finally, the subcellular localization of GRF1 was ascertained in the tobacco leaf through a vector. Therefore, this study provides a comprehensive overview of the GRF TF family in N. nucifera.

A phylotranscriptome study using silica gel-dried leaf tissues produces an updated robust phylogeny of Ranunculaceae

The utility of transcriptome data in plant phylogenetics has gained popularity in recent years. However, because RNA degrades much more easily than DNA, the logistics of obtaining fresh tissues has become a major limiting factor for widely applying this method. Here, we used Ranunculaceae to test whether silica-dried plant tissues could be used for RNA extraction and subsequent phylogenomic studies. We sequenced 27 transcriptomes, 21 from silica gel-dried (SD-samples) and six from liquid nitrogen-preserved (LN-samples) leaf tissues, and downloaded 27 additional transcriptomes from GenBank. Our results showed that although the LN-samples produced slightly better reads than the SD-samples, there were no significant differences in RNA quality and quantity, assembled contig lengths and numbers, and BUSCO comparisons between two treatments. Using these data, we conducted phylogenomic analyses, including concatenated- and coalescent-based phylogenetic reconstruction, molecular dating, coalescent simulation, phylogenetic network estimation, and whole genome duplication (WGD) inference. The resulting phylogeny was consistent with previous studies with higher resolution and statistical support. The 11 core Ranunculaceae tribes grouped into two chromosome type clades (T- and R-types), with high support. Discordance among gene trees is likely due to hybridization and introgression, ancient genetic polymorphism and incomplete lineage sorting. Our results strongly support one ancient hybridization event within the R-type clade and three WGD events in Ranunculales. Evolution of the three Ranunculaceae chromosome types is likely not directly related to WGD events. By clearly resolving the Ranunculaceae phylogeny, we demonstrated that SD-samples can be used for RNA-seq and phylotranscriptomic studies of angiosperms.

Paleopolyploidies and Genomic Fractionation in Major Eudicot Clades

Eudicots account for ~75% of living angiosperms, containing important food and energy crops. Recently, high-quality genome sequences of several eudicots including Aquilegia coerulea and Nelumbo nucifera have become available, providing an opportunity to investigate the early evolutionary characteristics of eudicots. We performed genomic hierarchical and event-related alignments to infer homology within and between representative species of eudicots. The results provide strong evidence for multiple independent polyploidization events during the early diversification of eudicots, three of which are likely to be allopolyploids: The core eudicot-common hexaploidy (ECH), Nelumbo-specific tetraploidy (NST), and Ranunculales-common tetraploidy (RCT). Using different genomes as references, we constructed genomic alignment to list the orthologous and paralogous genes produced by polyploidization and speciation. This could provide a fundamental framework for studying other eudicot genomes and gene(s) evolution. Further, we revealed significantly divergent evolutionary rates among these species. By performing evolutionary rate correction, we dated RCT to be ~118-134 million years ago (Mya), after Ranunculales diverged with core eudicots at ~123-139 Mya. Moreover, we characterized genomic fractionation resulting from gene loss and retention after polyploidizations. Notably, we revealed a high degree of divergence between subgenomes. In particular, synonymous nucleotide substitutions at synonymous sites (Ks) and phylogenomic analyses implied that A. coerulea might provide the subgenome(s) for the gamma-hexaploid hybridization.