The complete chloroplast genome of Rhododendron delavayi (Ericaceae)

Assemble and comparative analysis of the mitochondrial genome of Rhododendron delavayi: Insights into phylogenetic relationships and genomic variations

Rhododendron delavayi, a notable ornamental plant primarily found in regions of China like Yunnan and Guizhou provinces, holds substantial horticultural value. To elucidate the systematic phylogenetic relationships and organelle genomic differences within R. delavayi and related Rhododendron species, we conducted sequencing and assembly of the complete mitochondrial genome of R. delavayi. The full-length mitochondrial genome of it was a singular circular molecule spanning 1,009,263 bp, comprising 53 protein-coding genes, including 18 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 32 protein-coding genes. A total of 1,182 simple sequence repeats (SSRs) loci were identified in the R. delavayi mitochondrial genome, primarily consisting of single nucleotide, dinucleotide, and trinucleotide repeats. Nucleotide diversity analysis highlighted five genes (atp6, atp9, cox2, nad1, and rpl10) with the highest diversity within the mitochondrial genomes of Rhododendron genus. Comparative analysis of the mitochondrial genome of R. delavayi with those of four other Rhododendron species indicated complex rearrangements in 21 genes, including rps4, nad6, rps3, atp6, cob, atp9, nad7, among others. The mitochondrial phylogenetic tree revealed a close relationship between R. delavayi and R. decorum, forming a sister clade to R. × pulchrum and R. simsii. Furthermore, 126 plastid-to-mitochondrial gene transfers in R. delavayi were identified, ranging from 30 bp to 19,385 bp. These fragments collectively constituted 47.54 % and 9.52 % of the chloroplast and mitochondrial genomes (202,169 bp), respectively. Complex mitochondrial-to-mitochondrial transfers were also observed, with 843 identified fragments totaling 312,036 bp (30.92 % of the mitochondrial genome). Segments exceeding 10 kb may mediate homologous recombination within the mitochondrial molecules. Remarkably, our study underscores that the mitochondrial genome of R. delavayi was the largest reported within the Rhododendron genus to date. The intricate rearrangements observed in the mitochondrial genomes of Rhododendron species, alone with the identification of five potential molecular marker sites, provided valuable insights for species classification and parentage identification within the Rhododendron genus.

Comparison of the Complete Chloroplast Genomes of Astilbe: Two Korean Endemic Plant Species

BACKGROUND

Astilbe, consisting of about 18 species, is distributed throughout East Asia and Northeastern America, and most Astilbe species are widely cultivated as ornamental plants. A total of four species of Astilbe have been confirmed to be distributed throughout Korea, two of which are endemic to Korea.

METHODS

In this study, we sequenced and assembled the complete chloroplast genomes of two endemic Korean plants using Illumina sequencing technology, identified simple sequence repeats (SSRs) and repetitive sequences, and compared them with three previously reported chloroplast genomes.

RESULTS

The chloroplast genomes of the two species were 156,968 and 57,142 bp in length and had a four-part circular structure. They consisted of a large single-copy region of 87,223 and 87,272 bp and a small single-copy region of 18,167 and 18,138 bp, separated by a pair of inverted repeats (IRa and IRb, 25,789 and 25,866 bp). The genomes contained 130 genes, 49 SSRs, and 49 long repetitive sequences. Comparative analysis with the chloroplast genomes of five Astilbe species indicated that A. uljinensis was closely related to A. chinensis and A. taquetii to A. koreana.

CONCLUSIONS

This study provides valuable references for the identification of two endemic Korean Astilbe species and contributes to a deeper understanding of the phylogeny and evolution of the genus Astilbe.

The complete chloroplast genome of rhododendron williamsianum (ericaceae)

Rhododendron williamsianum Rehder & E. H. Wilson 1913, is a plant with important horticultural value. Here we report its chloroplast genome. The total length of the chloroplast genome was 205,424 bp, with a GC content of 35.8%. It consisted of a 107,968 bp large single copy, a 2606 bp small single copy, and a pair of 47,425 bp inverted repeats separating them. Within the chloroplast genome, there were a total of 110 unique genes, which included 76 protein-coding genes, 4 rRNA genes, and 30 tRNA genes. Our phylogenetic analyses indicated that R. williamsianum was closely genetically related to R. sutchuenense and R. jingangshanicum. The findings from this study not only contribute to the genetic database of Rhododendron plants but also have implications for evolutionary research within the family Ericaceae.

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

Comparative analysis of the complete chloroplast genome of Papaveraceae to identify rearrangements within the Corydalis chloroplast genome

