Draft genome of the famous ornamental plant Paeonia suffruticosa

Molecular dissection of the parental contribution in Paeonia Itoh hybrids

Hybrid breeding between herbaceous peonies (the maternal parent) and tree peonies (the paternal parent) results in Paeonia Itoh hybrids (Itoh peonies), a triploid species that combines advantageous traits from both parental species, thus offering great economic value. However, the exact genetic contribution of the two parents is unclear. In this study, we introduce a straightforward approach utilizing heterozygous single-nucleotide polymorphisms (SNPs) and Sanger sequencing of targeted gene fragments to trace the original bases back to their parents in Itoh peonies. Our results indicate that in triploid Itoh peonies, only one set of genes is derived from herbaceous peonies, and two sets of genes are derived from the tree peonies. Notably, the presence of three distinct bases of heterozygous SNPs across multiple Itoh cultivars suggests that the gametes from the paternal parents carry two sets of heterozygous homologous chromosomes, which could be due to Meiosis I failure during gamete formation. To validate our method's effectiveness in parentage determination, we analyze two Itoh hybrids and their parents, confirming its practical utility. This research presents a method to reveal the parental genetic contribution in Itoh peonies, which could enhance the efficiency and precision of hybrid breeding programs of triploids in Paeonia and other plant species.

Synergistic actions of 3 MYB transcription factors underpin blotch formation in tree peony

Blotches in floral organs attract pollinators and promote pollination success. Tree peony (Paeonia suffruticosa Andr.) is an internationally renowned cut flower with extremely high ornamental and economic value. Blotch formation on P. suffruticosa petals is predominantly attributed to anthocyanin accumulation. However, the endogenous regulation of blotch formation in P. suffruticosa remains elusive. Here, we identified the regulatory modules governing anthocyanin-mediated blotch formation in P. suffruticosa petals, which involves the transcription factors PsMYB308, PsMYBPA2, and PsMYB21. PsMYBPA2 activated PsF3H expression to provide sufficient precursor substrate for anthocyanin biosynthesis. PsMYB21 activated both PsF3H and PsFLS expressions and promoted flavonol biosynthesis. The significantly high expression of PsMYB21 in nonblotch regions inhibited blotch formation by competing for anthocyanin biosynthesis substrates, while conversely, its low expression in the blotch region promoted blotch formation. PsMYB308 inhibited PsDFR and PsMYBPA2 expressions to directly prevent anthocyanin-mediated blotch formation. Notably, a smaller blotch area, decreased anthocyanin content, and inhibition of anthocyanin structural gene expression were observed in PsMYBPA2-silenced petals, while the opposite phenotypes were observed in PsMYB308-silenced and PsMYB21-silenced petals. Additionally, PsMYBPA2 and PsMYB308 interacted with PsbHLH1-3, and their regulatory intensity on target genes was synergistically regulated by the PsMYBPA2-PsbHLH1-3 and PsMYB308-PsbHLH1-3 complexes. PsMYB308 also competitively bound to PsbHLH1-3 with PsMYBPA2 to fine-tune the regulatory network to prevent overaccumulation of anthocyanin in blotch regions. Overall, our study uncovers a complex R2R3-MYB transcriptional regulatory network that governs anthocyanin-mediated blotch formation in P. suffruticosa petals, providing insights into the molecular mechanisms underlying blotch formation in P. suffruticosa.

Genome-Wide Identification of NAC Gene Family Members of Tree Peony (Paeonia suffruticosa Andrews) and Their Expression under Heat and Waterlogging Stress

An important family of transcription factors (TFs) in plants known as NAC (NAM, ATAF1/2, and CUC2) is crucial for the responses of plants to environmental stressors. In this study, we mined the NAC TF family members of tree peony (Paeonia suffruticosa Andrews) from genome-wide data and analyzed their response to heat and waterlogging stresses in conjunction with transcriptome data. Based on tree peony's genomic information, a total of 48 PsNAC genes were discovered. Based on how similar their protein sequences were, these PsNAC genes were divided into 14 branches. While the gene structures and conserved protein motifs of the PsNAC genes within each branch were largely the same, the cis-acting elements in the promoter region varied significantly. Transcriptome data revealed the presence of five PsNAC genes (PsNAC06, PsNAC23, PsNAC38, PsNAC41, PsNAC47) and one PsNAC gene (PsNAC37) in response to heat and waterlogging stresses, respectively. qRT-PCR analysis reconfirmed the response of these five PsNAC genes to heat stress and one PsNAC gene to waterlogging stress. This study lays a foundation for the study of the functions and regulatory mechanisms of NAC TFs in tree peony. Meanwhile, the NAC TFs of tree peony in response to heat and waterlogging stress were excavated, which is of great significance for the selection and breeding of new tree peony varieties with strong heat and waterlogging tolerance.

Selection of Stable Reference Genes for QRT-PCR in Tree Peony 'Doulv' and Functional Analysis of PsCUC3

