Genome-Wide Characterization of Ascorbate Peroxidase Gene Family in Peanut (Arachis hypogea L.) Revealed Their Crucial Role in Growth and Multiple Stress Tolerance

Ascorbate peroxidase (APX), an important antioxidant enzyme, plays a significant role in ROS scavenging by catalyzing the decrease of hydrogen peroxide under various environmental stresses. Nevertheless, information about the APX gene family and their evolutionary and functional attributes in peanut (Arachis hypogea L.) was not reported. Therefore, a comprehensive genome-wide study was performed to discover the APX genes in cultivated peanut genome. This study identified 166 AhAPX genes in the peanut genome, classified into 11 main groups. The gene duplication analysis showed that AhAPX genes had experienced segmental duplications and purifying selection pressure. Gene structure and motif investigation indicated that most of the AhAPX genes exhibited a comparatively well-preserved exon-intron pattern and motif configuration contained by the identical group. We discovered five phytohormones-, six abiotic stress-, and five growth and development-related cis-elements in the promoter regions of AhAPX. Fourteen putative ah-miRNAs from 12 families were identified, targeting 33 AhAPX genes. Furthermore, we identified 3,257 transcription factors from 38 families (including AP2, ARF, B3, bHLH, bZIP, ERF, MYB, NAC, WRKY, etc.) in 162 AhAPX genes. Gene ontology and KEGG enrichment analysis confirm the role of AhAPX genes in oxidoreductase activity, catalytic activity, cell junction, cellular response to stimulus and detoxification, biosynthesis of metabolites, and phenylpropanoid metabolism. Based on transcriptome datasets, some genes such as AhAPX4/7/17/77/82/86/130/133 and AhAPX160 showed significantly higher expression in diverse tissues/organs, i.e., flower, leaf, stem, roots, peg, testa, and cotyledon. Likewise, only a few genes, including AhAPX4/17/19/55/59/82/101/102/137 and AhAPX140, were significantly upregulated under abiotic (drought and cold), and phytohormones (ethylene, abscisic acid, paclobutrazol, brassinolide, and salicylic acid) treatments. qRT-PCR-based expression profiling presented the parallel expression trends as generated from transcriptome datasets. Our discoveries gave new visions into the evolution of APX genes and provided a base for further functional examinations of the AhAPX genes in peanut breeding programs.

