Exposure to acidic and alkaline conditions were found to cause the excess accumulation of reactive oxygen species in tree peony, thereby causing damage and inhibiting plant growth and development. The activities of antioxidant enzymes were also found to be significantly up-regulated, especially under alkaline conditions; this explained why tree peony is better adapted to alkaline than to acidic conditions. Through pairwise comparisons, 144 differentially expressed genes (DEGs) associated with plant growth, photosynthesis, and stress were identified. The DEGs related to stress were up-regulated, whereas the remaining DEGs were almost all down-regulated after acid and alkaline treatments. The nutrient assimilation was greatly inhibited. Chlorophyll synthesis genes were suppressed, and chlorophyll content was reduced. The development and structures of stomata and chloroplasts and the transcription of related genes were also influenced. Among photosynthesis-related DEGs, electron transport chains were the most sensitive. The suppressed expression of photosynthesis genes and the reduced light-harvesting capacity, together with the impairment of chloroplasts and stomata, finally led to a sharp decrease in the net photosynthetic rate. Carbohydrate accumulation and plant biomass were also reduced. The present study provides a theoretical basis for the response mechanisms of tree peony to adverse pH conditions and enriches knowledge of plant adaptation to alkaline conditions.

Tree peony (Paeonia suffruticosa Andrews) is a globally famous ornamental flower, with large and colorful flowers and abundant flower types. However, a relatively short and uniform flowering period hinders the applications and production of ornamental tree peony. Unfortunately, the molecular mechanism of regulating flowering time and floral organ development in tree peony has yet to be elucidated. Because of the absence of genomic information, 454-based transcriptome sequence technology for de novo transcriptomics was used to identify the critical flowering genes using re-blooming, non-re-blooming, and wild species of tree peonies. A total of 29,275 unigenes were obtained from the bud transcriptome, with an N50 of 776 bp. The average length of unigenes was 677.18 bp, and the longest sequence was 5815 bp. Functional annotation showed that 22,823, 17,321, 13,312, 20,041, and 9940 unigenes were annotated by NCBI-NR, Swiss-Prot, COG, GO, and KEGG, respectively. Within the differentially expressed genes (DEGs) 64 flowering-related genes were identified and some important flowering genes were also characterized by bioinformatics methods, reverse transcript polymerase chain reaction (RT-PCR), and rapid-amplification of cDNA ends (RACE). Then, the putative genetic network of flowering induction pathways and a floral organ development model were put forward, according to the comparisons of DEGs in any two samples and expression levels of the important flowering genes in differentiated buds, buds from different developmental stages, and with GA or vernalization treated. In tree peony, five pathways (long day, vernalization, autonomous, age, and gibberellin) regulated flowering, and the floral organ development followed an ABCE model. Moreover, it was also found that the genes PsAP1, PsCOL1, PsCRY1, PsCRY2, PsFT, PsLFY, PsLHY, PsGI, PsSOC1, and PsVIN3 probably regulated re-blooming of tree peony. This study provides a comprehensive report on the flowering-related genes in tree peony for the first time and investigated the expression levels of the critical flowering related genes in buds of different cultivars, developmental stages, differentiated primordium, and flower parts. These results could provide valuable insights into the molecular mechanisms of flowering time regulation and floral organ development.