Genetic Contribution of Paleopolyploidy to Adaptive Evolution in Angiosperms

Evolutionary history of magnoliid genomes and benzylisoquinoline alkaloid biosynthesis

Benzylisoquinoline alkaloids (BIAs) are important metabolites synthesized in early-diverging eudicots and magnoliids, yet the genetic basis of BIA biosynthesis in magnoliids remains unclear. Here, we decode the genomes of two magnoliid species, Saruma henryi and Aristolochia manshuriensis, and reconstruct the ancestral magnoliid karyotype and infer the chromosomal rearrangement history following magnoliid diversification. Metabolomic, transcriptomic, and phylogenetic analyses reveal the intermediate chemical components and genetic basis of BIA biosynthesis in A. manshuriensis. Although the core enzymes involved in BIA synthesis appear to be largely conserved between early-diverging eudicots and magnoliids, the biosynthetic pathways in magnoliids seem to exhibit greater flexibility. Significantly, our investigation of the evolutionary history of BIA biosynthetic genes revealed that almost all were duplicated before the emergence of extant angiosperms, with only early-diverging eudicots and magnoliids preferentially retaining these duplicated genes, thereby enabling the biosynthesis of BIAs in these groups.

Nanoparticles alter the nature and strength of intraploidy and interploidy interactions in plants

Engineered nanoparticles have profound impacts on organisms, yet there is limited understanding of how nanoparticle exposure shapes species interactions that are key to natural community dynamics. By growing plants of the same (intraploidy) and different ploidy levels (interploidy) of Fragaria in axenic microcosms, we examined the influence of nanoparticles on species interactions in polyploid and diploid plants. Under copper oxide (CuO) nanoparticle exposure, polyploids experienced reduced competition and a shift towards facilitation, when growing with both polyploids (the effect of polyploids on polyploids measured by the relative interaction index, RII8x,8x) and diploids (the effect of diploids on polyploids, RII8x,2x). This reduction in competitive interactions in polyploids, in line with the stress gradient hypothesis, was primarily driven by nanoscale effects. In contrast, the strength of competitive interactions (RII8x,8× and RII8x,2x) increased under CuO bulk particles compared to control conditions. Different from polyploids, diploids experienced neutral interactions (RII2x,2x and RII2x,8x) under both nanoparticles and bulk particles. These findings highlight ploidy-specific interaction dynamics and the need to consider species interactions when predicting organismal responses to nanoparticle pollution in ecological communities, providing critical insights for conservation strategies and sustainable nanotechnology applications.

Dynamics of Mixed-Ploidy Populations under Demographic and Environmental Stochasticities

AbstractThe population dynamics of autopolyploids-organisms with more than two genome copies of a single species-and their diploid progenitors have been extensively studied. The acquisition of multiple genome copies is heavily influenced by stochasticity, which strongly suggests the efficacy of a probabilistic approach to examine the long-term dynamics of a population with multiple cytotypes. Yet our current understanding of the dynamics of autopolyploid populations has not incorporated stochastic population dynamics and coexistence theory. To investigate the factors contributing to the probability and stability of coexisting cytotypes, we designed a new population dynamics model that incorporates demographic and environmental stochasticities to simulate the formation, establishment, and persistence of diploids, triploids, and autotetraploids in the face of gene flow among cytotypes. We found that increased selfing rates and pronounced reproductive isolation promote coexistence of multiple cytotypes. In stressful environments and with strong competitive effects among cytotypes, these dynamics are more complex; our stochastic modeling approach reveals the resulting intricacies that give autotetraploids competitive advantage over their diploid progenitors. Our work is foundational for a better understanding of the dynamics of coexistence of multiple cytotypes.

Comprehensive Genomic Dataset of Chinese Lizardtail Herb and Comparative Genomic Analysis Provide Insights Into Its Paleo-Polyploidization Event

The Chinese lizardtail herb, Saururus chinensis, holds a prominent position in traditional Chinese medicine. In this study, we present a comprehensive genomic dataset for S. chinensis. Furthermore, comparative genomic analysis indicates that the extant genome of S. chinensis retains extensive traces of a paleo-tetraploidization event. These traces are observable at both the macroscopic level of chromosomes and the microscopic level within specific gene families, such as the PEL (pseudo-etiolation in light) gene family. Additionally, our findings further suggest that this paleo-tetraploidization event drives an expansion of the PEL gene family in the S. chinensis genome, potentially facilitating its neo- and sub-functionalization, and thereby contributing to the evolutionary adaptability of this species.

Chromosomal-Level Genome Suggests Adaptive Constraints Leading to the Historical Population Decline in an Extremely Endangered Plant

Increased human activity and climate change have significantly impacted wild habitats and increased the number of endangered species. Exploring evolutionary history and predicting adaptive potential using genomic data will facilitate species conservation and biodiversity recovery. Here, we examined the genome evolution of a critically endangered tree Pellacalyx yunnanensis, a plant species with extremely small populations (PSESP) that is narrowly distributed in Xishuangbanna, China. The species has neared extinction due to economic exploitation in recent decades. We assembled a chromosome-level genome of 334 Mb, with the N50 length of 20.5 Mb. Using the genome, we discovered that P. yunnanensis has undergone several population size reductions, leading to excess deleterious mutations. The species may possess low adaptive potential due to reduced genetic diversity and the loss of stress-responsive genes. We estimate that P. yunnanensis is the basal species of its genus and diverged from its relatives during global cooling, suggesting it was stranded in unsuitable environments during periods of dramatic climate change. In particular, the loss of seed dormancy leads to germination under unfavourable conditions and reproduction challenges. This dormancy loss may have occurred through genetic changes that suppress ABA signalling and the loss of genes involved in seed maturation. The high-quality genome has also enabled us to reveal phenotypic trait evolution in Rhizophoraceae and identify divergent adaptation to intertidal and inland habitats. In summary, our study elucidates mechanisms underlying the decline and evaluates the adaptive potential of P. yunnanensis to future climate change, informing future conservation efforts.

Evolutionary landscape of plant chalcone isomerase-fold gene families

Flavonoids are crucial for plant survival and adaptive evolution, and chalcone isomerase (CHI) genes serve as key rate-limiting gene in the flavonoid biosynthesis pathway. It is important for plant adaptive evolution to comprehensively study the evolution and diversity of the CHI gene families. However, the CHI gene families in many plant lineages remain elusive. This study systematically identified CHI genes from 259 species including algae, bryophytes, ferns, gymnosperms, and angiosperms. A total of 1,738 CHI gene family members were discovered. We analyzed the diversity, distribution trajectory, and the driving forces of gene duplication during the evolution of the plant lineages. The present study is the first to identify potential type II and type IV CHI genes in the extant liverwort model species Marchantia polymorpha. The distribution pattern of CHI genes across the plant kingdom reveals that the origin of type II CHI can be traced back to the last common ancestor of bryophytes and vascular plants, and type III CHI may represent the ancestral form of the CHI gene family. The identification of conserved motifs showed significant differences in motif distribution among different CHI gene types. It was found that the drivers of gene duplication varied across plant lineages: dispersed duplications (DSD) were predominant in algae and bryophytes, whole-genome duplication (WGD) was the main driver in basal angiosperms and monocots, while tandem duplications (TD) predominating in eudicots. Structural clustering analysis demonstrated the 3-layer sandwich structure in the CHI-fold proteins remained conserved in the central region, while repeated loss of N-terminal sequences contributed to structural diversity. This study provides a deeper understanding of the evolution and diversity of the CHI-fold proteins and lays a theoretical foundation for further studies of their function and the identification of new functional CHI genes.

