Chromosome-level genome assembly of a parent species of widely cultivated azaleas

Genome-wide identification, expression analysis, and stress response analysis of the RdbZIP gene family in Rhododendron delavayi

BACKGROUND

Basic leucine zipper (bZIP) gene family members represent one of the most diverse and largest groups of transcription factors in eukaryotes. Research has demonstrated that bZIP transcription factors play crucial roles not only in plant growth and development but also in response to various abiotic stresses. However, studies focusing on bZIP factors in Rhododendron delavayi (RdbZIPs) remain limited.

RESULT

In this study, a total of 59 RdbZIPs were identified using bioinformatics approaches, and these could be classified into 13 subfamilies based on the genomic data of R. delavayi. Members of the same RdbZIP subfamily exhibited similar gene structures and conserved motifs, and were unevenly distributed across the 13 chromosomes of R. delavayi. Collinearity analysis revealed a total of 20 duplication events, comprising 3 pairs of tandem duplications and 17 pairs of segmental duplications. Additionally, cis-acting element analysis indicated that RdbZIP family members may be involved in various biological processes, including transcription, development, hormone regulation, and responses to biotic and abiotic stresses. Transcriptomic analysis revealed that RdbZIP family genes were highly expressed in reproductive tissues. RT-qPCR expression analysis revealed that many selected RdbZIP genes were significantly upregulated under high salinity and drought conditions, suggesting their potential involvement in stress-responsive regulatory networks.

CONCLUSION

This study provides the first comprehensive characterization of the bZIP transcription factor family in Rhododendron delavayi, laying a foundational framework for functional studies of individual RdbZIP genes. The results highlight the pivotal role of RdbZIP genes in abiotic stress tolerance, which is crucial for understanding the adaptive mechanisms of R. delavayi. Future research should focus on the functional validation of key RdbZIP genes and elucidation of their regulatory pathways, which may contribute to the genetic improvement of Rhododendron species under adverse environmental conditions.

CLINICAL TRIAL

Not applicable.

SSR_VibraProfiler: a Python package for accurate classification of varieties using SSRs with intra-variety specificity and inter-variety polymorphism

BACKGROUND

Simple sequence repeats (SSRs) are widely used as molecular markers; however, traditional development of SSR molecular markers heavily relies on experimental methods. The advancement of modern sequencing technology has provided the possibility of directly extracting SSR characteristics from sequencing data and using them for variety identification.

RESULTS

We have developed a computational framework for variety identification, treating the presence or absence of each SSR in sequencing data as a numerical characteristic while ignoring specific loci, flanking sequences, and occurrence counts. Therefore, subsequent variety identification does not rely on experimental validation but is directly performed based on the numerical characteristic matrix. Using a formula, we measure the variance of these numerical characteristics both within and among varieties, and select SSRs that exhibit intra-variety specificity and inter-variety polymorphism, forming a 0,1 matrix. We use t-SNE (t-distributed Stochastic Neighbor Embedding) to project the matrix onto a two-dimensional plane, followed by K-means clustering of the individuals. The classification performance of the matrix is preliminarily assessed by comparing the cluster labels with the true labels, providing an initial evaluation of its effectiveness in variety detection. Ultimately, we construct a recognition model based on the SSRs matrix and apply it for variety identification. The process has been encapsulated into the package SSR_VibraProfiler, which can serve as a tool for constructing an SSR variety DNA fingerprint database. We tested this package on a Rhododendron dataset that included 40 individuals from 8 varieties. The accuracy achieved through t-SNE dimensionality reduction and K-means clustering was 100%. Furthermore, we used the leave-one-out method to validate the accuracy of our method in predicting variety, and confirmed the reliability of our method in detecting varieties. The package is freely available at https://github.com/Olcat35412/SSR_VibraProfiler .

CONCLUSION

We introduced SSR_VibraProfiler, a Python package for distinguishing and predicting individual varieties without a reference genome by extracting SSR numerical characteristics from next-generation sequencing data. This tool will contribute to the development, identification, and protection of new varieties.

Phylogenetic and expression analysis of HSP20 gene family in Rhododendron species of different altitudes

Small heat shock proteins (sHSPs/HSP20s) play important roles in regulating plant growth, development and stress responses, especially heat stress. Rhododendron plants are major components in landscaping or potting, but poor high-temperature tolerance limits their wide application. To elucidate the adaptive differences between Rhododendron species inhabiting high- and low-altitude regions, this study identified a total of 265 HSP20 genes across seven Rhododendron species, with categorizing into 11 subfamilies. In the CI subfamily, low-altitude Rhododendron species had more HSP20 genes than high-altitude species. Ka/Ks analysis indicated that nearly all HSP20 genes in the seven Rhododendron species have undergone purifying selection, with only a few in the low-altitude species exhibiting positive selection. Analysis of cis-acting elements revealed that most HSP20 genes in R. ovatum and R. simsii could respond to a variety of plant hormones and stresses. Expression pattern analysis revealed that HSP20 members are implicated in flower development and response to heat stress, with the CI subfamily being the main branch responsible for the heat stress response in low-altitude Rhododendron species. Heat stress treatment of transgenic yeast further validated the crucial role of CI subfamily genes in heat stress tolerance. This study provides the first analysis of evolutionary differences in the HSP20 gene families between high- and low-altitude Rhododendron species. It offers insights into the evolutionary direction of HSP20 genes and identifies key genes related to heat tolerance. Additionally, it highlights the role of CI subfamily genes in heat stress tolerance, contributing to the development of heat-tolerant Rhododendron varieties and advancing flower development research.

High-Quality Genome Assembly and Transcriptome of Rhododendron platypodum Provide Insights into Its Evolution and Heat Stress Response

R. platypodum (Rhododendron platypodum) is an endangered alpine species with a highly restricted distribution in the southwestern region of China, which possesses significant ornamental and horticultural value. In this study, the high-quality genome assembly of R. platypodum at the chromosomal level is reported. The total genome size was determined to be 642.25 Mb, with a contig N50 of 25.64 Mb, and it contains 36,522 predicted genes. Comparative genomic analysis between R. platypodum and other species revealed the expansion of gene families, such as those related to transition metal ion binding and sodium ion transport, as well as the contraction of gene families involved in the recognition of pollen and pollen-pistil interaction. These findings might explain the adaptation of R. platypodum to rocky habitats and contribute to its endangered status. Furthermore, a heat stress experiment was conducted on R. platypodum, followed by transcriptome sequencing and physiological co-analysis to construct a co-expression network. This analysis identified the candidate gene TAR1-A and other transcription factors exhibiting differential expression under heat stress. The whole-genome sequencing, transcriptome analysis, and physiological co-analysis of R. platypodum provide valuable resources for its conservation and offer insights into its mechanisms of heat stress.

Sucrose synthase 3 improves fruit quality in grape

Sucrose synthase (SS) is one of the key functional enzymes involved in sugar metabolism in plants, and its activity can directly affect sugar accumulation, thereby influencing fruit quality. Our previous research found an interaction between VvSS3 and VvSnRK1β. In this study, we discovered that SS enzyme activity was mainly oriented towards decomposition, with the highest activity occurring before the veraison stage, and the overall trend of SS enzyme activity changes was positively correlated with exogenous ABA concentration. Site-directed mutagenesis revealed that VvSS3S176 and VvSS3S381 were crucial sites for the interaction between VvSS3 and VvSnRK1β. Subcellular localization results showed that VvSS3S176 and VvSS3S381, as well as VvSS3, were located on the cell membrane. VvSS3 was sensitive to ABA, promoting the accumulation of soluble sugars and anthocyanins in transgenic callus, increasing endogenous ABA content, and reducing organic acid components. VvSS3S176 and VvSS3S381 only altered the activity of SS in the synthesis direction, while the transcription level of VvSnRK1β in transgenic callus significantly decreased. After exogenous ABA treatment, Ser176 and Ser381 reduced the inhibition of VvSS3 on VvSnRK1β expression. Mutation of the binding sites of VvSS3 prevented the formation of complexes VvSnRK1β-VvSS3, thus interfering with downstream metabolism. This suggests that VvSnRK1β may form a protein complex by interacting with VvSS3, participating in the accumulation of soluble sugars in grape fruits mediated by ABA response.

A Comprehensive Evolutionary Analysis of the Dihydroflavonol 4-Reductase (DFR) Gene Family in Plants: Insights from 237 Species

BACKGROUND