(1) Background: Tree peonies display extensive cultivar diversity due to widespread hybridization, resulting in a complex genetic architecture. This complexity complicates the selection of universal reference genes across different cultivars for qRT-PCR analyses. Paeonia suffruticosa 'Doulv', notable for its unique green blooms in China, exhibits chlorosis post-flowering and features petaloid stamens and pistils. (2) Methods: Based on published literature and RNA-seq data from 'Doulv', nine candidate reference genes-ACT (Actin), TUB (β-Tubulin), UBC (Ubiquitin Conjugating Enzyme), UBQ (Ubiquitin), UPL (Ubiquitin Protein Ligase), PP2A (Protein Phosphatase 2A), PP2C (Protein Phosphatase 2C), MBF1A (Multiprotein Bridging Factor 1A), and GAPDH (Glyceraldehyde-3-Phosphate Dehydrogenase)-were selected. Their expression stability was assessed across various tissues and developmental stages of 'Doulv' flowers using qRT-PCR, with evaluations conducted via GeNorm_v3.5, NormFinder_v20, and BestKeeper_v1.0. Gene cloning and expression analyses of PsCUC3, including its subcellular localization, were performed. (3) Results: GAPDH and ACT were identified as the most stable reference genes in petaloid stamens across various developmental stages of 'Doulv', whereas UBC and MBF1A were optimal across different tissues. Notably, specific conserved amino acids in PsCUC3 from 'Doulv' diverged from those in NAM/CUC3 proteins of other species, impacting its protein structure. PsCUC3 expression analysis revealed no correlation with chlorophyll content in petaloid stamens but an association with petaloid organ development. Furthermore, PsCUC3 was predominantly localized in the nucleus. (4) Conclusions: This study comprehensively evaluated suitable reference genes using GeNorm_v3.5, NormFinder_v20, and BestKeeper_v1.0 software, establishing a robust qRT-PCR detection system for 'Doulv' peony. These results provide a solid experimental foundation for further research on 'Doulv' peony. Building on this experimental foundation, the functional analysis of the PsCUC3 gene was conducted. The findings suggest a potential association between the PsCUC3 gene and floral morphology alterations in 'Doulv', identifying PsCUC3 as crucial for understanding the molecular mechanisms influencing floral structure in tree peonies.

Molecular and genetic regulation of petal number variation

Floral forms with an increased number of petals, also known as double-flower phenotypes, have been selected and conserved in many domesticated plants, particularly in ornamentals, because of their great economic value. The molecular and genetic mechanisms that control this trait are therefore of great interest, not only for scientists, but also for breeders. In this review, we summarize current knowledge of the gene regulatory networks of flower initiation and development and known mutations that lead to variation of petal number in many species. In addition to the well-accepted miR172/AP2-like module, for which many questions remain unanswered, we also discuss other pathways in which mutations also lead to the formation of extra petals, such as those involved in meristem maintenance, hormone signalling, epigenetic regulation, and responses to environmental signals. We discuss how the concept of 'natural mutants' and recent advances in genomics and genome editing make it possible to explore the molecular mechanisms underlying double-flower formation, and how such knowledge could contribute to the future breeding and selection of this trait in more crops.

Identification of ARF Genes and Elucidation of the Regulatory Effects of PsARF16a on the Dormancy of Tree Peony Plantlets

The low survival rate of transplanted plantlets, which has limited the utility of tissue-culture-based methods for the rapid propagation of tree peonies, is due to plantlet dormancy after rooting. We previously determined that the auxin response factor PsARF may be a key regulator of tree peony dormancy. To clarify the mechanism mediating tree peony plantlet dormancy, PsARF genes were systematically identified and analyzed. Additionally, PsARF16a was transiently expressed in the leaves of tree peony plantlets to examine its regulatory effects on a downstream gene network. Nineteen PsARF genes were identified and divided into four classes. All PsARF genes encoded proteins with conserved B3 and ARF domains. The number of motifs, exons, and introns varied between PsARF genes in different classes. The overexpression of PsARF16a altered the expression of NCED, ZEP, PYL, GA2ox1, GID1, and other key genes in abscisic acid (ABA) and gibberellin (GA) signal transduction pathways, thereby promoting ABA synthesis and decreasing GA synthesis. Significant changes to the expression of some key genes contributing to starch and sugar metabolism (e.g., AMY2A, BAM3, BGLU, STP, and SUS2) may be associated with the gradual conversion of sugar into starch. This study provides important insights into PsARF functions in tree peonies.

The Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome

Paeonia lactiflora (P. lactiflora), a perennial plant renowned for its medicinal roots, provides a unique case for studying the phylogenetic relationships of species based on organelle genomes, as well as the transference of DNA across organelle genomes. In order to investigate this matter, we sequenced and characterized the mitochondrial genome (mitogenome) of P. lactiflora. Similar to the chloroplast genome (cpgenome), the mitogenome of P. lactiflora extends across 181,688 base pairs (bp). Its unique quadripartite structure results from a pair of extensive inverted repeats, each measuring 25,680 bp in length. The annotated mitogenome includes 27 protein-coding genes, 37 tRNAs, 8 rRNAs, and two pseudogenes (rpl5, rpl16). Phylogenetic analysis was performed to identify phylogenetic trees consistent with Paeonia species phylogeny in the APG Ⅳ system. Moreover, a total of 12 MTPT events were identified and 32 RNA editing sites were detected during mitogenome analysis of P. lactiflora. Our research successfully compiled and annotated the mitogenome of P. lactiflora. The study provides valuable insights regarding the taxonomic classification and molecular evolution within the Paeoniaceae family.

Insights into the Superrosids phylogeny and flavonoid synthesis from the telomere-to-telomere gap-free genome assembly of Penthorum chinense Pursh

The completion of the first telomere-to-telomere (T2T) genome assembly of Penthorum chinense Pursh (PC), a prominent medicinal plant in China, represents a significant achievement. This assembly spans a length of 257.5 Mb and consists of nine chromosomes. PC's notably smaller genome size in Saxifragales, compared to that of Paeonia ostii, can be attributed to the low abundance of transposable elements. By utilizing single-copy genes from 30 species, including 28 other Superrosids species, we successfully resolved a previously debated Superrosids phylogeny. Our findings unveiled Saxifragales as the sister group to the core rosids, with both being the sister group to Vitales. Utilizing previously characterized cytochrome P450 (CYP) genes, we predicted the compound classes that most CYP genes of PC are involved in synthesizing, providing insight into PC's potential metabolic diversity. Metabolomic and transcriptomic data revealed that the richest sources of the three most noteworthy medicinal components in PC are young leaves and flowers. We also observed higher activity of upstream genes in the flavonoid synthesis pathway in these plant parts. Additionally, through weighted gene co-expression network analysis, we identified gene regulatory networks associated with the three medicinal components. Overall, these findings deepen our understanding of PC, opening new avenues for further research and exploration.

