Rhododendron meddianum is a critically endangered species with important ornamental value and is also a plant species with extremely small populations. In this study, we used double digest restriction-site-associated DNA sequencing (ddRAD) technology to assess the genetic diversity, genetic structure and demographic history of the three extant populations of R. meddianum. Analysis of SNPs indicated that R. meddianum populations have a high genetic diversity (π = 0.0772 ± 0.0024, H_E = 0.0742 ± 0.002). Both F_ST values (0.1582–0.2388) and AMOVA showed a moderate genetic differentiation among the R. meddianum populations. Meanwhile, STRUCTURE, PCoA and NJ trees indicated that the R. meddianum samples were clustered into three distinct genetic groups. Using the stairway plot, we found that R. meddianum underwent a population bottleneck about 70,000 years ago. Furthermore, demographic models of R. meddianum and its relative, Rhododendron cyanocarpum, revealed that these species diverged about 3.05 (2.21–5.03) million years ago. This divergence may have been caused by environmental changes that occurred after the late Pliocene, e.g., the Asian winter monsoon intensified, leading to a drier climate. Based on these findings, we recommend that R. meddianum be conserved through in situ, ex situ approaches and that its seeds be collected for germplasm.

Luculia yunnanensis (Rubiaceae), an evergreen shrub or small tree, is endemic to China and confined to Nujiang Prefecture, Yunnan Province. This plant is of high ornamental value owing to its attractive pink flowers, sweet fragrance, and long flowering period. Due to the influence of climate change and human factors, the distribution range of L. yunnanensis has exhibited a significant shrinking trend, and it has become a vulnerable species that is in urgent need of conservation and rational utilization research. In this study, the flower transcriptome sequencing of L. yunnanensis was conducted using an Illumina HiSeq platform. We designed and developed a series of EST-SSR primers based on the flower transcriptome data of L. yunnanensis. The results showed that 98,389 unigenes were obtained from the L. yunnanensis flower transcriptome, all of which were aligned with sequences in public databases. Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KEGG, and GO annotated 31,859, 13,853, 22,684, 10,947, 21,416, 9722, and 23,390 unigenes, respectively. The MISA (Microsatellite) tool was used to identify SSR loci from all unigenes, and a total of 15,384 SSRs were identified. Repeat motifs were given priority with mononucleotides, dinucleotides, and trinucleotides. The 81 primer pairs were synthesized randomly, of which 44 pairs showed effective amplification. A total of 17 primers showed stable amplification, and rich polymorphism was observed in 6 populations. We concluded via genetic diversity analysis that the average effective number of alleles (Ne), Shannon’s information index (I), and polymorphism information content (PIC) were 1.925, 0.837, and 0.403, respectively. In conclusion, 17 EST-SSR primers can be used for subsequent population genetic diversity analysis and molecular-marker-assisted breeding, which is of great significance for formulating resource conservation and utilization strategies for L. yunnanensis.