Dihydroflavonol 4-reductase (DFR) is a key enzyme in the flavonoid biosynthetic pathway that regulates anthocyanin and proanthocyanidin accumulation in plants. Although DFR genes have been studied in various species, their origin of the DFR gene family, its distribution across the plant kingdom, and the reasons behind the emergence of different DFR subtypes Methods: This study performed a whole-genome analysis of DFR genes in 237 plant species, including algae, mosses, ferns, gymnosperms, and angiosperms, integrating phylogeny, conserved motifs, duplication mechanisms, positive selection, and expression pattern analyses.

RESULTS

These results indicate that the DFR gene family originated from the common ancestor of extant ferns and seed plants, and the emergence of asparagine (Asn)-type and aspartic (Asp)-type DFRs is associated with gymnosperms. Notably, we report for the first time the presence of Asn-type, Asp-type, and arginine (Arg)-type DFRs in some species, which breaks the previous notion that Arg-type DFRs are exclusive to ferns. Tandem duplication is considered the primary driving force behind the expansion of the DFR family and is associated with the formation of different DFR subtypes. Furthermore, Asn-type DFRs were highly expressed during the early stages of seed development, suggesting their important role in seed development.

CONCLUSIONS

Overall, this study revealed the dynamic evolutionary trajectory of the DFR gene family in plants, providing a theoretical foundation for future research on DFR genes.

Genome assembly of Stewartia sinensis reveals origin and evolution of orphan genes in Theaceae

Orphan genes play crucial roles in diverse biological processes, but the evolutionary trajectories and functional divergence remain largely unexplored. The Theaceae family, including the economically and culturally important tea plant, offers a distinctive model to examine these aspects. Here, we integrated Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods to decode a pseudo-chromosomal genome assembly of Stewartia sinensis, from the earliest-diverging tribe of Theaceae, spanning 2.95 Gb. Comparative genomic analysis revealed the absence of recent whole-genome duplication events in the Theaceae ancestor, highlighting tandem duplications as the predominant mechanism of gene expansion. We identified 31,331 orphan genes, some of which appear to have ancient origins, suggesting early emergence with frequent gains and losses, while others seem more specific and recent. Notably, orphan genes are distinguished by shorter lengths, fewer exons and functional domains compared to genes that originate much earlier, like transcription factors. Moreover, tandem duplication contributes significantly to the adaptive evolution and characteristic diversity of Theaceae, and it is also a major mechanism driving the origination of orphan genes. This study illuminates the evolutionary dynamics of orphan genes, providing a valuable resource for understanding the origin and evolution of tea plant flavor and enhancing genetic breeding efforts.

Genome assembly and multiomic analyses reveal insights into flower and bark colors of Lagerstroemia excelsa

Lagerstroemia excelsa is a unique plant species from China, holds a significant aesthetic and economic value, and plays a crucial role in landscape architecture and horticulture. Thus far, there is little genetic and genomic information available about this species, which limits its use in development of new cultivars. In this study, a high-quality genome map of L. excelsa was obtained via whole-genome sequencing. Results showed that its genome size is about 330.4 Mb and a scaffold mapping rate is approximately 97.20%, resulting in 24 pseudochromosomes. L. excelsa might have undergone a recent whole-genome triplication event and diverged from the pomegranate about 32.3 million years ago (MYA). Subsequently, the divergence time between L. indica and L. excelsa was around 5.9 MYA. The transcriptomic and metabolomic analyses of L. excelsa and L. indica indicated that the chalcone synthase pathway may play a key role in regulating flower color differentiation between the two species. Additionally, a transcription factor LeMYB103 may be involved in regulating anthocyanin synthesis by interacting with LeMYB66, resulting in the accumulation of anthocyanins in the stem bark. This study is the first step toward genomic analysis of L. excelsa, which may provide a foundation for further molecular investigation of this species and offer valuable insights into the molecular mechanisms underlying the flower and stem bark colors in L. excelsa, two important ornamental traits in Lagerstroemia breeding.

Chromosome-level genome assembly of Pontederia cordata L. provides insights into its rapid adaptation and variation of flower colours

Pontederia cordata L. is an aquatic ornamental plant native to the Americas but has been widely distributed in South Asia, Australia, and Europe. The genetic mechanisms behind its rapid adaptation and spread have not yet been well understood. To understand the mechanisms for its rapid adaptation, this study assembled the first chromosome-level genome of P. cordata. The genome assembly, which spans 527.5 Mb, is anchored on 8 pseudochromosomes with a scaffold N50 of 48 Mb and encompasses 29,389 protein-coding genes. Further analyses revealed that P. cordata had experienced 3 whole-genome duplications (WGDs) events. These WGDs are associated with gene family expansion and increased numbers of resistance gene analogs and transcription factors. Positive selection analysis indicated that genes derived from tandem duplication (TD) and proximal duplication were more likely to undergo positive selection, and were enriched in plant defense and disease resistance. These results implied that WGDs, TD, and positive selection enhanced the environmental adaptability of P. cordata. In addition, we found that down-regulation of F3'5'H, DFR, ANS, and UFGT likely caused the flower colour variation for P. cordata from violet to white. The first chromosome-level genome of P. cordata here provides a valuable genomic resource for investigating the rapid adaptation and flower colour variation of the species.

The chromosome-scale genome of black wolfberry (Lycium ruthenicum) provides useful genomic resources for identifying genes related to anthocyanin biosynthesis and disease resistance

The black wolfberry (L ycium ruthenicum; 2n = 2x = 24) is an important medicinal plant with ecological and economic value. Its fruits have numerous beneficial pharmacological activities, especially those of anthocyanins, polysaccharides, and alkaloids, and have high nutritional value. However, the lack of available genomic resources for this species has hindered research on its medicinal and evolutionary mechanisms. In this study, we developed the telomere-to-telomere (T2T) nearly gapless genome of L. ruthenicum (2.26 Gb) by integrating PacBio HiFi, Nanopore Ultra-Long, and Hi-C technologies. The assembled genome comprised 12 chromosomes with 37,149 protein-coding genes functionally annotated. Approximately 80% of the repetitive sequences were identified, of which long terminal repeats (LTRs) were the most abundant, accounting for 73.01%. The abundance of LTRs might be the main reason for the larger genome of this species compared to that of other Lycium species. The species-specific genes of L. ruthenicum were related to defense mechanisms, salt tolerance, drought resistance, and oxidative stress, further demonstrating their superior adaptability to arid environments. Based on the assembled genome and fruit transcriptome data, we further constructed an anthocyanin biosynthesis pathway and identified 19 candidate structural genes and seven transcription factors that regulate anthocyanin biosynthesis in the fruit developmental stage of L. ruthenicum, most of which were highly expressed at a later stage in fruit development. Furthermore, 154 potential disease resistance-related nucleotide-binding genes have been identified in the L. ruthenicum genome. The whole-genome and proximal, dispersed, and tandem duplication genes in the L. ruthenicum genome enriched the number of genes involved in anthocyanin synthesis and resistance-related pathways. These results provide an important genetic basis for understanding genome evolution and biosynthesis of pharmacologically active components in the Lycium genus.

A genome-wide-level insight into the HSF gene family of Rhodomyrtus tomentosa and the functional divergence of RtHSFA2a and RtHSFA2b in thermal adaptation

Heat shock transcription factor (HSF) is one of the most important regulatory elements in plant development and stress response. Rhohomyrtus tomentosa has many advantages in adapting to high temperature and high humidity climates, whereas its inherence has barely been elucidated. In this study, we aimed to characterize the HSF family and investigate the thermal adaptation mechanisms of R. tomentosa. We identified 25 HSF genes in the R. tomentosa genome. They could be classified into three classes: HSFA, HSFB, and HSFC. Gene duplication events are major motivations for the expansion of the RtHSF gene family. Most of the genes in the same subclass share similar conserved motifs and gene structures. The cis-acting elements of the promoter regions of RtHSF genes are related to development, phytohormone signaling, and stress responses, and they vary among the genes even in the same subclass, resulting in different expression patterns. Especially, there exists subfunctionalization in the RtHSFA2 subfamily in responding to various abiotic stresses, viz. RtHSFA2a is sensitive to drought, salt, and cold stresses, whilst RtHSFA2b is mainly induced by heat stress. We further proved that RtHSFA2b might be of more importance in R. tomentosa thermotolerance, for Arabidopsis plants with overexpressed RtHSFA2b outperformed those with RtHSFA2a under heat stress, and RtHSFA2b had much higher transcription activity than RtHSFA2a in regulating certain heat shock response (HSR) genes. RtHSFA2a plays a role in transactivating RtHSFA2b. All these results provide a general prospect of the RtHSF gene family and enclose a basal thermal adaptation mechanism of R. tomentosa.