Chloroplast genomes are valuable for inferring evolutionary relationships. We report the complete chloroplast genomes of 36 Corydalis spp. and one Fumaria species. We compared these genomes with 22 other taxa and investigated the genome structure, gene content, and evolutionary dynamics of the chloroplast genomes of 58 species, explored the structure, size, repeat sequences, and divergent hotspots of these genomes, conducted phylogenetic analysis, and identified nine types of chloroplast genome structures among Corydalis spp. The ndh gene family suffered inversion and rearrangement or was lost or pseudogenized throughout the chloroplast genomes of various Corydalis species. Analysis of five protein-coding genes revealed simple sequence repeats and repetitive sequences that can be potential molecular markers for species identification. Phylogenetic analysis revealed three subgenera in Corydalis. Subgenera Cremnocapnos and Sophorocapnos represented the Type 2 and 3 genome structures, respectively. Subgenus Corydalis included all types except type 3, suggesting that chloroplast genome structural diversity increased during its differentiation. Despite the explosive diversification of this subgenus, most endemic species collected from the Korean Peninsula shared only one type of genome structure, suggesting recent divergence. These findings will greatly improve our understanding of the chloroplast genome of Corydalis and may help develop effective molecular markers.

Phylogenetic relationship and characterization of the complete chloroplast genome of Rhododendron przewalskii subsp. przewalskii Maximowicz 1877, an ornamental and medicinal plant in China

As one of the top 10 famous flowers in Chinese tradition, Rhododendron przewalskii subsp. przewalskii known as 'beauty in flowers,' which has high ornamental and medicinal value. The complete chloroplast (cp) genome of R. przewalskii subsp. przewalskii was determined in this study. The complete chloroplast genome of R. przewalskii subsp. przewalskii was 201,233 bp in length and contained a large single-copy region (LSC, 108,077 bp), and a small single-copy region (SSC, 2624 bp) and a pair of inverted repeat regions (IRa and IRb, 45,266 bp). A total of 142 functional genes were observed in this cp genome, including 91 protein-coding genes (PCGs), 43 transfer RNA genes (tRNAs), and eight ribosomal RNA genes (rRNAs). The R. przewalskii subsp. przewalskii cp genome has an A + T content of 64.06% and presents a positive AT-skew (0.53%) and a negative GC-skew (-1.56%). The maximum likelihood phylogenetic analysis based on the concatenated nucleotide sequences of 13 PCGs strongly supported the monophyletic relationship of R. przewalskii subsp. przewalskii the clade of R. henanense subsp. lingbaoense. This study provides genomic evidence for the vegetation classification of Rhododendron.

Comprehensive Biochemical, Physiological, and Transcriptomic Analyses Provide Insights Into Floral Bud Dormancy in Rhododendron delavayi Franch

Due to a scarcity of relevant data, the ornamental woody flower Rhododendron delavayi Franch. is examined in the current study for its low temperature-induced floral bud dormancy (late October-end December) aspect. This study used transcriptome data profiling and co-expression network analyses to identify the interplay between endogenous hormones and bud dormancy phases such as pre-dormancy, para-dormancy, endo-dormancy, eco-dormancy, and dormancy release. The biochemical and physiological assays revealed the significance of the abundance of phytohormones (abscisic acid, auxin, zeatin, and gibberellins), carbohydrate metabolism, oxidative species, and proteins (soluble proteins, proline, and malondialdehyde) in the regulatory mechanism of floral bud dormancy. The transcriptome sequencing generated 65,531 transcripts, out of which 504, 514, 307, and 240 expressed transcripts were mapped uniquely to pre-, para-, endo-, and eco-phases of dormancy, showing their roles in the stimulation of dormancy. The transcripts related to LEA29, PGM, SAUR family, RPL9e, ATRX, FLOWERING LOCUS T, SERK1, ABFs, ASR2, and GID1 were identified as potential structural genes involved in floral bud dormancy. The transcription factors, including Zinc fingers, CAD, MADS-box family, MYB, and MYC2, revealed their potential regulatory roles concerning floral bud dormancy. The gene co-expression analysis highlighted essential hub genes involved in cold stress adaptations encoding proteins, viz, SERPIN, HMA, PMEI, LEA_2, TRX, PSBT, and AMAT. We exposed the connection among low temperature-induced dormancy in floral buds, differentially expressed genes, and hub genes via strict screening steps to escalate the confidence in selected genes as being truly putative in the pathways regulating bud dormancy mechanism. The identified candidate genes may prove worthy of further in-depth studies on molecular mechanisms involved in floral bud dormancy of Rhododendron species.

First Report of Colletotrichum fioriniae Causing Anthracnose on Rhododendron delavayi in China

