Genome-Wide Identiﬁcation and Expression Proﬁling of Germin-Like Proteins Reveal Their Role in Regulating Abiotic Stress Response in Potato

Germin and germin-like proteins (GLPs) perform a signiﬁcant role in plants against biotic and abiotic stress. To understand the role of GLPs in potato, a comprehensive genome-wide analysis was performed in the potato genome. This study identiﬁed a total of 70 StGLPs genes in the potato genome, distributed among 11 chromosomes. Phylogenetic analysis exhibited that StGLPs were categorized into six groups with high bootstrap values. StGLPs gene structure and motifs analysis showed a relatively well-maintained intron–exon and motif formation within the cognate group. Additionally, several cis -elements in the promoter regions of GLPs were hormones, and stress-responsive and different families of miRNAs target StGLPs . Gene duplication under selection pressure also exhibited positive and purifying selections in StGLPs . In our results, the StGLP5 gene showed the highest expression in response to salt stress among all expressed StGLPs . Totally 19 StGLPs genes were expressed in response to heat stress. Moreover, three genes, StGLP30 , StGLP17 and StGLP14 , exhibited a relatively higher expression level in the potato after heat treatment. In total, 22 genes expressed in response to abscisic acid ( ABA) treatment indicated that ABA performed an essential role in the plant defense or tolerance mechanism to environmental stress. RNA-Seq data validated by RT-qPCR also conﬁrm that the StGLP5 gene showed maximum expression among selected genes under salt stress. Concisely, our results provide a platform for further functional exploration of the StGLPs against salt and heat stress conditions.


INTRODUCTION
Germin-like proteins (GLPs) are pervasive proteins regarded as water-soluble glycoproteins present in either cell apoplast or its extracellular matrix (Lane et al., 1992;Dunwell et al., 2008). These proteins were first identified in wheat germinating embryos (Lu et al., 2010). At present, several proteins have been identified in different plant groups that exhibit 30-70% sequence homogeneity with wheat germin proteins (Finn et al., 2007).
Germin-like genes are expressed in embryos, leaves, stems, roots, flowers, and seeds of plants responding to various environmental stimuli based upon the species or genes under consideration (Wang et al., 2013). Abiotic and biotic stresses significantly affect plant growth and production worldwide (Raza et al., 2020;Qasim et al., 2021Qasim et al., , 2018aRaza, 2021). Various evidence has been given as an indication of GLPs involvement in defense and tolerance mechanisms of plants (Breen and Bellgard, 2010). Incidents like pathogen infections, insect infestation, or exposure to various chemicals such as salicylic acid, hydrogen peroxide, or ethylene may trigger the expression level of GLPs (Godfrey et al., 2007). Many internal elements have influenced the expression of germin, as in wheat, where the germin gene gf-2.8 is induced by hormone auxin (Schweizer et al., 1999). Transient overexpression and silencing of some Hordeum vulgare GLP genes have imparted enhanced resistance against fungal infection of powdery mildew (Schweizer et al., 1999). The promoter alternative of rice oxalate oxidase genes envisaged in the resistance against Magnaporthe oryzae. For some subfamilies, transitory and stable expression depicted the superoxide dismutase activity of the encoded protein (EC1.15.1.1) (Zimmermann et al., 2006). The silencing of GLPs in Nicotiana attenuate has proved fruitful in enhancing performance against herbivores (Lou and Baldwin, 2006). Fluctuation in the transcription of germin genes in Brassica napus and nearly associated Arabidopsis is found during the circadian cycle (Heintzen et al., 1994). These GLPs are associated with the protective mechanisms of cereals against the infection of fungal pathogens (Lane, 2000). Studies have revealed that the transfer of GLPs from soybean and sunflower to Nicotiana tabacum can increase its resistance against Sclerotinia sclerotiorum (Beracochea et al., 2015;Zhang et al., 2018). Transient-enhanced expression of TaGLP4 and HvGLP4 makes wheat and barley plants more resistant to Blumeria graminis. At the same time, the ephemeral gene that is silenced by utilizing RNA interference minimizes the basal resistance in cereals (Christensen et al., 2004, p. 4). A study is conducted to explore that the regulatory factors in various GLP gene promoters and several regulatory factors are associated with many essential roles with variable copy numbers and dispersed in locations of various positive or negative strands. Cotton GLP as GhABP19 is involved in various essential tasks in the exhibition of resistance to Verticillium and Fusarium wilt in crops (Pei et al., 2019, p. 19). A few GLP subfamily members are intricated in a vast range of fungal disease resistance in the tissues of leaves (Davidson et al., 2009).
Germin-like protein genes have variable expression and multiple tasks in several tissues and at different phases of plant growth processes (Barman and Banerjee, 2015). This gene family is one of those gene families that are crucial and indispensable for the defense and stress tolerance mechanisms of plants. Expression of these genes varies against various abiotic stress such as salt (Wang et al., 2013), heavy metal (Berna and Bernier, 1999), and drought stress. Moreover, some purposes of GLP genes have been characterized in Hordeum vulgare (Zimmermann et al., 2006), Arabidopsis thaliana (Rietz et al., 2012), Oryza sativa (Li et al., 2016), and Glycine max . The expression levels of many rice GLP genes were altered under salt stress (Das et al., 2019). OsGLP8-14 was found to be involved in salt stress response (Banerjee et al., 2017). Liao et al. (2021) reported the differential expression of CsGLP genes in response to salt, drought, and ABA treatments (Liao et al., 2021).
Potato (Solanum tuberosum L.) is an economically essential field crop and is widely grown as a basic food commodity worldwide. About 368,168,914 tons of potatoes were harvested in 2018 on 17,578,672 ha of land, and over 1 billion people consumed them (Rashid et al., 2021). Like other plant species, potato yield is also in jeopardy due to many biotic and abiotic environmental fluctuations (Zaynab et al., 2021d). To date, the expression analysis of GLP genes against abiotic stress has not been reported in potatoes. Therefore, this study used genome-wide analysis to find GLP genes in the potato genome. Furthermore, their evolutionary relationships, synteny analysis, gene structures, conserved motifs, cis-elements, and miRNA predictions were investigated. Expression profiling in many tissues/organs and response to various hormones and abiotic stresses have been substantially explored, enabling us to explore and understand these GLPs in extensive detail in the potato genome.

