Genome-Wide Identification and Analysis of Class III Peroxidases in Allotetraploid Cotton (Gossypium hirsutum L.) and their Responses to PK Deficiency

Class III peroxidases (PODs), commonly known as secretable class III plant peroxidases, are plant-specific enzymes that play critical roles in not only plant growth and development but also the responses to biotic and abiotic stress. In this study, we identified 198 nonredundant POD genes, designated GhPODs, with 180 PODs being predicted to secrete into apoplast. These POD genes were divided into 10 sub-groups based on their phylogenetic relationships. We performed systematic bioinformatic analysis of the POD genes, including analysis of gene structures, phylogenetic relationships, and gene expression profiles. The GhPODs are unevenly distributed on both upland cotton sub-genome A and D chromosomes. Additionally, these genes have undergone 15 segmental and 12 tandem duplication events, indicating that both segmental and tandem duplication contributed to the expansion of the POD gene family in upland cotton. Ka/Ks analysis suggested that most duplicated GhPODs experienced negative selection, with limited functional divergence during the duplication events. High-throughput RNA-seq data indicated that most highly expressed genes might play significant roles in root, stem, leaf, and fiber development. Under K or P deficiency conditions, PODs showed different expression patterns in cotton root and leaf. This study provides useful information for further functional analysis of the POD gene family in upland cotton.