Next Article in Journal
Identification of Biochemical Determinants for Diagnosis and Prediction of Severity in 5q Spinal Muscular Atrophy Using 1H-Nuclear Magnetic Resonance Metabolic Profiling in Patient-Derived Biofluids
Previous Article in Journal
Raloxifene Protects Oxygen-Glucose-Deprived Astrocyte Cells Used to Mimic Hypoxic-Ischemic Brain Injury
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integration of Genome-Wide Identification and Transcriptome Analysis of Class III Peroxidases in Paeonia ostii: Insight into Their Roles in Adventitious Roots, Heat Tolerance, and Petal Senescence

College of Landscape Architecture and Art, Henan Agricultural University, Zhengzhou 450002, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(22), 12122; https://doi.org/10.3390/ijms252212122
Submission received: 16 October 2024 / Revised: 4 November 2024 / Accepted: 6 November 2024 / Published: 12 November 2024
(This article belongs to the Section Molecular Plant Sciences)

Abstract

:
As a plant-specific gene family, class III peroxidases (PODs) play an important role in plant growth, development, and stress responses. However, the POD gene family has not been systematically studied in Paeonia ostii. In this study, a total of 57 PoPOD genes were identified in the P. ostii genome. Subsequently, phylogenetic analysis and chromosome mapping revealed that PoPODs were classified into six subgroups and were unevenly distributed across five chromosomes. The gene structure and conserved motifs indicated the potential for functional divergence among the different subgroups. Meanwhile, four PoPODs were identified as tandem duplicated genes, with no evidence of segmental duplication. Using RNA-seq data from eight different tissues, multiple PoPODs exhibited enhanced expression in apical and adventitious roots (ARs). Next, RNA-seq data from AR development combined with trend analysis showed that PoPOD30/34/43/46/47/57 are implicated in the formation of ARs in tree peony. Through WGCNA based on RNA-seq, two key genes, PoPOD5/15, might be involved in heat tolerance via ABA and MeJA signaling. In addition, real-time quantitative PCR (qRT-PCR) analysis indicated that PoPOD23 may play an important role in flower senescence. These findings deepened our understanding of POD-mediated AR development, heat tolerance, and petal senescence in tree peony.

1. Introduction

As a large category of enzymes, peroxidases (EC 1.11.1.X) are bifunctional enzymes present in living organisms, which can generate reactive oxygen species (ROS) or oxidize various substrates with H2O2 to maintain intracellular ROS levels [1]. Based on their functional attributes and protein structure, peroxidases can be classified into two groups, namely hemoglobin peroxidases and non-hemoglobin peroxidases [2]. In addition to animal peroxidases, hemoglobin peroxidases are further subdivided into three categories. Class I peroxidases are distributed in microorganisms and plants, and their principal function is to eliminate excess H2O2 to prevent cellular damage [3]. Examples include catalase (CAT) and ascorbate peroxidase (APX). Class II peroxidases are found in fungi and are involved in lignin degradation [4]. Examples include manganese and lignin peroxidases. Class III peroxidases (PRXs; EC 1.11.1.7) are specific to plants and are often referred to by abbreviations such as POX, PER, or POD [1,5]. Recently, this class has been referred to as POD.
The role of Class III POD genes has been extensively demonstrated in numerous studies, highlighting their crucial involvement in plant growth and development as well as stress responses. For instance, PODs participate in various biological processes, including ROS scavenging, lignification [6], cell wall metabolism [7], auxin metabolism [8], wound healing [9], seed germination [10], pollen–pistil interactions [11], fruit ripening [12], and root development [13]. Arabidopsis POD genes (AtPrx33 and Atprx34) were predominantly expressed in roots, and the overexpression of AtPrx34 promoted root elongation [14]. In cucumber, CsPrx73 overexpression enhanced the formation of AR and increased the elimination of ROS under waterlogged conditions [13]. Additionally, PODs play significant roles in the presence of abiotic stressors such as salt [15], drought [16], and cold stress [17], as well as with various biological stressors such as bacteria and fungi [18]. The rice POD gene (OsPrx114) was localized in the endoplasmic reticulum and plasma membrane. Compared with the control, OsPrx114 overexpression in rice increased the activity of POD and CAT, reduced the accumulation of ROS, and improved drought tolerance [16]. In sweet potato, the overexpression of IbPRX17 can maintain a high SOD and POD activity and proline content, reduce the accumulation of ROS, and enhance salt and drought tolerance [19]. Therefore, a genome-wide analysis of POD will facilitate our understanding of its role in biology, physiology, and stress responses. In recent years, the POD family has been identified in numerous species, such as Arabidopsis thaliana [20], Zea mays [21], Populus trichocarpa [7], Pyrus bretschneideri [22], and Capsicum annuum [23]. Nevertheless, to date, no previous systematic study of the POD family has been reported in tree peony.
The tree peony is a perennial woody plant belonging to Paeonia, Paeoniaceae, and serves as an important ornamental plant in China [24]. Simultaneously, it possesses ornamental and medicinal properties. Nevertheless, the tissue culture plantlets of tree peony exhibited poor rooting and low transplantation survival rates [25], severely limiting their large-scale production. Under natural circumstances, a single tree peony flower persisted for three to five days, and petal senescence significantly affected its ornamental period [26]. Exploring the response mechanism of petal senescence in tree peony holds great significance. During the production process, the hot weather in southern China caused leaf damage in tree peony [27], resulting in reduced planting effectiveness and a considerable decline in their ornamental value. Considering that the PODs play critical roles in root development, senescence, and stress responses, yet their function in tree peony remains largely unreported, a comprehensive identification and analysis of the PoPOD gene family is of considerable significance in revealing root development, flowering regulation, and heat tolerance in tree peony.
In this study, all members of the PoPOD gene family were identified utilizing P. ostii genome data [24]. Firstly, the physicochemical properties, phylogenetic relationships, chromosome localization, collinearity, conserved motifs, gene structure, and cis-elements were systematically analyzed. Secondly, the expression pattern of the PoPOD gene family in adventitious root (AR) formation was disclosed through trend analysis using transcriptome data of AR development in tree peony. Additionally, WGCNA analysis was conducted using heat-resistant transcriptome data to construct regulatory networks between transcription factors and PoPODs. Finally, qRT-PCR was employed to analyze the expression patterns of tree peony petals during flower senescence. The results of this study provide a foundation for further functional research on PoPODs.

2. Results

2.1. Identification and Characterization of PoPODs

Based on hmmsearch and BLASTP, a total of 57 PODs were identified in the P. ostii genome and designated as PoPOD1-PoPOD57 according to their chromosomal locations (Table 1). These 57 PoPOD proteins encoded 183 (PoPOD7) to 436 (PoPOD23) amino acids, with molecular weights ranging from 19733.31 (PoPOD7) to 47449.47 (PoPOD23). The theoretical pI values ranged from 4.58 (PoPOD20) to 9.51 (PoPOD7), of which 25 PoPOD proteins were acidic amino acids (pI < 7) and 32 were basic amino acids (pI > 7). The grand average of hydropathicity of most PoPODs was lower than zero, indicating that the majority of PoPODs were hydrophilic proteins. Most PoPODs were predicted to localize to the chloroplast, while seven were located extracellularly, six were within the cytoplasm, four had dual localization (extracellular/vacuole, mitochondrion/vacuole, chloroplast/vacuole, and chloroplast/nucleus), two were in the nucleus and endoplasmic reticulum and one was in the plasma membrane (PoPOD32).