Genome-wide identification and analysis of phosphate utilization related genes (PURs) reveal their roles involved in low phosphate responses in Brassica napus L

BACKGROUND

Phosphorus (P) is an essential macronutrient for Brassica napus L. growth and development, and is mainly acquired from the soil as phosphate (Pi). However, there is no research on the system analysis of Pi utilization related genes (PURs) in B. napus yet.

RESULTS

In this study, 285 PURs were identified in B. napus genome, including 4 transcription factor (TF) gene families (83 genes) and 17 structural gene families (202 genes). Subcellular localization analysis showed that the proteins encoded by B. napus PURs were mainly located in the nucleus (~ 46.0%) and cell membrane (~ 36.5%). Chromosome localization analysis suggested that B. napus PURs were distributed on An (131) and Cn (149) subgenomes without bias. Analysis of 35 representative species confirmed that PURs were widely present in plants ranging from Chlorophyta to angiosperms with a rapid expansion trend. Collinearity analysis revealed that allopolyploidization and small-scale duplication events resulted in the large expansion of B. napus PURs. For each gene pair of B. napus PURs, the sequence identity of promoter was significantly lower than that of CDS, proving the significant difference in promoter region that might be related to the divergence of PURs expression and function. Transcription factor (TF) binding site prediction, cis-element analysis, and microRNA prediction suggested that the expressions of B. napus PURs are regulated by multiple factors including 32 TF gene families (362), 108 types of CRE (29,770) and 25 types of miRNAs (66). Spatiotemporal expression analysis demonstrated that B. napus PURs were widely expressed during the whole developmental stages, and most synteny-gene pairs (76.42%) shared conserved expression patterns. RNA-seq analyses revealed that most B. napus PURs were induced by low Pi stress, and the hub genes were generally the Pi transporter (PHT) family members. qRT-PCR analysis proved that the expression levels of four B. napus PURs were positively correlated with the root system architecture of three B. napus varieties under low Pi supply at the seedling stage.

CONCLUSION

The 285 PURs were identified from B. napus with strong LP inducible expression profile. Our findings regarding the evolution, transcriptional regulation, and expression of B. napus PURs provide valuable information for further functional research.

Ancient large-scale gene duplications and diversification in bryophytes illuminate the plant terrestrialization

Large-scale gene duplications (LSGDs) are crucial for evolutionary adaptation and recurrent in vascular plants. However, the role of ancient LSGDs in the terrestrialization and diversification of bryophytes, the second most species-rich group of land plants, remains largely elusive due to limited sampling in bryophytes. Employing the most extensive nuclear gene dataset in bryophytes to date, we reconstructed a time-calibrated phylogenetic tree from 209 species, covering virtually all key bryophyte lineages, for phylogenomic analyses of LSGDs and diversification. We newly identified two ancient LSGDs: one in the most recent common ancestor (MRCA) of extant bryophytes and another in the MRCA of the majority of Jungermanniales s. lato. Duplicated genes from these two LSGDs show significant enrichment in photosynthesis-related processes and structures. Rhizoid-responsive ROOTHAIR DEFECTIVE SIX-LIKE (RSL) genes from ancient LSGDs are present in rhizoidless bryophytes, challenging assumptions about rhizoid absence mechanisms. We highlighted four major diversification rate upshifts, two of which slightly postdated LSGDs, potentially linked to the flourishing of gymnosperms and angiosperms and explaining over 80% of bryophyte diversity. Our findings, supported by extensive bryophyte sampling, highlight the significance of LSGDs in the early terrestrialization and diversification of bryophytes, offering new insights into land plant evolution.

S-RNase evolution in self-incompatibility: Phylogenomic insights into synteny with Class I T2 RNase genes

S-RNases are essential in the gametophytic self-incompatibility (GSI) system of many flowering plants, where they act as stylar-S determinants. Despite their prominence, the syntenic genomic origin and evolutionary trajectory of S-RNase genes in eudicots have remained largely unclear. Here, we performed large-scale phylogenetic and microsynteny network analyses of T2 RNase genes across 130 angiosperm genomes, encompassing 35 orders and 56 families. S-like RNase genes in Cucurbitaceae species phylogenetically grouped with functionally characterized S-RNases in various species. Additionally, Cucurbitaceae S-like RNase genes showed conserved synteny with Class I T2 RNase genes. From this, we inferred that the well-characterized S-RNase genes (belonging to Class III-A genes) and Class I T2 RNase genes (located on duplicated genomic blocks) likely derived from the gamma triplication event shared by core eudicots. Additionally, we identified frequent lineage-specific gene transpositions of S-RNases and S-like RNases across diverse angiosperm lineages, including Rosaceae, Solanaceae, and Rutaceae families, accompanied by a significant increase in transposable element activity near these genes. Our findings delineate the genomic origin and evolutionary path of eudicot S-RNase genes, enhancing our understanding of the evolution of the S-RNase-based GSI system.

The role of DNA content in shaping chromatin architecture and gene expression

Whole-genome duplication is an evolutionary force that drives speciation in all living kingdoms and is notably prevalent in plants. The evolutionary history of plants involved at least two genomic duplications that significantly expanded the plant morphology and physiology spectrum. Many important crops are polyploids, showing valuable features relative to morphological and stress response traits. After genome duplication, diploidization processes facilitate genomic adjustments to restore disomic inheritance. However, little is known about the chromatin changes triggered by nuclear DNA content alterations. Here, we report that synthetically induced genome duplication leads to chromatinization and significant changes in gene expression, resulting in a transcriptional landscape resembling a natural tetraploid. Interestingly, synthetic diploidization elicits only minor alterations in transcriptional activity and chromatin accessibility compared to the more pronounced effects of tetraploidization. We identified epigenetic factors, including specific histone variants, that showed increased expression following genome duplication and decreased expression after genome reduction. These changes may play a key role in the epigenetic mechanisms underlying the phenotypic complexity after tetraploidization in plants. Our findings shed light on the mechanisms that modulate chromatin accessibility remodeling and gene transcription regulation underlying plant genome adaptation in response to changes in genome size.

Introgression among subgroups is an important driving force for genetic improvement and evolution of the Asian cultivated rice Oryza sativa L

Anagenesis accumulates favorable mutations that enable crops to adapt to continually improving artificial production environments, while cladogenesis results in the deposition of beneficial variations across diverse ecotypes. Integrating advantageous genetic variations from diverse evolutionary sources establishes the foundation for the continued genetic improvement of crops. For a long time, rice breeding practices have been guided by the established belief that the Asian cultivated rice consists of two subspecies: Oryza sativa subsp. indica and subsp. japonica. Integrating elite genetic variants from both subspecies has been a major strategy for genetic improvement. This approach has proven successful through the achievements of temperate japonica breeding programs in China, Japan, and Korea over the past decades. The genetic differentiation within the Asian cultivated rice has been successfully harnessed for heterosis breeding, thereby enhancing rice yield productivity. Genomic investigations have revealed more genetic divergences in the Asian cultivated rice, prompting the proposal of six subgroups within it. This indicates that there is greater potential for uncovering additional genetic divergences and diversity in future breeding practices. Genetic introgression and gene flow among subgroups have led to improvements in agronomic traits within the indica, temperate japonica, and tropical japonica subgroups during the modern rice breeding process. The introgression process has widened the genetic diversity within subgroups and reduced the genetic distance between them, resulting in the creation of new genetic blocks and subpopulations. Artificial introgression has accelerated the evolution process in rice breeding history. Advancements in the study of genetic divergence and diversity in rice offer valuable insights to guide breeding practices. The mini subgroups aus, basmatic, and rayada possess untapped genetic potential but have been poorly studied worldwide; more samples should be further investigated. This information will be invaluable for harnessing these advantageous variations through introgression breeding. Further studying the nature of reproductive barriers among subgroups will enhance our understanding of genetic differentiation, allow us to overcome these barriers and facilitate effective genetic exchange, and even enable us to harness heterosis among subgroups.