Divergent evolutionary paces among eudicot plants revealed by simultaneously duplicated genes produced billions of years ago

Polyploidization often occurs more than once along an evolutionary lineage to form extant plants. Major core eudicot plants share a whole-genome triplication (ceWGT), through which thousands of simultaneously duplicated genes are retained in extant genomes, providing a valuable starting line to check the difference in their evolutionary paces. Here, by characterizing the synonymous nucleotide substitutions (Ks) between these duplicates from 28 representative plants from 21 families, we checked the various evolutionary rates among plants among plants subjected to different rounds of extra polyploidization events. We found up to 68.04% difference in evolutionary rates among the selected plants. A statistical correlation analysis (correlation coefficient =0.57, at significant level = 0.01) indicated that plants affected by extra polyploidies have evolved faster than plants without such extra polyploidies showing that (additional) polyploidization has resulted in elevated genetic diversity. Comparing the plants affected by additional polyploidization and plants without it, the duplicated genes produced by the ceWGT and retained in extant genomes have gathered 4.75% more nucleotide substitutions in the former plants. By identifying the fast- and slowly evolving genes, we showed that genes evolving at divergent rates were often related to different evolutionary paths. By performing correction to evolutionary rates using a genome-scale approach, we revised the estimated timing of key evolutionary events. The present effort exploited the simultaneously duplicated genes produced by the shared polyploidization and help deepen the understanding of the role of polyploidization, especially its long-term effect in plant evolution and biological innovation.

Chromatin Topological Domains Associate With the Rapid Formation of Tandem Duplicates in Plants

In eukaryotes, chromatin is compacted within nuclei under the principle of compartmentalization. On top of that, condensin II establishes eukaryotic chromosome territories, while cohesin organizes the vertebrate genome by extruding chromatin loops and forming topologically associating domains (TADs). Thus far, the formation and roles of these chromatin structures in plants remain poorly understood. This study integrates Hi-C data from diverse plant species, demonstrating that nuclear DNA content influences large-scale chromosome conformation and affects the finer details of compartmental patterns. These contrasting compartmental patterns are distinguished by gene-to-gene loops and validated through cytological observations. Additionally, a novel chromatin domain type associated with tandem duplicate gene clusters is identified. These domains are independent of H3K27me3-mediated chromatin compartmentalization and exhibit evolutionary conservation across species. Gene pairs within TAD-like domains are younger and show higher levels of coexpression. These domains potentially promote the formation of tandem duplicates, a property appears unique to the Actinidia family. Overall, this study reveals functional chromatin domains in plants and provides evidence for the role of three-dimensional chromatin architecture in gene regulation and genome evolution.

Maximising informativeness for target capture-based phylogenomics in Erica (Ericaceae)

Plant phylogenetics has been revolutionised in the genomic era, with target capture acting as the primary workhorse of most recent research in the new field of phylogenomics. Target capture (aka Hyb-Seq) allows researchers to sequence hundreds of genomic regions (loci) of their choosing, at relatively low cost per sample, from which to derive phylogenetically informative data. Although this highly flexible and widely applicable method has rightly earned its place as the field's de facto standard, it does not come without its challenges. In particular, users have to specify which loci to sequence-a surprisingly difficult task, especially when working with non-model groups, as it requires pre-existing genomic resources in the form of assembled genomes and/or transcriptomes. In the absence of taxon-specific genomic resources, target sets exist that are designed to work across broad taxonomic scales. However, the highly conserved loci that they target may lack informativeness for difficult phylogenetic problems, such as that presented by the rapid radiation of Erica in southern Africa. We designed a target set for Erica phylogenomics intended to maximise informativeness and minimise paralogy while maintaining universality by including genes from the widely used Angiosperms353 set. Comprising just over 300 genes, the targets had excellent recovery rates in roughly 90 Erica species as well as outgroups from Calluna, Daboecia, and Rhododendron, and had high information content as measured by parsimony informative sites and Quartet Internode Resolution Probability (QIRP) at shallow nodes. Notably, QIRP was positively correlated with intron content, while including introns in targets-rather than recovering them via exon-flanking "bycatch"-substantially improved intron recovery. Overall, our results show the value of building a custom target set, and we provide a suite of open-source tools that can be used to replicate our approach in other groups (https://github.com/SethMusker/TargetVet).

Novel Insights into Hg0 Oxidation in Rice Leaf: Catalase Functions and Transcriptome Responses

Rice leaves can assimilate atmospheric mercury (Hg0), which is accumulated by grains and causes health risks to rice consumers. However, the molecular mechanisms underlying Hg0 assimilation in rice leaves remain poorly understood. Here, we investigated catalase's (CAT) function in Hg0 oxidation within rice leaves, as well as the Hg speciation and transcriptomic profiles of rice leaves exposed to Hg0. The inactivation of catalase reduced Hg0 oxidation by 91% in the leaf homogenate and the Hg0 oxidation rate increased along with CAT activity, showing the CAT's function in Hg0 oxidation. Hg0 was converted to Hg(cysteine)2 complexes in the leaf. Transcriptomic results revealed that the expression levels of both OsCATA and OsCATB (catalase-encoding genes) increased with Hg concentration, suggesting the involvement of catalase-related molecular network in Hg0 oxidation. Upstream transcription factors, including NAC (NAM-no apical meristem, ATAF-Arabidopsis transcription activation factor, and CUC-cup-shaped cotyledon), and ethylene-responsive transcription factor, are likely involved in catalase expression. Genes related to cysteine metabolism and amino acid transport appeared to regulate Hg accumulation. Our findings demonstrate the important function of catalase in Hg0 oxidation within rice and are fundamental for developing genetically modified rice cultivars to minimize human Hg exposure health risks.

R2R3 MYB Transcription Factors Involved in Flower Petal Pigmentation via Regulating Anthocyanin Synthesis in Rhododendron simsii

Azaleas (Rhododendron simsii) are popular ornamental woody plants known for their bright colors; however, very limited studies have been reported on the process of flower petal pigmentation. In this study, we found significant differences in the anthocyanin contents of petals from different colored azaleas, and the results of quantitative real-time PCR indicated that the R2R3 MYB genes, RsMYB12, RsMYB90, and RsMYB123, showed significant expression changes during the petal coloration in azalea petals; therefore, we hypothesized that RsMYB12, RsMYB90, and RsMYB123 might involve in the coloring process of azalea petals by regulating anthocyanin synthesis. This work provides insights into the underlying mechanisms of petal pigmentation in R. simsii and provides candidate genes for flower color breeding of azaleas and other ornamental flowers.

Lonicera caerulea genome reveals molecular mechanisms of freezing tolerance and anthocyanin biosynthesis

INTRODUCTION

Lonicera caerulea L. (blue honeysuckle) is a noteworthy fleshy-fruited tree and a prominent medicinal plant, which possesses notable characteristics such as exceptional resilience to winter conditions and early maturation, and the richest source of functional anthocyanins, particularly cyanidin-3-glucoside. The molecular mechanisms responsible for its freezing tolerance and anthocyanin biosynthesis remain largely unknown.

OBJECTIVES

Here, a chromosome-scale genome of L. caerulea was presented, aiming to examine the genetic foundations that underlie these characteristics of blue honeysuckle.

METHODS

The PacBio HiFi reads and Hi-C data were used to construct high-quality genome of blue honeysuckle. Comparative genomic and transcriptomic analyses were conducted to elucidate the molecular mechanisms of freezing tolerance and anthocyanin biosynthesis.

RESULTS

Comparative genomics analysis between L. caerulea and L. japonica revealed that the dynamic changes of duplicated genes contributed to their phytochemical reconstruction and environmental adaptation. Moreover, the ABA and ICE-CBF-COR signaling pathways were closely correlated to the freezing tolerance of L. caerulea. Genome-wide identification and biochemical function indicated that three anthocyanin 3',5'-O-methyltransferases (LcOMT2, LcOMT14, and LcOMT20) and two 3'-O-glycosyltransferases (LcUGT78X1 and LcUGT95P1) were responsible for anthocyanin biosynthesis. In addition, LcUGT78X1 was regarded as the potent glycosyltransferase for the accumulation of cyanidin-3-glucoside in L. caerulea.

CONCLUSION