Identification and Characterization of GLP Genes in Potato
In this study, we used two methods to identify GLP genes from the genome of the Solanum tuberosum. The first method was the Hidden Markov Model (HMM), and the second one was the BLASTP search Zaynab et al., 2021d). Potato genome sequences were retrieved from JGI Phytozome 12.0 1 database (Goodstein et al., 2012). GLP amino acids sequences of A. thaliana were taken from TAIR 2 (Rhee, 2003) database and used as a query to run BLASTP search. The HMM file with Pfam ID PF00190 containing GLP genes was downloaded from Pfam 3 (El-Gebali et al., 2019) database. Moreover, a local HMMER 3.1 4 webserver was utilized to identify the GLP genes with default parameters (Finn et al., 2015). Eventually, a total of 70 GLP genes were identified from the genome of the potato. Non-redundant genes with conserved domains were serially named through StGLP1-StGLP70 and used for further analysis.  Physiochemical characteristics, including molecular weight and isoelectric points, were identified through ProtParam 5 (Gasteiger et al., 2005).

Phylogenetic, Gene Structure, Motif, and Synteny Analysis
In order to investigate the evolutionary history of the StGLP family genes, a phylogenetic tree was constructed among the amino acid sequences of Solanum lycopersicum, S. tuberosum, and A. thaliana. We used MEGA 7 6 (Kumar et al., 2018) for the construction of phylogenetic tree. The neighbor-joining method was applied to make a phylogenetic tree with 1,000 bootstrap values. TBtools software 7 was deployed to construct the gene structure of StGLPs . Conserved motifs present in StGLPs sequences were recognized through MEME 8 (Bailey et al., 2009) server. For comparative synteny analysis, Circoletto Tool 9 was used.