Genome-wide identification of starch metabolism gene families in Potentilla anserina and the expression pattern in response to abiotic stress factors

BACKGROUND

The tuberous roots of Potentilla anserina (Pan), which are called silverweed cinquefoil roots, serve as a source of starch for the inhabitants of the Qinghai-Tibet Plateau. They are also regarded as a valuable tonic food and herbal medicine in the ethnic medicine system. Starch plays a crucial part in the plant's life cycle, particularly during the growth stage and in response to abiotic stress. Moreover, numerous biological processes and regulatory networks are involved in the synthesis and accumulation of starch.

RESULTS

In this research, a combination of transcriptome and metabolomics approaches were employed to analyze the genes related to starch synthesis and degradation in Pan. The crucial genes involved in starch metabolism were identified, and the response characteristics of these genes to drought and low temperature treatments were investigated. Seven AMYs (Alpha-amylases) and 18 BAMs (Beta-amylases) were identified from the genome of Pan. Molecular phylogenetic analyses of AMYs and BAMs derived from 11 species of rosids were conducted respectively. AMYs of Pan and other species were clustered into 3 groups, whereas BAMs were classified into 4 groups. In the tuberous roots, when compared with the control conditions, 2 AMYs and 4 BAMs were upregulated, while 3 BAMs were downregulated, meanwhile, the contents of maltose and glucose-6-phosphate (G6P) were decreased, while the content of glucose-1-phosphate (G1P) was increased under drought treatment. In the tuberous roots under low temperature treatment, 1 BAM was upregulated, while 2 AMYs and 4 BAMs were downregulated compared with control. The response characteristics of AMYs and BAMs to drought and low temperature treatments were further verified through qPCR analysis.

CONCLUSIONS

In this research, the genes associated with starch synthesis in Pan were characterized, and the functions of AMYs and BAMs in abiotic stress treatments were elucidated. An overview of evolution of AMY and BAM gene families within rosids was also provided.

Identification of a central regulator of ginkgolide biosynthesis in Ginkgo biloba that integrates jasmonate and light signaling

Ginkgolides are secondary metabolites unique to Ginkgo biloba with the potential to prevent and treat cardiovascular and cerebrovascular diseases. Although the biosynthetic pathways of ginkgolides have been partly uncovered, the mechanism regulating their biosynthesis is still largely unknown. Here, using multiomic and genetic analyses, we report the identification of a transcription factor, named ETHYLENE RESPONSE FACTOR ASSOCIATED WITH GINKGOLIDE BIOSYNTHESIS (GbEAG), as a critical regulator of ginkgolide biosynthesis. GbEAG is highly expressed in the roots of G. biloba, and its expression is significantly induced by methyl jasmonate (MeJA). Ginkgolide content was significantly increased in roots by overexpressing GbEAG using a "cut-dip-regeneration" system. GbEAG positively regulates ginkgolide biosynthesis by directly binding to the GCC-boxes in the promoter regions of genes involved in the biosynthesis of ginkgolides, such as ISOPENTENYL DIPHOSPHATE ISOMERASE (GbIDI) and CYTOCHROME P450 7005C3 (GbCYP7005C3). GbEAG mediates the jasmonic acid (JA)-activated ginkgolide synthesis through its direct interaction with the JASMONATE ZINC-FINGER INFLORESCENCE MERISTEM DOMAIN 3 (GbJAZ3) repressor. Importantly, we also found that the central light-response regulator ELONGATED HYPOCOTYL 5 (GbHY5) mediates light induction of ginkgolide biosynthesis by binding to the G-box in the GbEAG promoter. Our findings provide mechanistic insights into the coordinated regulation of ginkgolide biosynthesis via JA and light signals, with GbEAG as a central regulator in G. biloba, and shed light on the potential to develop ginkgolide-rich varieties through molecular breeding and gene editing.

Convergent evolution in angiosperms adapted to cold climates

Convergent and parallel evolution occur more frequently than previously thought. Here, we focus on the evolutionary adaptations of angiosperms at sub-zero temperatures. We begin by introducing the history of research on convergent and parallel evolution, defining all independent similarities as convergent evolution. Our analysis reveals that frost zones (periodic or constant), which cover 49.1% of Earth's land surface, host 137 angiosperm families, with over 90% of their species thriving in these regions. In this context, we revisit the global biogeography and evolutionary trajectories of plant traits, such as herbaceous form and deciduous leaves, that are thought to be evasion strategies for frost adaptation. At the physiological and molecular levels, many angiosperms have independently evolved cold acclimation mechanisms through multiple pathways in addition to the well-characterized C-repeat binding factor/dehydration-responsive element binding protein 1 (CBF/DREB1) regulatory pathway. These convergent adaptations have occurred across various molecular levels, including amino acid substitutions and changes in gene duplication and expression within the same or similar functional pathways; however, identical amino acid changes are rare. Our results also highlight the prevalence of polyploidy in frost zones and the occurrence of paleopolyploidization events during global cooling. These patterns suggest repeated evolution in cold climates. Finally, we discuss plant domestication and predict climate zone shifts due to global warming and their effects on plant migration and in situ adaptation. Overall, the integration of ecological and molecular perspectives is essential for understanding and forecasting plant responses to climate change.

Genomic insights into drought adaptation of the forage shrub Caragana korshinskii (Fabaceae) widely planted in drylands