This research elucidates the crucial roles of the ABA and ICE-CBF-COR signaling pathways in enhancing freezing tolerance, while also identifying highly efficient anthocyanin biosynthetic enzymes in L. caerulea. These findings advance the understanding of environmental adaptation and phytochemical production in Lonicera species.

Analysis of bZIP transcription factors in Rhododendron simsii and functional study of RsbZIP6 in regulating anthocyanin biosynthesis

The basic leucine zipper (bZIP) transcription factors play a critical role in various plant biological processes, including anthocyanin biosynthesis. This study focuses on Rhododendron simsii, a notable ornamental species with insufficiently explored bZIP transcription factors. We identified 66 bZIP transcription factors in the R. simsii genome and conducted comprehensive bioinformatics analyses to determine their gene localization, phylogenetic relationships, grouping, gene/protein structure, duplication events, synteny, and expression profiles. Our analysis identified RsbZIP6, a homolog of HY5 known to influence anthocyanin biosynthesis in many plants, as a potential regulator of this pathway. We cloned the complete coding sequence of RsbZIP6, which encodes a 170-amino acid protein spanning 510 bp. Subcellular localization analysis verified the nuclear presence of the RsbZIP6 protein. RT-qPCR analysis revealed the highest expression of RsbZIP6 in petals, which correlated with anthocyanin accumulation. Transgenic experiments indicated that overexpressing RsbZIP6 in Arabidopsis enhanced anthocyanin accumulation by upregulating genes involved in anthocyanin biosynthesis (4CL, CHS, CHI, DFR, F3H, F3'H, ANS and UF3GT). Our findings enhance understanding of the bZIP transcription factor family in R. simsii and underscore the vital role of RsbZIP6 in anthocyanin biosynthesis, providing insights for future genetic enhancement strategies.

First De Novo genome assembly and characterization of Gaultheria prostrata

Gaultheria Kalm ex L. (Ericaceae), a type of evergreen shrub, known as a natural source of methyl salicylate, possesses rich germplasm resources, strong habitat adaptability, significant ornamental value, and noteworthy pharmacological activities. However, due to the paucity of whole genomic information, genetically deep research in these areas remains limited. Consequently, we intend to obtain genome data through high-throughput sequencing, gene annotation, flow cytometry, transcription factors prediction and genetic marker analysis for a representative species of this genus, with Gaultheria prostrata selected for our study. In this study, we preliminarily obtained the genome of G. prostrata through next-generation sequencing methods. Utilizing 47.94 Gb of high-quality sequence data (108.95× coverage), assembled into 114,436 scaffolds, with an N50 length of 33,667 bp. The genome size assembled by SOAPdenovo, approximately 417 Mb, corresponded closely to predictions by flow cytometry (440 Mb) and k-mer analysis (447 Mb). The genome integrity was evaluated using BUSCO with 91%. The heterozygosity ratio was 0.159%, the GC content was 38.85%, and the repetitive regions encompassed over 34.6% of the genome. A total of 26,497 protein-coding genes have been predicted and annotated across Nr, Swissprot, GO, KEGG, and Pfam databases. Among these, 14,377 and 2,387 genes received functional annotation in Nr and Swissprot, respectively; 21,895, 24,424, and 22,330 genes were similarly annotated in GO, KEGG, and Pfam. Moreover, A total of 279,785 SSRs were identified and 345,270 primers for these SSRs were designed. Within the various nucleotide types of SSRs, AG/CT and AAG/CTT constituted the predominant dinucleotide and trinucleotide repeat types in G. prostrata. In addition, 1,395 transcription factors (TFs) from 75 TF families, 462 transcription regulators (TRs) from 33 TR families and 840 protein kinase (PKs) from 118 PK families were identified in this genome. We also performed phylogenetic analyses of G. prostrata and related species, including estimation of divergence times and expansion and contraction analyses, followed by positive selection analyses of orthologous gene pairs of G. prostrata and its close relative Vaccinium corymbosum. These results provide a reference for in-depth study of genus Gaultheria, contributing to future functional and comparative genomics analyses and providing supporting data for the development of molecular markers.

Genomic insights into adaptive evolution of the species-rich cosmopolitan plant genus Rhododendron

The species-rich cosmopolitan genus Rhododendron offers a good system for exploring the genomic mechanisms underlying adaptation to diverse habitats. Here, we report high-quality chromosomal-level genome assemblies of nine species, representing all five subgenera, different altitudinal distributions, and all flower color types of this genus. Further comprehensive genomic analyses indicate diverse adaptive strategies employed by Rhododendron, particularly adaptation to alpine and subalpine habitats by expansion/contraction of gene families involved in pathogen defense and oxidative phosphorylation, genomic convergent evolution, and gene copy-number variation. The convergent adaptation to high altitudes is further shown by population genomic analysis of R. nivale from the Himalaya-Hengduan Mountains. Moreover, we identify the genes involved in the biosynthesis of anthocyanins and carotenoids, which play a crucial role in shaping flower color diversity and environmental adaptation. Our study is significant for comprehending plant adaptive evolution and the uneven distribution of species diversity across different geographical regions.

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.

Identification of key genes involved in lignin and flavonoid accumulation during Tilia tuan seed maturation

KEY MESSAGE

Transcriptomics and phenotypic data analysis identified 24 transcription factors (TFs) that play key roles in regulating the competitive accumulation of lignin and flavonoids. Tilia tuan Szyszyl. (T. tuan) is a timber tree species with important ecological and commercial value. However, its highly lignified pericarp results in a low seed germination rate and a long dormancy period. In addition, it is unknown whether there is an interaction between the biosynthesis of flavonoids and lignin as products of the phenylpropanoid pathway during seed development. To explore the molecular regulatory mechanism of lignin and flavonoid biosynthesis, T. tuan seeds were harvested at five stages (30, 60, 90, 120, and 150 days after pollination) for lignin and flavonoid analyses. The results showed that lignin accumulated rapidly in the early and middle stages (S1, S3, and S4), and rapid accumulation of flavonoids during the early and late stages (S1 and S5). High-throughput RNA sequencing analysis of developing seeds identified 50,553 transcripts, including 223 phenylpropanoid biosynthetic pathway genes involved in lignin accumulation grouped into 3 clusters, and 106 flavonoid biosynthetic pathway genes (FBPGs) grouped into 2 clusters. Subsequent WGCNA and time-ordered gene co-expression network (TO-GCN) analysis revealed that 24 TFs (e.g., TtARF2 and TtWRKY15) were involved in flavonoids and lignin biosynthesis regulation. The transcriptome data were validated by qRT-PCR to analyze the expression profiles of key enzyme-coding genes. This study revealed that there existed a competitive relationship between flavonoid and lignin biosynthesis pathway during the development of T. tuan seeds, that provide a foundation for the further exploration of molecular mechanisms underlying lignin and flavonoid accumulation in T. tuan seeds.

Haplotype-resolved chromosome-level genome assembly of Ehretia macrophylla

Ehretia macrophylla Wall, known as wild loquat, is an ecologically, economically, and medicinally significant tree species widely grown in China, Japan, Vietnam, and Nepal. In this study, we have successfully generated a haplotype-resolved chromosome-scale genome assembly of E. macrophylla by integrating PacBio HiFi long-reads, Illumina short-reads, and Hi-C data. The genome assembly consists of two haplotypes, with sizes of 1.82 Gb and 1.58 Gb respectively, and contig N50 lengths of 28.11 Mb and 21.57 Mb correspondingly. Additionally, 99.41% of the assembly was successfully anchored into 40 pseudo-chromosomes. We predicted 58,886 protein-coding genes, of which 99.60% were functionally annotated from databases. We furthermore detected 2.65 Gb repeat sequences, 659,290 rRNAs, 4,931 tRNAs and 4,688 other ncRNAs. The high-quality assembly of the genome offers a solid basis for furthering the fields of molecular breeding and functional genomics of E. macrophylla.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

Genomic basis of the distinct biosynthesis of β-glucogallin, a biochemical marker for hydrolyzable tannin production, in three oak species

Hydrolyzable tannins (HTs), predominant polyphenols in oaks, are widely used in grape wine aging, feed additives, and human healthcare. However, the limited availability of a high-quality reference genome of oaks greatly hampered the recognition of the mechanism of HT biosynthesis. Here, high-quality reference genomes of three Asian oak species (Quercus variabilis, Quercus aliena, and Quercus dentata) that have different HT contents were generated. Multi-omics studies were carried out to identify key genes regulating HT biosynthesis. In vitro enzyme activity assay was also conducted. Dual-luciferase and yeast one-hybrid assays were used to reveal the transcriptional regulation. Our results revealed that β-glucogallin was a biochemical marker for HT production in the cupules of the three Asian oaks. UGT84A13 was confirmed as the key enzyme for β-glucogallin biosynthesis. The differential expression of UGT84A13, rather than enzyme activity, was the main reason for different β-glucogallin and HT accumulation. Notably, sequence variations in UGT84A13 promoters led to different trans-activating activities of WRKY32/59, explaining the different expression patterns of UGT84A13 among the three species. Our findings provide three high-quality new reference genomes for oak trees and give new insights into different transcriptional regulation for understanding β-glucogallin and HT biosynthesis in closely related oak species.

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

Metabolite analysis reveals flavonoids accumulation during flower development in Rhododendron pulchrum sweet (Ericaceae)

The azalea (Rhododendron simsii Planch.) is an important ornamental woody plant with various medicinal properties due to its phytochemical compositions and components. However little information on the metabolite variation during flower development in Rhododendron has been provided. In our study, a comparative analysis of the flavonoid profile was performed in Rhododendron pulchrum sweet at three stages of flower development, bud (stage 1), partially open flower (stage 2), and full bloom (stage 3). A total of 199 flavonoids, including flavone, flavonol, flavone C-glycosides, flavanone, anthocyanin, and isoflavone were identified. In hierarchical clustering analysis (HCA) and principal component analysis (PCA), the accumulation of flavonoids displayed a clear development stage variation. During flower development, 78 differential accumulated metabolites (DAMs) were identified, and most were enriched to higher levels at the full bloom stage. A total of 11 DAMs including flavone (chrysin, chrysoeriol O-glucuronic acid, and chrysoeriol O-hexosyl-O-pentoside), isoflavone (biochanin A), and flavonol (3,7-di-O-methyl quercetin and isorhamnetin) were significantly altered at three stages. In particular, 3,7-di-O-methyl quercetin was the top increased metabolite during flower development. Furthermore, integrative analyses of metabolomic and transcriptomic were conducted, revealing that the contents of isoflavone, biochanin A, glycitin, and prunetin were correlated with the expression of 2-hydroxyisoflavanone dehydratase (HIDH), which provide insight into the regulatory mechanism that controls isoflavone biosynthesis in R. pulchrum. This study will provide a new reference for increasing desired metabolites effectively by more accurate or appropriate genetic engineering strategies.

Metabolic and molecular mechanisms of spine color formation in Chinese red chestnut

The spines of Chinese red chestnut are red and the depth of their color gradually increases with maturity. To identify the anthocyanin types and synthesis pathways in red chestnut and to identify the key genes regulating the anthocyanin biosynthesis pathway, we obtained and analyzed the transcriptome and anthocyanin metabolism of red chestnut and its control variety with green spines at 3 different periods. GO and KEGG analyses revealed that photosynthesis was more highly enriched in green spines compared with red spines, while processes related to defense and metabolism regulation were more highly enriched in red spines. The analysis showed that the change in spine color promoted photoprotection in red chestnut, especially at the early growth stage, which resulted in the accumulation of differentially expressed genes involved in the defense metabolic pathway. The metabolome results revealed 6 anthocyanins in red spines. Moreover, red spines exhibited high levels of cyanidin, peonidin and pelargonidin and low levels of delphinidin, petunidin and malvidin. Compared with those in the control group, the levels of cyanidin, peonidin, pelargonidin and malvidin in red spines were significantly increased, indicating that the cyanidin and pelargonidin pathways were enriched in the synthesis of anthocyanins in red spines, whereas the delphinidin pathways were inhibited and mostly transformed into malvidin. During the process of flower pigment synthesis, the expression of the CHS, CHI, F3H, CYP75A, CYP75B1, DFR and ANS genes clearly increased, that of CYP73A decreased obviously, and that of PAL, 4CL and LAR both increased and decreased. Notably, the findings revealed that the synthesized anthocyanin can be converted into anthocyanidin or epicatechin. In red spines, the upregulation of BZ1 gene expression increases the corresponding anthocyanidin content, and the upregulation of the ANR gene also promotes the conversion of anthocyanin to epicatechin. The transcription factors involved in color formation included 4 WRKYs.

Transcriptional Regulation of SugarCane Response to Sporisorium scitamineum: Insights from Time-Course Gene Coexpression and Ca2+ Signaling

Sugarcane response to Sporisorium scitamineum is determined by multiple major genes and numerous microeffector genes. Here, time-ordered gene coexpression networks were applied to explore the interaction between sugarcane and S. scitamineum. Totally, 2459 differentially expressed genes were identified and divided into 10 levels, and several stress-related subnetworks were established. Interestingly, the Ca2+ signaling pathway was activated to establish the response to sugarcane smut disease. Accordingly, two CAX genes (ScCAX2 and ScCAX3) were cloned and characterized from sugarcane. They were significantly upregulated under ABA stress but inhibited by MeJA treatment. Furthermore, overexpression of ScCAX2 and ScCAX3 enhanced the susceptibility of transgenic plants to the pathogen infection, suggesting its negative role in disease resistance. A regulatory model for ScCAX genes in disease response was thus depicted. This work helps to clarify the transcriptional regulation of sugarcane response to S. scitamineum stress and the function of the CAX gene in disease response.

Chromosome-level genome of putative autohexaploid Actinidia deliciosa provides insights into polyploidisation and evolution

Actinidia ('Mihoutao' in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2n = 174) with diploid A. chinensis (2n = 58) as the putative parent. A CCS-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for hexaploid A. deliciosa yielded a 3.91-Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele-defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.

Chromosome-level genome assembly provides insights into the genome evolution and functional importance of the phenylpropanoid-flavonoid pathway in Thymus mongolicus

BACKGROUND

Thymus mongolicus (family Lamiaceae) is a Thyme subshrub with strong aroma and remarkable environmental adaptability. Limited genomic information limits the use of this plant.

RESULTS

Chromosome-level 605.2 Mb genome of T. mongolicus was generated, with 96.28% anchored to 12 pseudochromosomes. The repetitive sequences were dominant, accounting for 70.98%, and 32,593 protein-coding genes were predicted. Synteny analysis revealed that Lamiaceae species generally underwent two rounds of whole genome duplication; moreover, species-specific genome duplication was identified. A recent LTR retrotransposon burst and tandem duplication might play important roles in the formation of the Thymus genome. Using comparative genomic analysis, phylogenetic tree of seven Lamiaceae species was constructed, which revealed that Thyme plants evolved recently in the family. Under the phylogenetic framework, we performed functional enrichment analysis of the genes on nodes that contained the most gene duplication events (> 50% support) and of relevant significant expanded gene families. These genes were highly associated with environmental adaptation and biosynthesis of secondary metabolites. Combined transcriptome and metabolome analyses revealed that Peroxidases, Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferases, and 4-coumarate-CoA ligases genes were the essential regulators of the phenylpropanoid-flavonoid pathway. Their catalytic products (e.g., apigenin, naringenin chalcone, and several apigenin-related compounds) might be responsible for the environmental tolerance and aromatic properties of T. mongolicus.

CONCLUSION

This study enhanced the understanding of the genomic evolution of T. mongolicus, enabling further exploration of its unique traits and applications, and contributed to the understanding of Lamiaceae genomics and evolutionary biology.

Genome and GWAS analysis identified genes significantly related to phenotypic state of Rhododendron bark

As an important horticultural plant, Rhododendron is often used in urban greening and landscape design. However, factors such as the high rate of genetic recombination, frequent outcrossing in the wild, weak linkage disequilibrium, and the susceptibility of gene expression to environmental factors limit further exploration of functional genes related to important horticultural traits, and make the breeding of new varieties require a longer time. Therefore, we choose bark as the target trait which is not easily affected by environmental factors, but also has ornamental properties. Genome-wide association study (GWAS) of Rhododendron delavayi (30 samples), R. irroratum (30 samples) and their F1 generation R. agastum (200 samples) was conducted on the roughness of bark phenotypes. Finally, we obtained 2416.31 Gbp of clean data and identified 5 328 800 high-quality SNPs. According to the P-value and the degree of linkage disequilibrium of SNPs, we further identified 4 out of 11 candidate genes that affect bark roughness. The results of gene differential expression analysis further indicated that the expression levels of Rhdel02G0243600 and Rhdel08G0220700 in different bark phenotypes were significantly different. Our study identified functional genes that influence important horticultural traits of Rhododendron, and illustrated the powerful utility and great potential of GWAS in understanding and exploiting wild germplasm genetic resources of Rhododendron.

Genomics insights into flowering and floral pattern formation: regional duplication and seasonal pattern of gene expression in Camellia

BACKGROUND

The formation and domestication of ornamental traits are influenced by various aspects, such as the recognition of esthetic values and cultural traditions. Camellia japonica is widely appreciated and domesticated around the world mainly due to its rich variations in ornamental traits. Ornamental camellias have a diverse range of resources, including different bud variations from Camellia spp. as well as inter- and intra- specific hybridization. Despite research on the formation of ornamental traits, a basic understanding of their genetics and genomics is still lacking.

RESULTS

Here, we report the chromosomal-level reference genome of C. japonica through combining multiple DNA-sequencing technologies and obtain a high-density genetic linkage map of 4255 markers by sequencing 98 interspecific F1 hybrids between C. japonica and C. chekiangoleosa. We identify two whole-genome duplication events in C. japonica: one is a shared ancient γ event, and the other is revealed to be specific to genus Camellia. Based on the micro-collinearity analysis, we find large-scale segmental duplication of chromosome 8, resulting to two copies of the AGAMOUS loci, which may play a key role in the domestication of floral shapes. To explore the regulatory mechanisms of seasonal flowering, we have analyzed year-round gene expression patterns of C. japonica and C. azalea-a sister plant of continuous flowering that has been widely used for cross breeding. Through comparative analyses of gene co-expression networks and annual gene expression patterns, we show that annual expression rhythms of some important regulators of seasonal growth and development, including GIGANTEA and CONSTANS of the photoperiod pathway, have been disrupted in C. azalea. Furthermore, we reveal that the distinctive expression patterns of FLOWERING LOCUS T can be correlated with the seasonal activities of flowering and flushing. We demonstrate that the regulatory module involved in GIGANTEA, CONSTANS, and FLOWERING LOCUS T is central to achieve seasonality.

CONCLUSIONS

Through the genomic and comparative genomics characterizations of ornamental Camellia spp., we propose that duplication of chromosomal segments as well as the establishment of gene expression patterns has played a key role in the formation of ornamental traits (e.g., flower shape, flowering time). This work provides a valuable genomic platform for understanding the molecular basis of ornamental traits.

Differential gene expression and potential regulatory network of fatty acid biosynthesis during fruit and leaf development in yellowhorn (Xanthoceras sorbifolium), an oil-producing tree with significant deployment values

Xanthoceras sorbifolium (yellowhorn) is a woody oil plant with super stress resistance and excellent oil characteristics. The yellowhorn oil can be used as biofuel and edible oil with high nutritional and medicinal value. However, genetic studies on yellowhorn are just in the beginning, and fundamental biological questions regarding its very long-chain fatty acid (VLCFA) biosynthesis pathway remain largely unknown. In this study, we reconstructed the VLCFA biosynthesis pathway and annotated 137 genes encoding relevant enzymes. We identified four oleosin genes that package triacylglycerols (TAGs) and are specifically expressed in fruits, likely playing key roles in yellowhorn oil production. Especially, by examining time-ordered gene co-expression network (TO-GCN) constructed from fruit and leaf developments, we identified key enzymatic genes and potential regulatory transcription factors involved in VLCFA synthesis. In fruits, we further inferred a hierarchical regulatory network with MYB-related (XS03G0296800) and B3 (XS02G0057600) transcription factors as top-tier regulators, providing clues into factors controlling carbon flux into fatty acids. Our results offer new insights into key genes and transcriptional regulators governing fatty acid production in yellowhorn, laying the foundation for efforts to optimize oil content and fatty acid composition. Moreover, the gene expression patterns and putative regulatory relationships identified here will inform metabolic engineering and molecular breeding approaches tailored to meet biofuel and bioproduct demands.

The first high-altitude autotetraploid haplotype-resolved genome assembled (Rhododendron nivale subsp. boreale) provides new insights into mountaintop adaptation

BACKGROUND

Rhododendron nivale subsp. boreale Philipson et M. N. Philipson is an alpine woody species with ornamental qualities that serve as the predominant species in mountainous scrub habitats found at an altitude of ∼4,200 m. As a high-altitude woody polyploid, this species may serve as a model to understand how plants adapt to alpine environments. Despite its ecological significance, the lack of genomic resources has hindered a comprehensive understanding of its evolutionary and adaptive characteristics in high-altitude mountainous environments.

FINDINGS

We sequenced and assembled the genome of R. nivale subsp. boreale, an assembly of the first subgenus Rhododendron and the first high-altitude woody flowering tetraploid, contributing an important genomic resource for alpine woody flora. The assembly included 52 pseudochromosomes (scaffold N50 = 42.93 Mb; BUSCO = 98.8%; QV = 45.51; S-AQI = 98.69), which belonged to 4 haplotypes, harboring 127,810 predicted protein-coding genes. Conjoint k-mer analysis, collinearity assessment, and phylogenetic investigation corroborated autotetraploid identity. Comparative genomic analysis revealed that R. nivale subsp. boreale originated as a neopolyploid of R. nivale and underwent 2 rounds of ancient polyploidy events. Transcriptional expression analysis showed that differences in expression between alleles were common and randomly distributed in the genome. We identified extended gene families and signatures of positive selection that are involved not only in adaptation to the mountaintop ecosystem (response to stress and developmental regulation) but also in autotetraploid reproduction (meiotic stabilization). Additionally, the expression levels of the (group VII ethylene response factor transcription factors) ERF VIIs were significantly higher than the mean global gene expression. We suspect that these changes have enabled the success of this species at high altitudes.

CONCLUSIONS

We assembled the first high-altitude autopolyploid genome and achieved chromosome-level assembly within the subgenus Rhododendron. In addition, a high-altitude adaptation strategy of R. nivale subsp. boreale was reasonably speculated. This study provides valuable data for the exploration of alpine mountaintop adaptations and the correlation between extreme environments and species polyploidization.

Populus cathayana genome and population resequencing provide insights into its evolution and adaptation

Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype-environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

A Genome-Wide Analysis and Expression Profile of Heat Shock Transcription Factor (Hsf) Gene Family in Rhododendron simsii

Heat shock transcription factors are key players in a number of transcriptional regulatory pathways that function during plant growth and development. However, their mode of action in Rhododendron simsii is still unclear. In this study, 22 RsHsf genes were identified from genomic data of R. simsii. The 22 genes were randomly distributed on 12 chromosomes, and were divided into three major groups according to their phylogenetic relationships. The structures and conserved motifs were predicted for the 22 genes. Analysis of cis-acting elements revealed stress-responsive and phytohormone-responsive elements in the gene promoter regions, but the types and number varied among the different groups of genes. Transcriptional profile analyses revealed that RsHsfs were expressed in a tissue-specific manner, with particularly high transcript levels in the roots. The transcriptional profiles under abiotic stress were detected by qRT-PCR, and the results further validated the critical function of RsHsfs. This study provides basic information about RsHsf family in R. simsii, and paves the way for further research to clarify their precise roles and to breed new stress-tolerant varieties.

Comparative transcriptome analysis of dioecious floral development in Trachycarpus fortunei using Illumina and PacBio SMRT sequencing

BACKGROUND

Trachycarpus fortunei is a plant with significant economic and ornamental value. Both male and female flowers of T. fortunei originate as bisexual flowers, and selective abortion occurs during floral development. However, the regulatory mechanisms underlying this process remain unclear in T. fortunei. In this study, transcriptome sequencing with Illumina and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) platforms were used to investigate gene expression differences between male and female T. fortunei plants.

