The water lily genome and the early evolution of flowering plants

Genome assembly and resequencing shed light on evolution, population selection, and sex identification in Vernicia montana

Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69 Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933 bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree)

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

Multi-omic dissection of the cold resistance traits of white water lily

The white water lily (Nymphaea candida), exemplifying nature's resilience, thrives in the high-altitude terrains of Xinjiang, China, serving as an ideal model for investigating cold adaptation mechanisms in aquatic plants. This study meticulously elucidates the complex cold adaptation mechanisms of the white water lily through a comprehensive and integrated methodological approach. We discovered that the water lily undergoes ecodormancy in winter, retaining high cellular viability and growth potential. During overwintering, the white water lily demonstrates effective resource reallocation, a process facilitated by morphological adjustments, thereby strengthening its resistance to cold temperatures. This enhancement is achieved particularly through the compartmentalization of large vacuoles, the accumulation of osmoregulatory substances, and an increased antioxidant capacity. We established the first exhaustive full-length transcriptome for the white water lily. A subsequent comprehensive analysis of the transcriptome, phytohormones, and metabolome uncovered a multifaceted regulatory network orchestrating cold adaptation. Our research spotlights phytohormone signaling, amino acid metabolism, and circadian rhythms as key elements in the water lily's defense against cold. The results emphasize the critical role of nitrogen metabolism, especially amino acid-related pathways, during cold stress. Metabolite profiling revealed the importance of compounds like myo-inositol and L-proline in enhancing cold tolerance. Remarkably, our study demonstrates that the white water lily notably diminishes the utilization of unsaturated fatty acids in its temperature regulation strategies. In conclusion, this research substantially enriches our understanding of the white water lily's intricate cold adaptation mechanisms, offering new perspectives on the adaptive strategies of aquatic plants and potential applications in agricultural advancement.

Climatic niche evolution and niche conservatism of Nymphaea species in Africa, South America, and Australia

BACKGROUND

Interest in the evolution of climatic niches, particularly in understanding the potential adaptive responses of species under climate change, has increased both theoretically and within macroecological studies. These studies have provided valuable insights into how climatic traits of species influence their niche evolution. In this study, we aim to investigate whether niche conservatism plays a role in the species diversification of Nymphaea, a group of aquatic plants with a cosmopolitan distribution that is facing severe habitat loss. We will use climatic models and phylogenetic data for 23 species to reconstruct Nymphaea's niche evolution, measure niche overlap, and assess disparity through time while testing for evolutionary models.

RESULTS

There was a lot of overlap in niches both within and between groups, especially for species that can be found in many places. The breadth and peaks of the niche profile varied depending on the bioclimatic variables, which suggested that the species evolved differently to cope with changes in climate. The analysis also showed that evolutionary changes happened across the phylogeny, with weak to moderate signals. The morphological disparity index (MDI) values indicated that there were disparities within subclades over time but not between or among them. Niche reconstruction and evolution analysis revealed both convergent and divergent evolution among various variables. For example, N. immutabilis, N. atrans, N. violancea, and N. nouchali evolved towards intermediate temperatures for bio2 and bio3 (isothermity) while moving towards extreme temperatures for bio8 and bio9 (wettest and driest average quarterly temperatures).

CONCLUSION

Our study will improve our understanding of how changes in climatic niches are potentially driving the evolution of Nymphaea. It has significant scientific implications for the limits, assemblages, evolution, and diversification of species. This information is crucial for the ongoing efforts of conservation and management, particularly considering the inevitable effects of climate change.

A high-quality haplotype genome of Michelia alba DC reveals differences in methylation patterns and flower characteristics

Michelia alba DC is a highly valuable ornamental plant of the Magnoliaceae family. This evergreen tropical tree commonly grows in Southeast Asia and is adored for its delightful fragrance. Our study assembled the M. alba haplotype genome MC and MM by utilizing Nanopore ultralong reads, Pacbio Hifi long reads and parental second-generation data. Moreover, the first methylation map of Magnoliaceae was constructed based on the methylation site data obtained using Nanopore data. Metabolomic datasets were generated from the flowers of three different species to assess variations in pigment and volatile compound accumulation. Finally, transcriptome data were generated to link genomic, methylation, and morphological patterns to reveal the reasons underlying the differences between M. alba and its parental lines in petal color, flower shape, and fragrance. We found that the AP1 and AP2 genes are crucial in M. alba petal formation, while the 4CL, PAL, and C4H genes control petal color. The data generated in this study serve as a foundation for future physiological and biochemical research on M. alba, facilitate the targeted improvement of M. alba varieties, and offer a theoretical basis for molecular research on Michelia L.

On the importance of sequence alignment inspections in plastid phylogenomics - an example from revisiting the relationships of the water-lilies

The water-lily clade represents the second earliest-diverging branch of angiosperms. Most of its species belong to Nymphaeaceae, of which the "core Nymphaeaceae"-comprising the genera Euryale, Nymphaea and Victoria-is the most diverse clade. Despite previous molecular phylogenetic studies on the core Nymphaeaceae, various aspects of their evolutionary relationships have remained unresolved. The length-variable introns and intergenic spacers are known to contain most of the sequence variability within the water-lily plastomes. Despite the challenges with multiple sequence alignment, any new molecular phylogenetic investigation on the core Nymphaeaceae should focus on these noncoding plastome regions. For example, a new plastid phylogenomic study on the core Nymphaeaceae should generate DNA sequence alignments of all plastid introns and intergenic spacers based on the principle of conserved sequence motifs. In this investigation, we revisit the phylogenetic history of the core Nymphaeaceae by employing such an approach. Specifically, we use a plastid phylogenomic analysis strategy in which all coding and noncoding partitions are separated and then undergo software-driven DNA sequence alignment, followed by a motif-based alignment inspection and adjustment. This approach allows us to increase the reliability of the character base compared to the default practice of aligning complete plastomes through software algorithms alone. Our approach produces significantly different phylogenetic tree reconstructions for several of the plastome regions under study. The results of these reconstructions underscore that Nymphaea is paraphyletic in its current circumscription, that each of the five subgenera of Nymphaea is monophyletic, and that the subgenus Nymphaea is sister to all other subgenera of Nymphaea. Our results also clarify many evolutionary relationships within the Nymphaea subgenera Brachyceras, Hydrocallis and Nymphaea. In closing, we discuss whether the phylogenetic reconstructions obtained through our motif-based alignment adjustments are in line with morphological evidence on water-lily evolution.

Genome-Wide Identification and Expression Analysis of Bx Involved in Benzoxazinoids Biosynthesis Revealed the Roles of DIMBOA during Early Somatic Embryogenesis in Dimocarpus longan Lour

Benzoxazinoids (BXs) are tryptophan-derived indole metabolites and play a role in various physiological processes, such as auxin metabolism. Auxin is essential in the process of somatic embryogenesis (SE) in plants. In this study, we used bioinformatics, transcriptome data, exogenous treatment experiments, and qPCR analysis to study the evolutionary pattern of Bx genes in green plants, the regulatory mechanism of DlBx genes during early SE, and the effect of 2,4-dihydroxy-7-methoxy-1,4-benzoxazine-3-one (DIMBOA) on the early SE in Dimocarpus longan Lour. The results showed that 27 putative DlBxs were identified in the longan genome; the Bx genes evolved independently in monocots and dicots, and the main way of gene duplication for the DlBx was tandem duplication (TD) and the DlBx were strongly constrained by purification selection during evolution. The transcriptome data indicated varying expression levels of DlBx during longan early SE, and most DlBxs responded to light, temperature, drought stress, and 2,4-dichlorophenoxyacetic acid (2,4-D) treatment; qRT-PCR results showed DlBx1, DlBx6g and DlBx6h were responsive to auxin, and treatment with 0.1mg/L DIMBOA for 9 days significantly upregulated the expression levels of DlBx1, DlBx3g, DlBx6c, DlBx6f, DlB6h, DlBx7d, DlBx8, and DlBx9b. The correlation analysis showed a significantly negative correlation between the expression level of DlBx1 and the endogenous IAA contents; DIMBOA significantly promoted the early SE and significantly changed the endogenous IAA content, and the IAA content increased significantly at the 9th day and decreased significantly at the 13th day. Therefore, the results suggested that DIMBOA indirectly promote the early SE by changing the endogenous IAA content via affecting the expression level of DlBx1 and hydrogen peroxide (H2O2) content in longan.