2.2. Phylogenetic Analysis of PoPODs

To elucidate the evolutionary relationship of PoPOD proteins, a Maximum Likelihood (ML) phylogenetic tree was constructed using the amino acid sequences of 57 PoPODs and 73 Arabidopsis PODs (Figure 1). Based on the classification of the Arabidopsis POD family, the POD proteins were divided into six subgroups (I-VI). Phylogenetic analysis revealed that the PoPOD genes exhibited an uneven distribution across each subgroup. For instance, subgroup VI contained the highest number of PoPODs (25), while subgroup IV contained none. The number of Arabidopsis PODs was higher than that of the tree peony in all other subgroups, except for subgroup VI.

2.3. Conserved Motif and Gene Structure Analysis of PoPOD Family

The conserved domains of PoPOD proteins were predicted using the MEME website. A total of ten motifs were identified, ranging from motif 1 to motif 10. As shown in Figure 2, a higher number (55) of consensus sequences were found in motif 4, while fewer numbers (46) were observed in motifs 5 and 9. Additionally, 54 PoPODs, except for PoPOD10/38/56, contained motif 3, and 54 PoPOD genes, except for PoPOD4/15/23, contained motif 4. Excluding PoPOD7/10/21/38, 53 PoPODs contained motif 7. Other PoPOD proteins contained motifs 3, 4, and 7, with the order consistently being 7, 3, and 4. It was worth noting that PoPOD10/38 did not contain motifs 3 or 7, whereas PoPOD21 contained two copies of motif 4 (Figure 2). The structural analysis of PoPOD genes revealed that the number of exons ranged from 2 to 12. Most had four exons and three introns, with 27 genes (47.36%) having three introns, and 15 genes (26.32%) having two introns. Notably, the PoPOD23 gene had the highest number of exons and introns, with twelve exons and eleven introns, whereas PoPOD17/24 only had two exons and one intron (Figure S1). Additionally, most of the PoPODs in subgroups I and V contained three introns, suggesting a similarity between these subgroups. Subgroups II, III, and VI displayed a diverse range of intron and exon numbers. These findings indicated the potential for functional divergence among the different subgroups.

2.4. Chromosomal Localization and Gene Collinearity of PoPODs

To gain deeper insights into the evolutionary relationship within the POD gene family, a gene collinearity analysis between P. ostii and A. thaliana was constructed. The results revealed that seven homologous gene pairs were found in P. ostii and A. thaliana (Figure 3A). Among these, two PoPODs were located on chromosomes 01 and 04, with one PoPOD gene present on each of the other chromosomes.
To reveal the chromosomal distributions of PoPODs, a graph of chromosomes was drawn using the TBtools software (v. 2.119). A total of 56 out of the 57 PoPOD genes were mapped on the five chromosomes of P. ostii. As shown in Figure 3B, the number of PoPODs varied across the chromosomes. For instance, chromosomes 02, 04, and 05 contained 13 PoPODs, chromosome 01 contained 12 PoPODs, while only 5 PoPODs were located on chromosome 03. Two tandem duplication events were observed on chromosomes 04 (PoPOD36/PoPOD37) and 05 (PoPOD52/PoPOD53). Notably, both duplicated genes belong to subgroup VI.

2.5. Analysis of Cis-Regulatory Elements in PoPODs

The 2.0 kb upstream regions of all PoPODs were extracted from the P. ostii genome sequence as putative promoters, and their cis-elements were predicted using the PlantCARE website. Fifty-seven PoPODs contained 7359 TATA-box and CAAT-box elements, which are defined as common elements (Figure 4A). Additionally, other cis-elements were categorized into the following six groups: light response, phytohormone, stress response, transcription factor binding site, and plant tissue (Figure 4B,C). A significant number of light response cis-acting elements (589) were discovered in the promoters of 57 PoPODs, averaging 10.33 elements per gene. Notably, PoPOD57 had the highest number (20), while PoPOD13/45 each contained only 2 elements (Figure 4B,C). The hormone response elements included 149 ABA, 142 MeJA, 56 GA, 38 SA, 32 auxin, 25 zein, and 21 GA-related elements (Figure 4B,C). Among these, there were also 48 PoPODs containing ABA response elements (ABRE), with PoPOD47 having the highest number (14) (Figure 4B,C). Additionally, 38 PoPODs contained MeJA elements (CGtCA-motifs and TGACg-motifs); PoPOD27/57 had the highest number (8) (Figure 4B,C). For stress responses, the elements primarily included anaerobic (129), mixed stress (27), and low-temperature response elements (26). Transcription factor binding elements were mainly MYB binding sites, with a total of 110 (Figure 4B,C). Plant tissue-specific elements were also identified in some PoPODs (Figure 4B,C). These results suggested that PoPOD genes in tree peony may play significant roles in hormone signaling and environmental stress responses.

2.6. Expression Profiling of PoPOD Genes in Different Organs

To investigate tissue-specific expression patterns of the PoPOD genes in tree peony, transcriptome sequencing data were obtained from eight different tissues, including apical, adventitious roots, stems, leaves, petals, stamens, pistils, and seeds. As shown in Figure 5, PoPOD4/15/23/55 showed high expression levels in all tissues, while other PoPODs, such as PoPOD9/34/40/47/54, exhibited higher expression levels in the apical and ARs. In contrast, these five genes showed a relatively low expression in other tissues. Additionally, three genes (PoPOD27/44/57) were highly expressed in the stamens.

2.7. Expression Profiling and Trend Analysis of PoPOD Genes in an In Vitro AR Formation

In order to study the response of POD family genes within the in vitro ARs of tree peony, transcriptome data at three key time points during rooting were analyzed (Figure 6A). Initially, the genes were classified through trend analysis and divided into eight clusters according to their temporal expression trends. A total of 20 PoPODs were divided into eight clusters, with 6 PoPODs belonging to the eighth cluster and 4 PoPODs belonging to the fourth cluster (Figure 6B). The KEGG enrichment analysis revealed that genes within the eighth cluster were associated with the biosynthesis of secondary metabolites, starch and sucrose metabolism, MAPK signaling pathway-plant, endocytosis, and phenylpropanoid biosynthesis (Figure 6C). Subsequently, a total of eighty-six transcription factors including bHLH, C2H2, MYB, NAC, B3, and LBD were predicted using Plant TFDB, showing significant up-regulation on days five and ten (Figure 6D,E).

2.8. Expression Profiling and Co-Expression Network Analysis of PoPOD Genes in Response to High-Temperature Stress