RESULTS

A total of 833,137 full-length non-chimeric (FLNC) reads were obtained, and 726,846 high-quality full-length transcripts were identified. A total of 159 genes were differentially expressed between male and female flowers at all development stages. Some of the differentially expressed genes (DEGs) showed male bias, including serine/threonine-protein kinase (STPK), THUMP1 homolog and other genes. Through single-nucleotide polymorphisms(SNPs) identification, 28 genes were considered as potential sex-associated SNPs. Time-Ordered Gene Co-expression Network (TO-GCN) analysis revealed that MADS2 and MADS26 may play important roles in the development of female and male flowers T. fortune plants, respectively.

CONCLUSIONS

These findings provide a genetic basis for flower development and differentiation in T. fortunei, and improve our understanding of the mechanisms underlying sexual differentiation in T. fortunei.

A chromosome-level genome assembly provides insights into Cornus wilsoniana evolution, oil biosynthesis, and floral bud development

Cornus wilsoniana W. is a woody oil plant with high oil content and strong hypolipidemic effects, making it a valuable species for medicinal, landscaping, and ecological purposes in China. To advance genetic research on this species, we employed PacBio together with Hi-C data to create a draft genome assembly for C. wilsoniana. Based on an 11-chromosome anchored chromosome-level assembly, the estimated genome size was determined to be 843.51 Mb. The N50 contig size and N50 scaffold size were calculated to be 4.49 and 78.00 Mb, respectively. Furthermore, 30 474 protein-coding genes were annotated. Comparative genomics analysis revealed that C. wilsoniana diverged from its closest species ~12.46 million years ago (Mya). Furthermore, the divergence between Cornaceae and Nyssaceae occurred >62.22 Mya. We also found evidence of whole-genome duplication events and whole-genome triplication γ, occurring at ~44.90 and 115.86 Mya. We further inferred the origins of chromosomes, which sheds light on the complex evolutionary history of the karyotype of C. wilsoniana. Through transcriptional and metabolic analysis, we identified two FAD2 homologous genes that may play a crucial role in controlling the oleic to linoleic acid ratio. We further investigated the correlation between metabolites and genes and identified 33 MADS-TF homologous genes that may affect flower morphology in C. wilsoniana. Overall, this study lays the groundwork for future research aimed at identifying the genetic basis of crucial traits in C. wilsoniana.