Chromosomal Locations, Selection Pressure, and Cis-Elements Analysis
Data related to the chromosomal location was taken from general feature format (GFF) files of S. tuberosum genome database (PGSC Ref seq v4.03) (The Potato Genome Sequencing Consortium, 2011). The localization of 70 identified genes on chromosomes was performed by utilizing TBtools . To explore the putative cis-elements in the promoter region of StGLPs, we extracted a sequence of 2 Kb upstream of start codons in the genome of S. tuberosum. Afterward, to analyze the promoter region of each gene, we deployed them to PlantCARE 10 (Lescot, 2002) and figures made by TBtools . Besides, the synonymous and non-synonymous substitution rate was computed by employing KaKs Calculator 2.0 software using the MYN procedure. Ks and Ka refer to the number of synonymous and nonsynonymous substitutions concerning the synonymous nonsynonymous site . Divergence time (t = Ks/2r) was estimated through exchange rate as (r = 2.6 × 10 −9 ) .

Prediction of miRNA and Protein-Protein Interactions
For identifying miRNA targeting StGLPs, psRNAtarget 11 was used with default parameters (Dai et al., 2018). For illustrating the protein-protein interaction network, the STRING 11.0 12 server was employed, and for reference, the genome of S. tuberosum was used. The exhibit specifications were propped as the confidence value is 0.9 for network; and 10 for edges evidence and maximum number of interaction.

Expression Profiling of StGLP Genes in Different Tissues
In order to investigate the expression of StGLP genes in various tissues and under various stress conditions, the transcriptome data were obtained from the available public database BioProject 10 http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ 11 https://www.zhaolab.org/psRNATarget/ 12 https://string-db.org/ ID: Project: PRJEB2430. For the expression analysis of StGLP genes in different tissues and under different stresses, fragments per kilobase million (FPKM) were observed (Zaynab et al., 2021d). The obtained data were arranged according to their expression in different tissues. The observed FPKM values were utilized to construct a heatmap through TBtools after the normalization of log2 values.

Potato Growth Conditions and Stress Treatment
The rooster variety of S. tuberosum was selected as the sample for these experiments, which were conducted at the experimental station of the National Institute for Genomics and Advanced Biotechnology (NIGAB), Islamabad, Pakistan. After 6 weeks of growth, the sample was treated with 150 mM NaCl for 0 and 24 h, respectively, for salt treatment. The first time point (0 h) served as a control. After NaCl treatment, leaves were taken, and samples were immediately frozen in liquid nitrogen and stored at −80 • C for total RNA extraction.

RNA Extraction and RT-qPCR Analysis
Total RNA was extracted using TransZol Up Plus RNA Kit (TransGen Biotech, Beijing, China). For cDNA synthesis, cDNA Synthesis SuperMix was purchased from the above-mentioned Chinese company. Instructions given by manufacturers were followed during the whole process. The comprehensive knowledge about reactions of RT-qPCR has been reported in our earlier published work (Zaynab et al., 2021a). In our study, the expression data were evaluated by using the 2 − CT methods (Zaynab et al., 2021c), and the elongation factor 1-alpha was chosen as a housekeeping gene (Zaynab et al., 2021b). The primers used in this study are given in Supplementary  Table 1.

Identification, Phylogenetic, Gene Structure, Motif Analysis of GLP Genes in Potato
A total of 70 GLP genes were identified in the potato genome by employing sequences of A. thaliana as a query (Table 1). In this study, Batch-CD and SMART were utilized to confirm the domain with Pfam ID (PF00190). All these 70 genes contain a domain with Pfam ID (PF00190). The protein sequences of A. thaliana, S. lycopersicum, and identified StGLPs are given in Supplementary Table 2. Amino acids length varied from 94 to 567, the molecular weight from 10.650 to 65.780 kDa, and the isoelectric point values from 4.46 to 9.74 (Table 1). A phylogenetic tree was constructed using the neighbor-joining method through MEGA7 to illustrate the phylogenetic relationship among A. thaliana, S. lycoperiscum, and S. tuberosum. This study results showed that the phylogenetic tree was divided into six groups. Group-6 is the largest group, and group-3 is the smallest among all groups (Figure 1). Exons and introns are two main determinants in the gene family evolution. The detailed analysis of the phylogenetic relationship and gene composition illustration corroborated our understanding of the structures of StGLPs genes (Figure 2). The Frontiers in Plant Science | www.frontiersin.org results revealed that numerical attributes of exon and intron in StGLP genes possessed high inconsistency in numbers of exons and introns (Figure 2). In this study, the StGLP46 gene possessed the most extended sequence. Some genes such as StGLP41 and StGLP38 exhibited similarities in exon numbers. Additionally, StGLP21 and StGLP57 genes were equal in exons numbers but differed in their sequence lengths. For determining structural diversification and functional characterization of these StGLPs, we analyzed full-length sequences of 70 StGLPs through MEME software to localize their conserved motifs. The motif results revealed that a total of ten motifs were identified (Figure 2). The sequences of predicted motifs are given in Supplementary Table 3. All of the StGLPs genes contained motif 1 except genes StGLP StGLP18, StGLP26, StGLP29, StGLP31, StGLP51, StGLP52, StGLP63, StGLP65, StGLP21, StGLP6, StGLP20, StGLP12, StGLP22, StGLP8, and StGLP27. The length of the motif was varied from 21 to 50 amino acids. The highest 50 amino acids were possessed by motif 3,7, and 10, while motif 2 and 5 contained the lowest 21 amino acids. In our results, motif 4 had 39 amino acids, motif 8 had 35 amino acids, while motifs 1 and 9 had 41 amino acids in their sequence. In short, the reliability of group organization was exceedingly encouraged by analyzing conserved motif conformation, genetic structures, and phylogenetic relationship, indicating that StGLP proteins possess enormously well-sustained amino acid members and remain inside a group. Thus, it can be inferred that proteins with identical motifs and structures may possess functional similarities.