Integrative analysis of transcriptome and target metabolites uncovering flavonoid biosynthesis regulation of changing petal colors in Nymphaea 'Feitian 2'

BACKGROUND

Nymphaea (waterlily) is known for its rich colors and role as an important aquatic ornamental plant globally. Nymphaea atrans and some hybrids, including N. 'Feitian 2,' are more appealing due to the gradual color change of their petals at different flower developmental stages. The petals of N. 'Feitian 2' gradually change color from light blue-purple to deep rose-red throughout flowering. The mechanism of the phenomenon remains unclear.

RESULTS

In this work, flavonoids in the petals of N. 'Feitian 2' at six flowering stages were examined to identify the influence of flavonoid components on flower color changes. Additionally, six cDNA libraries of N. 'Feitian 2' over two blooming stages were developed, and the transcriptome was sequenced to identify the molecular mechanism governing petal color changes. As a result, 18 flavonoid metabolites were identified, including five anthocyanins and 13 flavonols. Anthocyanin accumulation during flower development is the primary driver of petal color change. A total of 12 differentially expressed genes (DEGs) in the flavonoid biosynthesis pathway were uncovered, and these DEGs were significantly positively correlated with anthocyanin accumulation. Six structural genes were ultimately focused on, as their expression levels varied significantly across different flowering stages. Moreover, 104 differentially expressed transcription factors (TFs) were uncovered, and three MYBs associated with flavonoid biosynthesis were screened. The RT-qPCR results were generally aligned with high-throughput sequencing results.

CONCLUSIONS

This research offers a foundation to clarify the mechanisms underlying changes in the petal color of waterlilies.

Chromosome-level genome assembly of the threatened resource plant Cinnamomum chago

Cinnamomum chago is a tree species endemic to Yunnan province, China, with potential economic value, phylogenetic importance, and conservation priority. We assembled the genome of C. chago using multiple sequencing technologies, resulting in a high-quality, chromosomal-level genome with annotation information. The assembled genome size is approximately 1.06 Gb, with a contig N50 length of 92.10 Mb. About 99.92% of the assembled sequences could be anchored to 12 pseudo-chromosomes, with only one gap, and 63.73% of the assembled genome consists of repeat sequences. In total, 30,497 genes were recognized according to annotation, including 28,681 protein-coding genes. This high-quality chromosome-level assembly and annotation of C. chago will assist us in the conservation and utilization of this valuable resource, while also providing crucial data for studying the evolutionary relationships within the Cinnamomum genus, offering opportunities for further research and exploration of its diverse applications.

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.

Broad Chain-Length Specificity of the Alkane-Forming Enzymes NoCER1A and NoCER3A/B in Nymphaea odorata

Many terrestrial plants produce large quantities of alkanes for use in epicuticular wax and the pollen coat. However, their carbon chains must be long to be useful as fuel or as a petrochemical feedstock. Here, we focus on Nymphaea odorata, which produces relatively short alkanes in its anthers. We identified orthologs of the Arabidopsis alkane biosynthesis genes AtCER1 and AtCER3 in N. odorata and designated them NoCER1A, NoCER3A and NoCER3B. Expression analysis of NoCER1A and NoCER3A/B in Arabidopsis cer mutants revealed that the N. odorata enzymes cooperated with the Arabidopsis enzymes and that the NoCER1A produced shorter alkanes than AtCER1, regardless of which CER3 protein it interacted with. These results indicate that AtCER1 frequently uses a C30 substrate, whereas NoCER1A, NoCER3A/B and AtCER3 react with a broad range of substrate chain lengths. The incorporation of shorter alkanes disturbed the formation of wax crystals required for water-repellent activity in stems, suggesting that chain-length specificity is important for surface cleaning. Moreover, cultured tobacco cells expressing NoCER1A and NoCER3A/B effectively produced C19-C23 alkanes, indicating that the introduction of the two enzymes is sufficient to produce alkanes. Taken together, our findings suggest that these N. odorata enzymes may be useful for the biological production of alkanes of specific lengths. 3D modeling revealed that CER1s and CER3s share a similar structure that consists of N- and C-terminal domains, in which their predicted active sites are respectively located. We predicted the complex structure of both enzymes and found a cavity that connects their active sites.

Fossil-calibrated molecular clock data enable reconstruction of steps leading to differentiated multicellularity and anisogamy in the Volvocine algae

BACKGROUND

Throughout its nearly four-billion-year history, life has undergone evolutionary transitions in which simpler subunits have become integrated to form a more complex whole. Many of these transitions opened the door to innovations that resulted in increased biodiversity and/or organismal efficiency. The evolution of multicellularity from unicellular forms represents one such transition, one that paved the way for cellular differentiation, including differentiation of male and female gametes. A useful model for studying the evolution of multicellularity and cellular differentiation is the volvocine algae, a clade of freshwater green algae whose members range from unicellular to colonial, from undifferentiated to completely differentiated, and whose gamete types can be isogamous, anisogamous, or oogamous. To better understand how multicellularity, differentiation, and gametes evolved in this group, we used comparative genomics and fossil data to establish a geologically calibrated roadmap of when these innovations occurred.

RESULTS

Our ancestral-state reconstructions, show that multicellularity arose independently twice in the volvocine algae. Our chronograms indicate multicellularity evolved during the Carboniferous-Triassic periods in Goniaceae + Volvocaceae, and possibly as early as the Cretaceous in Tetrabaenaceae. Using divergence time estimates we inferred when, and in what order, specific developmental changes occurred that led to differentiated multicellularity and oogamy. We find that in the volvocine algae the temporal sequence of developmental changes leading to differentiated multicellularity is much as proposed by David Kirk, and that multicellularity is correlated with the acquisition of anisogamy and oogamy. Lastly, morphological, molecular, and divergence time data suggest the possibility of cryptic species in Tetrabaenaceae.

CONCLUSIONS

Large molecular datasets and robust phylogenetic methods are bringing the evolutionary history of the volvocine algae more sharply into focus. Mounting evidence suggests that extant species in this group are the result of two independent origins of multicellularity and multiple independent origins of cell differentiation. Also, the origin of the Tetrabaenaceae-Goniaceae-Volvocaceae clade may be much older than previously thought. Finally, the possibility of cryptic species in the Tetrabaenaceae provides an exciting opportunity to study the recent divergence of lineages adapted to live in very different thermal environments.

Flowering genes identification, network analysis, and database construction for 837 plants

Flowering is one of the most important biological phenomena in the plant kingdom, which not only has important ecological significance, but also has substantial horticultural ornamental value. In this study, we undertook an exhaustive review of the advancements in our understanding of plant flowering genes. We delved into the identification and conducted comparative analyses of flowering genes across virtually all sequenced angiosperm plant genomes. Furthermore, we established an extensive angiosperm flowering atlas, encompassing a staggering 183 720 genes across eight pathways, along with 10 155 ABCDE mode genes, which play a pivotal role in plant flowering regulation. Through the examination of expression patterns, we unveiled the specificities of these flowering genes. An interaction network between flowering genes of the ABCDE model and their corresponding upstream genes offered a blueprint for comprehending their regulatory mechanisms. Moreover, we predicted the miRNA and target genes linked to the flowering processes of each species. To culminate our efforts, we have built a user-friendly web interface, named the Plant Flowering-time Gene Database (PFGD), accessible at http://pfgd.bio2db.com/. We firmly believe that this database will serve as a cornerstone in the global research community, facilitating the in-depth exploration of flowering genes in the plant kingdom. In summation, this pioneering endeavor represents the first comprehensive collection and comparative analysis of flowering genes in plants, offering valuable resources for the study of plant flowering genetics.