Chromosome-scale genome assembly provides insights into flower coloration mechanisms of Canna indica

Canna indica, which produce conspicuous and colorful flowers, are widely appreciated as ornamental plants. We used Pacific Biosciences sequencing (PacBio) and chromosome conformation capture (Hi-C) genome scaffolding to build a high-quality chromosome-scale genome assembly of C. indica and the genome assembly was 821Mb with a contig N50 of 48Mb assembled into nine chromosomes. The genome of C. indica was predicted to contain 31,130 genes and 30,816 genes were functionally annotated. Genome annotation identified 522 Mb (63.59 %) as repetitive sequences. Genome evolution analysis showed that whole-genome duplication occurred 53.4 million years ago. Transcriptome analysis revealed that petal coloration was linked with the expression of genes encoding enzymes involved in anthocyanin biosynthesis, carotenoid biosynthesis, and the methylerythritol phosphate (MEP) pathway. Furthermore, modules of co-expressed genes and hub genes were identified via weighted gene co-expression network analysis. These results suggested that, in Canna indica, deep red petal coloration was regulated by CHS2 and yellow petal coloration was associated with expression of ARF6 and NAC14. Considered together, the current study revealed a high-quality reference genome which may provide new insights into the molecular basis of flower coloration in Canna indica and help enhance the conservation and breeding of ornamental plants in general.