Chromosomal Locations, Synteny, and Gene Duplication Analysis
Chromosomal distribution analysis of StGLPs depicted their unequal localization on all chromosomes (Figure 3). A maximum number of 29 genes were present on the 1st chromosome, while the 9th chromosome had 11 genes. There are 9 genes present on chromosome number 11. Only 2 genes are located on chromosome number 02 (Figure 3) in our results. Furthermore, there was no gene located on chromosome number 08. To determine the rate of molecular evolution of each gene pair that has been duplicated, the Ka to Ks (synonymous/nonsynonymous substitution) ratios were calculated. When the value of the ratio is more than 1, then the selection effect will be positive; if its value is less than 1, then purifying selection will be exhibited, while when this ratio's value is equal to 1, then there will be the presence of neutral selection in gene pairs (Yang and Bielawski, 2000). The cosmic majority of StGLPs genes exhibited a value of Ks more than 0.52, and the corresponding deviation time may be more than 100 million years ago (MYA). In this study, duplicated genes (StGLP4/StGLP35) showed the Ks value 3.54 while its corresponding time may be 682.30 MYA ( Table 2). Comparative synteny analysis was analyzed among A. thaliana, S. tuberosum, and S. lycoperiscum which exhibited a phenomenal relationship in the evolution expression, duplication, triplication, and function of genes. Genes sequence of SlGLP17 exhibited synteny with the sequence of StGLP2. At the same time, two potato genes, StGLP20 and StGLP22, showed synteny with the SlGLP41 gene of tomato. The comparative synteny analysis  revealed the conservation and duplication of inward and outward tangling events in ribbons, respectively (Figure 4).

Cis-Elements and miRNA Analysis
To observe the StGLPs gene function and their regulatory role, promoter cis-elements were analyzed using the PlantCARE database. The cis-elements that are associated with plant hormones were identified. As shown in Figure 5, the number of genes associated with phytohormones indicates the importance of these genes in phytohormone stress. Furthermore, the role of anaerobic, light, and low-temperature responsive factors in this study was discovered (Figure 5). Mainly several elements that are responsive to light were categorized to be extensively distributed in all genes, indicating the significant role of StGLPs in response to light stress. Generally, findings envisaged that diversification in the expression levels of various StGLPs might be attributed to different plant hormones and environmental stress. The Hormones-and stress-related cis-elements found in the promoter regions of StGLPs are given in Supplementary Table 4. Many researchers have reported that microRNA mediated a significant role in enhancing stress response in plants Xu et al., 2021). Thus, for an in-depth comprehension of miRNA-triggered post-transcriptional modulations of StGLPs, we identified miRNAs targeting StGLPs genes. Several sites that are targeted by miRNA are illustrated in Supplementary Table 1. Our findings revealed that seven members of the miR166 family target six genes, including StGLP8, StGLP54, StGLP23, StGLP24, StGLP38, and StGLP28. A gene, StGLP2, was targeted by a single member of the miR169 family. Targets StGLP51, StGLP48, and StGLP47 (Figure 6 and Supplementary Table 5) are all three members of the miRNA gene family. The findings of our study revealed that two members of the miR8040 aimed at two genes, StGLP51 and StGLP62. A gene StGLP55 was the target of two members of the miR1886 family. Our results mainly estimated two genes, StGLP40 and StGLP8, to target numerous miRNAs.
The protein-protein interaction network provided details about StGLPs in the genome of potatoes. The interaction network was observed among several proteins, including StGLP57, StGLP44, StGLP46, StGLP48, and StGLP2. These proteins depicted the concurrence, co-expression, fusion, and homology of genes. Likewise, another framework was also identified among StGLP53, StGLP63, and StGLP4 protein. Among all StGLPs, two proteins, StGLP48 and StGLP46, are responsible for framing a protein-protein interaction network, so these can be reputed as hub proteins of StGLPs (Figure 7).