The miR166 targets CsHDZ3 genes to negatively regulate drought tolerance in tea plant (Camellia sinensis)

Drought is the stressor with a significant adverse impact on the yield stability of tea plants. HD-ZIP III transcription factors (TFs) play important regulatory roles in plant growth, development, and stress responses. However, whether and how HD-ZIP III TFs are involved in drought response and tolerance in tea plants remains unclear. Here, we identified seven HD-ZIP III genes (CsHDZ3-1 to CsHDZ3-7) in tea plant genome. The evolutionary analysis demonstrated that CsHDZ3 members were subjected to purify selection. Subcellular localization analysis revealed that all seven CsHDZ3s located in the nucleus. Yeast self-activation and dual-luciferase reporter assays demonstrated that CsHDZ3-1 to CsHDZ3-4 have trans-activation ability whereas CsHDZ3-5 to CsHDZ3-7 served as transcriptional inhibitors. The qRT-PCR assay showed that all seven CsHDZ3 genes could respond to simulated natural drought stress and polyethylene glycol treatment. Further assays verified that all CsHDZ3 genes can be cleaved by csn-miR166. Overexpression of csn-miR166 inhibited the expression of seven CsHDZ3 genes and weakened drought tolerance of tea leaves. In contrast, suppression of csn-miR166 promoted the expression of seven CsHDZ3 genes and enhanced drought tolerance of tea leaves. These findings established the foundation for further understanding the mechanism of CsHDZ3-miR166 modules' participation in drought responses and tolerance.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

Molecular and phylogenetic evidence of parallel expansion of anion channels in plants

Aluminum-activated malate transporters (ALMTs) and slow anion channels (SLACs) are important in various physiological processes in plants, including stomatal regulation, nutrient uptake, and in response to abiotic stress such as aluminum toxicity. To understand their evolutionary history and functional divergence, we conducted phylogenetic and expression analyses of ALMTs and SLACs in green plants. Our findings from phylogenetic studies indicate that ALMTs and SLACs may have originated from green algae and red algae, respectively. The ALMTs of early land plants and charophytes formed a monophyletic clade consisting of three subgroups. A single duplication event of ALMTs was identified in vascular plants and subsequent duplications into six clades occurred in angiosperms, including an identified clade, 1-1. The ALMTs experienced gene number losses in clades 1-1 and 2-1 and expansions in clades 1-2 and 2-2b. Interestingly, the expansion of clade 1-2 was also associated with higher expression levels compared to genes in clades that experienced apparent loss. SLACs first diversified in bryophytes, followed by duplication in vascular plants, giving rise to three distinct clades (I, II, and III), and clade II potentially associated with stomatal control in seed plants. SLACs show losses in clades II and III without substantial expansion in clade I. Additionally, ALMT clade 2-2 and SLAC clade III contain genes specifically expressed in reproductive organs and roots in angiosperms, lycophytes, and mosses, indicating neofunctionalization. In summary, our study demonstrates the evolutionary complexity of ALMTs and SLACs, highlighting their crucial role in the adaptation and diversification of vascular plants.

Variation in the Floral Scent Chemistry of Nymphaea 'Eldorado', a Valuable Water Lily, with Different Flowering Stages and Flower Parts

Nymphaea 'Eldorado', a valuable water lily, is a well-known fragrant plant in China. Studying the temporal and spatial characteristics of the floral components of this plant can provide a reference for the further development and utilization of water lily germplasm resources. In this study, headspace solid-phase microextraction (HS-SPME) combined with gas chromatography-mass spectrometry (GC-MS) was used to explore the types and relative contents of floral components at different flowering stages (S1: bud stage; S2: initial-flowering stage; S3: full-flowering stage; S4: end-flowering stage) and in different floral organs of N. 'Elidorado', combined with the observation of the microscopic structure of petals. A total of 60 volatile organic compounds (VOCs) were detected at different flowering stages, and there were significant differences in floral VOCs at different flowering stages and in different flower organs. The volatile compounds of N. 'Eldorado' can be divided into seven chemical classes,, namely, alkenes, alcohols, esters, aldehydes, ketones, alkanes, and others; the most common were alkenes and alkanes. A total of 39, 44, 47, and 42 volatile compounds were detected at S1, S2, S3, and S4. The VOCs present in high concentrations include benzaldehyde, benzyl alcohol, benzyl acetate, trans-α-bergamotene, α-curcumene, cis-α-farnesene, and so on. The types and total contents of volatiles at the full-flowering stage were higher than at other flowering stages. Comparing the VOCs in different parts of flower organs, it was found that the contents of alcohols, esters, and aldehydes were greatest in the petals, the alkenes in stamens were abundant with a relative content of up to 54.93%, and alkanes in the pistil were higher than in other parts. The types and total contents of volatiles in the stamens of N. 'Eldorado' were higher than those in other flower organs; they were the main part releasing fragrance. The observation of petal microstructure revealed that the size and quantity of the papillae on the epidermises of petals, the number of intracellular plastids, and the aggregates of floral components (osmophilic matrix granules) were significantly higher at the full-flowering stage than at the other flowering stages. This study suggested the main flowering stage and location at which the floral VOCs are released by N. 'Eldorado' and provided a reference for guiding the breeding of this water lily, exploring genetic patterns and developing related products.

NADPH oxidase-mediated reactive oxygen species, antioxidant isozymes, and redox homeostasis regulate salt sensitivity in maize genotypes

The aim of the study is to examine the relationship between oxidative bursts, their regulation with ion homeostasis, and NADPH oxidase (NOX) in different salt-sensitive maize genotypes. For this, in the first study, four differently salt-sensitive maize genotypes (BIL214 × BIL218 as tolerant, BHM-5 as sensitive, and BHM-7 and BHM-9 as moderate-tolerant) were selected on the basis of phenotype, histochemical detection of reactive oxygen species (ROS), malondialdehyde (MDA) content, and specific and in-gel activity of NOX. In the next experiment, these genotypes were further examined in 200 mM NaCl solution in half-strength Hoagland media for nine days to study salt-induced changes in NOX activity, ROS accumulation, ion and redox homeostasis, the activity of antioxidants and their isozyme responses, and to find out potential relationships among the traits. Methylglyoxal (MG) and glyoxalse enzymes (Gly I and II) were also evaluated. Fully expanded leaf samplings were collected at 0 (control), 3, 6, 9-day, and after 7 days of recovery to assay different parameters. Na+/K+, NOX, ROS, and MDA contents increased significantly with the progression of stress duration in all maize genotypes, with a significantly higher value in BHM-5 as compared to tolerant and moderate-tolerant genotypes. A continual induction of Cu/Zn-SOD was observed in BIL214 × BIL218 due to salt stress. Substantial decreases in CAT2 and CAT3 isozymes in BHM-5 might be critical for the highest H2O2 burst in that sensitive genotype under salt stress. The highest intensified POD isozymes were visualized in BHM-5, BHM-7, and BHM-9, whereas BIL214 × BIL218 showed a continual induction of POD isozymes, although GPX activity decreased in all the genotypes at 9 days. Under salt stress, the tolerant genotype BIL214 × BIL218 showed superior ASA- and GSH-redox homeostasis by keeping GR and MDHAR activity high. This genotype also had a stronger MG detoxification system by having higher glyoxalase activity. Correlation, comparative heatmap, and PCA analyses revealed positive correlations among Na+/K+, NOX, O2•-, H2O2, MG, proline, GR, GST, and Gly I activities. Importantly, the relationship depends on the salt sensitivity of the genotypes. The reduced CAT activity as well as redox homeostasis were critical to the survival of the sensitive genotype.

Comprehensive Evolutionary Analysis of the SMXL Gene Family in Rosaceae: Further Insights into Its Origin, Expansion, Diversification, and Role in Regulating Pear Branching