To explore how PoPODs responded to heat resistance in tree peony, transcriptome sequencing data from high-temperature treatment (0 h, 2 h, 6 h, 12 h, 24 h) samples were analyzed. The results indicated that PoPOD15/22/23/39/50 initially decreased followed by an increase when compared with the control group, whereas PoPOD13/30/41/55 continued to decrease (Figure 7A). Furthermore, to elucidate the regulatory network involved in the PoPODs’ heat resistance response, a WGCNA co-expression network analysis was conducted. Initially, the WGCNA analysis identified nineteen modules, among which the blue module significantly correlated with heat resistance (Figure 7B). PoPOD5/15 were found to belong to this module. Subsequently, a KEGG analysis for genes in the blue module revealed significant enrichments in pathways, including the biosynthesis of secondary metabolites, protein processing in the endoplasmic reticulum, plant-pathogen interactions, amino sugar and nucleotide sugar metabolism, ascorbate and aldarate metabolism, fructose and mannose metabolism, alpha-Linolenic acid metabolism, and glycerolipid metabolism (Figure 7C).
Furthermore, predictions of transcription factors within the blue module mainly included ERF, Trihelix, TCP, bHLH, and MYB (Figure 7D). A noteworthy observation was that the heatmap of the transcription factors primarily divided into two modules; those that were downregulated or upregulated after heat treatment (Figure 7E). Several ERF gene expressions were upregulated after treatment, while some MYBs and GRASs expression were shown to be downregulated following high-temperature treatment (Figure 7E). Finally, a regulatory network between transcription factors and PoPOD5/15 was constructed based on their interactions in the blue module (Figure 7F).

2.9. qRT-PCR Analysis of PoPOD Genes in Flower Senescence

In order to investigate the role of PoPODs in the petal senescence of tree peony, five highly expressed PoPODs in petals were selected and analyzed using qRT-PCR following cut flower treatment at time points of one, two, three, four, and five days. During petal senescence, PoPOD4/15/41/55 exhibited a downregulation followed by an upregulation trend, while PoPOD23 showed an upregulation followed by a downregulation trend (Figure 8). The expression of PoPOD4/15/41 peaked at day four, while PoPOD23 reached its peak at day two. These results suggested that PoPOD genes were involved in petal senescence.

3. Discussion

PODs are plant-specific enzymes, which play a crucial role in plant growth and development, and stress responses. Given the significant biological function of PODs, the POD gene family has been extensively studied in model plants, crops [21], fruits [22], and trees [7]. Nevertheless, studies have yet to identify or characterize POD genes in tree peony. In this study, a total of 57 PoPOD genes were identified in the P. ostii genome, which is fewer than in Arabidopsis thaliana (73) [20], Populus trichocarpa (93) [7], and Pyrus bretschneideri (94) [22], but slightly more than in Vitis vinifera (47) [28]. This finding indicated that the POD gene family in P. ostii and Vitis vinifera has not expanded significantly compared with other plants. Subsequently, the physical and chemical properties, phylogenetic analysis, conserved motifs, gene structure, chromosome localization, collinearity analysis, cis-acting element analysis, and gene expression profile of PoPODs were comprehensively analyzed.
The PoPODs were categorized into six distinct subgroups (I-VI) through phylogenetic analysis, consistent with previous reports on Glycine max [29] and sugarcane [30]. Conserved motif analysis demonstrated that motif 4 was present in the majority of PoPOD genes. Meanwhile, some PoPOD proteins contained motifs 3, 4, and 7, and the order of these motifs remained consistent. These results indicated that PoPOD proteins possessed highly conserved attributes, similar to those reported in Sorghum mosaic [17] and Pyrus bretschneideri [22]. It is widely acknowledged that the evolution of multigene families is caused by the structural diversification of genes. The introns within PoPODs exhibited substantial variation. In rare instances, the number of introns reached a considerable amount, such as PoPOD23, which contained 11 introns. This finding indicated that PoPOD members have a certain diversity. The majority of PoPODs (73.68%) possessed only two or three introns, a pattern associated with stress-related genes having fewer introns, a phenomenon also observed in plants like Vitis vinifera [28], Glycine max [29], and Pyrus bretschneideri [22]. Additionally, most of the PoPODs in subgroups I and V contained three introns, while subgroups II, III, and VI had varying numbers of introns. These results suggested a structural conservation and diversification within and between PoPOD subgroups, supporting the potential functional diversity of PoPODs.
Whole-genome duplication, segmental duplication, and tandem duplication are the three principal modes of gene family expansion. Chromosome localization revealed that the PoPOD genes were widely distributed across five chromosomes, which was consistent with the distribution of PODs in Arabidopsis thaliana [20], Populus trichocarpa [7], and Pyrus bretschneideri [22]. Additionally, the results revealed that no fragment repetition events were identified in the tree peony genome, and only four genes were determined as tandem duplicates. However, a total of 78 segmental duplications and 25 tandem duplications were identified in Nicotiana tabacum [18], and more duplication events were also discovered in Glycine max [29], Populus trichocarpa [7], and maize [21], suggesting that the expansion of the POD gene family in tree peony differs from that of most plants, which might be related to the γ whole-genome duplication (WGD) event [24]. PoPODs were minimally influenced by segmental or tandem duplications.
In this study, eight distinct tissues were utilized to verify the tissue-specific expression patterns of 57 PoPODs. Here, some PoPOD genes were highly expressed in all tissues, while others exhibited low or no expression in certain tissues. In Arabidopsis Thaliana, some POD genes (PRX39/40/57) were highly expressed in the apical meristem and elongation zone boundaries and were implicated in root development [31]. Intriguingly, in this research, multiple PoPOD genes were observed to have high expression levels in the apical and adventitious roots, and a similar phenomenon was observed in maize [21] and rice [32]. Hence, it is hypothesized that these highly expressed PoPOD genes play a crucial role in root development. To further validate this hypothesis, RNA-seq data were used to analyze the expression level of the PoPOD genes during in vitro adventitious root development of tree peony. These findings indicated that cluster 8 encompassed the majority of PoPOD genes (PoPOD30/34/43/46/47/57), and this cluster was enriched in the biosynthesis of secondary metabolites and phenylpropanoid biosynthesis pathways. Additionally, PoPOD30/34/57 contained meristem expression response elements, and PoPOD46/57 contained auxin response elements. It is noteworthy that the PoPOD43 homologous gene AtPRX64 has been reported to be regulated by AtMYB36 to maintain ROS balance and to participate in lateral root primordia formation [33], and the PoPOD30 homologous gene, AtIRR1, might participate in adventitious root formation through ROS signals [34]. In conclusion, six PoPODs are implicated in the formation of ARs in tree peony.
An increasing number of studies have demonstrated that PODs are widely involved in various plant stress responses. Overexpression of the POD gene TaPRX-2A in wheat has been reported to activate antioxidant enzymes and the ABA pathway, reducing ROS accumulation and enhancing salt tolerance [15]. Several studies have shown that PODs are regulated by plant hormones such as SA, ABA, JA, and ET, and play central roles in plant defense mechanisms [1]. In this study, the blue module significantly associated with heat resistance was identified through WGCNA analysis, which contained two PoPOD genes (PoPOD5/15), and the regulatory network of their interaction with transcription factors was constructed. Notably, PoPOD5/15 were observed to contain five ABA and 11 MeJA response elements. These results suggested that PoPOD5/15 could be involved in the heat resistance of tree peony via ABA and MeJA signaling. Senescence represents the final stage of plant development, during which endogenous SA and ROS levels increase in senescent plant tissues, with ROS acting as an upstream regulator of the SA pathway. In chrysanthemum florets, POD, CAT, and other antioxidant enzymes initially increased and then decreased during the senescence process. This suggested that an imbalance between ROS production and antioxidant defense led to senescence-related damage in florets [35]. In Iris versicolor flowers, some antioxidant enzymes increased as the flowers opened and decreased significantly toward senescence [36]. In this study, the expression of PoPOD23 exhibited an upregulation followed by a downregulation trend, and its promoter contained multiple hormone-responsive elements, such as SA and ABA. These results implied that PoPOD23 may play a significant role in tree peony senescence.