Rhododendron delavayi Franch, a member of Ericaceae family, is globally famous for its garden flowers with significant ornamental value (Liu et al., 2020). In July 2020 and 2021, a disease survey of R. delavayi groves was conducted in Baili Azalea Forest Area (N27°10'-27°20', E 105°04'-106°04'). We arbitrarily selected an area with around 280 R. delavayi trees covering 2.5 hectares in R. delavayi grove where 20-35% of leaves showed symptoms of anthracnose. Typical symptoms included elliptical to irregularly shaped brown lesions on leaves and masses of black dots clustered on it. About 30 pieces of leaves with anthracnose lesions were collected. A few black dots were picked from the lesions with a sterilized needle, plated on water agar and incubated at 25℃ for 12 h to observe spore germination (Fang, 2007). Then the germinated spores were transferred onto PDA medium for further purification and morphological observation. Fourteen single-spore isolates with similar morphology were obtained. The surface of the colony was white or gray and spongy; the edge was smooth; and the back side was pinkish brown after 7 days of growth on PDA. Conidia were spindle-shaped, transparent, 11.1-16.6×3.6-4.9 μm (n=50). Appressorium from conidia was nearly ovate or proximate, brown or dark brown in color, 4.3-10.3 ×3.2-7.6 μm (n=50). These characteristics are consistent with Colletotrichum fioriniae reported by Shivas and Tan (2009). DNA was extracted from a representative isolate MYDJ12. The internal transcribed spacer region (ITS), the glyceraldehyde-3-phosphate dehydrogenase (GAPDH), β-tubulin (TUB), actin (ACT), and chitin synthase 1 (CHS-1) genes were amplified using primer pairs described by Damm et al. (2012). The sequences were deposited in GenBank with accession number MW692854 (ITS), MW727518 (GAPDH), MW727519 (TUB2), MW727520 (ACT), and MW727521 (CHS-1). BLASTN searches of the ITS, GADPH, TUB2, ACT and CHS-1 genes revealed 100% (540/540 nucleotides), 100% (254/254 nucleotides), 99.38% (4488/491 nucleotides), 98.77% (242/245 nucleotides) and 100% (282/282 nucleotides) homology with those of C. fioriniae CBS:128517T in GenBank (NR_111747, JQ948622, JQ949943, JQ949613 and JQ948953 respectively). The phylogenetic tree showed the isolate MYDJ12 to cluster with C. fioriniae CBS:128517T. Finally, two-year old R. delavayi plants (n=5) were inoculated by wounding with a syringe needle and placing 10 μL of spore suspension (106 spores per mL) of the isolate MYDJ12 on three leaves per plant. Control leaves were inoculated with sterile water. The experiment was conducted twice. Inoculated leaves were wrapped in parafilm tape and then the plants were placed in a greenhouse at 25°C with high relative humidity (90 to 95%). Seven days after incubation, brown lesions appeared, similar to those observed in the grove. Black dots clustered on the lesions after 15 days. Re-isolation was conducted 20 days after inoculation. From all the five inoculated plants, similar symptoms were observed, and the same pathogen was re-isolated. One of the isolates was selected for morphological observation and multi-gene (ITS, GAPDH, ACT, TUB2 and CHS-1) analysis indicated the reisolated fungus to be C. fioriniae. No fungal pathogens were isolated from mock inoculated plants. This study can provide effective management and useful information for the control of this disease on R. delavayi in Baili Azalea Forest Area. References: Damm, U., et al. 2012. Stud Mycol 73: 37. Fang, Z. D. 2007. Research Methods of Plant Diseases (Third edition). China Agriculture Press. Liu, J., et al. 2020. Mitochondrial DNA B 5:37. Shivas, R, G; Tan, Y, P. 2009. Fungal Divers 39:111.

Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron

Standard plant DNA barcodes based on 2-3 plastid regions, and nrDNA ITS show variable levels of resolution, and fail to discriminate among species in many plant groups. Genome skimming to recover complete plastid genome sequences and nrDNA arrays has been proposed as a solution to address these resolution limitations. However, few studies have empirically tested what gains are achieved in practice. Of particular interest is whether adding substantially more plastid and nrDNA characters will lead to an increase in discriminatory power, or whether the resolution limitations of standard plant barcodes are fundamentally due to plastid genomes and nrDNA not tracking species boundaries. To address this, we used genome skimming to recover near-complete plastid genomes and nuclear ribosomal DNA from Rhododendron species and compared discrimination success with standard plant barcodes. We sampled 218 individuals representing 145 species of this species-rich and taxonomically difficult genus, focusing on the global biodiversity hotspots of the Himalaya-Hengduan Mountains. Only 33% of species were distinguished using ITS+matK+rbcL+trnH-psbA. In contrast, 55% of species were distinguished using plastid genome and nrDNA sequences. The vast majority of this increase is due to the additional plastid characters. Thus, despite previous studies showing an asymptote in discrimination success beyond 3-4 plastid regions, these results show that a demonstrable increase in discriminatory power is possible with extensive plastid genome data. However, despite these gains, many species remain unresolved, and these results also reinforce the need to access multiple unlinked nuclear loci to obtain transformative gains in species discrimination in plants.

