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.

From Floral Induction to Blooming: The Molecular Mysteries of Flowering in Woody Plants

Flowering is a pivotal developmental process in response to the environment and determines the start of a new life cycle in plants. Woody plants usually possess a long juvenile nonflowering phase followed by an adult phase with repeated flowering cycles. The molecular mechanism underlying flowering regulation in woody plants is believed to be much more complex than that in annual herbs. In this review, we briefly describe the successive but distinct flowering processes in perennial trees, namely the vegetative phase change, the floral transition, floral organogenesis, and final blooming, and summarize in detail the most recent advances in understanding how woody plants regulate flowering through dynamic gene expression. Notably, the florigen gene FLOWERING LOCUS T(FT) and its antagonistic gene TERMINAL FLOWER 1 (TFL1) seem to play a central role in various flowering transition events. Flower development in different taxa requires interactions between floral homeotic genes together with AGL6 conferring floral organ identity. Finally, we illustrate the issues and corresponding measures of flowering regulation investigation. It is of great benefit to the future study of flowering in perennial trees.

The chromosome-scale genome assembly, annotation and evolution of Rhododendron henanense subsp. lingbaoense

Rhododendron henanense subsp. lingbaoense (hereafter referred to as R. henanense) is an endemic species naturally distributed in the Henan province, China, with high horticultural, ornamental and medicinal value. Herein, we report a de novo genome assembly for R. henanense using a combination of PacBio long read and Illumina short read sequencing technologies. In total, we assembled 634.07 Mb with a contig N50 of 2.5 Mb, representing ~96.93% of the estimated genome size. By applying Hi-C data, 13 pseudochromosomes of R. henanense genome were assembled, covering ~98.21% of the genome assembly. The genome was composed of ~65.76% repetitive sequences and 31,098 protein-coding genes, 88.77% of which could be functionally annotated. Rhododendron henanense displayed a high level of synteny with other Rhododendron species from the Hymenanthes subgenus. Our data also suggests that R. henanense genes related to stress responses have undergone expansion, which may underly the unique abiotic and biotic stress resistance of the species. This alpine Rhododendron chromosome-scale genome assembly provides fundamental molecular resources for germplasm conservation, breeding efforts, evolutionary studies, and elucidating the unique biological characteristics of R. henanense.

Genome-wide transcriptional analysis unveils the molecular basis of organ-specific expression of isosteroidal alkaloids biosynthesis in critically endangered Fritillaria roylei Hook

Fritillaria roylei Hook. is a critically endangered high altitude Himalayan medicinal plant species with rich source of pharmaceutically active structurally diverse steroidal alkaloids. Nevertheless, except few marker compounds, the chemistry of the plant remains unexplored. Therefore, in the current study, transcriptome sequencing efforts were made to elucidate isosteroidal alkaloids biosynthesis by creating first organ-specific genomic resource using bulb, stem, and leaf tissues derived from natural populations of Indian Himalayan region. Overall, 349.9 million high quality paired-end reads obtained using NovaSeq 6000 platform were assembled (de novo) into 82,848 unigenes and 31,061 isoforms. Functional annotation and organ specific differential expression (DE) analysis identified 2488 significant DE transcripts, grouped into three potential sub-clusters (sub-cluster I: 728 transcripts; sub-cluster II: 446 transcripts and Sub-cluster III: 1314 transcripts). Subsequently, pathway enrichment (GO, KEGG) and protein-protein network analysis revealed significantly higher enrichment of phenyl-propanoid and steroid backbone including terpenoid, sesquiterpenoid and triterpenoid biosynthesis in bulb. Additionally, upregulated expression of cytochrome P450, UDP-dependent Glucuronosyltransferase families and key transcription factor families (FAR1, bHLH, GRAS, C2H2, TCP and MYB) suggests 'bulb' as a primary site of MVA mediated isosteroidal alkaloids biosynthesis. The comprehensive elucidation of molecular insights in this study is a first step towards the understanding of isosteroidal alkaloid biosynthesis pathway in F. roylei. Furthermore, key genes and regulators identified here can facilitate metabolic engineering of potential bioactive compounds at industrial scale.

The Rhododendron Plant Genome Database (RPGD): a comprehensive online omics database for Rhododendron

Background

The genus Rhododendron L. has been widely cultivated for hundreds of years around the world. Members of this genus are known for great ornamental and medicinal value. Owing to advances in sequencing technology, genomes and transcriptomes of members of the Rhododendron genus have been sequenced and published by various laboratories. With increasing amounts of omics data available, a centralized platform is necessary for effective storage, analysis, and integration of these large-scale datasets to ensure consistency, independence, and maintainability.

Results

