The first high-altitude autotetraploid haplotype-resolved genome assembled (Rhododendron nivale subsp. boreale) provides new insights into mountaintop adaptation.

BACKGROUND: Rhododendron nivale subsp. boreale Philipson et M. N. Philipson is an alpine woody species with ornamental qualities that serve as the predominant species in mountainous scrub habitats found at an altitude of ∼4,200 m. As a high-altitude woody polyploid, this species may serve as a model to understand how plants adapt to alpine environments. Despite its ecological significance, the lack of genomic resources has hindered a comprehensive understanding of its evolutionary and adaptive characteristics in high-altitude mountainous environments. FINDINGS: We sequenced and assembled the genome of R. nivale subsp. boreale, an assembly of the first subgenus Rhododendron and the first high-altitude woody flowering tetraploid, contributing an important genomic resource for alpine woody flora. The assembly included 52 pseudochromosomes (scaffold N50 = 42.93 Mb; BUSCO = 98.8%; QV = 45.51; S-AQI = 98.69), which belonged to 4 haplotypes, harboring 127,810 predicted protein-coding genes. Conjoint k-mer analysis, collinearity assessment, and phylogenetic investigation corroborated autotetraploid identity. Comparative genomic analysis revealed that R. nivale subsp. boreale originated as a neopolyploid of R. nivale and underwent 2 rounds of ancient polyploidy events. Transcriptional expression analysis showed that differences in expression between alleles were common and randomly distributed in the genome. We identified extended gene families and signatures of positive selection that are involved not only in adaptation to the mountaintop ecosystem (response to stress and developmental regulation) but also in autotetraploid reproduction (meiotic stabilization). Additionally, the expression levels of the (group VII ethylene response factor transcription factors) ERF VIIs were significantly higher than the mean global gene expression. We suspect that these changes have enabled the success of this species at high altitudes. CONCLUSIONS: We assembled the first high-altitude autopolyploid genome and achieved chromosome-level assembly within the subgenus Rhododendron. In addition, a high-altitude adaptation strategy of R. nivale subsp. boreale was reasonably speculated. This study provides valuable data for the exploration of alpine mountaintop adaptations and the correlation between extreme environments and species polyploidization.

Integrative study of transcriptome and microbiome to reveal the response of Rhododendron decorum to cadmium stress.

The anomalies of cadmium (Cd) in karst region pose a severe threat to plant growth and development. In this study, the responses of Rhododendron decorum to Cd stress were investigated at physiological, molecular, and endophytic microbial levels, and the potential correlation among these responses was assessed. The Cd stress impeded R. decorum growth and led to an increase in malondialdehyde (MDA) and hydrogen peroxide (H2O2) levels, as well as enhanced superoxide dismutase (SOD) and catalase (CAT) activities. Meanwhile, Cd stress increased the Cd (up to 80 times compared to the control), sodium (Na), aluminum (Al), and zinc (Zn) contents, while decreased the magnesium (Mg) and manganese (Mn) contents in R. decorum leaves. Transcriptome suggested that Cd significantly regulated the pathways including "protein repair", "hormone-mediated signaling pathway", and "ATP-binding cassette (ABC) transporters". Additionally, q-PCR analysis showed that Cd stress significantly up-regulated the expressions of ABCB19-like and pleiotropic drug resistance, while down-regulated the expressions of indole-3-acetic acid-amido synthetase and cytokinin dehydrogenase. The Cd stress influenced the composition of endophytic microbial communities in R. decorum leaves and enhanced the interspecific bacterial associations. Furthermore, the bacterial genera Achromobacter, Aureimonas and fungal genus Vishniacozyma exhibited a high degree of connectivity with other nodes in networks constructed by the metal element contents, differentially expressed genes (DEGs), and microbial communities, respectively. These findings provide a comprehensive insight into the response of R. decorum to Cd-induced stress, which might facilitate the breeding of the Cd-tolerant R. decorum.

Chalcone-Synthase-Encoding RdCHS1 Is Involved in Flavonoid Biosynthesis in Rhododendron delavayi.