The Korshinsk peashrub (Caragana korshinskii), known for its exceptional drought tolerance, is widely cultivated in arid and semi-arid regions for vegetation restoration and as a vital forage plant. To elucidate the genomic basis of its drought tolerance, we generated a chromosomal-scale genome sequence of C. korshinskii. Our synteny analysis disputes the previously hypothesized genus-specific whole-genome duplication event, as suggested by earlier transcriptome study of this species and its congeners. We identified that tandem duplications were critical for the expansion of gene families, such as early light-induced protein, heat shock protein 100, and Dehydrin, which are involved in cellular protection processes. These expansions are likely pivotal to the superior drought tolerance observed in C. korshinskii, as evidenced by the elevated gene expression of these genes under drought conditions. Furthermore, overexpression studies of seven tandemly duplicated DHN genes revealed a substantial enhancement in drought survival rates of seedlings, likely attributable to increased gene dosage effects. Conversely, gene silencing via virus-induced gene silencing demonstrated opposing effects. Additionally, we have established the CakorDB, a genomic resource database for C. korshinskii (https://bis.zju.edu.cn/cakordb/), accessible freely to the scientific community. Collectively, our study not only provides a valuable genomic resource for the Korshinsk peashrub but also highlights the genetic adaptations that enable C. korshinskii to thrive in desert environments, positioning its stress-responsive genes as a valuable genetic reservoir for breeding drought-resistant crops.

Multi-omics analysis reveals the role of UGT72 family genes in arbutin biosynthesis in Pyrus and evolution driven by whole genome duplication

The UGT72 gene family encodes proteins that glycosylate phenylpropanoids, and thus contribute to the synthesis of various phenolic substances. However, their functional role and evolutionary history in Pyrus spp. remains poorly understood. Here we explored the evolution, amplification, coding region structural variation, and functional divergence of the UGT72 gene family and its subfamilies. Further, we identified functional genes involved in arbutin synthesis and functionally validated the key genes. 15 UGT72 genes were identified in the complete genome sequence and classified into two subfamilies of Pyrus betulifolia. Significant expansion of the UGT72 gene family occurred after genome duplication in P. betulifolia. 53.33 % of all UGT72 family genes were found to have undergone expansion via WGD/segmental duplication. A noteworthy discovery was that the amplification of functional genes such as PbUGT72B1714 during polyploidization, combined with the loss of vital motifs and variations at important sites within these genes, significantly impacted the diversification of arbutin metabolism. These findings offer novel insights into how gene gains and losses caused by WGDs have contributed to metabolic diversification and evolutionary adaptation in Pyrus, as well as a groundwork for more detailed investigations into the mechanisms of arbutin metabolism.

A pivotal switch in β-domain determines the substrate selectivity of terpene synthases involved in Gardenia jasminoides floral scent synthesis

Plants have evolved a diverse array of secondary metabolites to enhance their adaptability to environmental stresses, with volatile terpenoids being a notable example. Gardenia (Gardenia jasminoides), celebrated for its unique fragrance, is a key natural source of volatile terpenoids. Using our chromosome-level genome and transcriptome data of G. jasminoides, we previously identified six terpene synthases (TPSs) involved in the production of its floral scent. Here, we functionally characterized these six key TPS enzymes, aligning their product profiles with volatile organic compound (VOC) analysis. Notably, we identified two highly similar TPSs, GjTPS1 and GjTPS2, which share high sequence homology but differ in substrate selectivity. Through AI-based predictions and site-directed mutagenesis, we pinpointed a single amino acid in the β-domain that acts as a "switch," modulating substrate selectivity-an unusual finding, as this residue is outside the active site region. Comparative genomic analyses with Coffea canephora and other eudicots revealed that G. jasminoides TPS genes primarily expanded through dispersed duplications (DSD) and tandem duplications (TD). Our study is the first to reveal a "switch" in a previously deemed "non-functional" region, regulating substrate selectivity in TPS proteins. These insights could facilitate breeding elite Gardenia varieties and support the rational engineering of terpene synthases.

DNA methylation dynamics in gymnosperm duplicate genes: implications for genome evolution and stress adaptation

Duplicate genes are pivotal in driving evolutionary innovation, often exhibiting expression divergence that offers a system to investigate the role of DNA methylation in transcriptional regulation. However, previous studies have predominantly focused on angiosperms, leaving the methylation patterns in major lineages of land plants still unclear. This study explores DNA methylation evolution in duplicate genes across representative gymnosperm species with large genomes, spanning over 300 million years, using genomic, transcriptomic, and high-depth DNA methylomic data. We observed variations in DNA methylation levels along gene bodies, flanking regions, and methylation statuses of coding regions across different duplication types. Biased divergences in DNA methylation and gene expression frequently occurred between duplicate copies. Specifically, methylation divergences in the 2-kb downstream regions negatively correlated with gene expression. Both CG and CHG DNA methylation in gene bodies were positively correlated with gene length, suggesting these methylation types may function as an epigenomic buffer to mitigate the adverse impact of gene length on expression. Duplicate genes exhibiting both methylation and expression divergences were notably enriched in adaptation-related biological processes, suggesting that DNA methylation may aid adaptive evolution in gymnosperms by regulating stress response genes. Changes in expression levels correlated with switches in methylation status within coding regions of transposed duplicates. Specifically, depletion for CG methylation or enrichment for non-CG methylation significantly reduced the expression of translocated copies. This correlation suggests that DNA methylation may reduce genetic redundancy by silencing translocated copies. Our study highlights the significance of DNA methylation in plant genome evolution and stress adaptation.

Bona Fide Plant Steroid Receptors are Innovated in Seed Plants and Angiosperms through Successive Whole-Genome Duplication Events

Whole-genome duplication (WGD) events are widespread in plants and animals, thus their long-term evolutionary contribution has long been speculated, yet a specific contribution is difficult to verify. Here, we show that ɛ-WGD and ζ-WGD contribute to the origin and evolution of bona fide brassinosteroid (BR) signaling through the innovation of active BR biosynthetic enzymes and active BR receptors from their respective ancestors. We found that BR receptors BRI1 (BR INSENSITIVE 1) and BRL1/3 (BRI1-LIKES 1/3) derived by ɛ-WGD and ζ-WGD, which occurred in the common ancestor of angiosperms and seed plants, respectively, while orphan BR receptor BRL2 first appeared in stomatophytes. Additionally, CYP85A enzymes synthesizing the bioactive BRs derived from a common ancestor of seed plants, while its sister enzymes CYP90 synthesizing BR precursors presented in all land plants, implying possible ligand-receptor coevolution. Consistently, the island domains (IDs) responsible for BR perception in BR receptors were most divergent among different receptor branches, supporting ligand-driven evolution. As a result, BRI1 was the most diversified BR receptor in angiosperms. Importantly, relative to the BR biosynthetic DET2 gene presented in all land plants, BRL2, BRL1/3 and BRI1 had high expression in vascular plants ferns, gymnosperms and angiosperms, respectively. Notably, BRI1 is the most diversified BR receptor with the most abundant expression in angiosperms, suggesting potential positive selection. Therefore, WGDs initiate a neofunctionalization process diverged by ligand-perception and transcriptional expression, which might optimize both BR biosynthetic enzymes and BR receptors, likely contributing to the evolution of land plants, especially seed plants and angiosperms.

Recent and Recurrent Autopolyploidization Fueled Diversification of Snow Carp on the Tibetan Plateau

Whole-genome duplication (WGD), or polyploidization, is a major contributor to biodiversity. However, the establishment and survival of WGDs are often considered to be stochastic, since elucidating the processes of WGD establishment remains challenging. In the current study, we explored the processes leading to polyploidy establishment in snow carp (Cyprinidae: Schizothoracinae), a predominant component of the ichthyofauna of the Tibetan Plateau and its surrounding areas. Using large-scale genomic data from isoform sequencing, we analyzed ohnolog genealogies and divergence in hundreds to thousands of gene families across major snow carp lineages. Our findings demonstrated that independent autopolyploidization subsequent to speciation was prevalent, while autopolyploidization followed by speciation also occurred in the diversification of snow carp. This was further supported by matrilineal divergence and drainage evolution evidence. Contrary to the long-standing hypothesis that ancient polyploidization preceded the diversification of snow carp, we determined that polyploidy in extant snow carp was established by recurrent autopolyploidization events during the Pleistocene. These findings indicate that the diversification of extant snow carp resembles a coordinated duet: first, the uplift of the Tibetan Plateau orchestrated the biogeography and diversification of their diploid progenitors; then, the extensive Pliocene-Pleistocene climate changes acted as relay runners, further fueling diversification through recurrent autopolyploidization. Overall, this study not only reveals a hitherto unrecognized recent WGD lineage in vertebrates but also advances current understanding of WGD processes, emphasizing that WGD establishment is a nonstochastic event, emerging from numerous adaptations to environmental challenges and recurring throughout evolutionary history rather than merely in plants.

Genomic convergence in terrestrial root plants through tandem duplication in response to soil microbial pressures

Despite increasing reports of convergent adaptation, evidence for genomic convergence across diverse species worldwide is lacking. Here, our study of 205 Archaeplastida genomes reveals evidence of genomic convergence through tandem duplication (TD) across different lineages of root plants despite their genomic diversity. TD-derived genes, notably prevalent in trees with developed root systems embedded in soil, are enriched in enzymatic catalysis and biotic stress responses, suggesting adaptations to environmental pressures. Correlation analyses suggest that many factors, particularly those related to soil microbial pressures, are significantly associated with TD dynamics. Conversely, flora transitioned to aquatic, parasitic, halophytic, or carnivorous lifestyles-reducing their interaction with soil microbes-exhibit a consistent decline in TD frequency. This trend is further corroborated in mangroves that independently adapted to hypersaline intertidal soils, characterized by diminished microbial activity. Our findings propose TD-driven genomic convergence as a widespread adaptation to soil microbial pressures among terrestrial root plants.

THI1 Gene Evolutionary Trends: A Comprehensive Plant-Focused Assessment via Data Mining and Large-Scale Analysis

Molecular evolution analysis typically involves identifying selection pressure and reconstructing evolutionary trends. This process usually requires access to specific data related to a target gene or gene family within a particular group of organisms. While recent advancements in high-throughput sequencing techniques have resulted in the rapid accumulation of extensive genomics and transcriptomics data and the creation of new databases in public repositories, extracting valuable insights from such vast data sets remains a significant challenge for researchers. Here, we elucidated the evolutionary history of THI1, a gene responsible for encoding thiamine thiazole synthase. The thiazole ring is a precursor for vitamin B1 and a crucial cofactor in primary metabolic pathways. A thorough search of complete genomes available within public repositories reveals 702 THI1 homologs of Archaea and Eukarya. Throughout its diversification, the plant lineage has preserved the THI1 gene by incorporating the N-terminus and targeting the chloroplasts. Likewise, evolutionary pressures and lifestyle appear to be associated with retention of TPP riboswitch sites and consequent dual posttranscriptional regulation of the de novo biosynthesis pathway in basal groups. Multicopy retention of THI1 is not a typical plant pattern, even after successive genome duplications. Examining cis-regulatory sites in plants uncovers two shared motifs across all plant lineages. A data mining of 484 transcriptome data sets supports the THI1 homolog expression under a light/dark cycle response and a tissue-specific pattern. Finally, the work presented brings a new look at public repositories as an opportunity to explore evolutionary trends to THI1.

Dosage sensitivity shapes balanced expression and gene longevity of homoeologs after whole-genome duplications in angiosperms

Following whole-genome duplication (WGD), duplicate gene pairs (homoeologs) can evolve varying degrees of expression divergence. However, the determinants influencing these relative expression level differences (RFPKM) between homoeologs remain elusive. In this study, we analyzed the RFPKM between homoeologs in 3 angiosperms, Nymphaea colorata, Nelumbo nucifera, and Acorus tatarinowii, all having undergone a single WGD since the origin of angiosperms. Our results show significant positive correlations in RFPKM of homoeologs among tissues within the same species, and among orthologs across these 3 species, indicating convergent expression balance/bias between homoeologous gene copies following independent WGDs. We linked RFPKM between homoeologs to gene attributes associated with dosage-balance constraints, such as protein-protein interactions, lethal-phenotype scores in Arabidopsis (Arabidopsis thaliana) orthologs, domain numbers, and expression breadth. Notably, homoeologs with lower RFPKM often had more interactions and higher lethal-phenotype scores, indicating selective pressures favoring balanced expression. Also, homoeologs with lower RFPKM were more likely to be retained after WGDs in angiosperms. Within Nelumbo, greater RFPKM between homoeologs correlated with increased cis- and trans-regulatory differentiation between species, highlighting the ongoing escalation of gene expression divergence. We further found that expression degeneration in 1 copy of homoeologs is inclined toward nonfunctionalization. Our research highlights the importance of balanced expression, shaped by dosage-balance constraints, in the evolutionary retention of homoeologs in plants.

Disentangling sources of gene tree discordance for Hordeum species via target-enriched sequencing assays

Hordeum is an economically and evolutionarily important genus within the Triticeae tribe of the family Poaceae, and contains 33 widely distributed and diverse species which cytologically represent four subgenomes (H, Xa, Xu and I). These wild species (except Hordeum spontaneum, which is the primary gene pool of barley) are secondary or tertiary gene-pool germplasms for barley and wheat improvement, and uncovering their complicated evolutionary relationships would benefit for future breeding programs. Here, we developed a complexity-reduced pipeline via capturing genome-wide distributed fragments via two novel target-enriched assays (HorCap v1.0 and BarPlex v1.0) in conjugation with high-throughput sequencing of the enrichments. Both assays were tested for genotyping 40 species from three genera (Hordeum, Triticum, and Aegilops) containing 82 samples 67 accessions. Either of both assays worked efficiently in genotyping, while integration of both assays can significantly improve the robustness and resolution of the Hordeum phylogenetic trees. Interestingly, the incomplete lineage sorting (ILS) was inferred for the first time as the major factor causing phylogenetic discordance among the four subgenomes, whereas in New World species (carrying I genome) post-speciation introgression events were revealed. Through revising the evolutionary relationships of the Hordeum species based on an ancestral state reconstruction for the diploids and parental donor inference for the polyploids, our results raised new queries about the Hordeum phylogeny. Moreover, both newly-developed assays are applicable in genotyping and phylogenetic analysis of Hordeum and other Triticeae wild species.

Insights from the first chromosome-level genome assembly of the alpine gentian Gentiana straminea Maxim

Gentiana straminea Maxim. is a perennial herb and mainly distributed in the Qinghai-Tibetan Plateau. To adapt to the extreme environment, it has developed particular morphological, physiological, and genetic structures. Also, rich in iridoids, it is one of the original plants of traditional Chinese herb 'Qinjiao'. Herein, we present its first chromosome-level genome sequence assembly and compare it with the genomes of other Gentiana species to facilitate the analysis of genomic characteristics. The assembled genome size of G. straminea was 1.25 Gb, with a contig N50 of 7.5 Mb. A total of 96.08% of the genome sequences was anchored on 13 pseudochromosomes, with a scaffold N50 of 92.70 Mb. A total of 54,310 protein-coding genes were predicted, 80.25% of which were functionally annotated. Comparative genomic analyses indicated that G. straminea experienced two whole-genome duplication events after the γ whole-genome triplication with other eudicots, and it diverged from other Gentiana species at ~3.2 Mya. A total of 142 enzyme-coding genes related to iridoid biosynthesis were identified in its genome. Additionally, we identified differences in the number and expression patterns of iridoid biosynthetic pathway genes in G. straminea compared with two other Gentiana species by integrating whole-genome sequence and transcriptomic analyses.

High-quality chromosome-level genomic insights into molecular adaptation to low-temperature stress in Madhuca longifolia in southern subtropical China

BACKGROUND

Madhuca longifolia, the energy-producing and medicinal tropical tree originally from southern India, faces difficulties in adapting to the low temperatures of late autumn and early winter in subtropical southern China, impacting its usability. Therefore, understanding the molecular mechanisms controlling the ability of this species to adapt to environmental challenges is essential for optimising horticulture efforts. Accordingly, this study aimed to elucidate the molecular responses of M. longifolia to low-temperature stress through genomic and transcriptomic analyses to inform strategies for its effective cultivation and utilisation in colder climates.

RESULTS

Herein, the high-quality reference genome and genomic assembly for M. longifolia are presented for the first time. Using Illumina sequencing, Hi-C technology, and PacBio HiFi sequencing, we assembled a chromosome-level genome approximately 737.92 Mb in size, investigated its genomic features, and conducted an evolutionary analysis of the genus Madhuca. Additionally, using transcriptome sequencing, we identified 17,941 differentially expressed genes related to low-temperature response. Through bioinformatics analysis of the WRKY gene family, 15 genes crucial for M. longifolia low-temperature resistance were identified.

CONCLUSIONS

This research not only lays the groundwork for the successful ecological adaptation and cultivation of M. longifolia in China's southern subtropical regions but also offers valuable insights for the genetic enhancement of cold tolerance in tropical species, contributing to their sustainable horticulture and broader industrial, medicinal, and agricultural use.

Chromosome-level assembly of Lindenbergia philippensis and comparative genomic analyses shed light on genome evolution in Lamiales

Lamiales, comprising over 23,755 species across 24 families, stands as a highly diverse and prolific plant group, playing a significant role in the cultivation of horticultural, ornamental, and medicinal plant varieties. Whole-genome duplication (WGD) and its subsequent post-polyploid diploidization (PPD) process represent the most drastic type of karyotype evolution, injecting significant potential for promoting the diversity of this lineage. However, polyploidization histories, as well as genome and subgenome fractionation following WGD events in Lamiales species, are still not well investigated. In this study, we constructed a chromosome-level genome assembly of Lindenbergia philippensis (Orobanchaceae) and conducted comparative genomic analyses with 14 other Lamiales species. L. philippensis is positioned closest to the parasitic lineage within Orobanchaceae and has a conserved karyotype. Through a combination of Ks analysis and syntenic depth analysis, we reconstructed and validated polyploidization histories of Lamiales species. Our results indicated that Primulina huaijiensis underwent three rounds of diploidization events following the γ-WGT event, rather than two rounds as reported. Besides, we reconfirmed that most Lamiales species shared a common diploidization event (L-WGD). Subsequently, we constructed the Lamiales Ancestral Karyotype (LAK), comprising 11 proto-chromosomes, and elucidated its evolutionary trajectory, highlighting the highly flexible reshuffling of the Lamiales paleogenome. We identified biased fractionation of subgenomes following the L-WGD event across eight species, and highlighted the positive impacts of non-WGD genes on gene family expansion. This study provides novel genomic resources and insights into polyploidy and karyotype remodeling of Lamiales species, essential for advancing our understanding of species diversification and genome evolution.

Tomato yellow leaf curl virus: Characteristics, influence, and regulation mechanism

Tomato yellow leaf curl virus (TYLCV), a DNA virus belonging to the genus Begomovirus, significantly impedes the growth and development of numerous host plants, including tomatoes and peppers. Due to its rapid mutation rate and frequent recombination events, achieving complete control of TYLCV proves exceptionally challenging. Consequently, identifying resistance mechanisms become crucial for safeguarding host plants from TYLCV-induced damage. This review article delves into the global distribution, dispersal patterns, and defining characteristics of TYLCV. Moreover, the intricate interplay between TYLCV and various influencing factors, such as insect vectors, susceptible host plants, and abiotic stresses, plays a pivotal role in plant-TYLCV interactions. The review offers an updated perspective on recent investigations focused on plant response mechanisms to TYLCV infection, including the intricate relationship between TYLCV, whiteflies, and regulatory factors. This comprehensive analysis aims to establish a foundation for future research endeavors exploring the molecular mechanisms underlying TYLCV infection and the development of plant resistance through breeding programs.

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.

Doubling down on polyploid discoveries: Global advances in genomics and ecological impacts of polyploidy

All flowering plants are now recognized as diploidized paleopolyploids (Jiao et al., 2011; One Thousand Plant Transcriptomes Initiative, 2019), and polyploid species comprise approximately 30% of contemporary plant species (Wood et al., 2009; Barker et al., 2016a). A major implication of these discoveries is that, to appreciate the evolution of plant diversity, we need to understand the fundamental biology of polyploids and diploidization. This need is broadly recognized by our community as there is a continued, growing interest in polyploidy as a research topic. Over the past 25 years, the sequencing and analysis of plant genomes has revolutionized our understanding of the importance of polyploid speciation to the evolution of land plants.

Species-tree topology impacts the inference of ancient whole-genome duplications across the angiosperm phylogeny

PREMISE

The history of angiosperms is marked by repeated rounds of ancient whole-genome duplications (WGDs). Here we used state-of-the-art methods to provide an up-to-date view of the distribution of WGDs in the history of angiosperms that considers both uncertainty introduced by different WGD inference methods and different underlying species-tree hypotheses.

METHODS

We used the distribution synonymous divergences (Ks) of paralogs and orthologs from transcriptomic and genomic data to infer and place WGDs across two hypothesized angiosperm phylogenies. We further tested these WGD hypotheses with syntenic inferences and Bayesian models of duplicate gene gain and loss.

RESULTS

The predicted number of WGDs in the history of angiosperms (~170) based on the current taxon sampling is largely similar across different inference methods, but varies in the precise placement of WGDs on the phylogeny. Ks-based methods often yield alternative hypothesized WGD placements due to variation in substitution rates among lineages. Phylogenetic models of duplicate gene gain and loss are more robust to topological variation. However, errors in species-tree inference can still produce spurious WGD hypotheses, regardless of method used.

CONCLUSIONS

Here we showed that different WGD inference methods largely agree on an average of 3.5 WGD in the history of individual angiosperm species. However, the precise placement of WGDs on the phylogeny is subject to the WGD inference method and tree topology. As researchers continue to test hypotheses regarding the impacts ancient WGDs have on angiosperm evolution, it is important to consider the uncertainty of the phylogeny as well as WGD inference methods.

Evolution and subfunctionalization of CIPK6 homologous genes in regulating cotton drought resistance

The occurrence of whole-genome duplication or polyploidy may promote plant adaptability to harsh environments. Here, we clarify the evolutionary relationship of eight GhCIPK6 homologous genes in upland cotton (Gossypium hirsutum). Gene expression and interaction analyses indicate that GhCIPK6 homologous genes show significant functional changes after polyploidy. Among these, GhCIPK6D1 and GhCIPK6D3 are significantly up-regulated by drought stress. Functional studies reveal that high GhCIPK6D1 expression promotes cotton drought sensitivity, while GhCIPK6D3 expression promotes drought tolerance, indicating clear functional differentiation. Genetic and biochemical analyses confirm the synergistic negative and positive regulation of cotton drought resistance through GhCBL1A1-GhCIPK6D1 and GhCBL2A1-GhCIPK6D3, respectively, to regulate stomatal movement by controlling the directional flow of K+ in guard cells. These results reveal differentiated roles of GhCIPK6 homologous genes in response to drought stress in upland cotton following polyploidy. The work provides a different perspective for exploring the functionalization and subfunctionalization of duplicated genes in response to polyploidization.

Genome-wide identification of starch phosphorylase gene family in Rosa chinensis and expression in response to abiotic stress

Chinese rose (Rosa chinensis) is an important ornamental plant, with economic, cultural, and symbolic significance. During the application of outdoor greening, adverse environments such as high temperature and drought are often encountered, which affect its application scope and ornamental quality. The starch phosphorylase (Pho) gene family participate in the synthesis and decomposition of starch, not only related to plant energy metabolism, but also plays an important role in plant stress resistance. The role of Pho in combating salinity and high temperature stress in R. chinensis remains unknown. In this work, 4 Phos from R. chinensis were detected with Pfam number of Pho (PF00343.23) and predicted by homolog-based prediction (HBP). The Phos are characterized by sequence lengths of 821 to 997 bp, and the proteins are predicted to subcellularly located in the plastid and cytoplasm. The regulatory regions of the Phos contain abundant stress and phytohormone-responsive cis-acting elements. Based on transcriptome analysis, the Phos were found to respond to abiotic stress factors such as drought, salinity, high temperature, and plant phytohormone of jasmonic acid and salicylic acid. The response of Phos to abiotic stress factors such as salinity and high temperature was confirmed by qRT-PCR analysis. To evaluate the genetic characteristics of Phos, a total of 69 Phos from 17 species were analyzed and then classified into 3 groups in phylogenetic tree. The collinearity analysis of Phos in R. chinensis and other species was conducted for the first time. This work provides a view of evolution for the Pho gene family and indicates that Phos play an important role in abiotic stress response of R. chinensis.

The evolution history of an allotetraploid mangrove tree analysed with a new tool Allo4D

Mangrove species are broadly classified as true mangroves and mangrove associates. The latter are amphibious plants that can survive in the intertidal zone and reproduce naturally in terrestrial environments. Their widespread distribution and extensive adaptability make them ideal research materials for exploring adaptive evolution. In this study, we de novo assembled two genomes of mangrove associates (the allotetraploid Barringtonia racemosa (2n = 4x = 52) and diploid Barringtonia asiatica (2n = 2x = 26)) to investigate the role of allopolyploidy in the evolutionary history of mangrove species. We developed a new allotetraploid-dividing tool Allo4D to distinguish between allotetraploid scaffold-scale subgenomes and verified its accuracy and reliability using real and simulated data. According to the two subgenomes of allotetraploid B. racemosa divided using Allo4D, the allopolyploidization event was estimated to have occurred approximately one million years ago (Mya). We found that B. racemosa, B. asiatica, and Diospyros lotus shared a whole genome duplication (WGD) event during the K-Pg (Cretaceous-Paleozoic) period. K-Pg WGD and recent allopolyploidization events contributed to the speciation of B. racemosa and its adaptation to coastal habitats. We found that genes in the glucosinolates (GSLs) pathway, an essential pathway in response to various biotic and abiotic stresses, expanded rapidly in B. racemosa during polyploidization. In summary, this study provides a typical example of the adaptation of allopolyploid plants to extreme environmental conditions. The newly developed tool, Allo4D, can effectively divide allotetraploid subgenomes and explore the evolutionary history of polyploid plants, especially for species whose ancestors are unknown or extinct.

Genome-wide characterization, transcriptome profiling, and functional analysis of the ALMT gene family in Medicago for aluminum resistance

Aluminum (Al) is the major limiting factor affecting plant productivity in acidic soils. Al3+ ions exhibit increased solubility at a pH below 5, leading to plant root tip toxicity. Alternatively, plants can perceive very low concentrations of Al3+, and Al triggers downstream signaling even at pH 5.7 without causing Al toxicity. The ALUMINUM-ACTIVATED-MALATE-TRANSPORTER (ALMT) family members act as anion channels, with some regulating the secretion of malate from root apices to chelate Al, which is a crucial mechanism for plant Al resistance. To date, the role of the ALMT gene family within the legume Medicago species has not been fully characterized. In this study, we investigated the ALMT gene family in M. sativa and M. truncatula and identified 68 MsALMTs and 18 MtALMTs, respectively. Phylogenetic analysis classified these genes into five clades, and synteny analysis uncovered genuine paralogs and orthologs. The real-time quantitative reverse transcription PCR (qRT-PCR) analysis revealed that MtALMT8, MtALMT9, and MtALMT15 in clade 2-2b are expressed in both roots and root nodules, and MtALMT8 and MtALMT9 are significantly upregulated by Al in root tips. We also observed that MtALMT8 and MtALMT9 can partially restore the Al sensitivity of Atalmt1 in Arabidopsis. Moreover, transcriptome analysis examined the expression patterns of these genes in M. sativa in response to Al at both pH 5.7 and pH 4.6, as well as to protons, and found that Al and protons can independently induce some Al-resistance genes. Overall, our findings indicate that MtALMT8 and MtALMT9 may play a role in Al resistance, and highlight the resemblance between the ALMT genes in Medicago species and those in Arabidopsis.

Integrated anatomical structure, physiological, and transcriptomic analyses to identify differential cold tolerance responses of Ziziphus jujuba mill. 'Yueguang' and its autotetraploid 'Hongguang'

Cold stress is a limiting stress factor that limits plant distribution and development; however, polyploid plants have specific characteristics such as higher resistance to abiotic stress, especially cold stress, that allow them to overcome this challenge. The cultivated cultivar Ziziphus jujuba Mill. 'Yueguang' (YG) and its autotetraploid counterpart 'Hongguang' (HG) exhibit differential cold tolerance. However, the underlying molecular mechanism and methods to enhance their cold tolerance remain unknown. Anatomical structure and physiological analysis indicated YG had a higher wood bark ratio, and xylem ratio under cold treatment compared to HG. However, the half-lethal temperature (LT50), cortex ratio, and malondialdehyde (MDA) content were significantly decreased in YG than HG, which indicated YG was cold tolerant than HG. Transcriptome analysis showed that 2084, 1725, 2888, and 2934 differentially expressed genes (DEGs) were identified in HC vs YC, H20 vs Y20, Y20 vs YC, and H20 vs HC treatment, respectively. Meanwhile, KEGG enrichment analysis of DEGs showed that several metabolic pathways, primarily plant hormone signal transduction and the MAPK signaling pathway, were involved in the differential regulation of cold tolerance between YG and HG. Furthermore, exogenous abscisic acid (ABA) and brassinolide (BR) treatments could improve their cold tolerance through increased SOD and POD activities, decreased relative electrical conductivity, and MDA content. All of these findings suggested that plant hormone signal transduction, particularly ABA and BR, might have an important role in the regulation of differential cold tolerance between YG and HG, laying the foundation for further improving cold tolerance in jujube and examining the molecular mechanisms underlying differences in cold tolerance among different ploidy cultivars.

Chromosome-level Alstonia scholaris genome unveils evolutionary insights into biosynthesis of monoterpenoid indole alkaloids

Alstonia scholaris of the Apocynaceae family is a medicinal plant with a rich source of bioactive monoterpenoid indole alkaloids (MIAs), which possess anti-cancer activity like vinca alkaloids. To gain genomic insights into MIA biosynthesis, we assembled a high-quality chromosome-level genome for A. scholaris using nanopore and Hi-C data. The 444.95 Mb genome contained 35,488 protein-coding genes. A total of 20 chromosomes were assembled with a scaffold N50 of 21.75 Mb. The genome contained a cluster of strictosidine synthases and tryptophan decarboxylases with synteny to other species and a saccharide-terpene cluster involved in the monoterpenoid biosynthesis pathway of the MIA upstream pathway. The multi-omics data of A. scholaris provide a valuable resource for understanding the evolutionary origins of MIAs and for discovering biosynthetic pathways and synthetic biology efforts for producing pharmaceutically useful alkaloids.

Functional Divergence in Orthologous Transcription Factors: Insights from AtCBF2/3/1 and OsDREB1C

Despite traditional beliefs of orthologous genes maintaining similar functions across species, growing evidence points to their potential for functional divergence. C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) are critical in cold acclimation, with their overexpression enhancing stress tolerance but often constraining plant growth. In contrast, a recent study unveiled a distinctive role of rice OsDREB1C in elevating nitrogen use efficiency (NUE), photosynthesis, and grain yield, implying functional divergence within the CBF/DREB1 orthologs across species. Here, we delve into divergent molecular mechanisms of OsDREB1C and AtCBF2/3/1 by exploring their evolutionary trajectories across rice and Arabidopsis genomes, regulatomes, and transcriptomes. Evolutionary scrutiny shows discrete clades for OsDREB1C and AtCBF2/3/1, with the Poaceae-specific DREB1C clade mediated by a transposon event. Genome-wide binding profiles highlight OsDREB1C's preference for GCCGAC compared to AtCBF2/3/1's preference for A/GCCGAC, a distinction determined by R12 in the OsDREB1C AP2/ERF domain. Cross-species multiomic analyses reveal shared gene orthogroups (OGs) and underscore numerous specific OGs uniquely bound and regulated by OsDREB1C, implicated in NUE, photosynthesis, and early flowering, or by AtCBF2/3/1, engaged in hormone and stress responses. This divergence arises from gene gains/losses (∼16.7% to 25.6%) and expression reprogramming (∼62.3% to 66.2%) of OsDREB1C- and AtCBF2/3/1-regulated OGs during the extensive evolution following the rice-Arabidopsis split. Our findings illustrate the regulatory evolution of OsDREB1C and AtCBF2/3/1 at a genomic scale, providing insights on the functional divergence of orthologous transcription factors following gene duplications across species.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