Tissue-Specific Expression Profiling
Expression profiling of StGLPs revealed that many genes exhibited high expression in the tissues of roots than that in the tissues of stems and leaves. Some genes showed expression in some tissues, while some genes have no expression in others. For example, the StGLP17 gene showed a higher expression level in the tissues of roots, leaves, and stems, while the StGLP26 and StGLP45 genes showed higher expression levels in the tissues of roots and stems. Furthermore, in the tissues of leaves, StGLP5 genes showed the maximum expression in our results (Figure 8).

Expression Patterns of StGLP Genes in Response to Heat and Salt Stress
Transcriptional profiling of GLPs was performed when potato plants were subjected to heat stress. StGLPs may help potato plants cope better with high temperatures, according to some research. A total of 19 StGLPs genes were expressed in response to heat stress. Moreover, three genes, StGLP30, StGLP17, and StGLP14, exhibited a relatively higher expression level in potatoes after heat treatment (Figure 9). In our results, StGLP30 gene expression was higher than StGLP17 and StGLP14.
Furthermore, StGLP16 genes showed the lowest expression compared with other expressed StGLPs after heat treatment. Likewise, 29 GLPs were expressed in response to salt treatment. This study showed that StGLP54, StGLP12, and StGLP5 genes were highly expressed in potatoes after salt treatment (Figure 9), but the expression level of the StGLP5 gene was relatively higher as compared to StGLP12 and StGLP54 genes in the potato genome.

Expression Patterns of StGLP Genes in Response to Phytohormones
Three hormones, abscisic acid (ABA), gibberellic acid (GA3), and indole acetic acid (IAA), were selected to evaluate the hormonal response of StGLPs.
This study results showed that 22 StGLP genes were expressed after the treatment of ABA. Similarly, in the case of IAA treatment, 24 StGLP genes were expressed. Furthermore, the FIGURE 4 | Synteny analysis of GLPs among potato, tomato, and Arabidopsis genomes. In comparative synteny analysis, the color and intensity of inward and outward tangling of ribbons showed conservation and duplication events, respectively.
StGLP12 gene was highly expressed compared with others after the treatment of ABA, while StGLP5 showed the highest expression after IAA treatment (Figure 10). Additionally, 27 StGLPs were expressed when treated with GA3, and the expression of gene StGLP14 was highest among all 27 expressed genes (Figure 10). The upregulation of StGLPs after ABA, IAA, and GA3 treatment revealed the significant role of StGLP genes in mitigating the lethal effects of different hormonal and environmental stress.

RT-qPCR Analysis
Real-time-qPCR was performed to confirm the expression of five genes, namely, StGLP63, StGLP17, StGLP51, StGLP14, and StGLP47, in the tissues of roots, stems, and leaves (Figure 11). Two genes, namely, StGLP17 and StGLP14, exhibited a comparatively higher expression level in the tissues of stems than that in the tissues of roots and leaves. In contrast, the expression of StGLP47 was relatively more in the tissues of leaves than that in the tissues of roots and stems. Generally, findings depicted that StGLP genes may perform particular functions in potato plant growth and developmental process. Similarly, RT-qPCR was performed to confirm the expression of five genes, namely, StGLP5, StGLP12, StGLP30, StGLP 36, and StGLP54, in leaves treated with 150 mM NaCl. The results showed that the StGLP5 genes showed the maximum FIGURE 6 | miRNAs targeting StGLP genes. Network figure of anticipated miRNAs targeting the StGLP genes. Green color represents miRNAs, while orange color corresponds to the target StGLP gene. expression after 24 h of NaCl stress. In addition, genes such as StGLP12, StGLP 36, and StGLP54 were highly expressed with the extension of stress time. Moreover, StGLP30 showed expression but low expression compared with the other four genes (Figure 12).