Genome-Wide Analysis of MBF1 Family Genes in Five Solanaceous Plants and Functional Analysis of SlER24 in Salt Stress

Multiprotein bridging factor 1 (MBF1) is an ancient family of transcription coactivators that play a crucial role in the response of plants to abiotic stress. In this study, we analyzed the genomic data of five Solanaceae plants and identified a total of 21 MBF1 genes. The expansion of MBF1a and MBF1b subfamilies was attributed to whole-genome duplication (WGD), and the expansion of the MBF1c subfamily occurred through transposed duplication (TRD). Collinearity analysis within Solanaceae species revealed collinearity between members of the MBF1a and MBF1b subfamilies, whereas the MBF1c subfamily showed relative independence. The gene expression of SlER24 was induced by sodium chloride (NaCl), polyethylene glycol (PEG), ABA (abscisic acid), and ethrel treatments, with the highest expression observed under NaCl treatment. The overexpression of SlER24 significantly enhanced the salt tolerance of tomato, and the functional deficiency of SlER24 decreased the tolerance of tomato to salt stress. SlER24 enhanced antioxidant enzyme activity to reduce the accumulation of reactive oxygen species (ROS) and alleviated plasma membrane damage under salt stress. SlER24 upregulated the expression levels of salt stress-related genes to enhance salt tolerance in tomato. In conclusion, this study provides basic information for the study of the MBF1 family of Solanaceae under abiotic stress, as well as a reference for the study of other plants.

Genome assembly and resequencing analyses provide new insights into the evolution, domestication and ornamental traits of crape myrtle

Crape myrtle (Lagerstroemia indica) is a globally used ornamental woody plant and is the representative species of Lagerstroemia. However, studies on the evolution and genomic breeding of L. indica have been hindered by the lack of a reference genome. Here we assembled the first high-quality genome of L. indica using PacBio combined with Hi-C scaffolding to anchor the 329.14-Mb genome assembly into 24 pseudochromosomes. We detected a previously undescribed independent whole-genome triplication event occurring 35.5 million years ago in L. indica following its divergence from Punica granatum. After resequencing 73 accessions of Lagerstroemia, the main parents of modern crape myrtle cultivars were found to be L. indica and L. fauriei. During the process of domestication, genetic diversity tended to decrease in many plants, but this was not observed in L. indica. We constructed a high-density genetic linkage map with an average map distance of 0.33 cM. Furthermore, we integrated the results of quantitative trait locus (QTL) using genetic mapping and bulk segregant analysis (BSA), revealing that the major-effect interval controlling internode length (IL) is located on chr1, which contains CDL15, CRG98, and GID1b1 associated with the phytohormone pathways. Analysis of gene expression of the red, purple, and white flower-colour flavonoid pathways revealed that differential expression of multiple genes determined the flower colour of L. indica, with white flowers having the lowest gene expression. In addition, BSA of purple- and green-leaved individuals of populations of L. indica was performed, and the leaf colour loci were mapped to chr12 and chr17. Within these intervals, we identified MYB35, NCED, and KAS1. Our genome assembly provided a foundation for investigating the evolution, population structure, and differentiation of Myrtaceae species and accelerating the molecular breeding of L. indica.

Genome-wide identification and bioinformatics analysis of the WD40 transcription factor family and candidate gene screening for anthocyanin biosynthesis in Rhododendron simsii

WD40 transcription factors (TFs) constitute a large gene family in eukaryotes, playing diverse roles in cellular processes. However, their functions in the major ornamental plant, Rhododendron simsii, remain poorly understood. In this study, we identified 258 WD40 proteins in the R. simsii genome, which exhibited an uneven distribution across chromosomes. Based on domain compositions and phylogenetic analysis, we classified these 258 RsWD40 proteins into 42 subfamilies and 47 clusters. Comparative genomic analysis suggested that the expansion of the WD40 gene family predates the divergence of green algae and higher plants, indicating an ancient origin. Furthermore, by analyzing the duplication patterns of RsWD40 genes, we found that transposed duplication played a major role in their expansion. Notably, the majority of RsWD40 gene duplication pairs underwent purifying selection during evolution. Synteny analysis identified significant orthologous gene pairs between R. simsii and Arabidopsis thaliana, Oryza sativa, Vitis vinifera, and Malus domestica. We also investigated potential candidate genes involved in anthocyanin biosynthesis during different flower development stages in R. simsii using RNA-seq data. Specifically, we identified 10 candidate genes during the bud stage and 7 candidate genes during the full bloom stage. GO enrichment analysis of these candidate genes revealed the potential involvement of the ubiquitination process in anthocyanin biosynthesis. Overall, our findings provide a valuable foundation for further investigation and functional analysis of WD40 genes, as well as research on the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species.

Insight into the Molecular Mechanism of Flower Color Regulation in Rhododendron latoucheae Franch: A Multi-Omics Approach

Rhododendron latoucheae Franch. (R. latoucheae) is a valuable woody plant known for its high ornamental value. While purple flowers are a distinct and attractive variant phenotype of R. latoucheae, the underlying mechanism regulating its flower color is still poorly understood. To investigate the molecular regulatory mechanism responsible for the variation in flower color, we selected plants with white-pink and purple petals as the object and conducted analyses of metabolites, key genes, and transcription factors associated with flower color. A combined metabolome-transcriptome analysis was performed, and the expression of key genes was subsequently verified through qRT-PCR experiments. The results of our study demonstrated a significant enrichment of differential metabolites in the flavonoid metabolic pathway. Changes in anthocyanin content followed the same trend as the observed flower color variations, specifically showing significant correlations with the contents of malvidin-3-O-glucoside, dihydromyricetin, gallocatechin, and peonidin-3-O-glucoside. Furthermore, we identified three key structural genes (F3GT1, LAR, ANR) and four transcription factors (bHLH130, bHLH41, bHLH123, MYB4) that are potentially associated with the biosynthesis of flavonoid compounds, thereby influencing the appearance of purple flower color in R. latoucheae.