High-quality assembly and methylome of a Tibetan wild tree peony genome (Paeonia ludlowii) reveal the evolution of giant genome architecture

Tree peony belongs to one of the Saxifragales families, Paeoniaceae. It is one of the most famous ornamental plants, and is also a promising woody oil plant. Although two Paeoniaceae genomes have been released, their assembly qualities are still to be improved. Additionally, more genomes from wild peonies are needed to accelerate genomic-assisted breeding. Here we assemble a high-quality and chromosome-scale 10.3-Gb genome of a wild Tibetan tree peony, Paeonia ludlowii, which features substantial sequence divergence, including around 75% specific sequences and gene-level differentials compared with other peony genomes. Our phylogenetic analyses suggest that Saxifragales and Vitales are sister taxa and, together with rosids, they are the sister taxon to asterids. The P. ludlowii genome is characterized by frequent chromosome reductions, centromere rearrangements, broadly distributed heterochromatin, and recent continuous bursts of transposable element (TE) movement in peony, although it lacks recent whole-genome duplication. These recent TE bursts appeared during the uplift and glacial period of the Qinghai-Tibet Plateau, perhaps contributing to adaptation to rapid climate changes. Further integrated analyses with methylome data revealed that genome expansion in peony might be dynamically affected by complex interactions among TE proliferation, TE removal, and DNA methylation silencing. Such interactions also impact numerous recently duplicated genes, particularly those related to oil biosynthesis and flower traits. This genome resource will not only provide the genomic basis for tree peony breeding but also shed light on the study of the evolution of huge genome structures as well as their protein-coding genes.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

RAD-sequencing improves the genetic characterization of a threatened tree peony (Paeonia ludlowii) endemic to China: Implications for conservation

Compared with traditional genetic markers, genomic approaches have proved valuable to the conservation of endangered species. Paeonia ludlowii having rarely and pure yellow flowers, is one of the world's most famous tree peonies. However, only several wild populations remain in the Yarlung Zangbo Valley (Nyingchi and Shannan regions, Xizang) in China due to increasing anthropogenic impact on the natural habitats. We used genome-wide single nucleotide polymorphisms to elucidate the spatial pattern of genetic variation, population structure and demographic history of P. ludlowii from the fragmented region comprising the entire range of this species, aiming to provide a basis for conserving the genetic resources of this species. Unlike genetic uniformity among populations revealed in previous studies, we found low but varied levels of intra-population genetic diversity, in which lower genetic diversity was detected in the population in Shannan region compared to those in Nyingzhi region. These spatial patterns may be likely associated with different population sizes caused by micro-environment differences in these two regions. Additionally, low genetic differentiation among populations (Fst = 0.0037) were detected at the species level. This line of evidence, combined with the result of significant genetic differentiation between the two closest populations and lack of isolation by distance, suggested that shared ancestry among now remnant populations rather than contemporary genetic connectivity resulted in subtle population structure. Demographic inference suggested that P. ludlowii probably experienced a temporal history of sharp population decline during the period of Last Glacial Maximum, and a subsequent bottleneck event resulting from prehistoric human activities on the Qinghai-Tibet Plateau. All these events, together with current habitat fragment and excavation might contribute to the endangered status of P. ludlowii. Our study improved the genetic characterization of the endangered tree peony (P. ludlowii) in China, and these genetic inferences should be considered when making different in situ and ex situ conservation actions for P. ludlowii in this evolutionary hotspot region.

Combined genome-wide association studies and expression quantitative trait locus analysis uncovers a genetic regulatory network of floral organ number in a tree peony (Paeonia suffruticosa Andrews) breeding population

Great progress has been made in our understanding of floral organ identity determination and its regulatory network in many species; however, the quantitative genetic basis of floral organ number variation is far less well understood for species-specific traits from the perspective of population variation. Here, using a tree peony (Paeonia suffruticosa Andrews, Paeoniaceae) cultivar population as a model, the phenotypic polymorphism and genetic variation based on genome-wide association studies (GWAS) and expression quantitative trait locus (eQTL) analysis were analyzed. Based on 24 phenotypic traits of 271 representative cultivars, the transcript profiles of 119 cultivars were obtained, which indicated abundant genetic variation in tree peony. In total, 86 GWAS-related cis-eQTLs and 3188 trans-eQTL gene pairs were found to be associated with the numbers of petals, stamens, and carpels. In addition, 19 floral organ number-related hub genes with 121 cis-eQTLs were obtained by weighted gene co-expression network analysis, among which five hub genes belonging to the ABCE genes of the MADS-box family and their spatial-temporal co-expression and regulatory network were constructed. These results not only help our understanding of the genetic basis of floral organ number variation during domestication, but also pave the way to studying the quantitative genetics and evolution of flower organ number and their regulatory network within populations.

Tree peony PsMYB44 negatively regulates petal blotch distribution by inhibiting dihydroflavonol-4-reductase gene expression

BACKGROUND AND AIMS

The tree peony (Paeonia suffruticosa Andr.) has been widely cultivated as a field plant, and petal blotch is one of its important traits, which not only promotes proliferation but also confers high ornamental value. However, the regulatory network controlling blotch formation remains elusive owing to the functional differences and limited conservation of transcriptional regulators in dicots.

METHODS