SMXL genes constitute a conserved gene family that is ubiquitous in angiosperms and involved in regulating various plant processes, including branching, leaf elongation, and anthocyanin biosynthesis, but little is known about their molecular functions in pear branching. Here, we performed genome-wide identification and investigation of the SMXL genes in 16 angiosperms and analyzed their phylogenetics, structural features, conserved motifs, and expression patterns. In total, 121 SMXLs genes were identified and were classified into four groups. The number of non-redundant SMXL genes in each species varied from 3 (Amborella trichopoda Baill.) to 18 (Glycine max Merr.) and revealed clear gene expansion events over evolutionary history. All the SMXL genes showed conserved structures, containing no more than two introns. Three-dimensional protein structure prediction revealed distinct structures between but similar structures within groups. A quantitative real-time PCR analysis revealed different expressions of 10 SMXL genes from pear branching induced by fruit-thinning treatment. Overall, our study provides a comprehensive investigation of SMXL genes in the Rosaceae family, especially pear. The results offer a reference for understanding the evolutionary history of SMXL genes and provide excellent candidates for studying fruit tree branching regulation, and in facilitating pear pruning and planting strategies.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

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

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 Analysis of APC E3 Ubiquitin Ligase Genes Family in Triticum aestivum

E3 ubiquitin ligases play a pivotal role in ubiquitination, a crucial post-translational modification process. Anaphase-promoting complex (APC), a large cullin-RING E3 ubiquitin ligase, regulates the unidirectional progression of the cell cycle by ubiquitinating specific target proteins and triggering plant immune responses. Several E3 ubiquitin ligases have been identified owing to advancements in sequencing and annotation of the wheat genome. However, the types and functions of APC E3 ubiquitin ligases in wheat have not been reported. This study identified 14 members of the APC gene family in the wheat genome and divided them into three subgroups (CCS52B, CCS52A, and CDC20) to better understand their functions. Promoter sequence analysis revealed the presence of several cis-acting elements related to hormone and stress responses in the APC E3 ubiquitin ligases in wheat. All identified APC E3 ubiquitin ligase family members were highly expressed in the leaves, and the expression of most genes was induced by the application of methyl jasmonate (MeJA). In addition, the APC gene family in wheat may play a role in plant defense mechanisms. This study comprehensively analyzes APC genes in wheat, laying the groundwork for future research on the function of APC genes in response to viral infections and expanding our understanding of wheat immunity mechanisms.

DataColor: unveiling biological data relationships through distinctive color mapping

In the era of rapid advancements in high-throughput omics technologies, the visualization of diverse data types with varying orders of magnitude presents a pressing challenge. To bridge this gap, we introduce DataColor, an all-encompassing software solution meticulously crafted to address this challenge. Our aim is to empower users with the ability to handle a wide array of data types through an assortment of tools, while simultaneously streamlining parameter selection for rapid insights and detailed enhancements. DataColor stands as a robust toolkit, encompassing 23 distinct tools coupled with over 600 parameters. The defining characteristic of this toolkit is its adept utilization of the color spectrum, allowing for the representation of data spanning diverse types and magnitudes. Through the integration of advanced algorithms encompassing data clustering, normalization, squarified layouts, and customizable parameters, DataColor unveils an abundance of insights that lay hidden within the intricate relationships embedded in the data. Whether you find yourself navigating the analysis of expansive datasets or embarking on the quest to visualize intricate patterns, DataColor stands as the comprehensive and potent solution. We extend the availability of DataColor to all users at no cost, accessible through the following link: https://github.com/frankgenome/DataColor.

Chromosomal-scale genomes of two Rosa species provide insights into genome evolution and ascorbate accumulation

Rosa roxburghii and Rosa sterilis, two species belonging to the Rosaceae family, are widespread in the southwest of China. These species have gained recognition for their remarkable abundance of ascorbate in their fresh fruits, making them an ideal vitamin C resource. In this study, we generated two high-quality chromosome-scale genome assemblies for R. roxburghii and R. sterilis, with genome sizes of 504 and 981.2 Mb, respectively. Notably, we present a haplotype-resolved, chromosome-scale assembly for diploid R. sterilis. Our results indicated that R. sterilis originated from the hybridization of R. roxburghii and R. longicuspis. Genome analysis revealed the absence of recent whole-genome duplications in both species and identified a series of duplicated genes that possibly contributing to the accumulation of flavonoids. We identified two genes in the ascorbate synthesis pathway, GGP and GalLDH, that show signs of positive selection, along with high expression levels of GDP-d-mannose 3', 5'-epimerase (GME) and GDP-l-galactose phosphorylase (GGP) during fruit development. Furthermore, through co-expression network analysis, we identified key hub genes (MYB5 and bZIP) that likely regulate genes in the ascorbate synthesis pathway, promoting ascorbate biosynthesis. Additionally, we observed the expansion of terpene synthase genes in these two species and tissue expression patterns, suggesting their involvement in terpenoid biosynthesis. Our research provides valuable insights into genome evolution and the molecular basis of the high concentration of ascorbate in these two Rosa species.

Reciprocal expression of MADS-box genes and DNA methylation reconfiguration initiate bisexual cones in spruce

The naturally occurring bisexual cone of gymnosperms has long been considered a possible intermediate stage in the origin of flowers, but the mechanisms governing bisexual cone formation remain largely elusive. Here, we employed transcriptomic and DNA methylomic analyses, together with hormone measurement, to investigate the molecular mechanisms underlying bisexual cone development in the conifer Picea crassifolia. Our study reveals a "bisexual" expression profile in bisexual cones, especially in expression patterns of B-class, C-class and LEAFY genes, supporting the out of male model. GGM7 could be essential for initiating bisexual cones. DNA methylation reconfiguration in bisexual cones affects the expression of key genes in cone development, including PcDAL12, PcDAL10, PcNEEDLY, and PcHDG5. Auxin likely plays an important role in the development of female structures of bisexual cones. This study unveils the potential mechanisms responsible for bisexual cone formation in conifers and may shed light on the evolution of bisexuality.

Genomic evolution and insights into agronomic trait innovations of Sesamum species

Sesame is an ancient oilseed crop with high oil content and quality. However, the evolutionary history and genetic mechanisms of its valuable agronomic traits remain unclear. Here, we report chromosome-scale genomes of cultivated sesame (Sesamum indicum L.) and six wild Sesamum species, representing all three karyotypes within this genus. Karyotyping and genome-based phylogenic analysis revealed the evolutionary route of Sesamum species from n = 13 to n = 16 and revealed that allotetraploidization occurred in the wild species Sesamum radiatum. Early divergence of the Sesamum genus (48.5-19.7 million years ago) during the Tertiary period and its ancient phylogenic position within eudicots were observed. Pan-genome analysis revealed 9164 core gene families in the 7 Sesamum species. These families are significantly enriched in various metabolic pathways, including fatty acid (FA) metabolism and FA biosynthesis. Structural variations in SiPT1 and SiDT1 within the phosphatidyl ethanolamine-binding protein gene family lead to the genomic evolution of plant-architecture and inflorescence-development phenotypes in Sesamum. A genome-wide association study (GWAS) of an interspecific population and genome comparisons revealed a long terminal repeat insertion and a sequence deletion in DIR genes of wild Sesamum angustifolium and cultivated sesame, respectively; both variations independently cause high susceptibility to Fusarium wilt disease. A GWAS of 560 sesame accessions combined with an overexpression study confirmed that the NAC1 and PPO genes play an important role in upregulating oil content of sesame. Our study provides high-quality genomic resources for cultivated and wild Sesamum species and insights that can improve molecular breeding strategies for sesame and other oilseed crops.

Stomata damage, photosynthesis, and transpiration evaluation of aquatic lirium after ultrasound irradiation

PURPOSE

Evaluate the structural damage and the changes in the photosynthesis and transpiration rates of aquatic lirium leaves caused by ultrasound (US) irradiation in search of environmentally friendly methodologies for the control of this weed.