Here, we report our development of the Rhododendron Plant Genome Database (RPGD; http://bioinfor.kib.ac.cn/RPGD/), which represents the first comprehensive database of Rhododendron genomics information. It includes large amounts of omics data, including genome sequence assemblies for R. delavayi, R. williamsianum, and R. simsii, gene expression profiles derived from public RNA-Seq data, functional annotations, gene families, transcription factor identification, gene homology, simple sequence repeats, and chloroplast genome. Additionally, many useful tools, including BLAST, JBrowse, Orthologous Groups, Genome Synteny Browser, Flanking Sequence Finder, Expression Heatmap, and Batch Download were integrated into the platform.

Conclusions

RPGD is designed to be a comprehensive and helpful platform for all Rhododendron researchers. Believe that RPGD will be an indispensable hub for Rhododendron studies.

Transcriptome characterization and generation of marker resource for Himalayan vulnerable species, Ulmus wallichiana

Ulmus wallichiana is a traditional medicinal plant listed as a vulnerable in the IUCN red list data. Genomic and transcriptomic resources for this species are lacking, hindering its genetic exploration. Further, no polymorphic marker resource is available for this species, thus limiting the elucidation of its underlying genetic diversity, which is a pre-requisite for its conservation. This study was therefore aimed to generate a functionally annotated transcriptomic resource and screen it for SSR regions. We used paired-end Illumina based RNAseq technology and trinity based de novo assembly approach to generate full length transcripts, which were screened for SSR regions and functionally annotated. Around 6.6 million raw reads were de novo assembled transcripts, which were clustered into 146,083 unigenes. 19,909 transcripts were provided with 3986 unique KEGG ids, 70,519 transcripts with 6621 unique Pfam domains, and 45,125 transcripts with 7302 unique INTERPRO domains. 1456 transcripts were identified as transcriptions factors (TFs). Further, 8868 unique GO terms were obtained for the unigenes. The transcripts mapped to 23,056 known pre-determined orthology clusters in the eggNOG database. A total of 16,570 SSRs were identified from the unigenes. Out of the 90 SSRs selected for characterization on 20 genotypes, 28 were polymorphic. Mean effective alleles (Ne) of 2.53, mean observed heterozygosity (Ho) of 0.77, and average polymorphic information content (PIC) of 0.57 were found. This study may facilitate the genetic exploration of this species. The polymorphic SSRs would prove useful to explore its genetic diversity patterns, required for its conservation.

Pharmaceutical resource discovery from traditional medicinal plants: Pharmacophylogeny and pharmacophylogenomics

The worldwide botanical and medicinal culture diversity are astonishing and constitute a Pierian spring for innovative drug R&D. Here, the latest awareness and the perspectives of pharmacophylogeny and pharmacophylogenomics, as well as their expanding utility in botanical drug R&D, are systematically summarized and highlighted. Chemotaxonomy is based on the fact that closely related plants contain the same or similar chemical profiles. Correspondingly, it is better to combine morphological characters, DNA markers and chemical markers in the inference of medicinal plant phylogeny. Medicinal plants within the same phylogenetic groups may have the same or similar therapeutic effects, thus forming the core of pharmacophylogeny. Here we systematically review and comment on the versatile applications of pharmacophylogeny in (1) looking for domestic resources of imported drugs, (2) expanding medicinal plant resources, (3) quality control, identification and expansion of herbal medicines, (4) predicting the chemical constituents or active ingredients of herbal medicine and assisting in the identification and determination of chemical constituents, (5) the search for new drugs sorting out, and (6) summarizing and improving herbal medicine experiences, etc. Such studies should be enhanced within the context of deeper investigations of molecular biology and genomics of traditional medicinal plants, phytometabolites and metabolomics, and ethnomedicine-based pharmacological activity, thus enabling the sustainable conservation and utilization of traditional medicinal resources.

The Rhododendron Genome and Chromosomal Organization Provide Insight into Shared Whole-Genome Duplications across the Heath Family (Ericaceae)

The genus Rhododendron (Ericaceae), which includes horticulturally important plants such as azaleas, is a highly diverse and widely distributed genus of >1,000 species. Here, we report the chromosome-scale de novo assembly and genome annotation of Rhododendron williamsianum as a basis for continued study of this large genus. We created multiple short fragment genomic libraries, which were assembled using ALLPATHS-LG. This was followed by contiguity preserving transposase sequencing (CPT-seq) and fragScaff scaffolding of a large fragment library, which improved the assembly by decreasing the number of scaffolds and increasing scaffold length. Chromosome-scale scaffolding was performed by proximity-guided assembly (LACHESIS) using chromatin conformation capture (Hi-C) data. Chromosome-scale scaffolding was further refined and linkage groups defined by restriction-site associated DNA (RAD) sequencing of the parents and progeny of a genetic cross. The resulting linkage map confirmed the LACHESIS clustering and ordering of scaffolds onto chromosomes and rectified large-scale inversions. Assessments of the R. williamsianum genome assembly and gene annotation estimate them to be 89% and 79% complete, respectively. Predicted coding sequences from genome annotation were used in syntenic analyses and for generating age distributions of synonymous substitutions/site between paralgous gene pairs, which identified whole-genome duplications (WGDs) in R. williamsianum. We then analyzed other publicly available Ericaceae genomes for shared WGDs. Based on our spatial and temporal analyses of paralogous gene pairs, we find evidence for two shared, ancient WGDs in Rhododendron and Vaccinium (cranberry/blueberry) members that predate the Ericaceae family and, in one case, the Ericales order.

Development of SSR markers via de novo transcriptome assembly in Akebia trifoliata (Thunb.) Koidz

Owing to its high nutritive, economic, and medicinal values, Akebia trifoliata has received increased attention, making worthy of being used as a new fruit crop for further domestication and commercialization in China. However, molecular research of A. trifoliata has lagged as investigations of its genomic resources and molecular markers are rare. In this study, a cDNA library of A. trifoliata leaves was sequenced using the Illumina NovaSeq. 6000 sequencing system. In total, 101 417 transcripts, 63 757 unigenes, and 9494 simple sequence repeats were assembled and identified from the transcriptome datasets. The majority of the SSRs were di- and trinucleotide repeats. Length and number of SSR motifs ranged from 15 to 66, and 5 to 48 bp, respectively. Of which, the A/T mononucleotide motif and AG/TC and CT/GA dinucleotide motifs were the most abundant. Furthermore, 100 SSR primers were randomly selected to validate amplification and polymorphism, and 88 A. trifoliata accessions were definitively distinguished by 49 primers. With the Qinling mountains and Huaihe River line as the boundaries, the northern and southern accessions were clustered into different groups, but no clear geographical patterns (city or origin) were observed in the southern accessions. These newly identified molecular markers may provide a foundation for the genetic identification and diversity analysis and marker-assisted selection breeding in species of Akebia.

Comparative analysis of microsatellite, SNP, and InDel markers in four Rhododendron species based on RNA-seq

Rhododendron possesses valuable horticultural and medicinal properties. However, the genetic studies have been hindered due to the lack of genetic markers. Based on RNA-seq, large-scale molecular markers were developed from four Rhododendron species endemic to Dabie Mountains (central China): R. fortunei (5.25 Gb; SSRs, 12,756, one/2.37 kb, 147 types; SNPs, 38,313; InDels, 3,174), R. simsii (5.80 Gb; SSRs, 13,294, one/2.58 kb, 167 types; SNPs, 136,590; InDels, 6,258), R. mariesii (6.53 Gb; SSRs, 15,724, one/2.51 kb, 170 types; SNPs, 44,942; InDels, 4,126), and R. molle (4.35 Gb; SSRs, 10,214, one/2.49 kb, 110 types; SNPs, 77,829; InDels, 3,416). Di-nucleotide repeats were the main type (59.126%-64.314%), and AG/CT repeat (55.18%-61.22%) was the most. In particular, 89 species-specific types had been found. Furthermore, C:G→T:A mutation was the main SNP type (30.475%-34.99%). However, C:G→G:C mutation was the least type in R. fortunei, while T:A→G:C mutation was the least in the other three species. InDels with length of 3 nt was most in R. fortunei, but 1 nt InDels were the main type in the other three species. Twelve microsatellite markers developed from R. simsii reveled high genetic diversity in the four populations, and heterozygote excess was observed. This research would benefit the genetic study, molecular marker-assisted selection, and breeding studies in Rhododendron species.

Transcriptome analysis and identification of genes associated with flower development in Rhododendron pulchrum Sweet (Ericaceae)

Flowering process is essential for plant development. However, the molecular mechanisms driving flower development of ornamental woody Rhododendron pulchrum Sweet are difficult to elucidate due to the lack of genomic data. In this research, high-throughput sequencing and comparative transcriptome analyses of R. pulchrum flowers collected at three key stages were performed: floral bud stage, early flowering stage, and full-flowering stage. Furthermore, expression of genes involved in flower development was also validated with quantitative real-time PCR (qRT-PCR). RNA-seq yielded 96,350,697 bp of clean reads, which were assembled into 98,610 unigenes with an average length of 717 bp. 58,279 (59.10%) unigenes could be annotated, including 324 major unigenes associated with floral development. In addition, ten modules (20,443 mRNAs) were dissected in the co-expression network. Especially, Flowering Locus (FLC) and Flowering Locus T (FT) were co-expressed. 9493 differentially expressed genes (DEGs) were scanned among three stages, and most DEGs existed between flower bud stage and early flowering stage. In particular, 79 DGEs associated with flowering process were enriched in 28 GO terms. Moreover, the expression levels of MYC2, EIN3, and ARR-B were all lowest at early flowering stage, while transcripts of MYC2, TIR1, CYCD3, COL-1, and EIN3 were all peaked at flower bud stage. Transcriptome profile presented here will benefit deep insights into molecular mechanism underlying R. pulchrum flowering process.