We performed phylogenetic analysis to identify MYB44-like transcription factors in P. suffruticosa blotched cultivar 'High noon' petals. A candidate MYB44-like transcription factor, PsMYB44, was analysed via expression pattern analysis, subcellular localization, target gene identification, gene silencing in P. suffruticosa petals and heterologous overexpression in tobacco.

KEY RESULTS

A blotch formation-related MYB44-like transcription factor, PsMYB44, was cloned. The C-terminal of the PsMYB44 amino acid sequence had a complete C2 motif that affects anthocyanin biosynthesis, and PsMYB44 was clustered in the MYB44-like transcriptional repressor branch. PsMYB44 was located in the nucleus, and its spatial and temporal expression patterns were negatively correlated with blotch formation. Furthermore, a yeast one-hybrid assay showed that PsMYB44 could target the promoter of the late anthocyanin biosynthesis-related dihydroflavonol-4-reductase (DFR) gene, and a dual-luciferase assay demonstrated that PsMYB44 could repress PsDFR promoter activity. On the one hand, overexpression of PsMYB44 significantly faded the red colour of tobacco flowers and decreased the anthocyanin content by 42.3 % by downregulating the expression level of the tobacco NtDFR gene. On the other hand, PsMYB44-silenced P. suffruticosa petals had a redder blotch colour, which was attributed to the fact that silencing PsMYB44 redirected metabolic flux to the anthocyanin biosynthesis branch, thereby promoting more anthocyanin accumulation in the petal base.

CONCLUSION

These results demonstrated that PsMYB44 negatively regulated the biosynthesis of anthocyanin by directly binding to the PsDFR promoter and subsequently inhibiting blotch formation, which helped to elucidate the molecular regulatory network of anthocyanin-mediated blotch formation in plants.

An R2R3-MYB network modulates stem strength by regulating lignin biosynthesis and secondary cell wall thickening in herbaceous peony

Stem strength is an important agronomic trait affecting plant lodging, and plays an essential role in the quality and yield of plants. Thickened secondary cell walls in stems provide mechanical strength that allows plants to stand upright, but the regulatory mechanism of secondary cell wall thickening and stem strength in cut flowers remains unclear. In this study, first, a total of 11 non-redundant Paeonia lactiflora R2R3-MYBs related to stem strength were identified and isolated from cut-flower herbaceous peony, among which PlMYB43, PlMYB83 and PlMYB103 were the most upregulated differentially expressed genes. Then, the expression characteristics revealed that these three R2R3-MYBs were specifically expressed in stems and acted as transcriptional activators. Next, biological function verification showed that these P. lactiflora R2R3-MYBs positively regulated stem strength, secondary cell wall thickness and lignin deposition. Furthermore, yeast-one-hybrid and dual luciferase reporter assays demonstrated that they could bind to the promoter of caffeic acid O-methyltransferase gene (PlCOMT2) and/or laccase gene (PlLAC4), two key genes involved in lignin biosynthesis. In addition, the function of PlLAC4 in increasing lignin deposition was confirmed by virus-induced gene silencing and overexpression. Moreover, PlMYB83 could also act as a transcriptional activator of PlMYB43. The findings of the study propose a regulatory network of R2R3-MYBs modulating lignin biosynthesis and secondary cell wall thickening for improving stem lodging resistance, and provide a resource for molecular genetic engineering breeding of cut flowers.

Genome-wide association study of 23 flowering phenology traits and 4 floral agronomic traits in tree peony (Paeonia section Moutan DC.) reveals five genes known to regulate flowering time

Tree peony is a unique traditional flower in China, with large, fragrant, and colorful flowers. However, a relatively short and concentrated flowering period limits the applications and production of tree peony. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of flowering phenology traits and ornamental phenotypes in tree peony. A diverse panel of 451 tree peony accessions was phenotyped for 23 flowering phenology traits and 4 floral agronomic traits over 3 years. Genotyping by sequencing (GBS) was used to obtain a large number of genome-wide single-nucleotide polymorphisms (SNPs) (107 050) for the panel genotypes, and 1047 candidate genes were identified by association mapping. Eighty-two related genes were observed during at least 2 years for flowering, and seven SNPs repeatedly identified for multiple flowering phenology traits over multiple years were highly significantly associated with five genes known to regulate flowering time. We validated the temporal expression profiles of these candidate genes and highlighted their possible roles in the regulation of flower bud differentiation and flowering time in tree peony. This study shows that GWAS based on GBS can be used to identify the genetic determinants of complex traits in tree peony. The results expand our understanding of flowering time control in perennial woody plants. Identification of markers closely related to these flowering phenology traits can be used in tree peony breeding programs for important agronomic traits.

Creating novel ornamentals via new strategies in the era of genome editing

Ornamental breeding has traditionally focused on improving novelty, yield, quality, and resistance to biotic or abiotic stress. However, achieving these goals has often required laborious crossbreeding, while precise breeding techniques have been underutilized. Fortunately, recent advancements in plant genome sequencing and editing technology have opened up exciting new frontiers for revolutionizing ornamental breeding. In this review, we provide an overview of the current state of ornamental transgenic breeding and propose four promising breeding strategies that have already proven successful in crop breeding and could be adapted for ornamental breeding with the help of genome editing. These strategies include recombination manipulation, haploid inducer creation, clonal seed production, and reverse breeding. We also discuss in detail the research progress, application status, and feasibility of each of these tactics.

Identification of differentially expressed miRNAs and their target genes in response to brassinolide treatment on flowering of tree peony (Paeonia ostii)