Cultivated peanut/groundnut (A. hypogaea L.), an allotetraploid crop, is one of the most valuable and economic oilseed food crops globally (Agarwal et al., 2018;Bertioli et al., 2019;Chen X. et al., 2019;Zhuang et al., 2019). This crop is being widely cultivated in the tropical and subtropical regions globally; however, several abiotic and biotic factors significantly affect its growth and production, including many important agronomic traits (Agarwal et al., 2018;Gangurde et al., 2020Gangurde et al., , 2021Kumar et al., 2020;Pandey et al., 2020;Shasidhar et al., 2020;Sinha et al., 2020;Jadhav et al., 2021;Soni et al., 2021;Aravind et al., 2022;Bomireddy et al., 2022;Liu et al., 2022;Patel et al., 2022). Therefore, it is vital to identify new potential genes associated with multiple stress tolerance and trait improvement in peanut for better protein-rich food supply, particularly in Asian and African countries. In this regard, the recently sequenced peanut genome and recent advances in genomics-assisted breeding make it easier for us to carry out a comprehensive systematic analysis of new gene families (Varshney et al., 2019(Varshney et al., , 2020(Varshney et al., , 2021a. To our best knowledge, APX gene family was yet to be comprehensively characterized in peanut. Thus, the current study performed a genome-wide identification and characterized the APX gene family in peanut (AhAPX). Several in silico analysis, such as characterization, genomic evolution, gene structure, conserved motifs, cis-regulatory elements, putative miRNA and transcription factors, functional annotations, etc., were utilized to get insights into the novel roles of AhAPX genes. Furthermore, their expression profiling in diverse tissues/organs, under phytohormones and abiotic stress conditions were also performed using transcriptome and qRT-PCR techniques. In short, this report offered evolutionary and functional roles of AhAPX genes which could open new windows for further functional studies on the novel roles of AhAPX genes in peanut breeding programs under stress conditions.

Discovery and Physicochemical Features of APX Genes
As explained earlier Raza et al., 2021b;Su et al., 2021), two approaches, i.e., BLASTP and the Hidden Markov Model (HMM), were applied to identify APX genes in the peanut (A. hypogea) genome. The peanut genome sequence was taken from peanut Genome Resource (PGR) database 1 (Zhuang et al., 2019). In the first approach, the sequences of eight Arabidopsis thaliana APX genes were gained from TAIR Arabidopsis genome database 2 (Rhee et al., 2003). Then, these sequences were utilized as a query to perform the BLASTP against peanut genome. In the second approach, HMMER 3.1 3 (Finn et al., 2015) software was employed to seek out the APX genes with default controls. Later, the HMM file of the ascorbic acid peroxidase domain (PF00141) was retrieved from the Pfam database 4 (El-Gebali et al., 2019). Lastly, the sequences comprising the PF00141 domain were chosen as putative APX genes, and finally, 166 AhAPX genes were discovered by uniting the results obtained from both approaches in the peanut genome. Following the same approaches, APX genes were also discovered in diploid parents, i.e., A. duranensis (90 genes; AdAPX1-AdAPX90) and A. ipaensis (102 genes; AiAPX1-AiAPX102). Their genome sequences were downloaded from PeanutBase database 5 (Dash et al., 2016). The detailed information (including gene name, gene ID, and protein sequences) of all identified APX genes is given in Supplementary Table 1.

Evaluation of Chromosomal Location, Phylogenetic Relationships, and Synteny Analysis of APX Genes
The data about the chromosomal location of AhAPX was attained from the PGR database, and the TBTools was utilized to map the genes on chromosomes. To discover the evolutionary link of the APX proteins, a phylogenetic tree among A. hypogea (AhAPXs), A. duranensis (AdAPXs), A. ipaensis (AiAPXs), and A. thaliana (AtAPXs) was created. Multiple sequence alignment was implemented using MEGA7 software 10 (Kumar et al., 2018).
The neighbor-joining (NJ) method was undertaken to design a phylogenetic tree with 1,000 bootstrap replicates and iTOL was used to beautify the tree 11 (Letunic and Bork, 2021).
The syntenic associations of APX genes between A. hypogea, A. duranensis, A. ipaensis, and A. thaliana were executed through the MCScanX toolkit and were pictured by the Advance Circos package in the TBTools software (Chen et al., 2020). Additionally, the multiple collinearity analysis of APX genes was completed via multiple synteny Plot packages in TBTools software. The Ka/Ks ratios of all AhAPX were predicted via simple Ka/Ks calculator in TBTools software.

Prediction of cis-Regulatory Elements in the AhAPX Promoters
To predict the putative cis-regulatory elements in the AhAPX promoters, the 2 Kb sequences upstream of start codons were separated from the peanut genome. The promoter sequences of all AhAPX genes were observed with PlantCARE website 12 (Lescot et al., 2002), and the picture was illustrated using TBtools software.

Prediction of Putative miRNAs Targeting AhAPX Genes and Functional Annotation Evaluation
The CDS of all AhAPX was used to predict the miRNA target sites with psRNATarget website 13 (Dai et al., 2018) with default considerations. The interactive network figure among the putative miRNAs and AhAPX genes was made via Cytoscape software (v3.9) 14 (Shannon et al., 2003). Gene ontology (GO) and Kyoto encyclopedia of genes and genomics (KEGG) annotation evaluation was undertaken by submitting all AhAPX protein sequences to the eggNOG v4.0 15 (Powell et al., 2014). At the same time, GO and KEGG enrichment evaluations were performed with TBtools software.

Prediction of Transcription Factor Regulatory Network of AhAPX Genes
To predict the putative transcription factors (TFs) and regulatory network, the 500 bp nucleotide sequences from upstream regions of AhAPX genes were removed and complied to the PlantRegMap (Transcriptional Regulatory Map) 16 with p-value ≤ 1e −6 . The regulatory network of predicted TFs and AhAPX genes was created with Cytoscape v3.9 software.

Expression Profiling of AhAPX Genes
The expression levels of all AhAPX genes at diverse developmental tissues/organs (embryo, cotyledon, testa, pericarp, peg, root and stem, root nodule, root tip, root, step tip, stem, leaf, and flower), under various hormones (ethylene, abscisic acid, paclobutrazol, brassinolide, and salicylic acid), and abiotic stress (drought and cold) conditions were evaluated using openly available transcriptome dataset of cultivated peanut (cultivar Shitouqi) at PGR database (see text footnote 1; BioProject PRJNA480120) (Zhuang et al., 2019). The detailed procedure for sample harvesting and data analysis is presented in our recent paper (Zhuang et al., 2019). Owing to the great differences in the expression trends, we normalize the log2 of fragments per kilobase of transcript per million (FPKM) values. Finally, the circular heat maps were designed by TBtools software.

Plant Material and Stress Conditions
In this study, a widely cultivated peanut variety in southeast China, "Minhua-6" was used for stress treatments. The same variety was also used for transcriptome analysis in our recent paper (Zhuang et al., 2019). The seeds of the "Minhua-6" cultivar were obtained from the FAFU, Fuzhou, China. The vigor seeds were cultured on small pots having a mix of vermicompost. For stress treatment, germinated seedlings at the four-leaf stage were exposed to cold stress at 4 • C and ABA (10 µg mL −1 ) for 0 (CK), 3, 6, 9, and 12 h with three biological repetitions. All of the samples were instantly frozen in liquid nitrogen and were kept at -80 • C until RNA extraction.

RNA Extraction and qRT-PCR-based Expression Analysis
Total RNA was isolated utilizing the CTAB method as described in our recent work (Sharif et al., 2021), and cDNA was prepared with the help of Evo M-MLV RT Kit with gDNA Clean for qPCR II (Code No. AG11711; Hunan Aikerui Biological Engineering Co., Ltd., China) following the developer guidelines. The comprehensive information on qRT-PCR reaction has been described in our recent work (Sharif et al., 2021). The peanut Actin gene was used as a housekeeping gene to stabilize the expression (Chi et al., 2012). The expression data of three biological repeats were normalized using the 2 − CT method (Livak and Schmittgen, 2001). All the primers used for qRT-PCR are given in Supplementary Table 2. The graphs were made with GraphPad Prism v9.0.0 software 17 (Swift, 1997).
On the other hand, eight genes (AtAPXs) from A. thaliana, 90 genes (AdAPX1-AdAPX90) from A. duranensis, and 102 genes (AiAPX1-AiAPX102) from A. ipaensis genomes were also recognized to study the evolution of APX genes between tetraploid and diploid parents (Supplementary Table 1).

Insights From Phylogenetic Relationships of APX Proteins
To determine the in-depth evolutionary and phylogenetic history between the AhAPX (166 members), AdAPX (90 members), AiAPX (102 members), and AtAPX proteins (8 members), an unrooted phylogenetic tree was built by a multiple sequence alignment, which was divided into 11 main groups (group1-group11) (Figure 2). The discoveries exposed that group1 comprised of seven APX members (2 AhAPX, 2 AiAPX, and 3 AdAPX) followed by group4/5 (eight APX members), and group3 (14 APX members). Notably, the maximum number of APX members (66 AhAPX, 39 AiAPX, and 33 AdAPX) were found in group 11 followed by group7 (38 APX members), group2/6 (37 APX members), group8/9 (28 APX members), and group10 (24 APX members) (Figure 2). All AtAPX members were clustered only in one group, i.e., group 2. In general, APXs grouped into the indistinguishable sub-group may retain corresponding functions. It is worth stating that A. hypogea APX (AhAPXs) were distributed in each group with homologs from A. duranensis, A. ipaensis, and A. thaliana., and group11 was detected to have more AhAPX members than the other 10 groups (Figure 2). Furthermore, it was observed that the AhAPXs      In the genomic position, the positive (+) and negative (−) sign shows the presence of a gene on the positive and negative strand of that specific marker correspondingly. MW, molecular weight; PI, isoelectric points; bp, base pair; aa, amino acids.
showed a greater phylogenetic network with the AdAPXs and AiAPXs in each group.

Insights Into Synteny and Collinearity of APX Genes
Gene duplications (i.e., tandem and segmental) are thought to be the main factors in supporting the expansion and evolution of new gene families in plants (Cannon et al., 2004). Hence, gene duplication procedures were assessed between AhAPXs, AdAPXs, AiAPXs, and AtAPXs (Supplementary Table 3). The results of gene duplication study showed that there were 92 AhAPX gene pairs, and these pairs were unevenly mapped on different chromosomes (Figure 3 and Supplementary Table 3). Mainly, chromosome 13 had a maximum number (i.e., 16) of AhAPX gene pairs, followed by chromosome 5 with 12 AhAPX gene pairs. The least number of gene pairs (i.e., two) was discovered on chromosome 12, and no gene pair was found on chromosome 2 (Figure 3 and Supplementary Table 3). The results reveal that segmental duplications have contributed to the expansion of AhAPX genes in the cultivated peanut genome (Supplementary Table 3). Notably, no tandem duplicated gene pairs were identified. Similarly, 10 duplicated gene pairs were detected between AhAPX and AtAPX (Supplementary  Table 3). All these gene pairs were irregularly mapped on different chromosomes. Taken together, these conclusions explained that the duplication activities played a vital role in enlarging the APX genes between diploid and tetraploid parents. Further, it can also be concluded that A. hypogea might have lost some genes during genome evolution.
Collinearity analysis was carried out to review the evolutionary association of the APX genes between A. hypogea, A. duranensis, A. ipaensis, and A. thaliana (Figure 4 and Supplementary Table 3). The results discovered a strong orthologous of APX genes among these four species (Figure 4). On the whole, several A. hypogea genes presented syntenic networks with different AdAPX, AiAPX, and AtAPX genes. Particularly, only one gene (AhAPX14) at chromosome Ah2     A iA P X 1 0 A h A P X 7 6 A d A P X 1 0 A d A P X 7 8 A h A P X 1 5 0 A iA P X 5 7 A h A P X 8 7 A iA P X 9 0 A h A P X 7 1 A h A P X 1 6 0 A iA P X 1 9 A h A P X 7 7 A d A P X 7 A h A P X 9 7 A d A P X 1 7 A h A P X 1 5 A iA P X 1 4 A h A P X 11 8 A iA P X 9 5 A h A P X 5 A d A P X Ah AP X3 7 Ai AP X4 3 A hA P X 86 A iA P X 37 A hA P X 7 A dA P X 37 A dA P X 42 A iA P X 3 9 A h A P X 1 3 3 A iA P X 3 6 A h A P X 4 2 A d A P X 6 9 A h A P X 1 3 2 A d A P X 3 4 A h A P X 4 1 A iA P X 3 0 A d A P X 5 6 A iA P X 5 0 A h A P X 6 5 A h A P X 1 5 1 A h A P X 4 3 A h A P X 1 3 4 A iA P X 9 8 A d A P X 5 1 A iA P X 6 4 A h A P X 1 2 5 A iA P X 5 5 A h A P X 9 9 A iA P X 3 5 A h A P X 7 3 A d A P X 3 A h A P X 9 8 A iA P X 7 0 A h A P X 7 4 A h A P X 5 0 A h A P X 1 3 1 A d A P X 8 4 A h A P X 1 4 6 A iA P X 2 7 A d A P X 3 0 A h A P X 6 1 A h A P X 1 3 0 A iA P X 3 4 A h A P X 4 9 A d A P X 4 0 A d A P X 6 1 A iA P X 6 0 A d A P X 6 5 A iA P X Ai AP X7 9 A hA P X 17 A dA P X 80 A tA P X 2 A tA P X 1 A hA P X 60 A d A P X 8 1 A h A P X 2 A iA P X 9 1 A d A P X 8 9 A d A P X 7 5 A iA P X 7 5 A h A P X 1 0 2 A h A P X 1 9 A tA P X 3 A tA P X 5 A d A P X 7 7 A iA P X 9 7 A h A P X 3 4 A h A P X 5 9 A h A P X 1 3 7 A d A P X 7 9 A iA P X 8 7 A tA P X S A tA P X T A h A P X 1 4 0 A iA P X 9 2 A h A P X 5 5 A d A P X 8 6 A tA P X 4 A h A P X 6 4 A d A P X 8 2 A h A P X 1 A iA P X 8 3 A tA P X 6 A d A P X 9 0 A h A P X 1 1 4 A iA P X 6 1 A h A P X 2 3 A d A P X 6 3 A d A P X 6 2 A d A P X 8 7 A iA P X 2 exhibited a syntenic connection with AdAPX85 gene at chromosome Ad02 (Figure 4 and Supplementary Table 3), while other homologous genes present on other A. hypogea chromosomes also showed a syntenic relationship with many AdAPX, AiAPX and AtAPX genes (Figure 4 and Supplementary Table 3). These findings indicate that whole-genome or segmental duplication procedures are considered a main evolutionary force in the evolution of AhAPX genes in the peanut genome (Figure 4 and Supplementary Table 3).
The Ka/Ks ratio is considered as a huge diagnostic marker in evaluating the sequence evolution in terms of selection pressures and duplication types (Hurst, 2002). Thus, to understand the evolutionary story of the AhAPX, the Ka, Ks, and Ka/Ks ratio was revealed (Supplementary Table 3 Ah AP X4 9 Ah AP X4 9 A hA P X 49 A hA P X 51 A h A P X 5 1 A h A P X 5 2 A h A P X 5 4 A h A P X 5 5 A h A P X 5 6 Ah AP X5 7 A h A P X 5 3 A h A P X 5 2 A h A P X 5 9 A h A P X 6 0 A h A P X 6 0 A h A P X 6 1 A h A P X 6 2 A h A P X 6 2 A h A P X 6 3 A h A P X 6 4 A h A P X 6 2 A h A P X 6 0 A h A P X 6 5 A h A P X 6 6 A h A P X 6 7 A h A P X 6 7 A h A P X 6 8 A h A P X 6 9 A h A P X 6 7 FIGURE 3 | Chromosomal positions and inter-chromosomal groups of duplicated AhAPX gene pairs. Gray lines in the background demonstrate all syntenic blocks in the A. hypogea genome, and the blue lines exhibit the segmental or tandem duplication network zones among AhAPX genes. The near location of AhAPX genes is marked with black lines outside with chromosome names. Table 3), demonstrating that the AhAPX genes may have experienced strong purifying selective pressure and segmental duplications throughout the evolution procedure (Supplementary Table 3).

Insights Into Gene Structures and Conserved Motifs of AhAPX Genes
The exon-intron arrangements and conserved motifs of the AhAPX genes were analyzed to get insights into the advancement of the APX family genes in peanut genome (Figure 5 and Supplementary Table 4). The outcomes revealed that the number of exons and introns varied from 16 to 1 and 0 to 15, respectively ( Figure 5B and Supplementary Table 10). In short, 6 genes have 1 exon and zero intron; 13 genes have 2 exons and 1 intron; 5 genes have 3 exons and 2 introns; 18 genes have 2 exons and 1 intron; 41 genes have 3 exons and 2 introns; 81 genes have 4 exons and 3 introns; 3 genes have 5/7 exons and 4/6 introns; 5 genes have 9 exons and 8 introns; only 1 gene has 10 exons and 9 introns; 3 genes have 11/12 exons and 10/11 introns; and only 2 genes have a maximum number of exons (16) and introns (15) ( Figure 5B and Supplementary Table 10). Above all, genes belonging to the same sub-tree almost had parallel structures apart from a few genes ( Figure 5B). Among all genes, AhAPX149 possess the longest structure, and only a few genes have a complex structure, such as AhAPX17, AhAPX19, AhAPX34, AhAPX45, AhAPX55, AhAPX59, AhAPX64, AhAPX92, AhAPX101, AhAPX102, and AhAPX145 ( Figure 5B). Exon loss or gain has been found during the evolution of APX family genes. The results recommended that APX genes held a somewhat frequent exon-intron composition throughout the evolution of peanut genome. Furthermore, AhAPX gene participants inside a sub-tree had exceptionally corresponding gene structures, steady with their phylogenetic clusters. The conserved motif of the AhAPX genes ranged from one (AhAPX55/140) to three (AhAPX2/71/117/131) ( Figure 5A). In total, 10 conserved motifs were recognized, and their complete dataset, including motif names, sequences, width, and E-value, is given in Supplementary Table 4. Similar to gene structure, the motif distributions were also similar within the sub-trees (Figure 5A), while some motifs were found to be specific to some genes. For instance, some genes such as AhAPX2/17/101/34/59/137 were limited to motifs 2, 3, and 10. While AhAPX60 gene only contained motifs 3 and 10; AhAPX71 contained motifs 3, 8, and 9; AhAPX117 contained motifs 4, 6, and 10; AhAPX131 contained motifs 4, 9, and 10; and AhAPX71 contained motifs 3, 8, and 9 ( Figure 5A). Almost all other motifs were present on all genes except in a few cases ( Figure 5A). In summary, the consistency of gene organizations within sub-trees was credibly constant by appraising the conserved motif structures, gene structures, and phylogenetic relations, representing that the APX proteins have enormously wellsustained amino acid deposits and APX members belonging to the same tree may hold corresponding roles.

Cis-Elements: Key Players in the Promoter Regions of AhAPX Genes
To better understand the regulatory role of AhAPX genes toward peanut growth and development, and tolerance to abiotic stress and phytohormones treatment, cis-regulatory elements in the promoter of AhAPX were explored. The complete dataset of cis-elements is presented in Supplementary Table 5. We emphasized and recognized three categories of cis-elements, including abiotic stress-responsive, phytohormones responsive, and growth and development responsive elements (Figures 6, 7 and Supplementary Table 5). Mainly, six abiotic stress-responsive (drought, light, low temperature, wound, defense and stress, and anaerobic) elements were detected. These elements consist of I-box, ATCT-motif, Box 4, GT1-motif, GA-motif, etc. (lightresponsive, 77%), ARE (13%), MBS (3%), TC-rich repeats (3%), LTR (3%), and WUN-motif (0.15%) (Figures 7A,B and Supplementary Table 5). Overall, results showed that most of the abiotic stress-related elements were predicted to be specific to some genes and unevenly distributed (Figure 6 and Supplementary Table 5), indicating their defensive role against stress conditions. Likewise, five phytohormone-responsive elements [methyl jasmonate (MeJA), abscisic acid, gibberellin, salicylic acid, and auxin] consist of CGTCA-motif/TGACG-motif (36%), ABRE (35%), P-box/TATC-box/GARE-motif (11%), TCAelement/SARE (11%), and AuxRR-core/TGA-element/TGA-box (7%) (Figures 7C,D and Supplementary Table 5). Some of the elements were found to be specific to some genes and unevenly distributed (Figure 6 and Supplementary Table 5). These outcomes suggest that element-specific genes could be considered as candidate players for further functional studies to reveal their protective role under hormone treatments.
Moreover, five growth and development-related (zein metabolism, meristem expression, endosperm expression, circadian control, and cell cycle regulation) elements were discovered. These key elements include O 2 -site (34%), CAT-box (31%), GCN4_motif/AACA_motif (17%), circadian (14%), and MSA-like (5%) (Figures 7E,F and Supplementary Table 5), suggesting their dynamic role in different growth and developmental stages of peanut. In a nutshell, these discoveries suggested that some of the key elements are widely and randomly distributed in some genes, while some of the elements are found to be specific to some genes. It can be concluded that the Frontiers in Plant Science | www.frontiersin.org 13 September 2022 | Volume 13 | Article 962182

Gene Ontology and Kyoto Encyclopedia of Genes and Genomics Enrichment Analysis of AhAPX Genes
To advance our knowledge about the dynamic roles of AhAPX genes at molecular level, GO and KEGG enrichment analysis were performed (Figure 10 and Supplementary Table 8).
Furthermore, KEGG pathway enrichment study discovered six pathways participating in diverse functions of AhAPX genes ( Figure 10B and Supplementary Table 8). The highly enriched pathways include biosynthesis of other secondary metabolites (B09110), metabolism (A09100), phenylpropanoid biosynthesis (00940), followed by ascorbate and aldarate metabolism (00053), glutathione metabolism (00480), and metabolism of other amino acids (B09106) (Figure 10B and Supplementary Table 8). Briefly, it can be concluded that GO and KEGG enrichment study validates the functional contribution of AhAPX genes in several cellular, molecular, and biological processes, that are associated with antioxidant defense systems, ROS scavenging, response to stresses, and biosynthesis of different metabolites.

Expression Profiling of AhAPX Genes at Diverse Developmental Tissues
The expression profiling of 166 AhAPX genes was observed in various tissues and organs, including embryo, cotyledon, testa, pericarp, peg, root and stem, root nodule, root tip, root, step tip, stem, leaf, and flower using openly available transcriptome dataset (Supplementary Table 9). Overall, the expression heatmap indicated that only a few genes were highly expressed in certain organs/tissues (Figure 11 and Supplementary Table 9). For example, some genes including AhAPX4, AhAPX7, AhAPX17, AhAPX19, AhAPX28, AhAPX42, AhAPX51, AhAPX76, AhAPX77, AhAPX82, AhAPX86, AhAPX101, AhAPX102, AhAPX130, AhAPX133, and AhAPX160 were highly expressed in almost all the organs/tissues (Figure 11). While some genes were found to be specific to some tissues like AhAPX12 showed considerable expression in cotyledon, root and stem, root tip, and stem; AhAPX109, AhAPX111, and AhAPX13 expressed in stem, roots and peg; AhAPX135 expressed in pericarp; and AhAPX138 expressed in cotyledon (Figure 11). Particularly, a few genes also exhibited modest expressions in a variety of tissues. On the whole, expression dataset shows that some particular genes may substantially participate in peanut growth and development. Hence, the functional characterization of these genes may perhaps be carried out in future studies.

Expression Profiling of AhAPX Genes Under Abiotic Stress and Hormones Treatments
To further study the contribution of all AhAPX genes toward abiotic and hormones stress tolerance in peanut, an openly available transcriptome dataset was used to evaluate the expression levels (Figure 12 and Supplementary Table 9). Similar to tissue-specific trend, only a few genes showed higher expressions in both cold and drought stresses. For instance, AhAPX4, AhAPX17, AhAPX19, AhAPX82, AhAPX101, and AhAPX102 were highly expressed under stress (cold and drought) and CK conditions. Likewise, some genes also showed moderate expression levels, such as AhAPX27, AhAPX34, AhAPX51, AhAPX55, AhAPX59, AhAPX113, AhAPX137, AhAPX138, AhAPX140, and AhAPX157 under stress (cold and drought) and normal conditions. On the other hand, AhAPX720, AhAPX21, AhAPX51, AhAPX77, AhAPX106, AhAPX130, AhAPX158, and AhAPX160 displayed considerable expression under cold stress compared to CK conditions ( Figure 12A).
Under phytohormones treatments, AhAPX4, AhAPX17, AhAPX55, AhAPX59, AhAPX82, AhAPX101, AhAPX102, AhAPX137, and AhAPX140 displayed significantly higher expression patterns throughout the treatments. In comparison   organelle membrane cellular response to stimulus cellular oxidant detoxification cellular detoxification symbiotic interaction external encapsulating structure organization response to stimulus response to chemical hydrogen peroxide catabolic process response to zinc ion modulation by symbiont of host defense response modulation by symbiont of host immune response obsolete oxidation-reduction process response to host immune response negative regulation of growth response to host cell wall organization reactive oxygen species metabolic process plant-type cell wall organization detoxification modulation by symbiont of host innate immune response response to defenses of other organism cellular response to toxic substance modulation by symbiont of host process response to host defenses plant-type cell wall organization or biogenesis biological process involved in interaction with host cellular component organization or biogenesis induction by symbiont of host innate immune response 0.00 3.00 6.00 9.00 12.00 15.00 18.00 to CK, some genes are specifically expressed under certain hormones, such as AhAPX21 under paclobutrazol, AhAPX27 under ethylene and abscisic acid, AhAPX51 under abscisic acid, and AhAPX88 under abscisic acid and paclobutrazol ( Figure 12B). Notably, most of the genes did not show any expression under any type of stress conditions. The candidate genes with higher expression could be genetically engineered to improve the tolerance against multiple hormones and abiotic stress (cold and drought) conditions.

qRT-PCR-Based Expression Profiling of AhAPX Genes Under Cold and ABA Treatment
For qRT-PCR-based expression profiling, 10 highly upregulated AhAPX genes were selected based on transcriptome datasets to validate their transcript levels under ABA and cold treatment at various time points (Figure 13). Under ABA treatment, almost all genes demonstrated higher expression levels at all time points compared to CK, excluding a few cases. Such as, AhAPX55 and AhAPX140 showed relatively low expression at 9 and 12 h compared to CK and other time points (Figure 13A). In response to cold stress, although all the genes were upregulated; nevertheless, some genes showed relatively low expression levels compared to CK, such as AhAPX4, AhAPX19, AhAPX55, AhAPX82, AhAPX102, AhAPX137, and AhAPX140. Whereas AhAPX17 and AhAPX59 showed considerably higher expression than CK ( Figure 13B). In short, all the preferred genes display parallel expression trends (i.e., upregulated) to those developed from transcriptome datasets (Supplementary Figure 4), therefore representing the reliability of the transcriptome datasets.

Characterization and Evolution of APX Gene Family in Plants
Cultivated peanut is widely known as an essential oilseed, protein-enrich food crop worldwide and retains important breeding traits during domestication (Zhuang et al., 2019;Bohra et al., 2022). Even so, peanut production is still substantially influenced by numerous biotic and abiotic factors (Agarwal et al., 2018;Gangurde et al., 2020Gangurde et al., , 2021Kumar et al., 2020;Shasidhar et al., 2020;Sinha et al., 2020;Jadhav et al., 2021;Pandey et al., 2021;Soni et al., 2021;Aravind et al., 2022;Bomireddy et al., 2022;Liu et al., 2022;Patel et al., 2022). When plants are exposed to diverse abiotic and biotic factors, APX enzyme as a primary marker can quickly eliminate unnecessary H 2 O 2 (i.e., ROS scavenging) from plant cells by adjusting several physiological and biochemical activities to safeguard cells from the noxiousness of overproduction of ROS (Das and Roychoudhury, 2014;Mittler, 2017;Hasanuzzaman et al., 2020Hasanuzzaman et al., , 2021. During the past few years, excessive advancement has been achieved in studying the mode of action of APX genes; however, their vital role still requires more examination. Recent peanut genome sequencing data allow us to comprehensively discover new gene family members and recognize their functional and defensive mechanisms against stress conditions. Usually, APX gene family of plants comprises a few genes. In this study, 166 AhAPX genes have been discovered in peanut genome (Supplementary Table 10), a larger APX family than previously reported APX gene families in watermelon (Malambane et al., 2018), shrub (Wang et al., 2022), rice (Teixeira et al., 2004), A. thaliana (Panchuk et al., 2002, Panchuk et al., 2005, sorghum (Akbudak et al., 2018), kiwifruit (Liao et al., 2020), tomato (Najami et al., 2008), wheat (Tyagi et al., 2020), and cotton (Tao et al., 2018). Deviations in the APX members amongst diverse plant species may perhaps be attributed to gene duplication events involving tandem and segmental repeats and play a role in expanding APXs for deviation. Repetition of APX genes was also discovered in several plant species (Teixeira et al., 2004;Panchuk et al., 2005;Akbudak et al., 2018;Liao et al., 2020;Wang et al., 2022). Our outcomes confirmed that AhAPXs had suffered segmental duplications (Supplementary Table 3). Subsequently, these reports recommended that AhAPXs duplicate cases may possibly play an essential role in gene evolution.
Previous studies showed that APX family genes are usually clustered into four major groups based on their subcellular localization or tree topologies (Pandey et al., 2011; Frontiers in Plant Science | www.frontiersin.org 20 September 2022 | Volume 13 | Article 962182 Malambane et al., 2018;Tyagi et al., 2020). In the present study, all APX genes from four plant species were grouped into 11 main groups based on tree topologies and sequence similarities (Figure 2). This grouping was also recently supported by a new study on brassica crops (B. napus and B. rapa), where all APX genes were grouped into 13 subfamilies (Ma et al., 2021). Further, gene structure analysis also showed that genes belonging to the same subtree possess almost similar exon-intron patterns, ranging from 16 to 1 (exons) and 0 to 15 (introns) (Figure 5B and Supplementary Table 10). A significant difference was observed in gene structures where some genes have many exons-introns while some lack introns. Similar gene structure patterns were also reported in previous reports, such as in wheat number of exons extended from 7 to 12 (Tyagi et al., 2020). In Actinidia chinensis, the gene structure possesses 4-22 introns (Liao et al., 2020), which are higher than our observations. The exon-intron organization disparity was practiced by three important methods (exon/intron gain/loss, exonization/pseudoexonization, and insertion/deletion), and they are clearly supported by structural inconsistency (Xu et al., 2012). Notably, the AhAPX genes in each group almost exhibited comparable exonintron group and conserved motifs (Figure 5), suggesting that these genes may possibly contribute to the similar tasks associated with several abiotic stressors. These outcomes are in agreement with earlier studies of kiwifruit (Liao et al., 2020), and wheat (Tyagi et al., 2020), where genes inside the same group comprise distinct structures and motifs organizations.

The Contribution of APX Genes Toward Stress Responses and Tolerance Mechanisms
To boost our understanding into the involvement of AhAPX genes contrary to numerous environmental factors, cis-elements were predicted in the promoter of AhAPX genes. The discoveries showed that three types of cis-elements were recognized, i.e., abiotic stress, phytohormones, and growth and developmentrelated elements (Figures 6, 7). Recent studies show that the cis-elements in APX genes contribute to the plant abiotic and phytohormones stress responses. Similar types of abiotic and phytohormone-related cis-regulatory elements have been identified in previous studies (Akbudak et al., 2018;Malambane et al., 2018;Tao et al., 2018;Liao et al., 2020;Tyagi et al., 2020;Wang et al., 2022). Furthermore, AhAPXs gene functions were further predicted by GO enrichment analysis (Figure 10), which also supported the role of these genes in ROS scavenging and stress response mechanisms. To get further insights into the role of AhAPX genes, their expression levels were examined under various hormones and abiotic stress treatments (Figures 12, 13). Our results showed that a few genes significantly contribute to specific stress responses like  cold, drought, and ABA. These results are in agreement with the previous reports of Akbudak et al. (2018), where some SbAPX genes were significantly induced by drought stress in the leaves and roots of two genotypes. Similarly, many genes showed higher expression levels in A. nanus under cold and osmotic stress (Wang et al., 2022). Many BrAPX and BnAPX genes showed higher expression trends in cold-tolerant varieties in response to cold stress (Ma et al., 2021). Under drought stress, most of the ClAPX genes were significantly upregulated and displayed elevated expression in watermelon (Malambane et al., 2018). These conclusions can enhance our perception of AhAPX genes under various stress conditions, especially cold and drought.
Recent reports also suggest that manipulating APX genes could contribute to stress tolerance in plants. For instance, a novel ScAPX6 gene from sugarcane was overexpressed in tobacco (Nicotiana benthamiana), and transgenic plants showed improved resistance to the biotic stress (Pseudomonas solanacearum and Fusarium solani) by positively regulating the phytohormones contents (Liu et al., 2018). The overexpression of PcAPX from Chinese white poplar (Populus tomentosa) improves tolerance to multiple stresses, including salinity, drought, and oxidative stress in transgenic tobacco plants by improving biochemical mechanisms (Cao et al., 2017). Likewise, the overexpression of Populus peroxisomal PpAPX gene enhances drought stress tolerance in transgenic tobacco plants (Li et al., 2009). Transgenic tobacco overexpressing cytosolic APX gene alleviated the drought stress tolerance (Faize et al., 2011). Ectopic overexpression of the peroxisomal SbpAPX gene improves salinity tolerance in transgenic peanut (Singh et al., 2014). So far, this is the only APX gene that has been functionally characterized in peanut. These studies recommend that the genetic engineering of APX genes is of great importance in conferring various stress tolerance in crop plants, including peanut.
Among various identified TFs, ERF TFs have been functionally characterized from peanut. The results exhibited that overexpression of AhERF008 impaired the root magnitude of A. thaliana; however, overexpression of AhERF019 improved tolerance to heat, salinity and drought stresses in A. thaliana (Wan et al., 2014). Ectopic overexpression of MYB repressor gene (GmMYB3a) increases drought tolerance and physiological mechanisms in transgenic peanut under drought stress . Another NAC TF gene (AhANC4) from peanut enhances drought tolerance in transgenic tobacco plants due to improved stomatal closure and advanced water use efficiency (Tang et al., 2017). A novel WRKY TF gene (AhWRKY75) improved salinity tolerance in transgenic peanut plants by improving antioxidant mechanisms, ROS scavenging, stomatal conductance, and photosynthesis under salinity stress (Zhu et al., 2021). All these studies suggest that the genetic engineering of TF is a promising approach to improve peanut performance under stressed conditions.

The Contribution of APX Genes in Numerous Organs/Tissues
Here, the tissue-specific expression profiling of 166 AhAPX genes was carried out in various organs/tissues using publically available transcriptome datasets. Overall, the results showed that only a few AhAPX genes showed higher expression levels, particularly in roots, stem, leaf, peg, pericarp, testa, and flowers (Figure 11). In A. nanus, the RNA-seq data was used to observe the expression levels in leaves. The results displayed that only one gene showed substantially higher expression in leaf (Wang et al., 2022). In wheat, most of the genes showed higher expression patterns in root, stem, leaf, spike, and grain. Especially, almost half of the APX genes were found to be leaf-specific due to significantly higher expression (Tyagi et al., 2020). In A. chinensis, qRT-PCR-based expression profiling of 13 AhAPX genes was performed in various fruit developmental stages. The outcomes demonstrated that eight AcAPX genes had the utmost expression patterns during the color turn-off phase (Liao et al., 2020). It can be concluded that the tissue-specific APX genes (such as AhAPX4, AhAPX17, AhAPX77, AhAPX82, AhAPX101, and AhAPX130) could be considered as target candidates for further molecular studies to fully reveal their role and mechanisms in peanut growth and development.

MicroRNA: Emerging Players for Crop Improvement and Stress Tolerance
MicroRNAs (miRNAs) are a group of tiny-non-coding RNAs formed from individual-strand hairpin RNA precursors. These miRNAs switch gene expression by attaching to corresponding sequences surrounded by target mRNAs Patil et al., 2021). Extensive progress has been put together in finding the targets of peanut miRNAs that contribute to various stresses and developmental activities (Zhao et al., 2010(Zhao et al., , 2015Chi et al., 2011;Zhang et al., 2017;Figueredo et al., 2020;Tong et al., 2021). The present predicted 14 miRNAs belonging to 12 different families targeting 33 AhAPX genes (Figure 8 and Supplementary Table 6). Notably, none of the previous studies predicted the miRNAs that can target APX genes, expect one study. A recent study supports our findings where 51 miRNAs have been identified targeting 29 TaAPX genes in wheat (Tyagi et al., 2020). However, these target genes are yet to be characterized in wheat. In another study, a new miRNA (ath-miR447a-3p) was found to be targeting APX3 gene, and its expression analysis showed that it negatively regulated the expression of APX3, which is directly involved in the APX synthesis under drought stress in Zanthoxylum bungeanum (Fei et al., 2020).
However, some of the identified miRNAs have been reported to take part in stress tolerance and developmental processes. For instance, spatio-temporal expression patterns of miRNA159 family representatives have been found targeting MYB genes in grapevine (Vitis vinifera L.). The results showed that miRNA159c-VvGAMYB module is involved in gibberellin-tempted parthenocarpy in grapevine . Another study discovered that miR167A is the main member of miR167 family that regulates the A. thaliana reproduction. Further, miR167A acts as a parental gene that works mostly via ARF6 and ARF8 genes in maternal management of embryonic and seed growth (Yao et al., 2019). A member of miRNA156 family has been reported to be involved in the interaction between ABA and miRNA156, which regulates the expression profile of anthocyanin biogenesis genes in drought-stressed plants (González-Villagra et al., 2017). Notably, several miRNA families such as miR3513, miR3518, miR3520, miR3513, miR3516, etc., have not been functionally characterized; therefore, the future work could also be focused on these unique miRNAs to reveal their potential in plant growth and development. Moreover, the expression profiling of prophesied miRNAs and their targeted genes demands validation to direct their biological roles in the peanut breeding programs.

CONCLUSION
Altogether, we recognized 166 putative AhAPX genes in the cultivated peanut genome, which are mapped on all chromosomes, including unassembled ones. Comprehensive in silico examination of AhAPX genes, i.e., characterization, evolution, gene structure, conserved motifs, cis-elements, putative miRNA and TFs prediction, GO and KEGG enrichment were executed to increase our understanding of AhAPX genes in peanut. Their expression trends were also evaluated in various developmental organs/tissues, phytohormones, and abiotic stress conditions. In brief, this report set the foundation for further functional experiments (such as overexpression, gene editing via CRISPR/Cas system, etc.) of some candidate genes such as AhAPX4/17/19/55/59/82/101/102/137 and AhAPX140, which can advance the peanut breeding programs under undesirable stress conditions.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI BioProject-PRJNA480120.

AUTHOR CONTRIBUTIONS
AR and WZ conceived the idea. AR analyzed the data and wrote the manuscript. YS and KC helped with qRT-PCR