Document Type : Original Article


1 Department of Biotechnology Research, Research Institute of Forests and Rangelands, Agricultural Research, Education and Extension Organization (AREEO), National Botanical Garden, Tehran Karaj Freeway, P.O. Box 13185-116, Tehran, Iran

2 Farzanegan Campus, Semnan University, Semnan, Iran


Biotic and abiotic stresses are the main problems affecting the quantity and quality losses of antioxidants in wild cherry (Prunus avium L), which provide health benefits to humans, such as anti-allergic, anti-inflammatory, antimicrobial, and anti-cancer benefits. Glutathione peroxidase (GPX) isoenzymes act a major role in responding to oxidative stresses through the multiple metabolic pathways. Protein sequences of GPXs family form Arabidopsis thaliana and Oryza sativa were used to identify homologous proteins by the BlastP in the sweet cherry’s genome. After the alignment, the existence of the specific domain, phylogenetic relationships, subcellular localization analysis, as well as exon-intron structures and conserved motives were evaluated in this study. MicroRNA target site and codon usage bias were also investigated. Furthermore, the existence of cis-regulatory element and SSR markers were identified. The three-dimensional structure of the GPXs protein prediction and the accuracy of the predicted model was also investigated. Regarding molecular docking, the best pose for considering the accuracy of binding site was evaluated. Finally, by analyzing biological process and molecular function, gene ontology of the predicted protein was considered. We identified the five GPX genes in P. avium that classified into 3 separate groups based on their phylogenetic relationships. The predicted subcellular localization showed the GPX proteins that localized both chloroplast and mitochondria. Also, results showed ENC values were ranked 45.92 to 59.21 in PavGPX3 and PavGPX2. The results of cis‐regulatory elements showed PavGPX8 contained ten cis‐regulatory elements which was the lowest number among PavGPXs. The best pose and protein interaction in the verified 3D model were the first pose with the -28.5 minimization energy in docking process. It was concluded that the precise information on putative roles of P. avium GPXs genes could be helpful in determining their functional characterizations. The presence of several developmental process, stresses, light and hormones- responsive elements in the upstream of these genes, and GPXs‐targeted miRNAs suggested that GPX family members may play various roles in different biological responses. Finally, the GPXs gene was proposed as a good candidate for genetic engineering and breeding programs in order to introduce an abiotic stress- tolerant cultivar.

Graphical Abstract

In Silico Genome-Wide Identification and Characterization of Glutathione Peroxidase Gene Family in Wild Cherries (Prunus avium L)


Main Subjects