DISCUSSION
Potato is an integral part of the food chain worldwide (Zaynab et al., 2017b). Potato plants are subjected to different biotic and abiotic stresses (Zaynab et al., 2021b). GLPs have a significant role in plant signaling and response to various environmental stresses (Yuan et al., 2021). Members of GLPs show variations in their expression and possess considerable functions in several tissues at different phases of growth cycle of the plant (Barman and Banerjee, 2015). This crucial gene family carried out vital defense mechanisms and is an essential part of plant defense processes. The expression of the essential genes increases exponentially when plants encounter stress like salinity (Wang et al., 2013;Zaynab et al., 2021d), water shortage (Komatsu et al., 2010), heavy metals (Berna and Bernier, 1999), and heat stress (Qasim et al., 2018b). Various functions of GLPs have also been determined in different crops such as Hordeum vulgare (Zimmermann et al., 2006), A. thaliana (Rietz et al., 2012), Oryza sativa (Li et al., 2016), and Glycine max . But till now, there has been no considerable in-depth study of GLPs in the potato. These days, the whole-genomic sequence of potato is comfortably accessible, allowing us to conduct genomewide analysis of the GLP genes family that may prove worthy for further research.
Structural analysis of StGLPs will prove beneficial for their functional characterization. The utilization of the evolutionary relic of GLPs has established that the evolution of genes has been affected by the arrangement of their exons and introns (Moore and Purugganan, 2005;Flagel and Wendel, 2009). This is in correlation with previous scientists finding that genes with short introns or no introns are retained in the genomic systems of plants during the evolutionary process. At the same time, the expression levels of such genes are also low in plants (Mattick and Gagen, 2001). Moreover, the compact arrangement of genes favors the expression response of the GLPs against different exo or endogenous stimuli (Jeffares et al., 2008). Our findings through structural analysis of these genes inferred that those sequences of StGLP depicted a similar number of exons-introns with identical functional properties because they may have arisen through duplication events in the course of evolution (Waqas et al., 2019).
Studying the defense mechanisms employed by plants against various types of environmental stress, cis-elements were investigated in the promoter region of these genes. According to our findings, different cis-elements were identified that show their response against light and hormone stress. Several identified cis-elements were associated with ABA, SA, GA, MeJA, low temperature, anaerobic induction, and light. Earlier reports have attributed the stress response of plants to these cis-acting elements (Osakabe et al., 2014). Cis-factors are integral parts of plant stress response machinery as they regulate stress-responsive genes (Wu et al., 2014). So, these crucial Cis-acting sites of StGLPs indicate their behavior when different stresses are imposed.
All identified StGLP genes were classified into six clades. Similarly, phylogenetic grouping was also observed in other plants. The comparative phylogenetic analysis revealed that the organization of S. lycoperiscum, S. tuberosum, and A. thaliana proteins was relatively similar in six clades, indicating that all StGLP genes in these groups may have descended from a common ancestor. Previous studies (Li et al., 2016) reported the classification of GLPs into six clades in rice and Arabidopsis, supporting our results.
The size of the genome, duplication events in the genome, and the distribution of genes are major determinants of genetic variability in land plants. Genetic duplication property has long been recognized throughout the origins of evolutionary novelty, expression, and complexity of the gene families. We also noticed some duplications in StGLPs that may carry out vital roles in their amplification. As gene duplication is an imperative character in neofunctionalization, expansion, and variability among gene families (Lavin et al., 2005), similarly, the arrangement and localization of StGLP genes at chromosomes will enable breeders of potato to breed better varieties with desirable features.
Several types of research have depicted the response of plants to salinity and other abiotic stresses, and it is now well understood that miRNAs and their objects have direct impacts on the stress tolerance processes of plants (Li, 2015;Villanueva et al., 2016). This study recognized micro RNAs from various families targeting StGLP genes. Likewise, the unique role of miR156 has been widely documented during various abiotic stress in several plants (Cui et al., 2014;Kohli et al., 2014). Similarly, a study on cotton established a critical role of miR827 in response to saline stress (Covarrubias and Reyes, 2010); miR167 has also been FIGURE 8 | Expression profiling of StGLP genes in leaves, roots, and stems. The orange and blue colors display high-to low-expression levels. The expression heat map was created by taking log10 normalized values of fragments per kilobase million (FPKM). recognized as one of the significant role players in response to various stressed environmental conditions (Khraiwesh et al., 2012). Concisely, these results are in concordance with our findings that revealed that miR156, miR167, and miR827 are significant determinants of tolerance in plants against several stresses as they alter the expression of StGLP genes in the genome of potatoes.
Germin-like gene miRNAs are generally present in embryos, cotyledons, leaves, stems, flowers, roots, and seeds. Some are expressed in response to different environmental conditions; this depends upon genes under analysis or species in which they are expressed. Many authentications have indicated that various GLPs may have a critical role in the defense mechanisms of plants (Breen and Bellgard, 2010). So, the expression profiling of GLPs in the potato is helpful for their in-depth comprehension. Many genes exhibited a high level of expression in tissues, suggesting the putative role of these genes in potato growth and developmental processes. Our RT-qPCR findings revealed that the upregulation of genes in the tissues of roots and stems indicated the role of GLP in potato plant growth.
Furthermore, GLPs expression was observed in response to hormones and abiotic stress. These findings are following earlier reports that GLPs are an integral part of the defense system in saline tolerance (Wang et al., 2013), drought stress FIGURE 9 | Expression profiling of StGLP genes in response to salt and heat treatments. The orange and blue colors display high-to low-expression levels. The expression heat map was created by taking log10 of fragments per kilobase million (FPKM). (Komatsu et al., 2010), cold tolerance, and heavy metals tolerance (Berna and Bernier, 1999). Up till now, numerous researches have depicted that GLPs are significant players of plant stress response mechanisms. In our results, StGLP14, StGLP17, and StGLP30 are upregulated in response to heat stress, and StGLP5, StGLP12, and StGLP54 are highly expressed FIGURE 10 | Expression profiling of StGLP genes in response to ABA, IAA, and GA3. The orange and blue colors display high-to low-expression levels. The expression heat map was created by taking log10 of fragments per kilobase million (FPKM). in response to salt stress supported by the findings of Wang et al. (2013). Hormones may alter the physio-biochemical pathways or metabolisms of plants by different pathways of signal transduction (Zaynab et al., 2017a;Fatima et al., 2021). ABA and IAA are major hormones of the immune system in plants. Several studies have established that besides being part of plant response to different stresses, GLPs also are involved in developing systems and hormonal signaling (Wang et al., 2020). To study how hormone signaling affects the production of StGLPs, leaves of potato were treated with ABA, GA3, and IAA, and the expression of genes was investigated. There was the induction of 22 genes on the treatment of  ABA, 24 genes on the treatment of IAA, and 27 genes on the treatment of GA3, suggesting that various members of StGLPs carry out crucial roles in the plant defense system induced by these hormones. When the treatment of ABA and IAA was made, increasing the expression level of potato GLPs exhibited that these hormones possess an essential part in this defense mechanism, similar to the results obtained by Wang et al. (2020). The strong correlation between clusters of genes and their expression was noted under various stresses and tissues.

CONCLUSION
In this study, we identified 70 GLPs in the genome of the potato by genome-wide analysis. The gene structure, synteny, phylogenetic motifs, promotor, and miRNA analysis were performed to enhance our comprehension. The expression profiling after various stresses has been analyzed. The findings exhibited that GLPs genes responded significantly against hormonal and abiotic stresses, strengthening our comprehension. In our results, the StGLP5 gene showed the highest expression in response to salt stress. So, further analysis is needed to corroborate the persistent role of GLPs genes in response to salt stress.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
MZ, JP, and YS gave the idea. MF, MA, and RA-Y performed the experiments. KK and SA wrote the manuscript. SL revised the manuscript. All authors contributed to the article and approved the submitted version.