Tree peony is a famous flower plant in China, but the short and concentrated flowering period limits its ornamental value and economic value. Brassinolide (BR) plays an important role in plant growth and development including flowering. There have been a large number of reports on the molecular aspects of the flowering process, but the genetic mechanism that was responsible for miRNA-guided regulation of tree peony is almost unclear. In this study, the leaves of tree peony cultivar, 'Feng Dan', were sprayed with different concentrations of BR, and the obvious bloom delay was found at the treatment with BR 50 μg/L. The small RNA sequencing and transcriptome sequencing were performed on the petals of tree peony under an untreated control (CK) and the treatment with BR 50 μg/L during four consecutive flowering development stages. A total of 22 known miRNAs belonging to 12 families were identified and 84 novel miRNAs were predicted. Combined with transcriptome data, a total of 376 target genes were predicted for the 18 differentially expressed known miRNAs and 177 target genes were predicted for the 23 differentially expressed novel miRNAs. Additionally, the potential miRNAs and their target genes were identified, including miR156b targeting SPL, miR172a_4 targeting AP2 and four novel miRNAs targeting SPA1, and revealed that they might affect the flowering time in tree peony. Collectively, these results would provide a theoretical basis for further analysis of miRNA-guided regulation on flowering period in tree peony.

Responses of the tree peony (Paeonia suffruticosa, Paeoniaceae) cultivar 'Yu Hong' to heat stress revealed by iTRAQ-based quantitative proteomics

Horticulture productivity has been increasingly restricted by heat stress from growing global warming, making it far below the optimum production capacity. As a popular ornamental cultivar of tree peony, Paeonia suffruticosa 'Yu Hong' has also been suffering from heat stress not suitable for its optimal growth. To better understand the response mechanisms against heat stress of tree peony, investigations of phenotypic changes, physiological responses, and quantitative proteomics were conducted. Phenotypic and physiological changes indicated that 24 h of exposure to heat stress (40 °C) was the critical duration of heat stress in tree peony. The proteomic analyses revealed a total of 100 heat-responsive proteins (HRPs). According to bioinformatic analysis of HRPs, the heat tolerance of tree peony might be related to signal transduction, synthesis/degradation, heat kinetic proteins, antioxidants, photosynthesis, energy conversion, and metabolism. Our research will provide some new insights into the molecular mechanism under the response against the heat stress of tree peony, which will benefit the future breeding of heat-resistant ornamental plants.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Integrated analysis of miRNAome transcriptome and degradome reveals miRNA-target modules governing floral florescence development and senescence across early- and late-flowering genotypes in tree peony

As a candidate national flower of China, tree peony has extremely high ornamental, medicinal and oil value. However, the short florescence and rarity of early-flowering and late-flowering varieties restrict further improvement of the economic value of tree peony. Specific miRNAs and their target genes engaged in tree peony floral florescence, development and senescence remain unknown. This report presents the integrated analysis of the miRNAome, transcriptome and degradome of tree peony petals collected from blooming, initial flowering, full blooming and decay stages in early-flowering variety Paeonia ostii 'Fengdan', an early-flowering mutant line of Paeonia ostii 'Fengdan' and late-flowering variety Paeonia suffruticosa 'Lianhe'. Transcriptome analysis revealed a transcript ('psu.G.00014095') which was annotated as a xyloglucan endotransglycosylase/hydrolase precursor XTH-25 and found to be differentially expressed across flower developmental stages in Paeonia ostii 'Fengdan' and Paeonia suffruticosa 'Lianhe'. The miRNA-mRNA modules were presented significant enrichment in various pathways such as plant hormone signal transduction, indole alkaloid biosynthesis, arachidonic acid metabolism, folate biosynthesis, fatty acid elongation, and the MAPK signaling pathway. Multiple miRNA-mRNA-TF modules demonstrated the potential functions of MYB-related, bHLH, Trihelix, NAC, GRAS and HD-ZIP TF families in floral florescence, development, and senescence of tree peony. Comparative spatio-temporal expression investigation of eight floral-favored miRNA-target modules suggested that transcript 'psu.T.00024044' and microRNA mtr-miR166g-5p are involved in the floral florescence, development and senescence associated agronomic traits of tree peony. The results might accelerate the understanding of the potential regulation mechanism in regards to floral florescence, development and abscission, and supply guidance for tree peony breeding of varieties with later and longer florescence characteristics.

Association study of SNP locus for color related traits in herbaceous peony (Paeonia lactiflora Pall.) using SLAF-seq

Paeonia lactiflora Pall. (P. lactiflora) is a famous ornamental plant with showy and colorful flowers that has been domesticated in China for 4,000 years. However, the genetic basis of phenotypic variation and genealogical relationships in P. lactiflora population is poorly understood due to limited genetic information, which brings about bottlenecks in the application of effective and efficient breeding strategies. Understanding the genetic basis of color-related traits is essential for improving flower color by marker-assisted selection (MAS). In this study, a high throughput sequencing of 99 diploid P. lactiflora accessions via specific-locus amplified fragment sequencing (SLAF-seq) technology was performed. In total, 4,383,645 SLAF tags were developed from 99 P. lactiflora accessions with an average sequencing depth of 20.81 for each SLAF tag. A total of 2,954,574 single nucleotide polymorphisms (SNPs) were identified from all SLAF tags. The population structure and phylogenetic analysis showed that P. lactiflora population used in this study could be divided into six divergent groups. Through association study using Mixed linear model (MLM), we further identified 40 SNPs that were significantly positively associated with petal color. Moreover, a derived cleaved amplified polymorphism (dCAPS) marker that was designed based on the SLAF tag 270512F co-segregated with flower colors in P. lactiflora population. Taken together, our results provide valuable insights into the application of MAS in P. lactiflora breeding programs.

Genomic basis of the giga-chromosomes and giga-genome of tree peony Paeonia ostii