Flower color is an important ornamental feature that is often modulated by the contents of flavonoids. Chalcone synthase is the first key enzyme in the biosynthesis of flavonoids, but little is known about the role of R. delavayi CHS in flavonoid biosynthesis. In this paper, three CHS genes (RdCHS1-3) were successfully cloned from R. delavayi flowers. According to multiple sequence alignment and a phylogenetic analysis, only RdCHS1 contained all the highly conserved and important residues, which was classified into the cluster of bona fide CHSs. RdCHS1 was then subjected to further functional analysis. Real-time PCR analysis revealed that the transcripts of RdCHS1 were the highest in the leaves and lowest in the roots; this did not match the anthocyanin accumulation patterns during flower development. Biochemical characterization displayed that RdCHS1 could catalyze p-coumaroyl-CoA and malonyl-CoA molecules to produce naringenin chalcone. The physiological function of RdCHS1 was checked in Arabidopsis mutants and tobacco, and the results showed that RdCHS1 transgenes could recover the color phenotypes of the tt4 mutant and caused the tobacco flower color to change from pink to dark pink through modulating the expressions of endogenous structural and regulatory genes in the tobacco. All these results demonstrate that RdCHS1 fulfills the function of a bona fide CHS and contributes to flavonoid biosynthesis in R. delavayi.

Spirosoma rhododendri sp. nov., isolated from a flower of royal azalea (Rhododendron schlippenbachii).

A Gram-negative, non-motile, strictly aerobic and rod- or filamentous-shaped strain, CJU-R4T, was isolated from a flower of royal azalea (Rhododendron schlippenbachii) collected in the Republic of Korea. Strain CJU-R4T was catalase-positive and oxidase-negative, and grew at 15-33 °C (optimum, 28-20 °C), at pH 5.0-8.0 (optimum, pH 7.0-8.0), and in the presence of 0-1 % NaCl (w/v; optimum, 0 %). Strain CJU-R4T had the highest 16S rRNA gene sequence similarity to Spirosoma oryzae RHs22T (96.6 %), revealing less than 93 % sequence similarity to other type strains. Phylogenetic and phylogenomic analysis also revealed strain CJU-R4T formed a robust cluster with S. oryzae RHs22T. The major fatty acids were summed feature 3 (comprising C16 : 1 ω7c and/or C16 : 1 ω6c; 33.0 %), C16 : 1 ω5c (22.1 %), iso-C15 : 0 (12.6 %) and C16 : 0 (10.7 %). The polar lipids were composed of phosphatidylethanolamine, three unidentified aminophospholipids, one unidentified phospholipid and four unidentified lipids. Menaquinone-7 was detected as the sole respiratory quinone. The genomic DNA G+C content was 55.2 mol%. The average nucleotide identity and digital DNA-DNA hybridization values between strain CJU-R4T and Spirosoma oryzae DSM 28354T were 81.5 and 23.9 %, respectively. Based on the results of the phenotypic and genotypic analyses, strain CJU-R4T is considered to represent a novel species of the genus Spirosoma, for which the name Spirosoma rhododendri sp. nov. is proposed. The type strain is CJU-R4T (=KACC 21264T=NBRC 114513T).

Characterization and functional analysis of RdDFR1 regulation on flower color formation in Rhododendron delavayi.

Rhododendron delavayi is a popular ornamental plant with globular flowers noted for their bright red color, but very limited studies have been reported on its flower color formation. In this study, we successfully isolated a novel DFR gene (RdDFR1) from red flowers of Rhododendron delavayi. Multiple sequence alignments revealed that RdDFR1 had the conserved NADP and substrate binding domain, and was classified into Asn-type DFR. Meanwhile, quantitative real-time PCR analysis showed that transcript levels of RdDFR1 matched the accumulation patterns of anthocyanins during flower development, hinting its potential role involved in anthocyanin biosynthesis. Then in vitro enzymatic analysis indicated that recombinant RdDFR1 protein could catalyze the production of leucoanthocyanidins from dihydroquercetin and dihydromyricetin. Furthermore, the in planta assay, using Arabidopsis thaliana dfr mutant (tt3-1) and tobacco, displayed that RdDFR1 transgenes recovered the defective proanthocyanidin and anthocyanin biosynthesis at seed coats, hypocotyl as well as cotyledon, and altered the flowers color of tobacco from pale pink to dark pink which demonstrated its function as dihydroflavonol 4-reductase in vivo. In summary, our findings suggest that RdDFR1 plays a crucial role in the biosynthesis of anthocyanin and will also make a contribution to understand the mechanisms of flower color formation in Rhododendron delavayi.