The high-quality genome of Cryptotaenia japonica and comparative genomics analysis reveals anthocyanin biosynthesis in Apiaceae

Cryptotaenia japonica, a traditional medicinal and edible vegetable crops, is well-known for its attractive flavors and health care functions. As a member of the Apiaceae family, the evolutionary trajectory and biological properties of C. japonica are not clearly understood. Here, we first reported a high-quality genome of C. japonica with a total length of 427 Mb and N50 length 50.76 Mb, was anchored into 10 chromosomes, which confirmed by chromosome (cytogenetic) analysis. Comparative genomic analysis revealed C. japonica exhibited low genetic redundancy, contained a higher percentage of single-cope gene families. The homoeologous blocks, Ks, and collinearity were analyzed among Apiaceae species contributed to the evidence that C. japonica lacked recent species-specific WGD. Through comparative genomic and transcriptomic analyses of Apiaceae species, we revealed the genetic basis of the production of anthocyanins. Several structural genes encoding enzymes and transcription factor genes of the anthocyanin biosynthesis pathway in different species were also identified. The CjANSa, CjDFRb, and CjF3H gene might be the target of Cjaponica_2.2062 (bHLH) and Cjaponica_1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae)

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

Genomic analysis reveals phylogeny of Zygophyllales and mechanism for water retention of a succulent xerophyte