Tree peony (Paeonia ostii) is an economically important ornamental plant native to China. It is also notable for its seed oil, which is abundant in unsaturated fatty acids such as α-linolenic acid (ALA). Here, we report chromosome-level genome assembly (12.28 Gb) of P. ostii. In contrast to monocots with giant genomes, tree peony does not appear to have undergone lineage-specific whole-genome duplication. Instead, explosive LTR expansion in the intergenic regions within a short period (~ two million years) may have contributed to the formation of its giga-genome. In addition, expansion of five types of histone encoding genes may have helped maintain the giga-chromosomes. Further, we conduct genome-wide association studies (GWAS) on 448 accessions and show expansion and high expression of several genes in the key nodes of fatty acid biosynthetic pathway, including SAD, FAD2 and FAD3, may function in high level of ALAs synthesis in tree peony seeds. Moreover, by comparing with cultivated tree peony (P. suffruticosa), we show that ectopic expression of class A gene AP1 and reduced expression of class C gene AG may contribute to the formation of petaloid stamens. Genomic resources reported in this study will be valuable for studying chromosome/genome evolution and tree peony breeding.

Characteristics of PoVIN3, a Key Gene of Vernalization Pathway, Affects Flowering Time

The tree peony (Paeonia section Moutan DC.) is the candidate flower in China, with abundant germplasm resources and high ornamental value. However, the short and concentrated flowering period severely restricted the improvement of the economic value of tree peonies. Based on the full-length transcriptome database of tree peonies, the PoVIN3 (GenBank ID: OP341879), involved in the flowering regulation of tree peonies were identified and cloned for the first time. The PoVIN3 was also characterized by bioinformatics methods, quantitative real-time PCR (qRT-PCR), and the establishment of a transgenic system. The expression levels of PoVIN3 in seven different petals developmental stages were the highest at the initial flowering stage of the variant cultivar of Paeonia ostii 'Fengdan,' the initial decay stage of the normal flowering Paeonia ostii 'Fengdan,' and the half opening stage of the late flowering Paeonia suffruticosa 'Lianhe.' Tissue-specific expression analysis showed that the relative expression levels of PoVIN3 were the highest in sepals of both normal flowering Paeonia ostii 'Fengdan' and the late flowering Paeonia suffruticosa 'Lianhe,' and the highest expression was in stamens of early flowering mutant Paeonia ostii 'Fengdan.' In addition, the flowering time of pCAMBIA2300-PoVIN3 transgenic plants was significantly earlier than that of the wild-type, indicating that PoVIN3 could promote plant flowering. The results provide a theoretical basis for exploring the role of PoVIN3 in the regulation of flowering in tree peonies.

Tree Peony R2R3-MYB Transcription Factor PsMYB30 Promotes Petal Blotch Formation by Activating the Transcription of the Anthocyanin Synthase Gene

Petal blotches are commonly observed in many angiosperm families and not only influence plant-pollinator interactions but also confer high ornamental value. Tree peony (Paeonia suffruticosa Andr.) is an important cut flower worldwide, but few studies have focused on its blotch formation. In this study, anthocyanins were found to be the pigment basis for blotch formation of P. suffruticosa, and peonidin-3,5-di-O-glucoside (Pn3G5G) was the most important component of anthocyanins, while the dihydroflavonol-4-reductase gene was the key factor contributing to blotch formation. Then, the R2R3-myeloblastosis (MYB) transcription factor PsMYB30 belonging to subgroup 1 was proven as a positive anthocyanin regulator with transcriptional activation and nuclear expression. Furthermore, silencing PsMYB30 in P. suffruticosa petals reduced blotch size by 37.9%, faded blotch color and decreased anthocyanin and Pn3G5G content by 23.6% and 32.9%, respectively. Overexpressing PsMYB30 increased anthocyanin content by 14.5-fold in tobacco petals. In addition, yeast one-hybrid assays, dual-luciferase assays and electrophoretic mobility shift assays confirmed that PsMYB30 could bind to the promoter of the anthocyanin synthase (ANS) gene and enhance its expression. Altogether, a novel MYB transcription factor, PsMYB30, was identified to promote petal blotch formation by activating the expression of PsANS involved in anthocyanin biosynthesis, which provide new insights for petal blotch formation in plants.

Nutrient Supply Is Essential for Shifting Tree Peony Reflowering Ahead in Autumn and Sugar Signaling Is Involved

The flowering time of tree peony is short and concentrated in spring, which limits the development of its industry. We previously achieved tree peony reflowering in autumn. Here, we further shifted its reflowering time ahead through proper gibberellin (GA) treatment plus nutrient supply. GA treatment alone initiated bud differentiation, but it aborted later, whereas GA plus nutrient (G + N) treatment completed the opening process 38 days before the control group. Through microstructural observation of bud differentiation and starch grains, we concluded that GA plays a triggering role in flowering induction, whereas the nutriment supply ensured the continuous developing for final opening, and both are necessary. We further determined the expression of five floral induction pathway genes and found that PsSOC1 and PsLFY probably played key integral roles in flowering induction and nutrient supply, respectively. Considering the GA signaling, PsGA2ox may be mainly involved in GA regulation, whereas PsGAI may regulate further flower formation after nutrient application. Furthermore, G + N treatment, but not GA alone, inhibited the expression of PsTPS1, a key restricting enzyme in sugar signaling, at the early stage, indicating that sugar signaling is also involved in this process; in addition, GA treatment induced high expression of PsSnRK1, a major nutrient insufficiency indicator, and the induction of PsHXK1, a rate-limiting enzyme for synthesis of sugar signaling substances, further confirmed the nutrient shortage. In short, besides GA application, exogenous nutrient supply is essential to shift tree peony reflowering ahead in autumn under current forcing culture technologies.

DNA Demethylation Induces Tree Peony Flowering with a Low Deformity Rate Compared to Gibberellin by Inducing PsFT Expression under Forcing Culture Conditions