Incongruent phylogenies and their implications for the study of diversification, taxonomy, and genome size evolution of Rhododendron.

PREMISE: Classification of taxa depends on the quality of inferred phylogenies. Rhododendron, a highly species-rich genus (>1156 species) of woody plants, has a highly debated infrageneric classification, due to its huge diversity, homoplasy in key characters, and incongruence among data sets. We provide a broad coverage of representative species to resolve Rhododendron infrageneric phylogeny and highlight the areas of incongruence. We further investigate the effect of polyploidy and genome size evolution on diversification of Rhododendron. METHODS: We generated two plastid and two nuclear loci for 260 Rhododendron species. We analyzed the loci separately as well as concatenated, utilizing both likelihood and Bayesian methods. We tested incongruence both among the data sets and with previous studies. We estimated genome sizes for 125 species through flow cytometry. RESULTS: Our results suggest stronger support for larger subgenera; however, the smaller subgenera pose several problems; for example, R. tomentosum (former genus Ledum) occupies incongruent positions based on different DNA regions. The main shift to higher diversification in the genus occurs in the Himalayan/Southeast Asian clade of R. subg. Hymenanthes. We found that polyploidy occurs in almost all subgenera but most frequently within R. subg. Rhododendron sections Rhododendron and Schistanthe. CONCLUSIONS: We endorse the recognition of five major clades at the subgeneric level, but a number of species cannot be confidently assigned to these clades due to incongruency. With regard to genome size evolution, results support previous reports that genome sizes of tropical plants are lower than those of colder and temperate regions and that genome downsizing promotes diversification.

An Aromatic Farnesyltransferase Functions in Biosynthesis of the Anti-HIV Meroterpenoid Daurichromenic Acid.

Rhododendron dauricum produces daurichromenic acid, an anti-HIV meroterpenoid, via oxidative cyclization of the farnesyl group of grifolic acid. The prenyltransferase (PT) that synthesizes grifolic acid is a farnesyltransferase in plant specialized metabolism. In this study, we demonstrated that the isoprenoid moiety of grifolic acid is derived from the 2-C-methyl-d-erythritol-4-phosphate pathway that takes place in plastids. We explored candidate sequences of plastid-localized PT homologs and identified a cDNA for this PT, RdPT1, which shares moderate sequence similarity with known aromatic PTs. RdPT1 is expressed exclusively in the glandular scales, where daurichromenic acid accumulates. In addition, the gene product was targeted to plastids in plant cells. The recombinant RdPT1 regiospecifically synthesized grifolic acid from orsellinic acid and farnesyl diphosphate, demonstrating that RdPT1 is the farnesyltransferase involved in daurichromenic acid biosynthesis. This enzyme strictly preferred orsellinic acid as a prenyl acceptor, whereas it had a relaxed specificity for prenyl donor structures, also accepting geranyl and geranylgeranyl diphosphates with modest efficiency to synthesize prenyl chain analogs of grifolic acid. Such a broad specificity is a unique catalytic feature of RdPT1 that is not shared among secondary metabolic aromatic PTs in plants. We discuss the unusual substrate preference of RdPT1 using a molecular modeling approach. The biochemical properties as well as the localization of RdPT1 suggest that this enzyme produces meroterpenoids in glandular scales cooperatively with previously identified daurichromenic acid synthase, probably for chemical defense on the surface of R. dauricum plants.

Identification and Characterization of Daurichromenic Acid Synthase Active in Anti-HIV Biosynthesis.