1. Choudhury F K, Rivero R M, Blumwald E, Mittler R. (2017). Reactive oxygen species, abiotic stress and stress combination. The Plant Journal, 90(5): 856-867. https://doi: 10.1111/tpj.13299  
2. Passaia G, Margis-Pinheiro M. (2015). Glutathione peroxidases as redox sensor proteins in plant cells. Plant Science, 234: 22-26.
3. Faltin Z, Holland D, Velcheva M, Tsapovetsky M, Roeckel-Drevet P, Handa A K, Abu-Abied M, Friedman-Einat M, Eshdat Y, Perl A. (2010). Glutathione peroxidase regulation of reactive oxygen species level is crucial for in vitro plant differentiation. Plant and cell physiology, 51(7): 1151-1162. https://doi: 10.1093/pcp/pcq082 
4. Gueta-Dahan Y, Yaniv Z, Zilinskas B A, Ben-Hayyim G. (1997). Salt and oxidative stress: similar and specific responses and their relation to salt tolerance in citrus. Planta, 203(4): 460-469.
5. Diao Y, Xu H, Li G, Yu A, Yu X, Hu W, Zheng X, Li S, Wang Y, Hu Z. (2014). Cloning a glutathione peroxidase gene from Nelumbo nucifera and enhanced salt tolerance by overexpressing in rice. Molecular biology reports, 41(8): 4919-4927. https://doi: 10.1007/s11033-014-3358-4
6. Avsian-Kretchmer O, Gueta-Dahan Y, Lev-Yadun S, Gollop R, Ben-Hayyim G. (2004). The salt-stress signal transduction pathway that activates the gpx1 promoter is mediated by intracellular H2O2, different from the pathway induced by extracellular H2O2. Plant Physiology, 135(3): 1685-1696. https://doi: 10.1104/pp.104.041921
7. Selles B, Hugo M, Trujillo M, Srivastava V, Wingsle G, Jacquot J-P, Radi R, Rouhier N. (2012). Hydroperoxide and peroxynitrite reductase activity of poplar thioredoxin-dependent glutathione peroxidase 5: kinetics, catalytic mechanism and oxidative inactivation. Biochemical Journal, 442(2): 369-380. https://doi: 10.1042/BJ20111378 
8. Zhai C-Z, Zhao L, Yin L-J, Chen M, Wang Q-Y, Li L-C, Xu Z-S, Ma Y-Z. (2013). Two wheat glutathione peroxidase genes whose products are located in chloroplasts improve salt and H 2 O 2 tolerances in Arabidopsis. PloS one, 8(10): e73989. https://doi: 10.1371/journal.pone.0073989
9. Kim Y-J, Jang M-G, Noh H-Y, Lee H-J, Sukweenadhi J, Kim J-H, Kim S-Y, Kwon W-S, Yang D-C. (2014). Molecular characterization of two glutathione peroxidase genes of Panax ginseng and their expression analysis against environmental stresses. Gene, 535(1): 33-41.
10. Gao F, Chen J, Ma T, Li H, Wang N, Li Z, Zhang Z, Zhou Y. (2014). The glutathione peroxidase gene family in Thellungiella salsuginea: genome-wide identification, classification, and gene and protein expression analysis under stress conditions. International journal of molecular sciences, 15(2): 3319-3335.
11. Islam T, Manna M, Kaul T, Pandey S, Reddy C S, Reddy M. (2015). Genome-wide dissection of Arabidopsis and rice for the identification and expression analysis of glutathione peroxidases reveals their stress-specific and overlapping response patterns. Plant molecular biology reporter, 33(5): 1413-1427.
12. Prvulović D, Malenčić D, Popović M, Ljubojević M, Ognjanov V. (2011). Antioxidant properties of sweet cherries (Prunus avium L.)–Role of phenolic compounds. World Academy of Science, Engineering and Technology, 59: 1149-1152
13. Lee T-H, Tang H, Wang X, Paterson A H. (2012). PGDD: a database of gene and genome duplication in plants. Nucleic acids research, 41(D1): D1152-D1158. https://doi: 10.1093/nar/gks1104
14. Milla M A R, Maurer A, Huete A R, Gustafson J P. (2003). Glutathione peroxidase genes in Arabidopsis are ubiquitous and regulated by abiotic stresses through diverse signaling pathways. The Plant Journal, 36(5): 602-615. https://doi: 10.1046/j.1365-313x.2003.01901.x
15. Letunic I, Doerks T, Bork P. (2012). SMART 7: recent updates to the protein domain annotation resource. Nucleic acids research, 40(D1): D302-D305. https://doi: 10.1093/nar/gkr931 
16. Thompson J D. (1994). Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res, 22: 4673-4680. https://doi: 10.1093/nar/22.22.4673
17. Felsenstein J. (1985). Confidence limits on phylogenies: an approach using the bootstrap. evolution, 39(4): 783-791. https://doi: 10.1111/j.1558-5646.1985.tb00420.x
18. Gasteiger E, Hoogland C, Gattiker A, Wilkins M R, Appel R D, Bairoch A. (2005). Protein identification and analysis tools on the ExPASy server. The proteomics protocols handbook: 571-607.
19. Krogh A, Larsson B, Von Heijne G, Sonnhammer E L. (2001). Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. Journal of molecular biology, 305(3): 567-580. https://doi: 10.1006/jmbi.2000.4315
20. Hall J á. (2002). Cellular mechanisms for heavy metal detoxification and tolerance. Journal of experimental botany, 53(366): 1-11.
21. Hu B, Jin J, Guo A-Y, Zhang H, Luo J, Gao G. (2015). GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics, 31(8): 1296-1297. https://doi: 10.1093/bioinformatics/btu817
22. Bailey T L, Boden M, Buske F A, Frith M, Grant C E, Clementi L, Ren J, Li W W, Noble W S. (2009). MEME SUITE: tools for motif discovery and searching. Nucleic acids research, 37(suppl_2): W202-W208. https://doi: 10.1093/nar/gkp335
23. Draghici S, Khatri P, Tarca A L, Amin K, Done A, Voichita C, Georgescu C, Romero R. (2007). A systems biology approach for pathway level analysis. Genome research, 17(10): 1537-1545. https://doi: 10.1101/gr.6202607 
24. Rozas J, Sánchez-DelBarrio J C, Messeguer X, Rozas R. (2003). DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics, 19(18): 2496-2497. https://doi: 10.1093/bioinformatics/btg359
25. Saeed A I, Bhagabati N K, Braisted J C, Liang W, Sharov V, Howe E A, Li J, Thiagarajan M, White J A, Quackenbush J. (2006). TM4 microarray software suite. Methods in enzymology, 411: 134-193. https://doi: 10.1016/S0076-6879(06)11009-5
26. You F M, Huo N, Gu Y Q, Luo M-c, Ma Y, Hane D, Lazo G R, Dvorak J, Anderson O D. (2008). BatchPrimer3: a high throughput web application for PCR and sequencing primer design. BMC bioinformatics, 9(1): 1-13.
27. Bitencourt-Ferreira G, de Azevedo W F. (2019). Molegro virtual docker for docking Docking Screens for Drug Discovery (pp. 149-167): Springer
28. Zhang L, Wu M, Yu D, Teng Y, Wei T, Chen C, Song W. (2018). Identification of glutathione peroxidase (GPX) gene family in Rhodiola crenulata and gene expression analysis under stress conditions. International journal of molecular sciences, 19(11): 3329. https://doi: 10.3390/ijms19113329
29. Chen M, Li K, Li H, Song C-P, Miao Y. (2017). The glutathione peroxidase gene family in Gossypium hirsutum: genome-wide identification, classification, gene expression and functional analysis. Scientific reports, 7(1): 1-15. https://doi: 10.1038/srep44743 (2017)
30. Zhou Y, Hu L, Ye S, Jiang L, Liu S. (2018). Genome-wide identification of glutathione peroxidase (GPX) gene family and their response to abiotic stress in cucumber. 3 Biotech, 8(3): 1-11. https://doi: 10.1007/s13205-018-1185-3
31. Wang Z, Qiao Y, Zhang J, Shi W, Zhang J. (2017). Genome wide identification of microRNAs involved in fatty acid and lipid metabolism of Brassica napus by small RNA and degradome sequencing. Gene, 619: 61-70. https://doi: 10.1016/j.gene.2017.03.040
32. Gielen H, Remans T, Vangronsveld J, Cuypers A. (2012). MicroRNAs in metal stress: specific roles or secondary responses? International journal of molecular sciences, 13(12): 15826-15847.
33. De Mandal S, Mazumder T H, Panda A K, Kumar N S, Jin F. (2020). Analysis of synonymous codon usage patterns of HPRT1 gene across twelve mammalian species. Genomics, 112(1): 304-311.
34. Lu B, Wang N, Xiao J, Xu Y, Murphy R W, Huang D. (2009). Expression and evolutionary divergence of the non-conventional olfactory receptor in four species of fig wasp associated with one species of fig. BMC Evolutionary Biology, 9(1): 1-16.
35. Wang P, Xiong Y, Gong R, Yang Y, Fan K, Yu S. (2019). A key variant in the cis-regulatory element of flowering gene Ghd8 associated with cold tolerance in rice. Scientific reports, 9(1): 1-14.
36. Haasl R J, Payseur B A. (2012). Microsatellites as targets of natural selection. Molecular biology and evolution, 30(2): 285-298.
37. Qin Z, Wang Y, Wang Q, Li A, Hou F, Zhang L. (2015). Evolution analysis of simple sequence repeats in plant genome. PLoS One, 10(12): e0144108. 10.1371/journal.pone.0144108