De novo RNA sequencing transcriptome of Rhododendron obtusum identified the early heat response genes involved in the transcriptional regulation of photosynthesis

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.

Transcriptome and photosynthetic analyses provide new insight into the molecular mechanisms underlying heat stress tolerance in Rhododendron × pulchrum Sweet

Rhododendron species provide excellent ornamental use worldwide, yet heat stress (HS) is one of the major threats to their cultivation. However, the intricate mechanisms underlying the photochemical and transcriptional regulations associated with the heat stress response in Rhododendron remain relatively unexplored. In this study, the analyses of morphological characteristics and chlorophyll fluorescence (ChlF) kinetics showed that HS (40 °C/35 °C) had a notable impact on both the donor's and acceptor's sides of photosystem II (PSII), resulting in reduced PSII activity and electron transfer capacity. The gradual recovery of plants observed following a 5-day period of culture under normal conditions indicates the reversible nature of the HS impact on Rhododendron × pulchrum. Analysis of transcriptome data unveiled noteworthy trends: four genes associated with photosynthesis-antenna protein synthesis (LHCb1, LHCb2 and LHCb3) and the antioxidant system (glutamate-cysteine ligase) experienced significant down-regulation in the leaves of R. × pulchrum during HS. Conversely, aseorbate peroxidase and glutathione S-transferase TAU 8 demonstrated an up-regulated pattern. Furthermore, six down-regulated genes (phos-phoenolpyruvate carboxylase 4, sedoheptulose-bisphosphatase, ribose-5-phosphate isomerase 2, high cyclic electron flow 1, beta glucosidase 32 and starch synthase 2) and two up-regulated genes (beta glucosidase 2 and UDP-glucose pyrophosphorylase 2) implicated in photosynthetic carbon fixation and starch/sucrose metabolism were identified during the recovery process. To augment these insights, a weighted gene co-expression network analysis yielded a co-expression network, pinpointing the hub genes correlated with ChlF dynamics' variation trends. The cumulative results showed that HS inhibited the synthesis of photosynthesis-antenna proteins in R. × pulchrum leaves. This disruption subsequently led to diminished photochemical activities in both PSII and PSI, albeit with PSI exhibiting heightened thermostability. Depending on the regulation of the reactive oxygen species scavenging system and heat dissipation, photoprotection sustained the recoverability of R. × pulchrum to HS.

Time-series RNA-Seq transcriptome profiling reveals novel insights about cold acclimation and de-acclimation processes in an evergreen shrub of high altitude

The high-altitude alpine regions are characterized by highly variable and harsh environmental conditions. However, relatively little is known about the diverse mechanisms adopted by alpine plants to adapt to these stressful conditions. Here, we studied variation in transcriptome and physiological adjustments occurring across the year at high elevation environments in the leaf tissue of Rhododendron anthopogon, an evergreen shrub of Himalaya. The samples were collected at 12 different time-points, from August until snowfall in November 2017, and then from June to September 2018. It was observed that with a drop in both ambient air temperature and photoperiod towards onset of winter, the freezing resistance of plants increased, resulting in 'cold acclimation'. Further, 'de-acclimation' was associated with a decrease in freezing resistance and increase in photosynthetic efficiency of leaves during spring. A considerable amount of variation was observed in the transcriptome in a time-dependent sequential manner, with a total of 9,881 differentially expressed genes. Based on gene expression profiles, the time-points could be segregated into four clusters directly correlating with the distinct phases of acclimation: non-acclimation (22-August-2017, 14-August-2018, 31-August-2018), early cold acclimation (12-September-2017, 29-September-2017), late cold acclimation (11-October-2017, 23-October-2017, 04-November-2017, 18-September-2018) and de-acclimation (15-June-2018, 28-June-2018, 14-July-2018). Cold acclimation was a gradual process, as indicated by presence of an intermediate stage (early acclimation). However, the plants can by-pass this stage when sudden decrease in temperature is encountered. The maximum variation in expression levels of genes occurred during the transition to de-acclimation, hence was 'transcriptionally' the most active phase. The similar or higher expression levels of genes during de-acclimation in comparison to non-acclimation suggested that molecular functionality is re-initiated after passing through the harsh winter conditions.

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.

SMRT sequencing of a full-length transcriptome reveals transcript variants involved in C18 unsaturated fatty acid biosynthesis and metabolism pathways at chilling temperature in Pennisetum giganteum

Background

Pennisetum giganteum, an abundant, fast-growing perennial C₄ grass that belongs to the genus Pennisetum, family Poaceae, has been developed as a source of biomass for mushroom cultivation and production, as a source of forage for cattle and sheep, and as a tool to remedy soil erosion. However, having a chilling-sensitive nature, P. giganteum seedlings need to be protected while overwintering in most temperate climate regions.

Results

To elucidate the cold stress responses of P. giganteum, we carried out comprehensive full-length transcriptomes from leaf and root tissues under room temperature (RT) and chilling temperature (CT) using PacBio Iso-Seq long reads. We identified 196,124 and 140,766 full-length consensus transcripts in the RT and CT samples, respectively. We then systematically performed functional annotation, transcription factor identification, long non-coding RNAs (lncRNAs) prediction, and simple sequence repeat (SSR) analysis of those full-length transcriptomes. Isoform analysis revealed that alternative splicing events may be induced by cold stress in P. giganteum, and transcript variants may be involved in C18 unsaturated fatty acid biosynthesis and metabolism pathways at chilling temperature in P. giganteum. Furthermore, the fatty acid composition determination and gene expression level analysis supported that C18 unsaturated fatty acid biosynthesis and metabolism pathways may play roles during cold stress in P. giganteum.

Conclusions

We provide the first comprehensive full-length transcriptomic resource for the abundant and fast-growing perennial grass Pennisetum giganteum. Our results provide a useful transcriptomic resource for exploring the biological pathways involved in the cold stress responses of P. giganteum.

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.

The Impact of Chronic Heat Stress on the Growth, Survival, Feeding, and Differential Gene Expression in the Sea Urchin Strongylocentrotus intermedius

To explore the impact of chronic heat stress on commercial echinoderms, the present study assessed the effects of chronic high temperature on the growth, survival, feeding, and differential gene expression in the sea urchin Strongylocentrotus intermedius cultured in northern Yellow Sea in China. One suitable seawater condition (20°C) and one laboratory-controlled high temperature condition (25°C) were set up. After 28 days incubation, our results showed that: (1) The specific growth, survival, and ingestion rates of S. intermedius reared under high temperature (25°C) decreased compared to those reared under optimal temperature (20°C) conditions; (2) comparative transcriptome analysis identified 2,125 differentially expressed genes (DEGs) in S. intermedius reared under high temperature (25°C) compared to those subjected to optimal temperature condition (20°C), which included 1,015 upregulated and 1,100 downregulated genes. The accuracy of the transcriptome profiles was verified by quantitative real-time PCR (qRT-PCR). Further Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses revealed that these DEGs mainly enriched the functional categories of ribosome, protein processing in endoplasmic reticulum, and prion diseases. A total of 732 temperature-induced expressed genes, such as ATP5, heat shock protein 70, and heat shock protein 90, were identified as candidates that were closely correlated with heat resistance in S. intermedius. Differentially expressed transcription factors (TFs), such as AP-1, Fos, CREB, and ZNF, were also identified as potential regulators that regulate the molecular network that was associated with responses to heat stress in sea urchins. Observations in the present study provide additional information that improves our understanding of the molecular mechanism of temperate echinoid species in response to heat stress, as well as theoretical basis for the molecular-assisted breeding of heat-resistant sea urchins.

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.