MATERIALS AND METHODS

Aquatic lirium plants were extracted from Xochimilco water canals in Mexico City. A part of the group of plants was selected for irradiation, and the rest formed the control group. The irradiation plants group was exposed to US irradiation of 17 kHz frequency and 30 W × 4 output power for 2 h, at noon and 25 °C room temperature. The structural analysis was done with a MOTICAM 1 digital camera, 800 × 600 pixels, incorporated into the MOTIC PSM-1000 optical microscope and edited with Motic Images Plus 2.0 ML software. The total stomata density and the damaged stomata density were determined by dividing the numbers of total and damaged stomata by the visual field area (67,917 mm2), respectively. The leaves' photosynthesis and transpiration rates were measured using an LI-6400XT Portable Photosynthesis System.

RESULTS

Significant damage was observed in the stomata and epidermal cells, finding that the average ratio between the damaged and total stomata densities as a function of time (days) showed an exponential increase described by a Box-Lucas equation with a saturation value near unity and a maximum rate of change of the density of damaged stomata on zero-day (immediately after irradiation), decreasing as the days go by. The transpiration rate showed a sudden increase during the first hour after irradiation, reaching a maximum of 36% of its value before irradiation. It then quickly fell during the next 6 days and more slowly until the 21st day, decreasing 79.9% of its value before irradiation. The photosynthetic rate showed similar behavior with a 37.7% maximum increment and a 73.6% minimum decrease of its value before irradiation.

CONCLUSIONS

The results of structural stomata damage on the ultrasound-irradiated aquatic lirium leaves are consistent with an excessive ultrasound stimulation on stomata's mechanical operation by guard cells that produce the measured significant increase of the photosynthetic and transpiration rates during the first hour after irradiation. The initial high evaporation could alter the water potential gradient, with a possible generation of tensions in the xylem that could cause embolism in their conduits. The loss of xylem conductivity or hydraulic failure would be consistent with the observed significant fall in the photosynthesis and transpiration rates of the aquatic lirium leaves after its sudden rise in the first hour after irradiation.

Divergent Retention of Sucrose Metabolism Genes after Whole Genome Triplication in the Tomato (Solanum lycopersicum)

Sucrose, the primary carbon transport mode and vital carbohydrate for higher plants, significantly impacts plant growth, development, yield, and quality formation. Its metabolism involves three key steps: synthesis, transport, and degradation. Two genome triplication events have occurred in Solanaceae, which have resulted in massive gene loss. In this study, a total of 48 and 65 genes from seven sucrose metabolism gene families in Vitis vinifera and Solanum lycopersicum were identified, respectively. The number of members comprising the different gene families varied widely. And there were significant variations in the pattern of gene duplication and loss in the tomato following two WGD events. Tandem duplication is a major factor in the expansion of the SWEET and Acid INV gene families. All the genes are irregularly distributed on the chromosomes, with the majority of the genes showing collinearity with the grape, particularly the CIN family. And the seven gene families were subjected to a purifying selection. The expression patterns of the different gene families exhibited notable variations. This study presents basic information about the sucrose metabolism genes in the tomato and grape, and paves the way for further investigations into the impact of SCT events on the phylogeny, gene retention duplication, and function of sucrose metabolism gene families in the tomato or Solanaceae, and the adaptive evolution of the tomato.

Diverse O-methyltransferases catalyze the biosynthesis of floral benzenoids that repel aphids from the flowers of waterlily Nymphaea prolifera

Nymphaea is a key genus of the ANA grade (Amborellales, Nymphaeales, and Austrobaileyales) of basal flowering plants, which serve as a key model to study the early evolution of floral traits. In this study, we comprehensively investigated the emission, biosynthesis, and biological function of the floral scent in a night-blossoming waterlily Nymphaea prolifera. The headspace volatile collection combined with GC-MS analysis showed that the floral scent of N. prolifera is predominately comprised by methylated benzenoids including anisole, veratrole, guaiacol, and methoxyanisole. Moreover, the emission of these floral benzenoids in N. prolifera exhibited temporal and spatial pattern with circadian rhythm and tissue specificity. By creating and mining transcriptomes of N. prolifera flowers, 12 oxygen methyltransferases (NpOMTs) were functionally identified. By in vitro enzymatic assay, NpOMT3, 6, and 7 could produce anisole and NpOMT5, 7, 9, produce guaiacol, whereas NpOMT3, 6, 9, 11 catalyzed the formation of veratrole. Methoxyanisole was identified as the universal product of all NpOMTs. Expression patterns of NpOMTs provided implication for their roles in the production of the respective benzenoids. Phylogenetic analysis of OMTs suggested a Nymphaea-specific expansion of the OMT family, indicating the evolution of lineage-specific functions. In bioassays, anisole, veratrole, and guaiacol in the floral benzenoids were revealed to play the critical role in repelling waterlily aphids. Overall, this study indicates that the basal flowering plant N. prolifera has evolved a diversity and complexity of OMT genes for the biosynthesis of methylated benzenoids that can repel insects from feeding the flowers. These findings provide new insights into the evolutional mechanism and ecological significance of the floral scent from early-diverged flowering plants.

Transcription factors KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4 regulate long intergenic noncoding RNAs expressed in Arabidopsis roots

Thousands of long intergenic noncoding RNAs (lincRNAs) have been identified in plant genomes. While some lincRNAs have been characterized as important regulators in different biological processes, little is known about the transcriptional regulation for most plant lincRNAs. Through the integration of 8 annotation resources, we defined 6,599 high-confidence lincRNA loci in Arabidopsis (Arabidopsis thaliana). For lincRNAs belonging to different evolutionary age categories, we identified major differences in sequence and chromatin features, as well as in the level of conservation and purifying selection acting during evolution. Spatiotemporal gene expression profiles combined with transcription factor (TF) chromatin immunoprecipitation (ChIP) data were used to construct a TF-lincRNA regulatory network containing 2,659 lincRNAs and 15,686 interactions. We found that properties characterizing lincRNA expression, conservation, and regulation differ between plants and animals. Experimental validation confirmed the role of 3 TFs, KANADI 1, MYB DOMAIN PROTEIN 44, and PHYTOCHROME INTERACTING FACTOR 4, as key regulators controlling root-specific lincRNA expression, demonstrating the predictive power of our network. Furthermore, we identified 58 lincRNAs, regulated by these TFs, showing strong root cell type-specific expression or chromatin accessibility, which are linked with genome-wide association studies genetic associations related to root system development and growth. The multilevel genome-wide characterization covering chromatin state information, promoter conservation, and chromatin immunoprecipitation-based TF binding, for all detectable lincRNAs across 769 expression samples, permits rapidly defining the biological context and relevance of Arabidopsis lincRNAs through regulatory networks.

The duplication of genomes and genetic networks and its potential for evolutionary adaptation and survival during environmental turmoil

The importance of whole-genome duplication (WGD) for evolution is controversial. Whereas some view WGD mainly as detrimental and an evolutionary dead end, there is growing evidence that polyploidization can help overcome environmental change, stressful conditions, or periods of extinction. However, despite much research, the mechanistic underpinnings of why and how polyploids might be able to outcompete or outlive nonpolyploids at times of environmental upheaval remain elusive, especially for autopolyploids, in which heterosis effects are limited. On the longer term, WGD might increase both mutational and environmental robustness due to redundancy and increased genetic variation, but on the short-or even immediate-term, selective advantages of WGDs are harder to explain. Here, by duplicating artificially generated Gene Regulatory Networks (GRNs), we show that duplicated GRNs-and thus duplicated genomes-show higher signal output variation than nonduplicated GRNs. This increased variation leads to niche expansion and can provide polyploid populations with substantial advantages to survive environmental turmoil. In contrast, under stable environments, GRNs might be maladaptive to changes, a phenomenon that is exacerbated in duplicated GRNs. We believe that these results provide insights into how genome duplication and (auto)polyploidy might help organisms to adapt quickly to novel conditions and to survive ecological uproar or even cataclysmic events.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

