Chromosome-level genome assembly and population genetic analysis of a critically endangered rhododendron provide insights into its conservation.

Rhododendrons are woody plants, famous throughout the world as having high horticultural value. However, many wild species are currently threatened with extinction. Here, we report for the first time a high-quality, chromosome-level genome of Rhododendron griersonianum, which has contributed to approximately 10% of all horticultural rhododendron varieties but which in its wild form has been evaluated as critically endangered. The final genome assembly, which has a contig N50 size of approximately 34 M and a total length of 677 M, is the highest-quality genome sequenced within the genus to date, in part due to its low heterozygosity (0.18%). Identified repeats constitute approximately 57% of the genome, and 38 280 protein-coding genes were predicted with high support. We further resequenced 31 individuals of R. griersonianum as well as 30 individuals of its widespread relative R. delavayi, and performed additional conservation genomic analysis. The results showed that R. griersonianum had lower genetic diversity (θ = 2.58e-3; π = 1.94e-3) when compared not only to R. delavayi (θ = 11.61e-3, π = 12.97e-3), but also to most other woody plants. Furthermore, three severe genetic bottlenecks were detected using both the Stairway plot and fastsimcoal2 analysis, which are thought to have occurred in the late Middle Pleistocene and the Last Glacial Maximum (LGM) period. After these bottlenecks, R. griersonianum recovered and maintained a constant effective population size (>25 000) until now. Intriguingly, R. griersonianum has accumulated significantly more deleterious mutations in the homozygous state than R. delavayi, and several deleterious mutations (e.g., in genes involved in the response to heat stress) are likely to have harmed the adaptation of this plant to its surroundings. This high-quality, chromosome-level genome and the population genomic analysis of the critically endangered R. griersonianum will provide an invaluable resource as well as insights for future study in this species to facilitate conservation and in the genus Rhododendron in general.

Floral Scent Chemistry of Luculia yunnanensis (Rubiaceae), a Species Endemic to China with Sweetly Fragrant Flowers.

Luculia plants are famed ornamentals with sweetly fragrant flowers. Luculia yunnanensis Hu is an endemic plant from Yunnan Province, China. Headspace-solid phase microextraction-gas chromatography-mass spectrometry (HS-SPME-GC-MS) was used to identify the volatile organic compounds (VOCs) of the different flower development stages of L. yunnanensis for the evaluation of floral volatile polymorphism. The results showed that a total of 40 compounds were identified at four different stages. The main aroma-active compounds were 3-carene, α-cubebene, α-copaene, δ-cadinene, and isoledene. Floral scent emission had the tendency to ascend first and descend in succession, reaching its peak level at the initial-flowering stage. The richest diversity of floral volatiles was detected at the full-flowering stage. Principal component analysis (PCA) indicated that the composition and its relative content of floral scent differed at the whole flower development stage. In comparison with the other two species of Luculia (L. pinceana and L. gratissima), the composition and its relative content of floral scent were also different among the tree species.

Volatile Organic Compounds Emissions from Luculia pinceana Flower and Its Changes at Different Stages of Flower Development.

Luculia plants are famed ornamental plants with sweetly fragrant flowers, of which L. pinceana Hooker, found primarily in Yunnan Province, China, has the widest distribution. Solid phase microextraction-gas chromatography-mass spectrometry (SPME-GC-MS) was employed to identify the volatile organic compounds (VOCs) emitted from different flower development stages of L. pinceana for the evaluation of floral volatile polymorphism. Peak areas were normalized as percentages and used to determine the relative amounts of the volatiles. The results showed that a total of 39 compounds were identified at four different stages of L. pinceana flower development, including 26 at the bud stage, 26 at the initial-flowering stage, 32 at the full-flowering stage, and 32 at the end-flowering stage. The most abundant compound was paeonol (51%-83%) followed by (E,E)-α-farnesene, cyclosativene, and δ-cadinene. All these volatile compounds create the unique fragrance of L. pinceana flower. Floral scent emission offered tendency of ascending first and descending in succession, meeting its peak level at the initial-flowering stage. The richest diversity of floral volatile was detected at the third and later periods of flower development. Principal component analysis (PCA) indicated that the composition and its relative content of floral scent differed throughout the whole flower development. The result has important implications for future floral fragrance breeding of Luculia. L. pinceana would be adequate for a beneficial houseplant and has a promising prospect for development as essential oil besides for a fragrant ornamental owing to the main compounds of floral scent with many medicinal properties.

A set of novel microsatellite markers developed for Luculia yunnanensis (Rubiaceae), an endangered plant endemic to Yunnan, China.

