Genome-wide characterization and functional analysis of class III peroxidase gene family in soybean reveal regulatory roles of GsPOD40 in drought tolerance

Class III peroxidases (PODs) are plant-specific glycoproteins, that play essential roles in various plant physiological processes and defence responses. To date, scarce information is available about the POD gene family in soybean. Hence, the present study is the first comprehensive report about the genome-wide characterization of GmPOD gene family in soybean (Glycine max L.). Here, we identified a total of 124 GmPOD genes in soybean, that are unevenly distributed across the genome. Phylogenetic analysis classified them into six distinct sub-groups (A-F), with one soybean specific subgroup. Exon-intron and motif analysis suggested the existence of structural and functional diversity among the sub-groups. Duplication analysis identified 58 paralogous gene pairs; segmental duplication and positive/Darwinian selection were observed as the major factors involved in the evolution of GmPODs. Furthermore, RNA-seq analysis revealed that 23 out of a total 124 GmPODs showed differential expression between drought-tolerant and drought-sensitive genotypes under stress conditions; however, two of them (GmPOD40 and GmPOD42) revealed the maximum deregulation in all contrasting genotypes. Overexpression (OE) lines of GsPOD40 showed considerably higher drought tolerance compared to wild type (WT) plants under stress treatment. Moreover, the OE lines showed enhanced photosynthesis and enzymatic antioxidant activities under drought stress, resulting in alleviation of ROS induced oxidative damage. Hence, the GsPOD40 enhanced drought tolerance in soybean by regulating the key physiological and biochemical pathways involved in the defence response. Lastly, the results of our study will greatly assist in further functional characterization of GsPODs in plant growth and stress tolerance in soybean.