Daurichromenic acid (DCA) synthase catalyzes the oxidative cyclization of grifolic acid to produce DCA, an anti-HIV meroterpenoid isolated from Rhododendron dauricum We identified a novel cDNA encoding DCA synthase by transcriptome-based screening from young leaves of R. dauricum The gene coded for a 533-amino acid polypeptide with moderate homologies to flavin adenine dinucleotide oxidases from other plants. The primary structure contained an amino-terminal signal peptide and conserved amino acid residues to form bicovalent linkage to the flavin adenine dinucleotide isoalloxazine ring at histidine-112 and cysteine-175. In addition, the recombinant DCA synthase, purified from the culture supernatant of transgenic Pichia pastoris, exhibited structural and functional properties as a flavoprotein. The reaction mechanism of DCA synthase characterized herein partly shares a similarity with those of cannabinoid synthases from Cannabis sativa, whereas DCA synthase catalyzes a novel cyclization reaction of the farnesyl moiety of a meroterpenoid natural product of plant origin. Moreover, in this study, we present evidence that DCA is biosynthesized and accumulated specifically in the glandular scales, on the surface of R. dauricum plants, based on various analytical studies at the chemical, biochemical, and molecular levels. The extracellular localization of DCA also was confirmed by a confocal microscopic analysis of its autofluorescence. These data highlight the unique feature of DCA: the final step of biosynthesis is completed in apoplastic space, and it is highly accumulated outside the scale cells.

In vitro and in vivo evidence for oxalate oxidase activity of a germin-like protein from azalea.

Germins and germin-like proteins (GLPs) comprise large families of extracellular plant glycoproteins that are structurally similar, yet they have been reported to have distinct biochemical activities: oxalate oxidase and superoxide dismutase activities, respectively. We expressed an azalea GLP (RmGLP2) in cultured cells of tobacco, and determined that the extracellular protein fraction and the recombinant RmGLP2 protein purified from these cells catalyzed the oxidation of oxalate. Notably, this activity is purportedly restricted to germin and has not been demonstrated for a GLP. Although the specific activity of the purified RmGLP2 protein was low compared with that of a previously characterized barley germin/oxalate oxidase, tobacco cells expressing RmGLP2 exhibited significantly reduced oxalate levels. Thus, RmGLP2 represents the first reported GLP with oxalate oxidase activity.

Molecular characterization of atmospheric NO2-responsive germin-like proteins in azalea leaves.

Atmospheric nitrogen dioxide (NO(2)) is an environmental oxidant that is removed through direct uptake by foliage, but plant responses to this highly reactive gas are not well understood at the molecular level. From NO(2)-exposed leaves of a woody azalea (Rhododendron mucronatum), we cloned two cDNAs (RmGLP1 and RmGLP2) for germin-like proteins (GLPs), a group of ubiquitous plant proteins that have been implicated in various plant physiological and developmental processes. Quantitative analysis of mRNA expression, together with immunoblotting data, showed that foliar exposure to NO(2) caused a robust induction of these GLP-encoding genes. When produced in tobacco cell culture, recombinant RmGLP2 was secreted into the apoplast, where it exhibited superoxide dismutase activity. RmGLP1 and RmGLP2 represent the first examples of plant genes that are responsive to airborne NO(2). These enzymes might have a potential role in extracellular defense mechanisms through attenuation of interactions between reactive nitrogen and oxygen species.

Ploidy analysis of azalea flower colour sports.

Flower colour variegation is not only a phenomenon of importance to horticulture, the phenotype involved is also often used as a scientific model system for the study of complex gene regulation processes. In the course of such studies on azalea, we observed a correlation between flower colour patterns, flower morphology and somatic polyploidy. Using high-resolution flow cytometry of nuclear DNA, the ploidy level was determined in flowers of different azalea sport families. Sports exhibiting variegated flowers with broad (> 7mm), differently coloured, petal edges (picotee type) proved to be tetraploid in the petal edge while diploid in the rest of the flower tissue. Neither flower colour pattern nor ploidy differences are chimeral in origin, but seem to be correlated with the topographic location of the cells within the flower tissue, i.e. the margin of the petals. The possible role of gene dosage effects and cell size involved in the remarkable correlation between somatic polyploidy, (flavonoid) gene expression and the flower morphology is discussed.

In vitro polyploidy induction in Rhododendron simsii-hybrids.

Three in vitro protocols for the induction of tetraploid Rhododendron simsii hybrids were evaluated. The most successful technique was a daily application of mitotic poison on the apical meristem of seedlings. However, this technique was rather labour intensive. Direct sowing on media enriched with these compounds required less work but was less efficient. Mitotic inhibitors, when added to multiplication media, induced ploidy chimeras of existing genotypes.