The genus Luculia Sweet contains about five species of small trees or shrubs and is a member of the family Rubiaceae (tribe Cinchoneae). Luculia yunnanensis is an endangered ornamental shrub endemic to southwest China. Only two natural populations of L. yunnanensis exist in the wild according to our field investigation. It can be inferred that L. yunnanensis is facing a very high risk of extinction in the wild and an urgent conservation strategy is required. By using a modified biotin-sterptavidin capture method, 24 primer sets were identified in two wild populations. Of these primers, 11 displayed polymorphisms and 13 were monomorphic. The number of alleles per locus ranged from two to four, values for observed and expected heterozygosities ranged from 0.000 to 0.833 and from 0.431 to 0.771, with averages of 0.389 and 0.614, respectively. These markers will be useful for further investigation of conservation of resources, selecting parental types in cross-breeding, evolution of this species at the molecular level and related research in Luculia species.

Morphological, Physiological, and Molecular Responses of Sweetly Fragrant Luculia gratissima During the Floral Transition Stage Induced by Short-Day Photoperiod.

Photoperiod-regulated floral transition is vital to the flowering plant. Luculia gratissima "Xiangfei" is a flowering ornamental plant with high development potential economically and is a short-day woody perennial. However, the genetic regulation of short-day-induced floral transition in L. gratissima is unclear. To systematically research the responses of L. gratissima during this process, dynamic changes in morphology, physiology, and transcript levels were observed and identified in different developmental stages of long-day- and short-day-treated L. gratissima plants. We found that floral transition in L. gratissima occurred 10 d after short-day induction, but flower bud differentiation did not occur at any stage under long-day conditions. A total of 1,226 differentially expressed genes were identified, of which 146 genes were associated with flowering pathways of sugar, phytohormones, photoperiod, ambient temperature, and aging signals, as well as floral integrator and meristem identity genes. The trehalose-6-phosphate signal positively modulated floral transition by interacting with SQUAMOSA PROMOTER-BINDING-LIKE PROTEIN 4 (SPL4) in the aging pathway. Endogenous gibberellin, abscisic acid, cytokinin, and jasmonic acid promoted floral transition, whereas strigolactone inhibited it. In the photoperiod pathway, FD, CONSTANS-LIKE 12, and nuclear factors Y positively controlled floral transition, whereas PSEUDO-RESPONSE REGULATOR 7, FLAVIN-BINDING KELCH REPEAT F-BOX PROTEIN 1, and LUX negatively regulated it. SPL4 and pEARLI1 positively affected floral transition. Suppressor of Overexpression of Constans 1 and AGAMOUSLIKE24 integrated multiple flowering signals to modulate the expression of FRUITFULL/AGL8, AP1, LEAFY, SEPALLATAs, SHORT VEGETATIVE PHASE, and TERMINAL FLOWER 1, thereby regulating floral transition. Finally, we propose a regulatory network model for short-day-induced floral transition in L. gratissima. This study improves our understanding of flowering time regulation in L. gratissima and provides knowledge for its production and commercialization.

Genetic Diversity of Phyllanthus emblica From Two Different Climate Type Areas.

Phyllanthus emblica L. is a well-known medicinal and edible plant species. Various medicinal compounds in the fruit make it an important medicinal and promising economic material. The plant is widely distributed in Southwestern and Southern China. However, due to massive deforestation and land reclamation as well as deterioration of its natural habitat in recent years, the wild resources of this species have been sharply reduced, and it is rare to see large-scale wild P. emblica forests so far. In order to effectively protect and rationally utilize this species, we investigated the genetic diversity, genetic structure, and population dynamics of 260 individuals from 10 populations of P. emblica sampled from the dry climate area in Yunnan and wet climate area in Guangxi using 20 polymorphic EST-SSR markers. We found high genetic diversity at the species level (He = 0.796) and within populations (He = 0.792), but low genetic differentiation among populations (F ST = 0.084). In addition, most genetic variation existed within populations (92.44%) compared with variation among the populations (7.56%). Meanwhile, the NJ tree, STRUCTURE, and hierarchical analysis suggested that the sampled individuals were clustered into two distinct genetic groups. In contrast, the genetic diversity of the dry climate group (He = 0.786, Na = 11.790, I = 1.962) was higher than that of the wet climate group (He = 0.673, Na = 9.060, I = 1.555), which might be attributed to the combined effects of altitude, precipitation, and geographic distance. Interestingly, only altitude and precipitation had significant pure effects on the genetic diversity, and the former was slightly stronger. In addition, DIYABC analysis suggested the effective population size of P. emblica might have contracted in the beginning of the Last Glacial Maximum. These genetic features provided vital information for the conservation and sustainable development of genetic resources of P. emblica, and they also provided new insights and guidelines for ecological restoration and economic development in dry-hot valleys of Yunnan and karst areas in Guangxi.

