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

Cell wall remodeling confers plant architecture with distinct wall structure in Nelumbo nucifera

Lotus (Nelumbo nucifera G.) is a perennial aquatic horticultural plant with diverse architectures. Distinct plant architecture (PA) has certain attractive and practical qualities, but its genetic morphogenesis in lotus remains elusive. In this study, we employ genome-wide association analysis (GWAS) for the seven traits of petiole length (PLL), leaf length (LL), leaf width (LW), peduncle length (PLF), flower diameter (FD), petal length (PeL), and petal width (PeW) in 301 lotus accessions. A total of 90 loci are identified to associate with these traits across 4 years of trials. Meanwhile, we perform RNA sequencing (RNA-seq) to analyze the differential expression of the gene (DEG) transcripts between large and small PA (LPA and SPA) of lotus stems (peduncles and petioles). As a result, eight key candidate genes are identified that are all primarily involved in plant cell wall remodeling significantly associated with PA traits by integrating the results of DEGs and GWAS. To verify this result, we compare the cell wall compositions and structures of LPA versus SPA in representative lotus germplasms. Intriguingly, compared with the SPA lotus, the LPA varieties have higher content of cellulose and hemicellulose, but less filling substrates of pectin and lignin. Additionally, we verified longer cellulose chains and higher cellulose crystallinity with less interference in LPA varieties. Taken together, our study illustrates how plant cell wall remodeling affects PA in lotus, shedding light on the genetic architecture of this significant ornamental trait and offering a priceless genetic resource for future genomic-enabled breeding.

Functional Analysis of the PoSERK-Interacting Protein PorbcL in the Embryogenic Callus Formation of Tree Peony (Paeonia ostii T. Hong et J. X. Zhang)

SERK is a marker gene for early somatic embryogenesis. We screened and functionally verified a SERK-interacting protein to gain insights into tree-peony somatic embryogenesis. Using PoSERK as bait, we identified PorbcL (i.e., the large subunit of Rubisco) as a SERK-interacting protein from a yeast two-hybrid (Y2H) library of cDNA from developing tree-peony somatic embryos. The interaction between PorbcL and PoSERK was verified by Y2H and bimolecular fluorescence complementation analyses. PorbcL encodes a 586-amino-acid acidic non-secreted hydrophobic non-transmembrane protein that is mainly localized in the chloroplast and plasma membrane. PorbcL was highly expressed in tree-peony roots and flowers and was up-regulated during zygotic embryo development. PorbcL overexpression caused the up-regulation of PoSERK (encoding somatic embryogenesis receptor-like kinase), PoAGL15 (encoding agamous-like 15), and PoGPT1 (encoding glucose-6-phosphate translocator), while it caused the down-regulation of PoLEC1 (encoding leafy cotyledon 1) in tree-peony callus. PorbcL overexpression led to increased indole-3-acetic acid (IAA) content but decreasing contents of abscisic acid (ABA) and 6-benzyladenosine (BAPR). The changes in gene expression, high IAA levels, and increased ratio of IAA to ABA, BAPR, 1-Aminocyclopropanecarboxylic acid (ACC), 5-Deoxystrigol (5DS), and brassinolide (BL) promoted embryogenesis. These results provide a foundation for establishing a tree-peony embryogenic callus system.

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.