Gibberellin (GA) is frequently used in tree peony forcing culture, but inappropriate application often causes flower deformity. Here, 5-azacytidine (5-azaC), an efficient DNA demethylating reagent, induced tree peony flowering with a low deformity rate by rapidly inducing PsFT expression, whereas GA treatment affected various flowering pathway genes with strong pleiotropy. The 5-azaC treatment, but not GA, significantly reduced the methylation level in the PsFT promoter with the demethylation of five CG contexts in a 369 bp CG-rich region, and eight light-responsive related cis-elements were also predicted in this region, accompanied by enhanced leaf photosynthetic efficiency. Through GO analysis, all methylation-closer differentially expressed genes (DEGs) were located in the thylakoid, the main site for photosynthesis, and were mainly involved in response to stimulus and single-organism process, whereas GA-closer DEGs had a wider distribution inside and outside of cells, associated with 12 categories of processes and regulations. We further mapped five candidate DEGs with potential flowering regulation, including three kinases (SnRK1, WAK2, and 5PTase7) and two bioactive enzymes (cytochrome P450 and SBH1). In summary, 5-azaC and GA may have individual roles in inducing tree peony flowering, and 5-azaC could be a preferable regulation approach; DNA demethylation is suggested to be more focused on flowering regulation with PsFT playing a core role through promoter demethylation. In addition, 5-azaC may partially undertake or replace the light-signal function, combined with other factors, such as SnRK1, in regulating flowering. This work provides new ideas for improving tree peony forcing culture technology.

A Preliminary Investigation on the Functional Validation and Interactions of PoWOX Genes in Peony (Paeonia ostii)

As a woody plant, peony (Paeonia suffruticosa) has a long growth cycle and inefficient traditional breeding techniques. There is an urgent need in peony molecular breeding to establish an efficient and stable in vitro regeneration and genetic transformation system, in order to overcome the recalcitrant characteristics of peony regeneration and shorten the breeding cycle. The development of plant somatic embryos is an important way to establish an efficient and stable in vitro regeneration and genetic transformation system. Plant-specific WUSCHEL-related homeobox (WOX) family transcription factors play important roles in plant development, from embryogenesis to lateral organ development. Therefore, in this research, four PoWOX genes of “Fengdan” (Paeonia ostii) were cloned from the peony genome and transcriptome data of preliminary peony somatic embryos. The sequence characteristics and evolutionary relationships of the PoWOX genes were analyzed. It was demonstrated that the four PoWOX genes, named PoWOX1, PoWOX4, PoWOX11, and PoWOX13, belonged to three branches of the WOX gene family. Their expression patterns were analyzed at different stages of development and in different tissues of peony seedlings. The expression localization of the PoWOX genes was determined to be the nucleus via subcellular localization assay. Finally, the interaction protein of the PoWOX genes was identified via yeast two-hybrid assay combined with bimolecular fluorescence complementation assay. It was shown that PoWOX1 and PoWOX13 proteins could form homodimers by themselves, and PoWOX11 interacted with PoWOX1 and PoWOX13 to form heterodimers. Peony stem cell activity may be regulated from PoWOX1 and PoWOX13 by forming dimers and moving to peony stem cells through plasmodesmata. Additionally, PoWOX11–PoWOX1 and PoWOX11–PoWOX13 may play important regulatory functions in promoting the proliferation of stem cells and maintaining the homeostasis of stem cells in the SAM of peony stems. Exploring the critical genes and regulatory factors in the development of the peony somatic embryo is beneficial not only to understand the molecular and regulatory mechanisms of peony somatic embryo development but also to achieve directed breeding and improvements in efficiency through genetic engineering breeding technology to accelerate the fundamental process of molecular breeding in peony.

A chromosome-level genome assembly of an alpine plant Crucihimalaya lasiocarpa provides insights into high-altitude adaptation

It remains largely unknown how plants adapt to high-altitude habitats. Crucihimalaya (Brassicaceae) is an alpine genus occurring in the Qinghai-Tibet Plateau characterized by cold temperatures and strong ultraviolet radiation. Here, we generated a chromosome-level genome for C. lasiocarpa with a total size of 255.8 Mb and a scaffold N50 size of 31.9 Mb. We first examined the karyotype origin of this species and found that the karyotype of five chromosomes resembled the ancestral karyotype of the Brassicaceae family, while the other three showed strong chromosomal structural variations. In combination with the rough genome sequence of another congener (C. himalaica), we found that the significantly expanded gene families and positively selected genes involved in alpine adaptation have occurred since the origin of this genus. Our new findings provide valuable information for the chromosomal karyotype evolution of Brassicaceae and investigations of high-altitude environment adaptation of the genus.

RNA-seq and phytohormone analysis reveals the culm color variation of Bambusa oldhamii Munro

BACKGROUND

The clumping bamboo Bambusa oldhamii Munro, known as "green bamboo", is famous for its edible bamboo shoots and fast-growing timber. The green and yellow striped-culm B. oldhamii variety, named B. oldhamii f. revoluta W.T. Lin & J. Y. Lin, is an attractive system for researching the culm color variation of B. oldhamii.

METHODS

Millions of clean reads were generated and assembled into 604,900 transcripts, and 383,278 unigenes were acquired with RNA-seq technology. The quantification of ABA, IAA, JA, GA₁, GA₃, GA₄, and GA₇ was performed using HPLC-MS/MS platforms.

RESULTS