The complete chloroplast genome of Rhododendron datiandingense (Ericaceae)

Rhododendron datiandingense is newly reported and endemic to China. The genome of R. datiandingense is 207,311 bp in length, including a large single-copy region of 190,689 bp and a small single-copy region of 2582 bp, a pair of inverted repeat regions (IRA) of 7020 bp each. The genome encodes 110 genes, comprising 77 protein-coding genes, four ribosomal RNA genes, and 29 transfer RNA genes. Repeat analysis revealed 62 simple sequence repeats (SSRs) in the genome. Phylogenetic analysis revealed that R. datiandingense is clearly separated from the other Rhododendron species and shown in the basal position.

The Rhododendron Plant Genome Database (RPGD): a comprehensive online omics database for Rhododendron

Background

The genus Rhododendron L. has been widely cultivated for hundreds of years around the world. Members of this genus are known for great ornamental and medicinal value. Owing to advances in sequencing technology, genomes and transcriptomes of members of the Rhododendron genus have been sequenced and published by various laboratories. With increasing amounts of omics data available, a centralized platform is necessary for effective storage, analysis, and integration of these large-scale datasets to ensure consistency, independence, and maintainability.

Results

Here, we report our development of the Rhododendron Plant Genome Database (RPGD; http://bioinfor.kib.ac.cn/RPGD/), which represents the first comprehensive database of Rhododendron genomics information. It includes large amounts of omics data, including genome sequence assemblies for R. delavayi, R. williamsianum, and R. simsii, gene expression profiles derived from public RNA-Seq data, functional annotations, gene families, transcription factor identification, gene homology, simple sequence repeats, and chloroplast genome. Additionally, many useful tools, including BLAST, JBrowse, Orthologous Groups, Genome Synteny Browser, Flanking Sequence Finder, Expression Heatmap, and Batch Download were integrated into the platform.

Conclusions

RPGD is designed to be a comprehensive and helpful platform for all Rhododendron researchers. Believe that RPGD will be an indispensable hub for Rhododendron studies.

Seven Complete Chloroplast Genomes from Symplocos: Genome Organization and Comparative Analysis

In the present study, chloroplast genome sequences of four species of Symplocos (S. chinensis for. pilosa, S. prunifolia, S. coreana, and S. tanakana) from South Korea were obtained by Ion Torrent sequencing and compared with the sequences of three previously reported Symplocos chloroplast genomes from different species. The length of the Symplocos chloroplast genome ranged from 156,961 to 157,365 bp. Overall, 132 genes including 87 functional genes, 37 tRNA genes, and eight rRNA genes were identified in all Symplocos chloroplast genomes. The gene order and contents were highly similar across the seven species. The coding regions were more conserved than the non-coding regions, and the large single-copy and small single-copy regions were less conserved than the inverted repeat regions. We identified five new hotspot regions (rbcL, ycf4, psaJ, rpl22, and ycf1) that can be used as barcodes or species-specific Symplocos molecular markers. These four novel chloroplast genomes provide basic information on the plastid genome of Symplocos and enable better taxonomic characterization of this genus.

Complete plastid genome of Rhododendron griersonianum, a critically endangered plant with extremely small populations (PSESP) from southwest China

The complete plastid genome of Rhododendron griersonianum, a critically endangered plant species with extremely small populations, was obtained using Illumina HiSeq X Ten and ONT PromethION sequencing. The full length of the plastid genome is 206,467 bp with an overall GC content of 35.8%, which encodes 118 unique genes, including 78 protein-coding genes, 36 tRNA and 4 rRNA genes. Phylogenetic analysis revealed that all Rhododendron species formed a monophyletic clade. This study provides a valuable reference and will facilitate future studies related to the general characteristics and evolution of plastid genomes in the genus Rhododendron.

Characterization of the complete chloroplast genome sequence of Koelreuteria bipinnata

Koelreuteria bipinnata is an important ornamental tree with attractive flowers and fruits. In this study, we used next-generation sequencing technology to obtain the complete chloroplast genome of K. bipinnata. The entire genome was determined to be 163,863 bp in size, harboring a typical quadripartite structure with a large single copy (LSC) region of 90,240 bp, a small single copy (SSC) region of 18,883 bp, and a pair of 27,370 bp inverted repeat (IR) regions. The genome was predicted to contain 132 genes, including 84 protein-coding genes, 40 tRNA genes, and 8 rRNA genes. The overall GC content of K. bipinnata chloroplast genome was 37.29%. Phylogenetic analysis based on complete chloroplast genome sequences indicated that K. bipinnata was closely related to K. paniculate. This study would be useful for future population genetics studies and phylogenetic analysis of K. bipinnata.