Development of novel EST-SSR markers for Rhododendron longipedicellatum (Ericaceae) and cross-amplification in two congeners.

PREMISE OF THE STUDY: To investigate the genetic background and population characteristics of Rhododendron longipedicellatum (Ericaceae), a newly discovered and critically endangered species, expressed sequence tag-simple sequence repeat markers were developed, and transferability was tested in two congeners, R. molle and R. simsii. METHODS AND RESULTS: Based on the transcriptome sequences of R. longipedicellatum, 102 primer sets were designed; 48 primer sets were successfully amplified, with 15 showing polymorphisms in 150 individuals from five extant populations of R. longipedicellatum. The number of alleles per locus ranged from four to 18, and the levels of observed and expected heterozygosity for the 15 loci varied from 0.255 to 0.913 and from 0.306 to 0.851, respectively. All 15 loci were found to amplify in R. molle and R. simsii. CONCLUSIONS: These polymorphic SSR markers can be used in conservation genetic and phylogeographic studies to elucidate the rarity and origin of R. longipedicellatum.

Study on Reproductive Biology of Rhododendron longipedicellatum: A Newly Discovered and Special Threatened Plant Surviving in Limestone Habitat in Southeast Yunnan, China.

Rhododendron longipedicellatum is a narrow endemic species and a subject of urgent demand in the domestic market and overseas. Its fascinating shapes, brilliantly gilvous flowers, and unusual flowering time endow this species with extremely high ornamental value. However, only five wild populations of R. longipedicellatum surviving in limestone habitat have been found through elaborate field investigation, and the number of the populations decreases further or is even confronted with risk of extinction due to the damage of human activities. To enhance the protection and utilization of R. longipedicellatum, this study systematically investigated several important aspects of reproductive biology, including floral syndrome, pollen viability and stigma receptivity, petal color reflectance, breeding system, and pollination biology. The results demonstrated that arched styles not only create obvious herkogamy that avoide self-pollination, but also effectively reduce rain damage to the intrinsic characteristics of the stigma surface secretions, promoting the female fitness of R. longipedicellatum in poor weather. Pollen viability maintained a high level over the flowering period. The reflectance spectrum of petals had two peaks at wavelengths of 360 and 580 nm. Tests of OCI, P/O and artificial pollination all indicated that R. longipedicellatum was self-compatible and that the breeding system was mixed mating. Geitonogamy mediated by Bombus braccatus was the primary pollination route in the natural environment, which suggested that the breeding system of R. longipedicellatum might be evolving from selfing to outcrossing. The pollination vector of R. longipedicellatum was very specific, in that only B. braccatus was confirmed to deliver pollen to the stigmas. Visitation frequency was influenced by the activity rhythms and resource requirements of the different castes (i.e., sex). B. braccatus workers were the most effective pollinators because of higher visitation frequency and more effective contribution to fruit production, whereas the presence of B. braccatus males might enhance pollen flow within the population to a certain extent. Finally, these findings not only provided a reliable theoretical basis for hybridization breeding of R. longipedicellatum as parents, but also laid a solid foundation for further molecular biology studies to more broadly reveal the mechanisms of its endangerment in the future.

Development of novel EST-SSR markers for Phyllanthus emblica (Phyllanthaceae) and cross-amplification in two related species.

PREMISE OF THE STUDY: A novel set of EST-SSR markers was developed for Phyllanthus emblica (Phyllanthaceae) to investigate the genetic structure and gene flow, identify novel genes of interest, and develop markers for assisted breeding. METHODS AND RESULTS: Based on the transcriptome data of P. emblica, 83 EST-SSR primer pairs were designed; 52 primer pairs were successfully amplified, with 20 showing polymorphisms in 90 individuals from three populations of P. emblica. The number of alleles per locus varied from 11 to 44. The observed and expected levels of heterozygosity for the 20 loci ranged from 0.240 to 0.868 and 0.754 to 0.933, respectively. Cross-species amplification was successful for all 20 loci in each of the two related species, P. reticulatus and Leptopus chinensis. CONCLUSIONS: These markers will be valuable for studying the population genetics and for mining genes of P. emblica, and may be useful for studies of related species.