High-quality Cymbidium mannii genome and multifaceted regulation of crassulacean acid metabolism in epiphytes

Epiphytes with crassulacean acid metabolism (CAM) photosynthesis are widespread among vascular plants, and repeated evolution of CAM photosynthesis is a key innovation for micro-ecosystem adaptation. However, we lack a complete understanding of the molecular regulation of CAM photosynthesis in epiphytes. Here, we report a high-quality chromosome-level genome assembly of a CAM epiphyte, Cymbidium mannii (Orchidaceae). The 2.88-Gb orchid genome with a contig N50 of 22.7 Mb and 27 192 annotated genes was organized into 20 pseudochromosomes, 82.8% of which consisted of repetitive elements. Recent expansions of long terminal repeat retrotransposon families have made a major contribution to the evolution of genome size in Cymbidium orchids. We reveal a holistic scenario of molecular regulation of metabolic physiology using high-resolution transcriptomics, proteomics, and metabolomics data collected across a CAM diel cycle. Patterns of rhythmically oscillating metabolites, especially CAM-related products, reveal circadian rhythmicity in metabolite accumulation in epiphytes. Genome-wide analysis of transcript and protein level regulation revealed phase shifts during the multifaceted regulation of circadian metabolism. Notably, we observed diurnal expression of several core CAM genes (especially βCA and PPC) that may be involved in temporal fixation of carbon sources. Our study provides a valuable resource for investigating post-transcription and translation scenarios in C. mannii, an Orchidaceae model for understanding the evolution of innovative traits in epiphytes.

Genome-wide identification and analysis of MIKC-type MADS-box genes expression in Chimonanthus salicifolius

BACKGROUND

MIKC type MADS-box transcription factors are one of the largest gene families and play a pivotal role in flowering time and flower development. Chimonanthus salicifolius belongs to the family Calycanthaceae and has a unique flowering time and flowering morphology compared to other Chimonanthus species, but the research on MIKC type MADS-box gene family of C. salicifolius has not been reported.

OBJECTIVE

Identification, comprehensive bioinformatic analysis, the expression pattern of MIKC-type MADS-box gene family from different tissues of C. salicifolius.

METHODS

Genome-wide investigation and expression pattern under different tissues of the MIKC-type MADS-box gene family in C. salicifolius, and their phylogenetic relationships, evolutionary characteristics, gene structure, motif distribution, promoter cis-acting element were performed.

RESULTS

A total of 29 MIKC-type MADS-box genes were identified from the whole genome sequencing. Interspecies synteny analysis revealed more significant collinearity between C. salicifolius and the magnoliids species compared to eudicots and monocots. MIKC-type MADS-box genes from the same subfamily share similar distribution patterns, gene structure, and expression patterns. Compared with Arabidopsis thaliana, Nymphaea colorata, and Chimonanthus praecox, the FLC genes were absent in C. salicifolius, while the AGL6 subfamily was expanded in C. salicifolius. The selectively expanded promoter (AGL6) and lack of repressor (FLC) genes may explain the earlier flowering in C. salicifolius. The loss of the AP3 homologous gene in C. salicifolius is probably the primary cause of the morphological distinction between C. salicifolius and C. praecox. The csAGL6a gene is specifically expressed in the flowering process and indicates the potential function of promoting flowering.

CONCLUSION

This study offers a genome-wide identification and expression profiling of the MIKC-types MADS-box genes in the C. salicifolius, and establishes the foundation for screening flowering development genes and understanding the potential function of the MIKC-types MADS-box genes in the C. salicifolius.

Genome-Wide Analysis of the Membrane Attack Complex and Perforin Genes and Their Expression Pattern under Stress in the Solanaceae

The Membrane Attack Complex and Perforin (MACPF) proteins play a crucial role in plant development and adaptation to environmental stresses. Heretofore, few MACPF genes have been functionally identified, leaving gaps in our understanding of MACPF genes in other plants, particularly in the Solanaceae family, which includes economically and culturally significant species, such as tomato, potato, and pepper. In this study, we have identified 26 MACPF genes in three Solanaceae species and in the water lily, which serves as the base group for angiosperms. Phylogenetic analysis indicates that angiosperm MACPF genes could be categorized into three distinct groups, with another moss and spikemoss lineage-specific group, which is further supported by the examination of gene structures and domain or motif organizations. Through inter-genome collinearity analysis, it is determined that there are 12 orthologous SolMACPF gene pairs. The expansion of SolMACPF genes is primarily attributed to dispersed duplications, with purifying selection identified as the principal driving force in their evolutionary process, as indicated by the ω values. Furthermore, the analysis of expression patterns revealed that Solanaceae genes are preferentially expressed in reproductive tissues and regulated by various environmental stimuli, particularly induced by submergence. Taken together, these findings offer valuable insights into and a fresh perspective on the evolution and function of SolMACPF genes, thereby establishing a foundation for further investigations into their phenotypic and functional characteristics.

Structural and Functional Analysis of the MADS-Box Genes Reveals Their Functions in Cold Stress Responses and Flower Development in Tea Plant (Camellia sinensis)

MADS-box genes comprise a large family of transcription factors that play crucial roles in all aspects of plant growth and development. However, no detailed information on the evolutionary relationship and functional characterization of MADS-box genes is currently available for some representative lineages, such as the Camellia plant. In this study, 136 MADS-box genes were detected from a reference genome of the tea plant (Camellia sinensis) by employing a 569 bp HMM (Hidden Markov Model) developed using nucleotide sequencing including 73 type I and 63 type II genes. An additional twenty-seven genes were identified, with five MIKC-type genes. Truncated and/or inaccurate gene models were manually verified and curated to improve their functional characterization. Subsequently, phylogenetic relationships, chromosome locations, conserved motifs, gene structures, and gene expression profiles were systematically investigated. Tea plant MIKC genes were divided into all 14 major eudicot subfamilies, and no gene was found in Mβ. The expansion of MADS-box genes in the tea plant was mainly contributed by WGD/fragment and tandem duplications. The expression profiles of tea plant MADS-box genes in different tissues and seasons were analyzed, revealing widespread evolutionary conservation and genetic redundancy. The expression profiles linked to cold stress treatments suggested the wide involvement of MADS-box genes from the tea plant in response to low temperatures. Moreover, a floral 'ABCE' model was proposed in the tea plant and proved to be both conserved and ancient. Our analyses offer a detailed overview of MADS-box genes in the tea plant, allowing us to hypothesize the potential functions of unknown genes and providing a foundation for further functional characterizations.

Small Auxin Up RNA 56 (SAUR56) regulates heading date in rice

Heading date is a critical agronomic trait that determines crop yield. Although numerous genes associated with heading date have been identified in rice, the mechanisms involving Small Auxin Up RNA (SAUR) family have not been elucidated. In this study, the biological function of several SAUR genes was initially investigated using the CRISPR-Cas9 technology in the Japonica cultivar Zhonghua11 (ZH11) background. Further analysis revealed that the loss-of-function of OsSAUR56 affected heading date in both NLD (natural long-day) and ASD (artificial short-day). OsSAUR56 exhibited predominant expression in the anther, with its protein localized in both the cytoplasm and nucleus. OsSAUR56 regulated flowering time and heading date by modulating the expression of the clock gene OsGI, as well as two repressors Ghd7 and DTH8. Furthermore, haplotype-phenotype association analysis revealed a strong correlation between OsSAUR56 and heading date, suggesting its role in selection during the domestication of rice. In summary, these findings highlights the importance of OsSAUR56 in the regulation of heading date for further potential facilitating genetic engineering for flowering time during rice breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01409-w.

Climatic stability and geological history shape global centers of neo- and paleoendemism in seed plants