4. Materials and Methods

4.1. Genome-Wide Identification of PoPOD Genes

The genome data of P. ostii and A. thaliana (TAIR 10) were obtained from the China National GeneBank DataBase (https://ftp.cngb.org/pub/CNSA/data5/CNP0003098/CNS0560369/CNA0050666/, accessed on 19 February 2023) [24] and the TAIR database (https://www.arabidopsis.org/, accessed on 12 August 2023) [37], respectively. The protein sequences of 57 POD proteins in Arabidopsis were downloaded from the TAIR database (https://www.arabidopsis.org/, accessed on 12 June 2024) [37]. These sequences were used as a query to perform the BLASTp search (E value < 10−5). On the other hand, the peroxidase domain (PF00141) was downloaded from the Pfam database (http://pfam.xfam.org/, accessed on 12 June 2024) [38], and we used the HMMER search to identify the possible PoPODs. Subsequently, the putative PoPODs were confirmed using CDD (https://www.ncbi.nlm.nih.gov/cdd, accessed on 14 June 2024) [39]. Finally, the physical and chemical properties and subcellular localization of the identified PoPOD proteins were calculated using the ExPASy-ProtParam tool (https://web.expasy.org/protparam/, accessed on 14 June 2024) [40] and the WoLF PSORT (https://wolfpsort.hgc.jp/, accessed on 14 June 2024), respectively.

4.2. Phylogenetic Analysis, Conserved Motif and Gene Structure

The Muscle method was used to construct the alignment of POD proteins from Arabidopsis and P. ostii. Then, the alignment results were subjected to IQTREE to perform a maximum likelihood (ML) phylogenetic tree with the bootstrap set to 5000 [41]. Using TBtools (v. 2.119) software, phylogenetic trees were created [42]. For conserved motifs, the amino acid sequences of PoPOD were submitted to the MEME online website (https://meme-suite.org/meme/tools/meme, accessed on 17 June 2024) [43] with a maximum of ten motifs, and images were generated by TBtools software [42]. The GFF3 file of the P. ostii was submitted to the TBtools software for the gene structure of PoPODs [42].

4.3. Chromosomal Localization, Gene Collinearity and CIS-Acting Element Analysis

PoPOD genes were mapped to the P. ostii chromosome using the TBtools software based on P. ostii genomic data [42]. For the chromosomal locations of PoPOD genes, we utilized the P. ostii genomic data and visualized the results using the TBtools software [42]. McscanX software was used to analyze the segment and tandem duplication events of PoPOD genes [44]. The intra-species synteny relationships of P. ostii genome PODs and the inter-species synteny relationships between Arabidopsis and P. ostii were identified by MCScanX software [44]. We used the TBtools to extract the PoPODs promoter sequences (upstream 2.0 kb) from the P. ostii genomic sequence data and then uploaded this to the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 12 July 2024) for cis-acting element analysis [45]. TBtools was used for visualization [42].

4.4. Expression Profiles of PoPOD and Data Analysis

Based on the RNA-seq data from eight different tissues of the tree peony, including apical, stems, leaves, petals, stamens, pistils, seeds, (https://ftp.cngb.org/pub/CNSA/data5/CNP0003098/, accessed on 13 June 2024) [24] and adventitious roots, the tissue-specific expression of PoPODs was examined. Additionally, the expression of PoPODs at zero, three, five, and ten days during AR formation was examined in tree peony (the RNA-seq data were acquired from the corresponding author). The expression trends of genes during AR formation were analyzed using the OmicShare tools (https://www.omicshare.com/tools/, accessed on 15 September 2024). To identify the responses of PoPODs to high-temperature stress (0 h, 2 h, 6 h, 12 h, and 24 h), the transcriptome raw data were downloaded from the NCBI SRA (accession numbers: PRJNA1079236) [27]. The WGCNA was used to analyze the modules which significantly correlated with heat resistance. The WGCNA-shinyApp (https://github.com/ShawnWx2019/WGCNA-shinyApp, accessed on 15 September 2024) was used for analysis in TBtools [42]. The PlantTFDB website (http://planttfdb.gao-lab.org/, accessed on 15 September 2024) was used to predict the transcription factors. The KEGG enrichment analysis of the identified clusters and the blue module was analyzed using KEGG Enrichment in TBtools [42] and images were demonstrated by the OmicShare tools (https://www.omicshare.com/tools/, accessed on 15 September 2024). Co-expression network diagrams of transcription factors and PoPOD5/15 in the blue module were constructed using the OmicShare tools (https://www.omicshare.com/tools/, accessed on 16 September 2024).

4.5. Plant Materials, RNA Extraction and qRT-PCR

The tree peony (Paeonia suffruticosa Andr.) cultivar ‘Xue Ying Tao Hua’ originated from Luoyang, Henan province, China. The cultivar was chosen because of its upright flowers and early flowering, making it an important cut flower for tree peony. In this study, flower petal samples at one, two, three, four, and five days after the bottle insertion were collected for qRT-PCR analysis. The total RNA from the petal samples was isolated using a Quick RNA extraction Kit (0416-50, Huayueyang, Beijing, China). Using a kit with gDNA Clean (AG11728, AG, Hunan, China), 1 μg of total RNA was converted to cDNA. The qRT-PCR was performed using a SYBR® Green Pro Taq HS Mix kit (AG11719, AG, Hunan, China). Then, the relative expression levels of PoPODs were calculated using the 2−ΔΔt method for relative quantification. The primers were designed using TBtools and were listed in the Supplementary Table S1.

5. Conclusions

In this study, a genome-wide characterization of PoPOD genes was conducted for the first time in tree peony, with a particular focus on their role in AR development, petal senescence, and heat tolerance. Fifty-seven PoPODs distributed across five chromosomes were systematically identified in tree peony. Conserved motifs and gene structure indicated structural conservation and diversification within and between PoPOD subgroups, supporting the potential functional diversity of PoPODs. The chromosomal location and gene collinearity analyses revealed that PoPODs were minimally influenced by segmental or tandem duplications. RNA-seq data from eight different tissues showed that multiple PoPOD genes were observed to have high expression levels in the apical and ARs. Expression profiling combined with cis-elements analysis revealed that PoPOD30/34/57 contained meristem expression response elements, and PoPOD46/57 contained auxin response elements. WGCNA analysis of the transcriptome data revealed that PoPOD5/15 were located in modules significantly associated with heat tolerance. The qRT-PCR analysis indicated that the expression of PoPOD23 exhibited an upregulation followed by a downregulation trend and its promoter contained multiple hormone-responsive elements, such as SA and ABA. These findings offer valuable insights into the function of PoPODs, especially in root development, heat tolerance, and petal senescence, and establish a foundation for future research to validate their precise regulatory mechanisms.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms252212122/s1.

Author Contributions

L.L., W.L., S.H. and Z.W. planned and designed the research. L.L., W.S. and W.L. conducted experiments and material collection. L.L. and W.L. analyzed data. L.L., P.Z., D.L. and Y.S. conducted fieldwork and material maintenance. L.L., W.L. and Z.W. wrote the manuscript. W.L., S.H. and Z.W. revised and finalized the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Beijing Laboratory of Urban and Rural Ecological Environment Open Project (Grant No. BFUKF202408), the China Postdoctoral Science Foundation (Grant No. GZC20240433), the Natural Science Foundation of China (Grant No. 32472788, 32271956), and the Central Plains Academics of Henan Province (Grant No. 244000510002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data contained within the article.

Acknowledgments

All individuals included in this section have consented to the acknowledgement.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cosio, C.; Dunand, C. Specific functions of individual class III peroxidase genes. J. Exp. Bot. 2009, 60, 391–408. [Google Scholar] [CrossRef] [PubMed]
  2. Hiraga, S.; Ichinose, C.; Onogi, T.; Niki, H.; Yamazoe, M. Bidirectional migration of SeqA-bound hemimethylated DNA clusters and pairing of oriC copies in Escherichia coli. Genes Cells 2000, 5, 327–341. [Google Scholar] [CrossRef] [PubMed]
  3. Erman, J.E.; Vitello, L.B. Yeast cytochrome c peroxidase: Mechanistic studies via protein engineering. Biochim. Biophys. Acta 2002, 1597, 193–220. [Google Scholar] [CrossRef] [PubMed]
  4. Piontek, K.; Smith, A.T.; Blodig, W. Lignin peroxidase structure and function. Biochem. Soc. Trans 2001, 29, 111–116. [Google Scholar] [CrossRef]
  5. Shigeto, J.; Tsutsumi, Y. Diverse functions and reactions of class III peroxidases. New Phytol 2016, 209, 1395–1402. [Google Scholar] [CrossRef]
  6. Herrero, J.; Fernandez-Perez, F.; Yebra, T.; Novo-Uzal, E.; Pomar, F.; Pedreno, M.A.; Cuello, J.; Guera, A.; Esteban-Carrasco, A.; Zapata, J.M. Bioinformatic and functional characterization of the basic peroxidase 72 from Arabidopsis thaliana involved in lignin biosynthesis. Planta 2013, 237, 1599–1612. [Google Scholar] [CrossRef]
  7. Ren, L.L.; Liu, Y.J.; Liu, H.J.; Qian, T.T.; Qi, L.W.; Wang, X.R.; Zeng, Q.Y. Subcellular Relocalization and Positive Selection Play Key Roles in the Retention of Duplicate Genes of Populus Class III Peroxidase Family. Plant Cell 2014, 26, 2404–2419. [Google Scholar] [CrossRef]
  8. Hoson, T. Regulation of polysaccharide breakdown during auxin-induced cell wall loosening. J. Plant Physiol. 1993, 106, 369–381. [Google Scholar] [CrossRef]
  9. Arrieta-Baez, D.; Stark, R.E. Modeling suberization with peroxidase-catalyzed polymerization of hydroxycinnamic acids: Cross-coupling and dimerization reactions. Phytochemistry 2006, 67, 743–753. [Google Scholar] [CrossRef]
  10. Jemmat, A.M.; Ranocha, P.; Le Ru, A.; Neel, M.; Jauneau, A.; Raggi, S.; Ferrari, S.; Burlat, V.; Dunand, C. Coordination of five class III peroxidase-encoding genes for early germination events of Arabidopsis thaliana. Plant Sci. 2020, 298, 110565. [Google Scholar] [CrossRef]
  11. Tung, C.W.; Dwyer, K.G.; Nasrallah, M.E.; Nasrallah, J.B. Genome-wide identification of genes expressed in Arabidopsis pistils specifically along the path of pollen tube growth. Plant Physiol. 2005, 138, 977–989. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, C.J.; Chan, Y.L.; Shien, C.H.; Yeh, K.W. Molecular characterization of fruit-specific class III peroxidase genes in tomato (Solanum lycopersicum). J. Plant Physiol. 2015, 177, 83–92. [Google Scholar] [CrossRef] [PubMed]
  13. Pan, J.; Song, J.; Sohail, H.; Sharif, R.; Yan, W.; Hu, Q.; Qi, X.; Yang, X.; Xu, X.; Chen, X. RNA-seq-based comparative transcriptome analysis reveals the role of CsPrx73 in waterlogging-triggered adventitious root formation in cucumber. Hortic. Res. 2024, 11, uhae062. [Google Scholar] [CrossRef] [PubMed]
  14. Passardi, F.; Tognolli, M.; De Meyer, M.; Penel, C.; Dunand, C. Two cell wall associated peroxidases from Arabidopsis influence root elongation. Planta 2006, 223, 965–974. [Google Scholar] [CrossRef]
  15. Su, P.; Yan, J.; Li, W.; Wang, L.; Zhao, J.; Ma, X.; Li, A.; Wang, H.; Kong, L. A member of wheat class III peroxidase gene family, TaPRX-2A, enhanced the tolerance of salt stress. BMC Plant Biol. 2020, 20, 392. [Google Scholar] [CrossRef]
  16. Zheng, C.; Wang, X.; Xu, Y.; Wang, S.; Jiang, X.; Liu, X.; Cui, W.; Wu, Y.; Yan, C.; Liu, H.; et al. The peroxidase gene OsPrx114 activated by OsWRKY50 enhances drought tolerance through ROS scavenging in rice. Plant Physiol. Bioch. 2023, 204, 10813824. [Google Scholar] [CrossRef]
  17. Wu, C.; Ding, X.; Ding, Z.; Tie, W.; Yan, Y.; Wang, Y.; Yang, H.; Hu, W. The Class III Peroxidase (POD) Gene Family in Cassava: Identification, Phylogeny, Duplication, and Expression. Int. J. Mol. Sci. 2019, 20, 2730. [Google Scholar] [CrossRef]
  18. Cheng, L.; Ma, L.; Meng, L.; Shang, H.; Cao, P.; Jin, J. Genome-Wide Identification and Analysis of the Class III Peroxidase Gene Family in Tobacco (Nicotiana tabacum). Front. Genet. 2022, 13, 916867. [Google Scholar] [CrossRef]
  19. Zhang, H.; Wang, Z.; Li, X.; Gao, X.; Dai, Z.; Cui, Y.; Zhi, Y.; Liu, Q.; Zhai, H.; Gao, S.; et al. The IbBBX24-IbTOE3-IbPRX17 module enhances abiotic stress tolerance by scavenging reactive oxygen species in sweet potato. New Phytol. 2022, 233, 1133–1152. [Google Scholar] [CrossRef]
  20. Tognolli, M.; Penel, C.; Greppin, H.; Simon, P. Analysis and expression of the class III peroxidase large gene family in Arabidopsis thaliana. Gene 2002, 288, 129–138. [Google Scholar] [CrossRef]
  21. Wang, Y.; Wang, Q.; Zhao, Y.; Han, G.; Zhu, S. Systematic analysis of maize class III peroxidase gene family reveals a conserved subfamily involved in abiotic stress response. Gene 2015, 566, 95–108. [Google Scholar] [CrossRef] [PubMed]
  22. Cao, Y.; Han, Y.; Meng, D.; Li, D.; Jin, Q.; Lin, Y.; Cai, Y. Structural, Evolutionary, and Functional Analysis of the Class III Peroxidase Gene Family in Chinese Pear (Pyrus bretschneideri). Front. Plant Sci. 2016, 7, 1874. [Google Scholar] [CrossRef] [PubMed]
  23. Gonzalez-Gordo, S.; Munoz-Vargas, M.A.; Palma, J.M.; Corpas, F.J. Class III Peroxidases (POD) in Pepper (Capsicum annuum L.): Genome-Wide Identification and Regulation during Nitric Oxide (NO)-Influenced Fruit Ripening. Antioxidants 2023, 12, 1013. [Google Scholar] [CrossRef] [PubMed]
  24. Yuan, J.; Jiang, S.; Jian, J.; Liu, M.; Yue, Z.; Xu, J.; Li, J.; Xu, C.; Lin, L.; Jing, Y.; et al. Genomic basis of the giga-chromosomes and giga-genome of tree peony Paeonia ostii. Nat. Commun. 2022, 13, 7328. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, J.; Song, Y.; Wang, G.; Shi, L.; Shen, Y.; Liu, W.; Xu, Y.; Lou, X.; Jia, W.; Zhang, M.; et al. PoARRO-1 regulates adventitious rooting through interaction with PoIAA27b in Paeonia ostii. Plant Sci. 2024, 347, 112204. [Google Scholar] [CrossRef]
  26. Yuan, Y.; Zhou, N.; Bai, S.; Zeng, F.; Liu, C.; Zhang, Y.; Gai, S.; Gai, W. Evolutionary and Integrative Analysis of the Gibberellin 20-oxidase, 3-oxidase, and 2-oxidase Gene Family in Paeonia ostii: Insight into Their Roles in Flower Senescence. Agronomy 2024, 14, 590. [Google Scholar] [CrossRef]
  27. Yang, X.; Huang, Y.; Yao, Y.; Bu, W.; Zhang, M.; Zheng, T.; Luo, X.; Wang, Z.; Lei, W.; Tian, J.; et al. Mining Heat-Resistant Key Genes of Peony Based on Weighted Gene Co-Expression Network Analysis. Genes 2024, 15, 383. [Google Scholar] [CrossRef]
  28. Xiao, H.; Wang, C.; Khan, N.; Chen, M.; Fu, W.; Guan, L.; Leng, X. Genome-wide identification of the class III POD gene family and their expression profiling in grapevine (Vitis vinifera L). BMC Genom. 2020, 21, 444. [Google Scholar] [CrossRef]
  29. Aleem, M.; Riaz, A.; Raza, Q.; Aleem, M.; Aslam, M.; Kong, K.; Atif, R.M.; Kashif, M.; Bhat, J.A.; Zhao, T. Genome-wide characterization and functional analysis of class III peroxidase gene family in soybean reveal regulatory roles of GsPOD40 in drought tolerance. Genomics 2022, 114, 45–60. [Google Scholar] [CrossRef]
  30. Shang, H.; Fang, L.; Qin, L.; Jiang, H.; Duan, Z.; Zhang, H.; Yang, Z.; Cheng, G.; Bao, Y.; Xu, J.; et al. Genome-wide identification of the class III peroxidase gene family of sugarcane and its expression profiles under stresses. Front. Plant Sci. 2023, 14, 1101665. [Google Scholar] [CrossRef]
  31. Eljebbawi, A.; Guerrero, Y.; Dunand, C.; Estevez, J.M. Highlighting reactive oxygen species as multitaskers in root development. iScience 2021, 24, 101978. [Google Scholar] [CrossRef] [PubMed]
  32. Passardi, F.; Longet, D.; Penel, C.; Dunand, C. The class III peroxidase multigenic family in rice and its evolution in land plants. Phytochemistry 2004, 65, 1879–1893. [Google Scholar] [CrossRef] [PubMed]
  33. Fernandez-Marcos, M.; Desvoyes, B.; Manzano, C.; Liberman, L.M.; Benfey, P.N.; Del Pozo, J.C.; Gutierrez, C. Control of Arabidopsis lateral root primordium boundaries by MYB36. New Phytol. 2017, 213, 105–112. [Google Scholar] [CrossRef]
  34. Yan, J.; Song, Y.; Li, M.; Hu, T.; Hsu, Y.F.; Zheng, M. IRR1 contributes to de novo root regeneration from Arabidopsis thaliana leaf explants. Physiol. Plant 2023, 175, e14047. [Google Scholar] [CrossRef]
  35. Chakrabarty, D.; Chatterjee, J.; Datta, S.K. Oxidative stress and antioxidant activity as the basis of senescence in chrysanthemum florets. Plant Growth Regul. 2007, 53, 107–115. [Google Scholar] [CrossRef]
  36. Ahmad, S.S.; Tahir, I. Increased oxidative stress, lipid peroxidation and protein degradation trigger senescence in Iris versicolor L. flowers. Physiol. Mol. Biol. Plants 2016, 22, 507–514. [Google Scholar] [CrossRef]
  37. Lamesch, P.; Berardini, T.Z.; Li, D.; Swarbreck, D.; Wilks, C.; Sasidharan, R.; Muller, R.; Dreher, K.; Alexander, D.L.; Garcia-Hernandez, M.; et al. The Arabidopsis Information Resource (TAIR): Improved gene annotation and new tools. Nucleic Acids Res. 2012, 40, D1202–D1210. [Google Scholar] [CrossRef]
  38. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam protein families database: Towards a more sustainable future. Nucleic Acids Res. 2016, 44, D279–D285. [Google Scholar] [CrossRef]
  39. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C.J.; Lu, S.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; et al. CDD/SPARCLE: Functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017, 45, D200–D203. [Google Scholar] [CrossRef]
  40. Wilkins, M.R.; Gasteiger, E.; Bairoch, A.; Sanchez, J.C.; Hochstrasser, D.F. Protein identification and analysis tools in the ExPASy server. Methods Mol. Biol. 1999, 112, 531–552. [Google Scholar]
  41. Nguyen, L.T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  42. Chen, C.; Wu, Y.; Li, J.; Wang, X.; Zeng, Z.; Xu, J.; Liu, Y.; Feng, J.; Chen, H.; He, Y.; et al. TBtools-II: A “one for all, all for one” bioinformatics platform for biological big-data mining. Mol. Plant 2023, 16, 1733–1742. [Google Scholar] [CrossRef] [PubMed]
  43. Bailey, T.L.; Williams, N.; Misleh, C.; Li, W.W. MEME: Discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34, W369–W373. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef]
  45. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2022, 30, 325–327. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic relationship of POD genes between Arabidopsis and P. ostii. The phylogenetic tree was constructed by TBtools (v. 2.119) software using the Maximum Likelihood Method (5000 bootstrap). The numbers at nodes indicated bootstrap values per 5000 replicates as determined by the Maximum Likelihood Method. These proteins were divided into six subgroups (I–VI)), and were represented by bright red, dark blue, strong yellow, strong pink, strong green, and deep purple, respectively.
Figure 1. Phylogenetic relationship of POD genes between Arabidopsis and P. ostii. The phylogenetic tree was constructed by TBtools (v. 2.119) software using the Maximum Likelihood Method (5000 bootstrap). The numbers at nodes indicated bootstrap values per 5000 replicates as determined by the Maximum Likelihood Method. These proteins were divided into six subgroups (I–VI)), and were represented by bright red, dark blue, strong yellow, strong pink, strong green, and deep purple, respectively.
Ijms 25 12122 g001
Figure 2. Phylogenetic evolutionary tree and motif composition of POD proteins in P. ostii. (A) The phylogenetic tree was constructed using full-length protein sequences of PoPOD proteins using the Maximum Likelihood Method (5000 bootstrap). (B) The amino acid motifs (numbered 1–10) in PoPOD proteins are displayed in ten colored boxes, and black lines indicate amino acid length.
Figure 2. Phylogenetic evolutionary tree and motif composition of POD proteins in P. ostii. (A) The phylogenetic tree was constructed using full-length protein sequences of PoPOD proteins using the Maximum Likelihood Method (5000 bootstrap). (B) The amino acid motifs (numbered 1–10) in PoPOD proteins are displayed in ten colored boxes, and black lines indicate amino acid length.
Ijms 25 12122 g002
Figure 3. A gene collinearity analysis between P. ostii and A. thaliana and the chromosomal location of PoPODs. (A) A collinearity analysis of PoPOD genes between P. ostii and A. thaliana was constructed, and the red lines represented seven homologous gene pairs. (B) The blue bars represent the chromosomes of P. ostii; Chr 01–05 represent each corresponding chromosome, and the 0–3000 Mb scale represents chromosome length. A total of 56 out of the 57 PoPOD genes were located on the five chromosomes of P. ostii in addition to PoPOD57.
Figure 3. A gene collinearity analysis between P. ostii and A. thaliana and the chromosomal location of PoPODs. (A) A collinearity analysis of PoPOD genes between P. ostii and A. thaliana was constructed, and the red lines represented seven homologous gene pairs. (B) The blue bars represent the chromosomes of P. ostii; Chr 01–05 represent each corresponding chromosome, and the 0–3000 Mb scale represents chromosome length. A total of 56 out of the 57 PoPOD genes were located on the five chromosomes of P. ostii in addition to PoPOD57.
Ijms 25 12122 g003
Figure 4. Cis-element analysis of PoPOD in P. ostii. (A) The proportion of cis-elements predicted in the promoters of PoPOD. (B) Numbers of the cis-elements involved in light response, phytohormone, stress response, transcription factor binding site, and plant tissue. (C) The distribution of the main cis-elements in PoPOD gene promoters.
Figure 4. Cis-element analysis of PoPOD in P. ostii. (A) The proportion of cis-elements predicted in the promoters of PoPOD. (B) Numbers of the cis-elements involved in light response, phytohormone, stress response, transcription factor binding site, and plant tissue. (C) The distribution of the main cis-elements in PoPOD gene promoters.
Ijms 25 12122 g004
Figure 5. Expression profiles of tissue-specific PoPOD genes in different tissues. Heatmap of the expression in eight different tissues, including apical, root, seed, stamen, petal, pistil, stem, and leaf, by columns.
Figure 5. Expression profiles of tissue-specific PoPOD genes in different tissues. Heatmap of the expression in eight different tissues, including apical, root, seed, stamen, petal, pistil, stem, and leaf, by columns.
Ijms 25 12122 g005
Figure 6. Transcriptome data analysis within the in vitro adventitious roots of the tree peony. (A) Expression profiles of PoPOD genes during in vitro adventitious roots of tree peony. (B) Trend analysis of the gene expression during in vitro adventitious roots of tree peony. The distinct colors in the background indicate different clusters. (C) KEGG enrichment analysis of cluster 8 genes. (D) Classification of transcription factors during in vitro adventitious roots in cluster 8. The distinct colors in the background indicate different transcription factors. (E) Heat map of transcription factor expression during in vitro adventitious roots in cluster 8.
Figure 6. Transcriptome data analysis within the in vitro adventitious roots of the tree peony. (A) Expression profiles of PoPOD genes during in vitro adventitious roots of tree peony. (B) Trend analysis of the gene expression during in vitro adventitious roots of tree peony. The distinct colors in the background indicate different clusters. (C) KEGG enrichment analysis of cluster 8 genes. (D) Classification of transcription factors during in vitro adventitious roots in cluster 8. The distinct colors in the background indicate different transcription factors. (E) Heat map of transcription factor expression during in vitro adventitious roots in cluster 8.
Ijms 25 12122 g006
Figure 7. Transcriptome data analysis of heat resistance in the tree peony. (A) Expression profiles of PoPOD genes during high temperature treatment. (B) Heat map of the correlation between modules and traits. (C) KEGG enrichment analysis of blue module genes. (D) Classification of transcription factors in blue module. The distinct colors in the background indicate different transcription factors. (E) Heat map of transcription factor expression during high temperature treatment in blue module. (F) Co-expression network diagram of transcription factors and PoPOD5/15 in blue module. The small blue circles indicate different transcription factors. The co-expression network was constructed using the OmicShare tools (https://www.omicshare.com/tools/, accessed on 16 September 2024).
Figure 7. Transcriptome data analysis of heat resistance in the tree peony. (A) Expression profiles of PoPOD genes during high temperature treatment. (B) Heat map of the correlation between modules and traits. (C) KEGG enrichment analysis of blue module genes. (D) Classification of transcription factors in blue module. The distinct colors in the background indicate different transcription factors. (E) Heat map of transcription factor expression during high temperature treatment in blue module. (F) Co-expression network diagram of transcription factors and PoPOD5/15 in blue module. The small blue circles indicate different transcription factors. The co-expression network was constructed using the OmicShare tools (https://www.omicshare.com/tools/, accessed on 16 September 2024).
Ijms 25 12122 g007
Figure 8. qRT-PCR analysis of five PoPoD genes in five different time periods. The mean ± SD from three biological replicates (n = 3) are shown. Different letters above the bars indicate a significant difference by Tukey’s test (p < 0.05, one-way ANOVA).
Figure 8. qRT-PCR analysis of five PoPoD genes in five different time periods. The mean ± SD from three biological replicates (n = 3) are shown. Different letters above the bars indicate a significant difference by Tukey’s test (p < 0.05, one-way ANOVA).
Ijms 25 12122 g008
Table 1. The physicochemical parameters of POD proteins in the P. ostii genome.
Table 1. The physicochemical parameters of POD proteins in the P. ostii genome.
Gene NameGene IDAAMWpIIIAIGAHSubcellular
Localization
PoPOD1Pos.gene8430433437,494.927.1242.5294.01−0.094Endoplasmic reticulum
PoPOD2Pos.gene1397732834,971.58.8738.6680.03−0.145Chloroplast
PoPOD3Pos.gene4053937840,274.496.1940.0873.57−0.131Nucleus
PoPOD4Pos.gene4855025027,425.215.7836.4576.56−0.367Cytoplasm
PoPOD5Pos.gene285631334,2726.3346.6681.34−0.156Chloroplast
PoPOD6Pos.gene4104224025,949.549.433383.75−0.15Chloroplast
PoPOD7Pos.gene2034018319,733.319.5120.1472.08−0.365Nucleus
PoPOD8Pos.gene3016031333,414.738.4333.7581.44−0.092Chloroplast
PoPOD9Pos.gene1541129931,983.028.535.0180.07−0.1Chloroplast
PoPOD10Pos.gene7629628330,991.789.0424.1587.240.016Extracellular
PoPOD11Pos.gene5872035038,414.899.3427.1391.940.039Chloroplast
PoPOD12Pos.gene3628935438,758.278.0847.0688.73−0.02Cytoplasm
PoPOD13Pos.gene5448336741,265.985.1642.0888.45−0.269Cytoplasm
PoPOD14Pos.gene5096831135,284.086.3740.7478.68−0.446Chloroplast/Nucleus
PoPOD15Pos.gene7527319821,8176.6536.7980.81−0.366Cytoplasm
PoPOD16Pos.gene5593833737,005.339.3938.7181.39−0.188Chloroplast
PoPOD17Pos.gene4137125527,465.187.0536.8681.88−0.16Cytoplasm
PoPOD18Pos.gene3900037640,087.664.938.9681.28−0.123Chloroplast
PoPOD19Pos.gene3900333535,654.326.234.5887.580.015Chloroplast
PoPOD20Pos.gene5997234536,541.64.5840.3179.83−0.1Vacuole
PoPOD21Pos.gene3467440245,305.819.3140.4274.55−0.398Mitochondrion/Vacuole
PoPOD22Pos.gene8364536640,172.876.6531.9783.93−0.13Chloroplast
PoPOD23Pos.gene4961243647,449.477.6951.5675.05−0.376Chloroplast
PoPOD24Pos.gene2843632134,356.948.3831.1984.240.034Chloroplast
PoPOD25Pos.gene5187035138,243.438.9940.6185.36−0.061Chloroplast
PoPOD26Pos.gene6604532635,222.765.3538.2585.890.051Vacuole
PoPOD27Pos.gene2739634437,363.274.6435.3691.280.041Extracellular
PoPOD28Pos.gene5419432435,276.158.954487.56−0.152Extracellular/Vacuole
PoPOD29Pos.gene1581533837,600.38.835.5486.54−0.224Extracellular/Vacuole
PoPOD30Pos.gene411634237,835.896.4429.3181.55−0.179Chloroplast
PoPOD31Pos.gene8127641245,147.998.7338.9791.380.085Chloroplast
PoPOD32Pos.gene7475035138,756.255.3237.5593.110.12Chloroplast
PoPOD33Pos.gene7482035338,242.516.834.7587.080.143Chloroplast
PoPOD34Pos.gene8349931734,306.229.0637.2986.25−0.096Chloroplast
PoPOD35Pos.gene4811731734,341.459.4934.5686.21−0.117Chloroplast
PoPOD36Pos.gene7159036338,969.516.9338.2987.41−0.021Chloroplast
PoPOD37Pos.gene4986031634,026.938.9133.6383.13−0.083Chloroplast
PoPOD38Pos.gene3170320922,213.275.8732.2487.850.063Extracellular
PoPOD39Pos.gene6548332635,923.895.9941.5976.35−0.145Chloroplast
PoPOD40Pos.gene3646733236,630.939.0743.6186.69−0.214Chloroplast/Vacuole
PoPOD41Pos.gene4999633136,771.25.7142.1183.9−0.183Extracellular
PoPOD42Pos.gene2319432435,145.024.936.6484.910.007Chloroplast
PoPOD43Pos.gene7935738842,578.377.1239.1881.24−0.256Endoplasmic reticulum
PoPOD44Pos.gene5467342247,085.995.0468.3769.38−0.39Chloroplast
PoPOD45Pos.gene4172935338,929.487.5739.8679.58−0.073Chloroplast
PoPOD46Pos.gene7676233236,370.448.9342.6587.32−0.173Extracellular
PoPOD47Pos.gene5250134136,930.059.442.7685.25−0.065Chloroplast
PoPOD48Pos.gene5090932736,281.486.4334.4985.29−0.103Extracellular
PoPOD49Pos.gene4582333836,567.498.0642.7681.92−0.124Extracellular
PoPOD50Pos.gene7647435138,792.669.143.3778.12−0.238Chloroplast
PoPOD51Pos.gene2611337440,966.098.9537.8892.62−0.009Chloroplast
PoPOD52Pos.gene2611032834,992.918.8436.5388.96−0.023Chloroplast
PoPOD53Pos.gene2837332834,992.918.8436.5388.96−0.023Chloroplast
PoPOD54Pos.gene2670332535,275.996.332.6684.62−0.066Vacuole
PoPOD55Pos.gene8200530534,282.116.840.9185.31−0.393Cytoplasm
PoPOD56Pos.gene8091927830,845.74.6750.2993.67−0.088Chloroplast
PoPOD57Pos.gene7239735138,650.28.3943.7589.2−0.093Vacuole
AA, number of amino acid; MW, molecular weight; pI, theoretical pI; II, instability index; AI, aliphatic index; and GAH, grand average of hydropathicity.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, L.; He, S.; Zhang, P.; Li, D.; Song, Y.; Shang, W.; Liu, W.; Wang, Z. Integration of Genome-Wide Identification and Transcriptome Analysis of Class III Peroxidases in Paeonia ostii: Insight into Their Roles in Adventitious Roots, Heat Tolerance, and Petal Senescence. Int. J. Mol. Sci. 2024, 25, 12122. https://doi.org/10.3390/ijms252212122

AMA Style

Li L, He S, Zhang P, Li D, Song Y, Shang W, Liu W, Wang Z. Integration of Genome-Wide Identification and Transcriptome Analysis of Class III Peroxidases in Paeonia ostii: Insight into Their Roles in Adventitious Roots, Heat Tolerance, and Petal Senescence. International Journal of Molecular Sciences. 2024; 25(22):12122. https://doi.org/10.3390/ijms252212122

Chicago/Turabian Style

Li, Li, Songlin He, Peidong Zhang, Dengpeng Li, Yinglong Song, Wenqian Shang, Weichao Liu, and Zheng Wang. 2024. "Integration of Genome-Wide Identification and Transcriptome Analysis of Class III Peroxidases in Paeonia ostii: Insight into Their Roles in Adventitious Roots, Heat Tolerance, and Petal Senescence" International Journal of Molecular Sciences 25, no. 22: 12122. https://doi.org/10.3390/ijms252212122

APA Style

Li, L., He, S., Zhang, P., Li, D., Song, Y., Shang, W., Liu, W., & Wang, Z. (2024). Integration of Genome-Wide Identification and Transcriptome Analysis of Class III Peroxidases in Paeonia ostii: Insight into Their Roles in Adventitious Roots, Heat Tolerance, and Petal Senescence. International Journal of Molecular Sciences, 25(22), 12122. https://doi.org/10.3390/ijms252212122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop