Rhododendron kuomeianum Y.H. Chang, J. Nielsen & Y.P. Ma, a new species of Rhododendron (Ericaceae) within subsect. Maddenia in sect. Rhododendron from Yiliang County, NE Yunnan, China, is described and illustrated. The new species is similar to R. valentinianum, but it can be easily distinguished by its sparse scales on the abaxial surface of the leaf blade, fewer flowers per inflorescence and white corolla with pale red margins. There are also differences in the widths of calyx lobes, leaf blade shape and indumentum characteristics of the petiole between the new species and Rhododendron linearilobum. We confirmed that R. kuomeianum is a new species closely related to R. valentinianum and R. changii with phylogenomic studies of 10 species within this subsection based on restriction site-associated DNA sequencing (RAD-seq) data. These phylogenomic analyses also clarified additional taxonomic problems in this subsection previously raised by morphological analysis. Our findings make a strong case for using next-generation sequencing to explore phylogenetic relationships and identify new species, especially in plants groups with complicated taxonomic problems.

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 F₁ 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.

Knowledge about how species richness varies along spatial and environmental gradients is important for the conservation and use of biodiversity. The Ericaceae is a major component of alpine and subalpine vegetation globally. However, little is known about the spatial pattern of species richness and the factors that drive that richness in Ericaceae. We investigated variation in species richness of Ericaceae along an elevational gradient in Yunnan, China, and used a variation partitioning analysis based on redundancy analysis ordination to examine how those changes might be influenced by the mid-domain effect, the species-area relationship, and climatic variables. Species richness varied significantly with elevation, peaking in the upper third of the elevational gradient. Of the factors examined, climate explained a larger proportion of the variance in species richness along the elevational gradient than either land area or geometric constraints. Species richness showed a unimodal relationship with mean annual temperature and mean annual precipitation. The elevational pattern of species richness for Ericaceae was shaped by the combined effects of climate and competition. Our findings contribute to a better understanding of the potential effects of climate change on species richness for Ericaceae.

Rhododendron delavayi Franch. is a famous ornamental plant. However, seasonal drought caused by a monsoon climate seriously limits its growth and flowering performance in urban gardens. Our previous study has found that brassinosteroids (BRs) can improve the tolerance of R. delavayi to drought stress. Here, we employed a data-independent acquisition (DIA) approach to compare the protein expression profiles under drought treatment (D) and pre-treatment with BR before drought treatment (BR). With an increase in drought stress, the net photosynthetic rate, stomatal conductance, and transpiration rate in the BR treatment showed more stable changes that were significantly higher than those in the D treatment. However, the contents of malondialdehyde, soluble sugar, soluble protein, and the activity of superoxide dismutase (SOD), peroxidase, and catalase showed opposite trends. The pre-treatment with BR alleviated the negative effect of drought stress on the photosynthetic performance of R. delavayi. A total of 3453 differentially expressed proteins (DEPs) were identified, and 683 DEPs were significantly expressed in the D and BR treatments. The DEPs uniquely expressed in the BR treatment participated in the pathways of “ribosome”, “ether lipid metabolism”, “photosynthesis”, and “oxidative phosphorylation”. The improvement effect of the BR treatment on the drought tolerance of R. delavayi was mainly attributed to improved photosynthesis by alleviating stomatal closure and oxidative stress, maintaining the integrity and stability of the ribosomal complex to mediate protein synthesis and the balance between energy metabolism and carbon metabolism. Our study presents a comprehensive understanding of the key proteins and metabolic pathways related to the response of R. delavayi to drought and will contribute to the breeding of drought-tolerant rhododendrons.

The cytokinin response regulator (RR) gene is essential for cytokinin signal transduction, which plays a crucial role in plant growth and development. Here, we applied bioinformatics to Rhododendron delavayi’s genome to identify its RR gene family and systematically analyzed their gene characteristics, phylogenetic evolution, chromosomal localization, collinearity analysis, promoter cis-elements, and expression patterns. Overall, 33 RdRR genes were distinguished and classified into three types. All these genes harbored motif 5 (YEVTTVNSGLEALELLRENKB), the most conserved one, along with the plant-conserved domain (REC domain), and could be mapped to 10 chromosomes with four gene pairs of segmental replication events but no tandem replication events; 13 RdRR genes showed collinearity with Arabidopsis thaliana genes. Promoter analysis revealed multiple hormone-related cis-elements in the RR genes. After a TDZ (thidiazuron) treatment, 13 genes had higher expression levels than the control, whose magnitude of change depended on the developmental stage of leaves’ adventitious buds. The expression levels of RdRR14, RdRR17, RdRR20, and RdRR24 agreed with the average number of adventitious buds post-TDZ treatment. We speculate that these four genes could figure prominently in bud regeneration from R. delavayi leaves in vitro. This study provides detailed knowledge of RdRRs for research on cytokinin signaling and RdRR functioning in R. delavayi.