Transcriptomic Analysis Reveals the Regulatory Mechanism of Color Diversity in Rhododendron pulchrum Sweet (Ericaceae)

Rhododendron pulchrum Sweet is a traditional ornamental plant cultivated in China and presents a great variation in petal coloration. However, few studies have been performed to reveal the genes involved and the regulatory mechanism of flower color formation in this plant. In this study, to explore the underlying genetic basis of flower color formation, transcriptome analysis was performed by high-throughput sequencing techniques on four petal samples of different colors: purple, pink, light pink, and white. Results show that a total of 35.55 to 40.56 million high-quality clean reads were obtained, of which 28.56 to 32.65 million reads were mapped to the reference genome. For their annotation, 28,273, 18,054, 24,301, 19,099, and 11,507 genes were allocated to Nr, Swiss-Prot, Pfam, GO, and KEGG databases, correspondingly. There were differentially expressed genes among the four different petal samples, including signal-transduction-related genes, anthocyanin biosynthesis genes, and transcription factors. We found that the higher expressed levels of genes associated with flavonol synthase (FLS) might be the key to white formation, and the formation of red color may be related to the higher expression of flavanone 4-reductase (DFR) families. Overall, our study provides some valuable information for exploring and understanding the flower color intensity variation in R. pulchrum.

A haplotype-resolved genome assembly of Rhododendron vialii based on PacBio HiFi reads and Hi-C data

Rhododendron vialii (subgen. Azaleastrum) is an evergreen shrub with high ornamental value. This species has been listed as a plant species with extremely small populations (PSESP) for urgent protection by China's Yunnan provincial government in 2021, due to anthropogenic habitat fragmentation. However, limited genomic resources hinder scientifically understanding of genetic threats that the species is currently facing. In this study, we assembled a high-quality haplotype-resolved genome of R. vialii based on PacBio HiFi long reads and Hi-C reads. The assembly contains two haploid genomes with sizes 532.73 Mb and 521.98 Mb, with contig N50 length of 35.67 Mb and 34.70 Mb, respectively. About 99.92% of the assembled sequences could be anchored to 26 pseudochromosomes, and 14 gapless assembled chromosomes were included in this assembly. Additionally, 60,926 protein-coding genes were identified, of which 93.82% were functionally annotated. This is the first reported genome of R. vialii, and hopefully it will lay the foundations for further research into the conservation genomics and horticultural domestication of this ornamentally important species.

Time-series transcriptome provides insights into the gene regulation network involved in the icariin-flavonoid metabolism during the leaf development of Epimedium pubescens

Herba Epimedii (Epimedium) leaves are rich in prenylated flavonol glycosides (PFGs) with high medicinal value. However, the dynamics and regulatory network of PFG biosynthesis remain largely unclear. Here, we combined metabolite profiling (targeted to PFGs) and a high-temporal-resolution transcriptome to elucidate PFGs' regulatory network in Epimedium pubescens and identified key candidate structural genes and transcription factors (TFs) involved in PFG accumulation. Chemical profile analysis revealed that PFG content was quite different between buds and leaves and displayed a continuous decline with leaf development. The structural genes are the determinant reasons, and they are strictly regulated by TFs under temporal cues. We further constructed seven time-ordered gene co-expression networks (TO-GCNs) of PFG biosynthesis genes (including EpPAL2, EpC4H, EpCHS2, EpCHI2, EpF3H, EpFLS3, and EpPT8), and three flavonol biosynthesis routines were then predicted. The TFs involved in TO-GCNs were further confirmed by WGCNA analysis. Fourteen hub genes, comprising 5 MYBs, 1 bHLH, 1 WD40, 2 bZIPs, 1 BES1, 1 C2H2, 1 Trihelix, 1 HD-ZIP, and 1 GATA were identified as candidate key TFs. The results were further validated by TF binding site (TFBS) analysis and qRT-PCR. Overall, these findings provide valuable information for understanding the molecular regulatory mechanism of PFGs biosynthesis, enriching the gene resources, which will guide further research on PFG accumulation in Epimedium.

Ploidy variation in Rhododendron subsection Maddenia and its implications for conservation

Polyploidy, which is common in plants, can confound taxon recognition and hence conservation assessments. In the taxonomically complex genus Rhododendron, 25 % of the over 1,300 taxa are considered under threat and 27 % Near Threatened or Data Deficient, with their taxonomy needing to be resolved urgently. Although ploidy levels of Rhododendron taxa range from diploid (2x) to dodecaploid (12x) according to previous reports, the extent of polyploidy across the genus has not been examined. We first summarized the taxonomic distribution of polyploids in the genus based on the literature. Then as a case study, we estimated ploidy levels of 47 taxa in subsection Maddenia (subgenus Rhododendron, section Rhododendron) using flow cytometry, together with verification of meiotic chromosome counts for representative taxa. The summary of reported ploidy in Rhododendron indicates that polyploidy is most common in subgenera Pentanthera and Rhododendron. In subsection Maddenia, all examined taxa are diploids except for the R. maddenii complex that shows a high ploidy variation (2-8x, 12x). We investigated ploidy level of 12 taxa in subsection Maddenia for the first time, and estimated genome sizes of two Rhododendron species. Knowledge of ploidy levels will inform phylogenetic analysis of unresolved species complexes. Overall, our study of subsection Maddenia provides a model for examining multiple issues including taxonomic complexity, ploidy variation and geographic distribution in relation to biodiversity conservation.

Chromosome-scale genome assembly and insights into the metabolome and gene regulation of leaf color transition in an important oak species, Quercus dentata

Quercus dentata Thunb., a dominant forest tree species in northern China, has significant ecological and ornamental value due to its adaptability and beautiful autumn coloration, with color changes from green to yellow into red resulting from the autumnal shifts in leaf pigmentation. However, the key genes and molecular regulatory mechanisms for leaf color transition remain to be investigated. First, we presented a high-quality chromosome-scale assembly for Q. dentata. This 893.54 Mb sized genome (contig N50 = 4.21 Mb, scaffold N50 = 75.55 Mb; 2n = 24) harbors 31 584 protein-coding genes. Second, our metabolome analyses uncovered pelargonidin-3-O-glucoside, cyanidin-3-O-arabinoside, and cyanidin-3-O-glucoside as the main pigments involved in leaf color transition. Third, gene co-expression further identified the MYB-bHLH-WD40 (MBW) transcription activation complex as central to anthocyanin biosynthesis regulation. Notably, transcription factor (TF) QdNAC (QD08G038820) was highly co-expressed with this MBW complex and may regulate anthocyanin accumulation and chlorophyll degradation during leaf senescence through direct interaction with another TF, QdMYB (QD01G020890), as revealed by our further protein-protein and DNA-protein interaction assays. Our high-quality genome assembly, metabolome, and transcriptome resources further enrich Quercus genomics and will facilitate upcoming exploration of ornamental values and environmental adaptability in this important genus.

Genome-wide identification and bioinformatics analysis of the WRKY transcription factors and screening of candidate genes for anthocyanin biosynthesis in azalea (Rhododendron simsii)

WRKY transcription factors have been demonstrated to influence the anthocyanin biosynthesis in many plant species. However, there is limited knowledge about the structure and function of WRKY genes in the major ornamental plant azalea (Rhododendron simsii). In this study, we identified 57 RsWRKY genes in the R. simsii genome and classified them into three main groups and several subgroups based on their structural and phylogenetic characteristics. Comparative genomic analysis suggested WRKY gene family has significantly expanded during plant evolution from lower to higher species. Gene duplication analysis indicated that the expansion of the RsWRKY gene family was primarily due to whole-genome duplication (WGD). Additionally, selective pressure analysis (Ka/Ks) suggested that all RsWRKY duplication gene pairs underwent purifying selection. Synteny analysis indicated that 63 and 24 pairs of RsWRKY genes were orthologous to Arabidopsis thaliana and Oryza sativa, respectively. Furthermore, RNA-seq data was used to investigate the expression patterns of RsWRKYs, revealing that 17 and 9 candidate genes may be associated with anthocyanin synthesis at the bud and full bloom stages, respectively. These findings provide valuable insights into the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species and lay the foundation for future functional studies of WRKY genes.