Assessing the distribution of geographically restricted and evolutionarily unique species and their underlying drivers is key to understanding biogeographical processes and critical for global conservation prioritization. Here, we quantified the geographic distribution and drivers of phylogenetic endemism for ~320,000 seed plants worldwide and identified centers and drivers of evolutionarily young (neoendemism) and evolutionarily old endemism (paleoendemism). Tropical and subtropical islands as well as tropical mountain regions displayed the world's highest phylogenetic endemism. Most tropical rainforest regions emerged as centers of paleoendemism, while most Mediterranean-climate regions showed high neoendemism. Centers where high neo- and paleoendemism coincide emerged on some oceanic and continental fragment islands, in Mediterranean-climate regions and parts of the Irano-Turanian floristic region. Global variation in phylogenetic endemism was well explained by a combination of past and present environmental factors (79.8 to 87.7% of variance explained) and most strongly related to environmental heterogeneity. Also, warm and wet climates, geographic isolation, and long-term climatic stability emerged as key drivers of phylogenetic endemism. Neo- and paleoendemism were jointly explained by climatic and geological history. Long-term climatic stability promoted the persistence of paleoendemics, while the isolation of oceanic islands and their unique geological histories promoted neoendemism. Mountainous regions promoted both neo- and paleoendemism, reflecting both diversification and persistence over time. Our study provides insights into the evolutionary underpinnings of biogeographical patterns in seed plants and identifies the areas on Earth with the highest evolutionary and biogeographical uniqueness-key information for setting global conservation priorities.

Comparative genomics reveals a unique nitrogen-carbon balance system in Asteraceae

The Asteraceae (daisy family) is one of the largest families of plants. The genetic basis for its high biodiversity and excellent adaptability has not been elucidated. Here, we compare the genomes of 29 terrestrial plant species, including two de novo chromosome-scale genome assemblies for stem lettuce, a member of Asteraceae, and Scaevola taccada, a member of Goodeniaceae that is one of the closest outgroups of Asteraceae. We show that Asteraceae originated ~80 million years ago and experienced repeated paleopolyploidization. PII, the universal regulator of nitrogen-carbon (N-C) assimilation present in almost all domains of life, has conspicuously lost across Asteraceae. Meanwhile, Asteraceae has stepwise upgraded the N-C balance system via paleopolyploidization and tandem duplications of key metabolic genes, resulting in enhanced nitrogen uptake and fatty acid biosynthesis. In addition to suggesting a molecular basis for their ecological success, the unique N-C balance system reported for Asteraceae offers a potential crop improvement strategy.

Genome-Wide Analysis of the MADS-box Gene Family and Expression Analysis during Anther Development in Salvia miltiorrhiza

MADS-box genes constitute a large family of transcription factors that play important roles in plant growth and development. However, our understanding of MADS-box genes involved in anther development and male sterility in Salvia miltiorrhiza is still limited. In this study, 63 MADS-box genes were identified from the genome of the male sterility ecotype Sichuan S. miltiorrhiza (S. miltiorrhiza_SC) unevenly distributed among eight chromosomes. Phylogenetic analysis classified them into two types and 17 subfamilies. They contained 1 to 12 exons and 10 conserved motifs. Evolution analysis showed that segmental duplication was the main force for the expansion of the SmMADS gene family, and duplication gene pairs were under purifying selection. Cis-acting elements analysis demonstrated that the promoter of SmMADS genes contain numerous elements associated with plant growth and development, plant hormones, and stress response. RNA-seq showed that the expression levels of B-class and C-class SmMADS genes were highly expressed during anther development, with SmMADS11 likely playing an important role in regulating anther development and male fertility in S. miltiorrhiza_SC. Overall, this study provides a comprehensive analysis of the MADS-box gene family in S. miltiorrhiza, shedding light on their potential role in anther development and male sterility.

Dim artificial light at night alters gene expression rhythms and growth in a key seagrass species (Posidonia oceanica)

Artificial light at night (ALAN) is a globally spreading anthropogenic stressor, affecting more than 20% of coastal habitats. The alteration of the natural light/darkness cycle is expected to impact the physiology of organisms by acting on the complex circuits termed as circadian rhythms. Our understanding of the impact of ALAN on marine organisms is lagging behind that of terrestrial ones, and effects on marine primary producers are almost unexplored. Here, we investigated the molecular and physiological response of the Mediterranean seagrass, Posidonia oceanica (L.) Delile, as model to evaluate the effect of ALAN on seagrass populations established in shallow waters, by taking advantage of a decreasing gradient of dim nocturnal light intensity (from < 0.01 to 4 lx) along the NW Mediterranean coastline. We first monitored the fluctuations of putative circadian-clock genes over a period of 24 h along the ALAN gradient. We then investigated whether key physiological processes, known to be synchronized with day length by the circadian rhythm, were also affected by ALAN. ALAN influenced the light signalling at dusk/night in P. oceanica, including that of shorter blue wavelengths, through the ELF3-LUX1-ZTL regulatory network, and suggested that the daily perturbation of internal clock orthologs in seagrass might have caused the recruitment of PoSEND33 and PoPSBS genes to mitigate the repercussions of a nocturnal stress on photosynthesis during the day. A long-lasting impairment of gene fluctuations in sites characterised by ALAN could explain the reduced growth of the seagrass leaves when these were transferred into controlled conditions and without lighting during the night. Our results highlight the potential contribution of ALAN to the global loss of seagrass meadows, posing questions about key interactions with a variety of other human-related stressors in urban areas, in order to develop more efficient strategies to globally preserve these coastal foundation species.

Identification of annotation artifacts concerning the chalcone synthase (CHS)

OBJECTIVE

Chalcone synthase (CHS) catalyzes the initial step of the flavonoid biosynthesis. The CHS encoding gene is well studied in numerous plant species. Rapidly growing sequence databases contain hundreds of CHS entries that are the result of automatic annotation. In this study, we evaluated apparent multiplication of CHS domains in CHS gene models of four plant species.

MAIN FINDINGS

CHS genes with an apparent triplication of the CHS domain encoding part were discovered through database searches. Such genes were found in Macadamia integrifolia, Musa balbisiana, Musa troglodytarum, and Nymphaea colorata. A manual inspection of the CHS gene models in these four species with massive RNA-seq data suggests that these gene models are the result of artificial fusions in the annotation process. While there are hundreds of seemingly correct CHS records in the databases, it is not clear why these annotation artifacts appeared.

Diploid and tetraploid genomes of Acorus and the evolution of monocots

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

Comparative genomics analysis of WAK/WAKL family in Rosaceae identify candidate WAKs involved in the resistance to Botrytis cinerea

BACKGROUND

Wall associated kinase (WAK) and WAK-like (WAKL) are typical pattern recognition receptors act as the first sentry of plant defense. But little of WAK/WAKL family is known in Rosaceae.

RESULTS

In this study, 131 WAK/WAKL genes from apple, peach and strawberry were identified using a bioinformatics approach. Together with 68 RcWAK/RcWAKL in rose, we performed a comparative analysis of 199 WAK/WAKL in four Rosaceae crops. The phylogenetic analysis divided all the WAK/WAKL into five clades. Among them, the cis-elements of Clade II and Clade V promoters were enriched in jasmonic acid (JA) signaling and abiotic stress, respectively. And this can also be verified by the rose transcriptome responding to different hormone treatments. WAK/WAKL families have experienced a considerable proportion of purifying selection during evolution, but still 26 amino acid sites evolved under positive selection, which focused on extracellular conserved domains. WAK/WAKL genes presented collinearity relationship within and between crops, throughout four crops we mined four orthologous groups (OGs). The WAK/WAKL genes in OG1 and OG4 were speculated to involve in plant-Botrytis cinerea interaction, which were validated in rose via VIGS as well as strawberry by qRT-PCR.

CONCLUSIONS

These results not only provide genetic resources and valuable information for the evolutionary relationship of WAK/WAKL gene family, but also offer a reference for future in-depth studies of Rosaceae WAK/WAKL genes.

Selection and Validation of Optimal RT-qPCR Reference Genes for the Normalization of Gene Expression under Different Experimental Conditions in Lindera megaphylla