Revealing the genetic basis for stress-resistant traits in extremophile plants will yield important information for crop improvement. Zygophyllum xanthoxylum, an extant species of the ancient Mediterranean, is a succulent xerophyte that can maintain a favorable water status under desert habitats; however, the genetic basis of this adaptive trait is poorly understood. Furthermore, the phylogenetic position of Zygophyllales, to which Z. xanthoxylum belongs, remains controversial. In this study, we sequenced and assembled the chromosome-level genome of Z. xanthoxylum. Phylogenetic analysis showed that Zygophyllales and Myrtales form a separated taxon as a sister to the clade comprising fabids and malvids, clarifying the phylogenetic position of Zygophyllales at whole-genome scale. Analysis of genomic and transcriptomic data revealed multiple critical mechanisms underlying the efficient osmotic adjustment using Na+ and K+ as "cheap" osmolytes that Z. xanthoxylum has evolved through the expansion and synchronized expression of genes encoding key transporters/channels and their regulators involved in Na+/K+ uptake, transport, and compartmentation. It is worth noting that ZxCNGC1;1 (cyclic nucleotide-gated channels) and ZxCNGC1;2 constituted a previously undiscovered energy-saving pathway for Na+ uptake. Meanwhile, the core genes involved in biosynthesis of cuticular wax also featured an expansion and upregulated expression, contributing to the water retention capacity of Z. xanthoxylum under desert environments. Overall, these findings boost the understanding of evolutionary relationships of eudicots, illustrate the unique water retention mechanism in the succulent xerophyte that is distinct from glycophyte, and thus provide valuable genetic resources for the improvement of stress tolerance in crops and insights into the remediation of sodic lands.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BACKGROUND

