Chromosome-scale genome assembly of Rhododendron molle provides insights into its evolution and terpenoid biosynthesis

Characterization and Engineering of Two Highly Paralogous Sesquiterpene Synthases Reveal a Regioselective Reprotonation Switch

Sesquiterpene synthases (STPSs) catalyze carbocation-driven cyclization reactions that can generate structurally diverse hydrocarbons. The deprotonation-reprotonation process is widely used in STPSs to promote structural diversity, largely attributable to the distinct regio/stereoselective reprotonations. However, the molecular basis for reprotonation regioselectivity remains largely understudied. Herein, we analyzed two highly paralogous STPSs, Artabotrys hexapetalus (-)-cyperene synthase (AhCS) and ishwarane synthase (AhIS), which catalyze reactions that are distinct from the regioselective protonation of germacrene A (GA), resulting in distinct skeletons of 5/5/6 tricyclic (-)-cyperene and 6/6/5/3 tetracyclic ishwarane, respectively. Isotopic labeling experiments demonstrated that these protonations occur at C3 and C6 of GA in AhCS and AhIS, respectively. The cryo-electron microscopy-derived AhCS complex structure provided the structural basis for identifying different key active site residues that may govern their functional disparity. The structure-guided mutagenesis of these residues resulted in successful functional interconversion between AhCS and AhIS, thus targeting the three active site residues [L311-S419-C458]/[M311-V419-A458] that may act as a C3/C6 reprotonation switch for GA. These findings facilitate the rational design or directed evolution of STPSs with structurally diverse skeletons.

A haplotype-resolved genome assembly of Rhododendron vialii based on PacBio HiFi reads and Hi-C data

Rhododendron vialii (subgen. Azaleastrum) is an evergreen shrub with high ornamental value. This species has been listed as a plant species with extremely small populations (PSESP) for urgent protection by China's Yunnan provincial government in 2021, due to anthropogenic habitat fragmentation. However, limited genomic resources hinder scientifically understanding of genetic threats that the species is currently facing. In this study, we assembled a high-quality haplotype-resolved genome of R. vialii based on PacBio HiFi long reads and Hi-C reads. The assembly contains two haploid genomes with sizes 532.73 Mb and 521.98 Mb, with contig N50 length of 35.67 Mb and 34.70 Mb, respectively. About 99.92% of the assembled sequences could be anchored to 26 pseudochromosomes, and 14 gapless assembled chromosomes were included in this assembly. Additionally, 60,926 protein-coding genes were identified, of which 93.82% were functionally annotated. This is the first reported genome of R. vialii, and hopefully it will lay the foundations for further research into the conservation genomics and horticultural domestication of this ornamentally important species.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.