Lindera megaphylla, a broad-leaved evergreen that is used as a landscape ornamental plant and medicinal plant, is an ecologically important and dominant tree species. However, little is known about the molecular mechanisms of its growth, development, and metabolism. The selection of suitable reference genes is critical for molecular biological analyses. To date, no research on reference genes as a foundation for gene expression analysis has been undertaken in L. megaphylla. In this study, 14 candidate genes were selected from the transcriptome database of L. megaphylla for RT-qPCR assay under different conditions. Results showed that helicase-15 and UBC28 were most stable in different tissues of seedlings and adult trees. For different leaf developmental stages, the best combination of reference genes was ACT7 and UBC36. UBC36 and TCTP were the best under cold treatment, while PAB2 and CYP20-2 were the best under heat treatment. Finally, a RT-qPCR assay of LmNAC83 and LmERF60 genes were used to further verify the reliability of selected reference genes above. This work is the first to select and evaluate the stability of reference genes for the normalization of gene expression analysis in L. megaphylla and will provide an important foundation for future genetic studies of this species.

Characterization and co-expression analysis of ATP-binding cassette transporters provide insight into genes related to cannabinoid transport in Cannabis sativa L

Plant ATP-binding cassette (ABC) transporters contribute the transport of diverse secondary metabolites. However, their roles in cannabinoid trafficking are still unsolved in Cannabis sativa. In this study, 113 ABC transporters were identified and characterized in C. sativa from their physicochemical properties, gene structure, and phylogenic relationship, as well as spatial gene expression patterns. Eventually, seven core transporters were proposed including one member in ABC subfamily B (CsABCB8) and six ABCG members (CsABCG4, CsABCG10, CsABCG11, CsABCG32, CsABCG37, and CsABCG41), harboring potential in participating cannabinoid transport, by combining phylogenetic and co-expression analysis from the gene and metabolite level. The candidate genes exhibited a high correlation with cannabinoid biosynthetic pathway genes and the cannabinoid content, and they were highly expressed where cannabinoids appropriately biosynthesized and accumulated. The findings underpin further research on the function of ABC transporters in C. sativa, especially in unveiling the mechanisms of cannabinoid transport to boost systematic and targeted metabolic engineering.

Comparative transcriptome sequencing analysis to postulate the scheme of regulated leaf coloration in Perilla frutescens

Perilla as herb, ornamental, oil and edible plant is widely used in East Asia. Until now, the mechanism of regulated leaf coloration is still unclear. In this study, four different kinds of leaf colors were used to measure pigment contents and do transcriptome sequence to postulate the mechanism of leaf coloration. The measurements of chlorophyll, carotenoid, flavonoid, and anthocyanin showed that higher contents of all the aforementioned four pigments were in full purple leaf 'M357', and they may be determined front and back leaf color formation with purple. Meanwhile, the content of anthocyanin was controlled back leaf coloration. The chromatic aberration analysis and correlative analysis between different pigments and L*a*b* values analysis also suggested front and back leaf color change was correlated with the above four pigments. The genes involved in leaf coloration were identified through transcriptome sequence. The expression levels of chlorophyll synthesis and degradation related genes, carotenoid synthesis related genes and anthocyanin synthesis genes showed up-/down-regulated expression in different color leaves and were consistent of accumulation of these pigments. It was suggested that they were the candidate genes regulated perilla leaf color formation, and genes including F3'H, F3H, F3',5'H, DFR, and ANS are probably important for regulating both front and back leaf purple formation. Transcription factors involved in anthocyanin accumulation, and regulating leaf coloration were also identified. Finally, the probable scheme of regulated both full green and full purple leaf coloration and back leaf coloration was postulated.

Green land: Multiple perspectives on green algal evolution and the earliest land plants

Green plants, broadly defined as green algae and the land plants (together, Viridiplantae), constitute the primary eukaryotic lineage that successfully colonized Earth's emergent landscape. Members of various clades of green plants have independently made the transition from fully aquatic to subaerial habitats many times throughout Earth's history. The transition, from unicells or simple filaments to complex multicellular plant bodies with functionally differentiated tissues and organs, was accompanied by innovations built upon a genetic and phenotypic toolkit that have served aquatic green phototrophs successfully for at least a billion years. These innovations opened an enormous array of new, drier places to live on the planet and resulted in a huge diversity of land plants that have dominated terrestrial ecosystems over the past 500 million years. This review examines the greening of the land from several perspectives, from paleontology to phylogenomics, to water stress responses and the genetic toolkit shared by green algae and plants, to the genomic evolution of the sporophyte generation. We summarize advances on disparate fronts in elucidating this important event in the evolution of the biosphere and the lacunae in our understanding of it. We present the process not as a step-by-step advancement from primitive green cells to an inevitable success of embryophytes, but rather as a process of adaptations and exaptations that allowed multiple clades of green plants, with various combinations of morphological and physiological terrestrialized traits, to become diverse and successful inhabitants of the land habitats of Earth.

The chromosome-level genome of Eucommia ulmoides provides insights into sex differentiation and α-linolenic acid biosynthesis

Eucommia ulmoides Oliver is a typical dioecious plant endemic to China that has great medicinal and economic value. Here, we report a high-quality chromosome-level female genome of E. ulmoides obtained by PacBio and Hi-C technologies. The size of the female genome assembly was 1.01 Gb with 17 pseudochromosomes and 31,665 protein coding genes. In addition, Hi-C technology was used to reassemble the male genome released in 2018. The reassembled male genome was 1.24 Gb with the superscaffold N50 (48.30 Mb), which was increased 25.69 times, and the number of predicted genes increased by 11,266. Genome evolution analysis indicated that E. ulmoides has undergone two whole-genome duplication events before the divergence of female and male, including core eudicot γ whole-genome triplication event (γ-WGT) and a recent whole genome duplication (WGD) at approximately 27.3 million years ago (Mya). Based on transcriptome analysis, EuAP3 and EuAG may be the key genes involved in regulating the sex differentiation of E. ulmoides. Pathway analysis showed that the high expression of ω-3 fatty acid desaturase coding gene EU0103017 was an important reason for the high α-linolenic acid content in E. ulmoides. The genome of female and male E. ulmoides presented here is a valuable resource for the molecular biological study of sex differentiation of E. ulmoides and also will provide assistance for the breeding of superior varieties.

Genome-Wide Identification of the MADS-Box Gene Family during Male and Female Flower Development in Chayote (Sechium edule)

The MADS-box gene plays an important role in plant growth and development. As an important vegetable of Cucurbitaceae, chayote has great edible and medicinal value. So far, there is little molecular research on chayote, and there are no reports on the MADS-box transcription factor of chayote. In this study, the MADS-box gene family of chayote was analyzed for the first time, and a total of 70 MADS-box genes were identified, including 14 type I and 56 type II MICK MADS genes. They were randomly distributed on 13 chromosomes except for chromosome 11. The light response element, hormone response element and abiotic stress response element were found in the promoter region of 70 MADS genes, indicating that the MADS gene can regulate the growth and development of chayote, resist abiotic stress, and participate in hormone response; GO and KEGG enrichment analysis also found that SeMADS genes were mainly enriched in biological regulation and signal regulation, which further proved the important role of MADS-box gene in plant growth and development. The results of collinearity showed that segmental duplication was the main driving force of MADS gene expansion in chayote. RNA-seq showed that the expression levels of SeMADS06, SeMADS13, SeMADS26, SeMADS28, SeMADS36 and SeMADS37 gradually increased with the growth of chayote, indicating that these genes may be related to the development of root tubers of 'Tuershao'. The gene expression patterns showed that 12 SeMADS genes were specifically expressed in the male flower in 'Tuershao' and chayote. In addition, SeMADS03 and SeMADS52 may be involved in regulating the maturation of male flowers of 'Tuershao' and chayote. SeMADS21 may be the crucial gene in the development stage of the female flower of 'Tuershao'. This study laid a theoretical foundation for the further study of the function of the MADS gene in chayote in the future.