Whole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.

RESULTS

Karyotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.

CONCLUSION

Chromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

Adaptation in Unstable Environments and Global Gene Losses: Small but Stable Gene Networks by the May-Wigner Theory

Although gene loss is common in evolution, it remains unclear whether it is an adaptive process. In a survey of seven major mangrove clades that are woody plants in the intertidal zones of daily environmental perturbations, we noticed that they generally evolved reduced gene numbers. We then focused on the largest clade of Rhizophoreae and observed the continual gene set reduction in each of the eight species. A great majority of gene losses are concentrated on environmental interaction processes, presumably to cope with the constant fluctuations in the tidal environments. Genes of the general processes for woody plants are largely retained. In particular, fewer gene losses are found in physiological traits such as viviparous seeds, high salinity, and high tannin content. Given the broad and continual genome reductions, we propose the May-Wigner theory (MWT) of system stability as a possible mechanism. In MWT, the most effective solution for buffering continual perturbations is to reduce the size of the system (or to weaken the total genic interactions). Mangroves are unique as immovable inhabitants of the compound environments in the land-sea interface, where environmental gradients (such as salinity) fluctuate constantly, often drastically. Extending MWT to gene regulatory network (GRN), computer simulations and transcriptome analyses support the stabilizing effects of smaller gene sets in mangroves vis-à-vis inland plants. In summary, we show the adaptive significance of gene losses in mangrove plants, including the specific role of promoting phenotype innovation and a general role in stabilizing GRN in unstable environments as predicted by MWT.

Ancient Duplication and Lineage-Specific Transposition Determine Evolutionary Trajectory of ERF Subfamily across Angiosperms

AP2/ERF transcription factor family plays an important role in plant development and stress responses. Previous studies have shed light on the evolutionary trajectory of the AP2 and DREB subfamilies. However, knowledge about the evolutionary history of the ERF subfamily in angiosperms still remains limited. In this study, we performed a comprehensive analysis of the ERF subfamily from 107 representative angiosperm species by combining phylogenomic and synteny network approaches. We observed that the expansion of the ERF subfamily was driven not only by whole-genome duplication (WGD) but also by tandem duplication (TD) and transposition duplication events. We also found multiple transposition events in Poaceae, Brassicaceae, Poales, Brassicales, and Commelinids. These events may have had notable impacts on copy number variation and subsequent functional divergence of the ERF subfamily. Moreover, we observed a number of ancient tandem duplications occurred in the ERF subfamily across angiosperms, e.g., in Subgroup IX, IXb originated from ancient tandem duplication events within IXa. These findings together provide novel insights into the evolution of this important transcription factor family.

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.