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

Functional Characterization of PoEP1 in Regulating the Flowering Stage of Tree Peony

The tree peony, a traditional flower in China, has a short and concentrated flowering period, restricting the development of the tree peony industry. To explore the molecular mechanism of tree peony flowering-stage regulation, PoEP1, which regulated the flowering period, was identified and cloned based on the transcriptome and degradome data of the early-flowering mutant Paeonia ostii 'Fengdan' (MU) and Paeonia ostii 'Fengdan' (FD). Through bioinformatics analysis, expression pattern analysis, and transgene function verification, the role of PoEP1 in the regulation of tree peony flowering was explored. The open-reading frame of PoEP1 is 1161 bp, encoding 386 amino acids, containing two conserved domains. PoEP1 was homologous to the EP1 of other species. Subcellular localization results showed that the protein was localized in the cell wall and that PoEP1 expression was highest in the initial decay stage of the tree peony. The overexpression of PoEP1 in transgenic plants advanced and shortened the flowering time, indicating that PoEP1 overexpression promotes flowering and senescence and shorten the flowering time of plants. The results of this study provide a theoretical basis for exploring the role of PoEP1 in the regulation of tree peony flowering.

Multi-locus genome-wide association studies reveal the dynamic genetic architecture of flowering time in chrysanthemum

KEY MESSAGE

The dynamic genetic architecture of flowering time in chrysanthemum was elucidated by GWAS. Thirty-six known genes and 14 candidate genes were identified around the stable QTNs and QEIs, among which ERF-1 was highlighted. Flowering time (FT) adaptation is one of the major breeding goals in chrysanthemum, a multipurpose ornamental plant. In order to reveal the dynamic genetic architecture of FT in chrysanthemum, phenotype investigation of ten FT-related traits was conducted on 169 entries in 2 environments. The broad-sense heritability of five non-conditional FT traits, i.e., budding (FBD), visible coloring (VC), early opening (EO), full-bloom (OF) and decay period (DP), ranged from 56.93 to 84.26%, which were higher than that of the five derived conditional FT traits (38.51-75.13%). The phenotypic variation coefficients of OF_EO and DP_OF were relatively large ranging from 30.59 to 36.17%. Based on 375,865 SNPs, the compressed variance component mixed linear model 3VmrMLM was applied for a multi-locus genome-wide association study (GWAS). As a result, 313 quantitative trait nucleotides (QTNs) were identified for the non-conditional FT traits in single-environment analysis, while 119 QTNs and 67 QTN-by-environment interactions (QEIs) were identified in multi-environment analysis. As for the conditional traits, 343 QTNs were detected in single-environment analysis, and 119 QTNs and 83 QEIs were identified in multi- environment analysis. Among the genes around stable QTNs and QEIs, 36 were orthologs of known FT genes in Arabidopsis and other plants; 14 candidates were mined by combining the transcriptomics data and functional annotation, including ERF-1, ACA10, and FOP1. Furthermore, the haplotype analysis of ERF-1 revealed six elite accessions with extreme FBD. Our findings contribute to the understanding of dynamic genetic architecture of FT and provide valuable resources for future chrysanthemum molecular breeding programs.

Multi-locus genome-wide association study and genomic prediction for flowering time in chrysanthemum

MAIN CONCLUSION

Multi-locus GWAS detected several known and candidate genes responsible for flowering time in chrysanthemum. The associations could greatly increase the predictive ability of genome selection that accelerates the possible application of GS in chrysanthemum breeding. Timely flowering is critical for successful reproduction and determines the economic value for ornamental plants. To investigate the genetic architecture of flowering time in chrysanthemum, a multi-locus genome-wide association study (GWAS) was performed using a collection of 200 accessions and 330,710 single-nucleotide polymorphisms (SNPs) via 3VmrMLM method. Five flowering time traits including budding (FBD), visible colouring (VC), early opening (EO), full-bloom (OF) and senescing (SF) stages, plus five derived conditional traits were recorded in two environments. Extensive phenotypic variations were observed for these flowering time traits with coefficients of variation ranging from 6.42 to 38.27%, and their broad-sense heritability ranged from 71.47 to 96.78%. GWAS revealed 88 stable quantitative trait nucleotides (QTNs) and 93 QTN-by-environment interactions (QEIs) associated with flowering time traits, accounting for 0.50-8.01% and 0.30-10.42% of the phenotypic variation, respectively. Amongst the genes around these stable QTNs and QEIs, 21 and 10 were homologous to known flowering genes in Arabidopsis; 20 and 11 candidate genes were mined by combining the functional annotation and transcriptomics data, respectively, such as MYB55, FRIGIDA-like, WRKY75 and ANT. Furthermore, genomic selection (GS) was assessed using three models and seven unique marker datasets. We found the prediction accuracy (PA) using significant SNPs identified by GWAS under SVM model exhibited the best performance with PA ranging from 0.90 to 0.95. Our findings provide new insights into the dynamic genetic architecture of flowering time and the identified significant SNPs and candidate genes will accelerate the future molecular improvement of chrysanthemum.

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.

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.

Identification of Floral Volatile Components and Expression Analysis of Controlling Gene in Paeonia ostii 'Fengdan' under Different Cultivation Conditions

In order to explore the release rule of floral volatile substances and the diurnal variation of different flower development stages of Paeonia ostii 'Fengdan' in potted and ground-planted conditions, dynamic headspace adsorption combined with gas chromatography-mass spectrometry(GC-MS) was used to analyze the dynamic changes in floral volatile components and contents. Quantitative real-time PCR (qRT-PCR) was used to analyze changes in flower fragrance-regulating genes PsPAL, PsTPSs, and PsbHLH at different flower development stages and a daily change process at the full-blooming stage. The results show that there were differences in aroma components and contents of Paeonia ostii 'Fengdan' at different flower development stages and different time quantum of every day. There were 25 and 28 aroma components identified in 7 flower development stages of tree peonies planted in pots and in the field, respectively, and 23 and 22 aroma components identified at different time quantum of the day, of which the largest and highest content was alkanes. The main characteristic aroma substances were (E)-β-ocimene, 1,3,5-trimethoxybenzene, 2,4-di-tert-butylphenol, methyl jasmonate, nerol, and cinnamyl alcohol; released amounts of the abovementioned substances varied depending on the development stage and the time of the day. The expression of flower fragrance-controlling genes (PsPAL, PsTPSs, and PsbHLH) in tree peonies varied greatly in different conditions. The results of this study provide a valuable resource to investigate floral fragrance formation in tree peonies.