Differential expression analysis showed that 449 unigenes were differentially expressed genes (DEGs), among which 190 DEGs were downregulated and 259 DEGs were upregulated in B. oldhamii f. revoluta. Phytohormone contents, especially GA₁ and GA_7, were higher in B. oldhamii. Approximately 21 transcription factors (TFs) were differentially expressed between the two groups: the bZIP, MYB, and NF-YA transcription factor families had the most DEGs, indicating that those TFs play important roles in B. oldhamii culm color variation. RNA-seq data were confirmed by quantitative RT-PCR analysis of the selected genes; moreover, phytohormone contents, especially those of ABA, GA₁ and GA₇, were differentially accumulated between the groups. Our study provides a basal gene expression and phytohormone analysis of B. oldhamii culm color variation, which could provide a solid fundamental theory for investigating bamboo culm color variation.

Seed development in Paeonia ostii (Paeoniaceae), with particular reference to embryogeny

BACKGROUND

Seeds of Paeonia ostii have been proposed as a source of raw material for the production of edible oil; however, lack of information about the developmental biology of the seeds hampers our ability to use them. Our aim was to investigate development of the seed coat, endosperm and embryo of P. ostii in relation to timing of accumulation of nutrient reserves from pollination to seed maturity. Ovules and developing seeds of P. ostii were collected at various stages of development from zygote to maturity. Seed fresh mass, dry mass, germination, moisture, soluble sugars, starch, protein and oil content were determined. Ontogeny of seeds including embryo, endosperm and seed coat were analyzed histologically.

RESULTS

The ovule of P. ostii is anatropous, crassinucellate and bitegmic. The zygote begins to divide at about 5 days after pollination (DAP), and the division is not accompanied by cell wall formation. By 25 DAP, the proembryo begins to cellularize. Thereafter, several embryo primordia appear at the surface of the cellularized proembryo, but only one matures. Endosperm development follows the typical nuclear type. The seed coat is derived from the outer integument. During seed development, soluble sugars, starch and crude fat content increased and then decreased, with maximum contents at 60, 80 and 100 DAP, respectively. Protein content was relatively low compared with soluble sugars and crude fat, but it increased throughout seed development.

CONCLUSIONS

During seed development in P. ostii, the seed coat acts as a temporary storage tissue. Embryo development of P. ostii can be divided into two stages: a coenocytic proembryo from zygote (n + n) that degenerates and a somatic embryo from peripheral cells of the proembryo (2n → 2n). This pattern of embryogeny differs from that of all other angiosperms, but it is similar to that of gymnosperms.

The R2R3-MYB gene PsMYB58 positively regulates anthocyanin biosynthesis in tree peony flowers

Tree peony (Paeonia suffruticosa Andr.) is a well-known ornamental flower in China with diverse colors. Flower color is one of the most important economic characteristics of tree peony and is mainly determined by anthocyanins. In this study, we cloned a PsMYB58 gene, which contained a 654 bp open reading frame (ORF), encoding a polypeptide of 218 amino acids. Sequence and phylogenetic analysis indicated that PsMYB58 is an anthocyanin regulatory R2R3-MYB gene. The transcription levels of PsMYB58 in different developmental stages of tree peony flowers were similar to those of the anthocyanin biosynthetic genes PsCHS, PsCHI, PsDFR, and PsANS. A bimolecular fluorescence complementation assay showed that PsMYB58 interacted with PsbHLH1 and PsbHLH3 in vivo. The overexpression of PsMYB58 in tobacco enhanced anthocyanin accumulation in various organs. Comparative transcriptome analysis showed that 943 genes were upregulated and 1203 downregulated in PsMYB58 transgenic tobacco, among which genes involved in the anthocyanin pathway were positively activated. Real-time quantitative PCR analysis verified that anthocyanin biosynthetic genes, including NtCHS, NtCHI, NtF3H, NtF3'H, NtDFR, and NtANS, and an anthocyanin regulatory bHLH gene, NtAN1b, were significantly upregulated in PsMYB58 transgenic tobacco. Our results indicated that PsMYB58 is a positive anthocyanin regulator in tree peony flowers. In summary, the functional identification of PsMYB58 furthers our understanding of the mechanism of peony flower color formation, thus providing a foundation for flower color improvement and molecular breeding.

Research advances in and prospects of ornamental plant genomics

The term 'ornamental plant' refers to all plants with ornamental value, which generally have beautiful flowers or special plant architectures. China is rich in ornamental plant resources and known as the "mother of gardens". Genomics is the science of studying genomes and is useful for carrying out research on genome evolution, genomic variations, gene regulation, and important biological mechanisms based on detailed genome sequence information. Due to the diversity of ornamental plants and high sequencing costs, the progress of genome research on ornamental plants has been slow for a long time. With the emergence of new sequencing technologies and a reduction in costs since the whole-genome sequencing of the first ornamental plant (Prunus mume) was completed in 2012, whole-genome sequencing of more than 69 ornamental plants has been completed in <10 years. In this review, whole-genome sequencing and resequencing of ornamental plants will be discussed. We provide analysis with regard to basic data from whole-genome studies of important ornamental plants, the regulation of important ornamental traits, and application prospects.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.

Germplasm resources and genetic breeding of Paeonia: a systematic review

Members of the genus Paeonia, which consists of globally renowned ornamentals and traditional medicinal plants with a rich history spanning over 1500 years, are widely distributed throughout the Northern Hemisphere. Since 1900, over 2200 new horticultural Paeonia cultivars have been created by the discovery and breeding of wild species. However, information pertaining to Paeonia breeding is considerably fragmented, with fundamental gaps in knowledge, creating a bottleneck in effective breeding strategies. This review systematically introduces Paeonia germplasm resources, including wild species and cultivars, summarizes the breeding strategy and results of each Paeonia cultivar group, and focuses on recent progress in the isolation and functional characterization of structural and regulatory genes related to important horticultural traits. Perspectives pertaining to the resource protection and utilization, breeding and industrialization of Paeonia in the future are also briefly discussed.