1. Introduction
Aspartic proteases (APs) (EC 3.4.23) are widely distributed among living organisms and viruses [1,2,3,4]. These enzymes have been extensively studied and constitute one of the four superfamilies of proteolytic enzymes. APs are characterized by the presence of two aspartic acid residues located within the conserved Asp-Thr/Ser-Gly motif, responsible for catalytic activity [2,5,6].
The MEROPs database classified APs into 15 families based on their amino acid sequence similarity and grouped them into 5 different clans based on their evolutionary relationship and tertiary structure [7]. Plant APs belong to families A1, A3, A11, and A12 of clan AA and family 22 of clan AD [8]. Most of the plant APs belong to the A1 family, are active at acidic pH, are specifically inhibited by pepstatin A, and present a great structural diversity [8,9]. In silico analysis of the Arabidopsis thaliana genome revealed the presence of 51 genes that encode possible APs (AtAP), 46 of which presented different primary structure characteristics compared to the canonical forms. Based on the sequence of the active sites and the organization of domains, AtAPs are currently sorted into three groups or categories: typical, atypical, and nucellin-like [10].
Typical plant AP precursors present a protein domain of 100 amino acids long known as the plant-specific insert (PSI), which is highly similar to saposin-like proteins [8]. In monomeric typical APs, this domain still remains in the mature enzyme [11], while PSI is removed from the precursors upon activation in most of the heterodimeric typical APs enzymes [8]. Atypical and nucellin-like APs share several common features such as the absence of the internal segment plant-specific insert in their sequence, variable N-terminal regions, differences in active site environment, unusually high cysteine content and localization, and have different primary structure organizations from typical APs [10,12]. Coincidentally with A. thaliana, most potential APs sequences from grape and rice genomes correspond to atypical and nucellin-like APs [9]. Regarding typical APs, their expression has been determined simultaneously in several different plant organs, including flowers, leaves, roots, stems, and seeds, as well as in different stages of development [9,13,14,15]. Typical AP functions have been associated with development, growth, lipid metabolism, and protein degradation [16,17].
Only nine atypical and one nucellin-like APs have been partially studied and characterized in A. thaliana [17,18,19,20,21,22,23,24,25]. In addition, atypical APs have been reported in other plant species such as Oryza sativa [16,26,27,28,29], Lotus japonicus [30], Fagopyrum tataricum [31,32], Solanum tuberosum [33], Nicotiana tabacum [34], several carnivorous plants [35,36,37], Phaseolus vulgaris [38], and Vitis vinifera [39]. Biological functions of atypical APs have been related to plant defense response mechanisms, hybrid sterility, reproductive development, abiotic and biotic stresses, chloroplast homeostasis, and lateral root development [12,13,14,15,16,17]. Current knowledge on the participation of APs in plants’ signaling pathways in response to water deficit and their relationship with plant stress tolerance was summarized recently [40]. These AP roles have been studied in several plant species [21,33,41,42,43,44]. In addition, we have demonstrated that leaves from A. thaliana plants overexpressing a typical AP (AtAPA1) showed changes in morphological and physiological features related to water loss deficit [45]. We concluded that AtAPA1 participates as intermediate in the ABA-induced stomatal closure, as well as in the stomatal density regulation conferring tolerance to mild water deficit.
For the last two decades, we focused our research on the biochemical and molecular characterization of Solanum tuberosum L. aspartic proteases (StAPs) [11,46,47,48,49]. Additionally, we reported the role/s of two typical StAPs (StAP1 and StAP3) in the plant’s mechanisms in response to water deficit and pathogen attack [45,46,47,48,49]. Potato is an ancestral crop domesticated at least 7000 years ago [50,51], and today it represents the third most important crop for human consumption and is the fourth most important crop in terms of production after wheat, rice, and maize (
The aim of this work was to identify and characterize the AP genes present on the S. tuberosum genome, focused on abiotic stress conditions. Data obtained provide new knowledge about gene structure and protein organization, tissue, and subcellular localization of these enzymes in potato plants. In addition, this study provides the baseline material for functional gene analysis and for potato breeding programs, including gene-editing mediated by CRISPR.
2. Results and Discussion
2.1. StAPs Potato Genome-Wide Identification
The Hidden Markov Model (HMMER3.1) profile, built from 59 aspartic protease amino acid sequences from A. thaliana (AtAPs), was tested against the TAIR 10 2010-12-2014 peptide database rendering 78 AtAPs amino acid sequences up to an e-value of 1.5 × 10−21 and a score of 77 on full sequence (Table S3). The HMMER 3.1 search allowed the detection of 149 amino acid sequences on the Potato Genome Sequence Consortium Doubled Monoploide 3.4 version (PGSC DM3.4) peptide database with a maximum e-value of 0.0024 and a score of 18.5 on full sequence (Table S4). These corresponded to 97 tentative StAP genes.
A total of 98 AP amino acid sequences codified by 62 AP genes were retrieved from the potato genome. The longest amino acid sequence was selected from each gene that was confirmed to have the Asp domain (PF00026), including two conserved aspartic acid catalytic residues (Figure 1). From these 149 amino acid sequences, there were 54 sequences with partial or absent Asp domain PF00026 (PFam database) and were thus not taken into consideration for further analyses. The corresponding genes, transcripts and protein names, lengths, and chromosome positions are presented in Table S5.
This AP genome representation is close to what has been previously reported: 51 and 67 AtAPs reported in A. thaliana [10,52], 50 APs in V. vinifera (VvAPs) [9], 67 APs in P. trichocarp (PtAPs) [53], and 96 APs from O. sativa (OsAPs) [16].
2.2. Phylogenetic Analyses
In order to know the phylogenetic relationship between the 62 identified StAPs and previously reported APs, we constructed a phylogenetic tree with 51 known AtAP amino acid sequences from the plant model species A. thaliana [10] and 6 reference AP amino acid sequences: Nucellin (U87148), Procardosin A (CAB40134), StAP3, OsAP2 (AAK81699), Nepenthesin I (AB114914), and Phytepsin (CAA39602) (cited in Materials and Methods section).
The phylogenetic analyses evidenced the presence of eight aspartic protease clusters, which were named from “a” to “h” (Figure 1). These clusters grouped typical, nucellin-like, and atypical APs in a similar distribution observed in other plant species [9,16,53]. Cluster “a” represented typical aspartic proteases, with three reference typical aspartic proteases previously cited (StAP3, Procardosin A, and Phytepsin) and four AtAPs reported for this category [10]. Cluster “a” included six StAPs: StAP1.6, StAP7.5, StAP2.2, StAP2.4, StAP3.1, and StAP9.3. Similar numbers of typical APs are present in grape (V. vinifera) [39], poplar (Populus sp.) [53], and rice (O. sativa) [16]. As expected, StAP1.6 clustered closely with its homolog from S.tuberosum ssp. tuberosum cv. Pampeana, StAP3 [11]. StAP2.2 and StAP 2.4 were related with two typical AtAPs, NP_176419, a seed-specific AP [54], and NP_172655, while StAP7.5 was closely related with Phytepsin from barley (Hordeum vulgare) [55] and Procardosin A from cardon (C. carduculus L.) [56].
Cluster “d” represents nucellin-like APs, since it included the two reference proteins Nucellin (U87148) from barley [57] and OsAP2 (AAK81699) from rice [27] and four AP proteins from A. thaliana [10]) for this category. This cluster grouped four StAPs: StAP6.1, StAP6.2, StAP6.3, and StAP6.5. Possible orthologous pairs present in this cluster were StAP6.2/NP_567922.2 without known function and StAP6.5/NP_564539.1_APCB1 involved in triggering autophagy and defense resistance [58]. Clusters “b”, “c”, “e”, “f”, “g”, and “h” grouped atypical aspartic proteases. Atypical aspartic proteases are divided into six categories [10] based on structural primary sequence features that could be reflected in the clustering.
Cluster “b” and “c” included mostly AtAP proteases (8 out of 10) from atypical category C4 (with C-terminal extension), 2 AtAPs from category Cx (not assigned type), and 1 from category C1. Cluster “b” included StAP3.2 and StAP2.6 related to NP_849967.1, StAP2.5 joined with NP_195313.2, all from category C4, which present a similar primary structure with a C-terminal extension, StAP4.3 (also with a C-terminal extension), and StAP0.1, closely related to NP_196570.1 (Category Cx). Cluster “c” presented StAPs related to AtAP category C4, such as StAP11.2/NP197676.2, StAP10.3/NP_563808.1, and StAP3.3/StAP4.2/StAP5.5/NP186923.1, while StAP10.2 and StAP9.2 grouped with NP181205.2 (category C1). All these last three sequences present a C-terminal extension characteristic of Category C4.
Cluster “e” included three AtAPs from category C1 related with StAP5.3, StAP6.7, StAP11.1, StAP11.3, and three AtAPs from category C2 (in [10] with a serine-rich pro sequence). These AtAPs were close related to a StAP each NP_565911.1/StAP2.3, NP_195839.1_PCS1/StA9.1, and NP19325.1/StAP8.11, however, none of the StAPs presented the characteristic serine pro sequence.
Cluster “f” included all three AtAPs from category C5 (“atypical without the DTS active site motif”), which were related with StAP1.4, StAP4.4, and StAP12.1; two AtAPs from category C3, which includes NP196638.2 homolog to CND41 [34], were closely related with StAP8.7, StAP8.8, and StAP8.9, and NP_565219.1 from category Cx was related with StAP6.4 and StAP3.4.
Cluster “g” grouped AtAPs from category C3 (known as “atypical with long prosegment”) and two sequences from category Cx and C1. This cluster contained NP_188478.1 ASPG1, which participates in drought tolerance and protein degradation [21], related with StAP3.5 and StAP6.6. In this cluster was also present NP_187876.2_NANA, closely related with StAP5.1, which plays regulatory functions in chloroplasts [22].
Cluster “h” includes AtAPs mostly from category C1, two AtAPs from category Cx, and one AtAP from category C2. Among category C1 there were present NP_198319.1_CDR1 and NP_176663.1, closely related to StAP1.1. NP_198319.1_CDR1 is a rice ortholog linked to disease resistance mechanisms in A. thaliana [26]. Cluster “h” also grouped Nepenthesin I from N. distillatoria [59] reference sequence closely related to StAP5.6.
To analyze StAPs phylogenetic relation together with protein and gene structures, an additional phylogenetic tree was built among the identified StAPs and eight reference APs from A. thaliana, H. vulgare, N. distillatoria, O. sativa, and C. cardunculus (Figure 2). StAPs were grouped in eight clusters, named from StAPI to StAPVIII, with a similar distribution to Figure 1. StAPI grouped sequences as described in cluster “a” (Figure 1) (StAP1.6, StAP7.5, StAP2.4, StAP2.2, StAP9.3, and StAP3.1). StAPIII presented two subgroups, one that was named StAPIIIa, which included nucellin-like StAPs (StAP6.5, StAP6.1, StAP6.2, and StAP6.3) as in cluster “d” (Figure 1), and another subgroup, StAPIIIb, representing atypical sequences from cluster “c” (StAP10.2, StAP9.2, StAP11.2, StAP10.3, StAP5.4, StAP9.4, StAP3.3, StAP4.2, StAP5.5). In both phylogenetic trees, some atypical sequences related to categories C4 and Cx were observed (Cluster “b” and “c” from Figure 1 and Cluster StAPII and subgroup StAPIIIb from Figure 2). A similar AtAP distribution was obtained in a phylogenetic tree where Populus and grape AP amino acid sequences were studied (Figure 1 in [53]).
To analyze StAPs phylogenetic relation together with protein and gene structures, an additional phylogenetic tree was built among the identified StAPs and eight reference APs from A. thaliana, H., N. distillatoria, O. sativa, and C. cardunculus (Figure 2). StAPs were grouped in eight clusters, named from StAPI to StAPVIII, with a similar distribution to Figure 1. StAPI grouped sequences as described in cluster “a” (Figure 1) (StAP1.6, StAP7.5, StAP2.4, StAP2.2, StAP9.3, and StAP3.1). StAPIII presented two subgroups, one that was named StAPIIIa, which included nucellin-like StAPs (StAP6.5, StAP6.1, StAP6.2, and StAP6.3) as in cluster “d” (Figure 1), and another subgroup, StAPIIIb, representing atypical sequences from cluster “c” (StAP10.2, StAP9.2, StAP11.2, StAP10.3, StAP5.4, StAP9.4, StAP3.3, StAP4.2, StAP5.5). In both phylogenetic trees, some atypical sequences related to categories C4 and Cx were observed (Cluster “b” and “c” from Figure 1 and Cluster StAPII and subgroup StAPIIIb from Figure 2). A similar AtAP distribution was obtained in a phylogenetic tree where Populus and grape AP amino acid sequences were studied (Figure 1 in [53]).
To analyze StAPs phylogenetic relation together with protein and gene structures, an additional phylogenetic tree was built among the identified StAPs and eight reference APs from A. thaliana, H., N. distillatoria, O. sativa, and C. cardunculus (Figure 2). StAPs were grouped in eight clusters, named from StAPI to StAPVIII, with a similar distribution to Figure 1. StAPI grouped sequences as described in cluster “a” (Figure 1) (StAP1.6, StAP7.5, StAP2.4, StAP2.2, StAP9.3, and StAP3.1). StAPIII presented two subgroups, one that was named StAPIIIa, which included nucellin-like StAPs (StAP6.5, StAP6.1, StAP6.2, and StAP6.3) as in cluster “d” (Figure 1), and another subgroup, StAPIIIb, representing atypical sequences from cluster “c” (StAP10.2, StAP9.2, StAP11.2, StAP10.3, StAP5.4, StAP9.4, StAP3.3, StAP4.2, StAP5.5). In both phylogenetic trees, some atypical sequences related to categories C4 and Cx were observed (Cluster “b” and “c” from Figure 1 and Cluster StAPII and subgroup StAPIIIb from Figure 2). A similar AtAP distribution was obtained in a phylogenetic tree where Populus and grape AP amino acid sequences were studied (Figure 1 in [53]).
The remaining clusters correspond to the rest of atypical APs: StAPIV includes five StAPs and the reference NP_565219.1 aspartic protease; StAPV includes seven StAPs and the reference NP_563851.1_AED3 aspartic protease; StAPVI includes six StAPs; StAPVII includes nine StAPs. Finally, StAPVIII includes 10 StAPs and the Nepenthensin I reference AP.
2.3. StAPs Domains and Conserved Motifs
In order to characterize StAPs at the primary structure level, the presence of domains and conserved motifs were analyzed on the 62 amino acid sequences (Figure 2A, B).
The main StAPs domains scanned with the Pfam database are shown in Figure 2A, while Figure 2B provides complementary primary sequence-structure information with motifs assessed with Multiple EM for Motif Elicitation (MEME) online software. The majority of StAPs presented a signal peptide and a pro sequence segment, in accordance with previous reports for typical, atypical, and nucellin-like aspartic proteases, [9,10,12,13,53,60]. Transmembrane signals at the C-terminal end were presented exclusively in Aps from StAPII and StAPIII groups. Figure 2A also showed aspartic domains (Asp) (PF00026) that contained the catalytic aspartate at the active site in all StAPs. Xylanase inhibitor domains (PF14543), called TAXI_N and TAXI_C [61,62,63], were found overlapping with the Asp domain, and they were thus not included in Figure 2A. Both TAXI domains were present in all StAPs except for the TAXI_C domain that was absent in typical StAPs. These domains were described as jointly necessary for creating a catalytic pocket for xylanase cleavage, a product generated by pathogens that destroys plant cells [62,63,64]. The role of these domains in aspartic proteases require further studies [64]. Another characteristic feature is the presence of the PSI domain or Saposin-like SapB 1 and SapB 2 domains [65] in all aspartic protease sequences from group StAPI, coincident with the reported for the close related typical AP reference sequences.
Regarding StAP motifs assessed with MEME online software (Figure 2B), the number of motifs per protein ranged from 4 to 12, with a length between 15 and 50 amino acids long.
A total of six StAP sequences clustered in Group StAPI along with the typical aspartic proteases, presented at least six group exclusive motifs, such as 7, 8, 11, 15, 16, and 19. Motifs 8 and 16 were detected as part of PSI domains SapB_1 and SapB_2 respectively. While motifs 1, 8, 11, and 16 were characterized by the presence of conserved cysteine residues.
Subgroup StAPIIIb exhibited motif 20 at the N-terminal region in most sequences, and a remarkable characteristic is the substitution of hydrophobic by polar amino acid (glutamine) at the surroundings of the first catalytic site (Table S1).
Group StAVIII related to NepenthensinI presented 2 distinctive conserved motifs, motif 14 at the N-terminal region and motif 18, both of which were present in 7 out of 10 StAPVIII sequences (Figure 2B).
Additionally, some motifs were shared by atypical and nucellin-like AP groups (Groups from StAPII to StAPVIII) but were absent in Typical StAPI group, such as motifs 6, 10, 12, and 13. Interestingly, motif 10 presented a segment that has been previously reported to be an important conserved sequence GCGYDQ in nucellin-like APs [10], but although the motif showed highly conserved Gly and Cys residues in the majority of StAPs, this six amino acid sequence was present in only two sequences of the nucellin-like, StAP6.1 and StAP6.2 (in StAPIIIb).
Motif 5 was present in StAPs from the majority of groups. Sequence inspection allowed the identification of a conserved Tyr residue of the flap region (Tyr75 of pepsin numbering [10]) in all 62 StAP sequences, although the motif was not present in groups StAPI and StAPII. It has been previously reported that nucellin-like APs had a conserved sequence QCDYE in the vicinity of Y75 [10]. Coincidentally, this sequence was found exclusively in StAPs that clustered with the reference nucellin-like APs in group StAPIIIa.
Finally, some motifs such as 1, 2, 3, and 4 were present in all eight StAPs groups, highlighting the importance of these conserved overall regions in aspartic proteases, while a deeper analysis showed group or subgroup distinctive amino acid sequences. Motif 2 was present at the C-terminal region of most 62 StAPs. A closer inspection of motif 2 confirmed the presence of HTVFD amino acid sequence exclusively in five out of six StAPs from Group StAPI. This short amino acid sequence was reported previously to be exclusive of typical StAPs [10]. Motif 4 comprises the GLIGL sequence reported to be in the loop that forms the active site [3,10].
Motif 1 and 3 contained the DT/SG from the first and second catalytic sites, respectively, and were both present in all 62 StAPs, except StAP5.3, which, although lacking motif 1, presented the DTG triad of the catalytic site.
The surroundings of the active sites of StAPs were distinctive and displayed the three characteristic variations reported previously for typical, nucellin-like, and atypical APs.
The main sequence features surrounding the characteristic first and second catalytic aspartate sequences and the C-terminal end that characterize typical, nucellin-like, and atypical aspartic proteases in potato were detected (as indicated in Table S1). The conserved amino acid signature assigned to typical aspartic proteases was found in StAPI at the first catalytic aspartate (hydrophobic-hydrophobic-hydrophobic (F)-DTG-serine-serine) and in the C-terminal region, which includes the characteristic saposin-like sequence HTVFD (except for StAP9.3). Most typical aspartic proteases presented a DSG sequence in the second catalytic aspartate, which was not present in StAP9.3 and StAP3.1, which had DTG instead. Moreover, the C-terminal end of amino acid sequences in this group confirmed the absence of a Cys residue, which was observed in most sequences of the rest of the groups.
StAPs from cluster “d” and subgroup from cluster StAPIII presented a similar first catalytic aspartate signature (DϕDTGS(D/N/T)) tha was reported in other species for nucellin-like aspartic proteases (DϕDTGS(D/E)) [10,53] (Table S3). While StAPs from StAPII, StAPIIIb, and StAPIV to StAPVIII presented the reported first catalytic aspartate signature for atypical APs (ϕϕDTGS(D/E) [10,12,53], 15 out of 51 atypical StAPs sequences were more diverse (ϕ(ϕ/A)D(T/A/S)GS(D/S/N/E/Q)), which could probably be reflected in their diversity of functions.
2.4. Cysteine Distribution and Glycosylation Sites
Regarding the presence of cysteines in StAPs, the cysteine content was similar within AP groups, as previously reported [10,12]. It was found that 13 cysteines were the most frequent number, with an upper and lower limit of 9 to 17.
Mature typical APs present six conserved Cys residues involved in three disulfide bonds that are crucial to stabilizing the three-dimensional structure of the enzymes [66]. The plant-specific insert (PSI) also has six cys residues that form three disulfide bridges with the main role of conferring rigidity for antimicrobial activity [65]. In agreement, conserved Cys residues were present in all sequences from group StAPI. Most of the group StAPI conserved cysteine residues were present within motifs 1, 8, 11, and 16. Cys residues from motifs 8 and 16 belong to the PSI sequence.
The Cys-rich region located between the first catalytic site and the conserved Tyr75 (of pepsin numbering, [10]) described in atypical and nucellin-like APs was coincident with motifs 5 and 6. Motif 5 was present in StAPs from all groups except StAPI and StAPII and exhibits a conserved Cys residue, while motif 6 is present in all groups except in StAPI and has two conserved Cys of the Cys rich region.
Conserved Cys residues described for atypical and nucellin-like AP were also found in motifs 10, 12, and 13. Motif 17 in groups StAPIV to StAPVIII included a conserved Cys that was present in all 62 StAPs. Although the structures of both typical and nucellin-like proteins have not yet been determined experimentally, molecular modeling predictions for a Nepethesin I protein indicate that the number of disulfide bonds might be similar to those of typical APs [35].
N-linked glycosylation plays an important role in AP biological functions [67]. Typical three-dimensional AP structures have been solved and present two or even three conserved glycosylation sites [66]. StAPI N-linked glycosylation sites were present within the conserved motifs 7, 11, and 15. It was reported that the PSI of populus APs presents a glycosylation site, which was also predicted on StAP in motif 8 [53].
Although the structures of atypical and nucellin-like APs are yet to be available, the presence of a glycosylation site has been determined [53]. Interestingly, opposite to what is observed for the presence of glycosylation sites in typical APs within conserved motifs, an atypical AP glycosylation site was found to be in a non-motif region right before motif 10.
2.5. StAPs Cellular and Subcellular Predicted Location
Primary protein structure was employed to predict cellular and subcellular location of StAPs. Most StAPs were predicted to be localized on vacuole (29 StAPs), plasma membrane (27 StAPs), extracellular space (32 StAPs), Golgi apparatus (21 StAPs), and cytoplasm (14 StAPs), while a few were located in the nucleus (9 StAPs), chloroplast (3 StAPs), peroxisome (2 StAPs), and mitochondria (1 StAPs) (Table S7). About 85% of the 62 StAPs were predicted to be multi-located. StAPs’ multilocation could be reflecting their possibility to act in different tissues and/or stress conditions. One of the functions of aspartic proteases is to process enzymes and includes turnover and enzyme activation [68].
Aspartic proteases from StAPI were preferentially predicted to be located on vacuole as other typical APs reported in barley [69], castor bean (Ricinus communis) [70], and A. thaliana [71]. To a lesser extent, typical APs were also predicted to be in the extracellular space coincident with those reported for tomato (Solanum lycopersicum) [72] and tobacco (Nicotiana tabacum) [73]. Finally, nucellin-like (Cluster StAPIIIa) and atypical aspartic proteases (Clusters StAPIIIb and StAPIV to StAPVIII) had been predicted to be in several cellular and subcellular locations, especially the plasma membrane, extracellular space, and the Golgi apparatus. Several examples of atypical and nucellin-like have widespread distribution in other plant species [17].
2.6. Chromosome Location
StAPs were distributed along all 12 potato chromosomes, the majority of them localized in chromosomes I, II, V, VI, VII, and VIII (Figure 3). Some AP genes were grouped closely, and others were located in tandem arrays as it could be observed with StAP1.1 and StAP1.2; StAP2.5 and StAP2.6; StAP6.5 and StAP 6.6; StAP 7.1 to StAP 7.3; StAP 8.1 to StAP 8.6; StAP 8.8, and StAP 8.9. However, the StAP0.1 chromosome location could not be determined on S. tuberosum group Phureja DM1-3 516 R44 genome (version 4.3 [74]). In order to find its probable location, a blast search on the NCBI database was performed. The results obtained showed that this sequence matched a sequence on chromosome 4 of S. tuberosum cultivar Solyntus (sequence ID: CP055237.1, Id: 4000/4130(97%); gaps: 32/4130) and on the same chromosome in S. pinnelli, and in S. lycopersicum.
Segmental and tandem duplications have played a fundamental role in gene family expansions and functional diversification [75]. Wang et al. (2018) investigated the variation in gene family sizes across species and the likely factors contributing to the variation, using the Solanaceae family as a model and Pfam domain families as a proxy for gene families [76]. They found that genes in high- and low-variability gene families tend to be duplicated by tandem and whole genome duplication, respectively. This could be the case of StAP genes, where gene duplications and tandem arrays were found (Figure 3, Table S2). Tandem array distribution has also been recently described in APs from A. thaliana [77], V. Vinifera [9], and P. trichocarpa [53].
As observed in the phylogenetic tree (Figure 2), clusters grouped highly similar proteins codified by genes located in different chromosomes. In some cases, it could be attributed to segmental duplications (blue lines in Figure 3), as it was confirmed on the potato genome duplication database (Table S2). In all these cases, a negative selection or purifying selection (Ka/Ks<1) was found with a Ka/Ks 0.211 to 0.24. We also tested other potential gene segmental duplications by using the criteria of a) the length of the aligned sequence covers (at least 75% of the longer gene) and b) a similarity of aligned regions equal or superior to 75% [78]. We found potential cases of segmental duplications for some pairs of genes: StAP5.4 and StAP9.4; StAP12.1 and StAP4.4; StAP1.3 and StAP10.1 (green lines in Figure 3), which presented a Ka/Ks<1 from 0.14 to 0.18 and approximately a duplication date of 24.9 MYA [79].
Segmental and tandem duplications observed in StAPs genes (approximately 19 of 62) were lower than what was reported in AP genes from P. trichocarpa (52 out of 62), however, both segmental duplication events were under purifying selection force. It had been proposed that this could mean the retention of their ancestral functions and redundancy, which will also be related to the sequence environment where they were located [16,53]. Thus, the sequence environment plays an important role and could lead to changes in gene expression patterns of the duplicated copies, as could be the cis-acting regulatory elements that would be analyzed in the next section.
2.7. StAP Gene Structure
Structural gene analyses were performed comparing exon–intron organization based on the predicted longest primary transcript structure of the 62 StAP nucleotide sequences. Twenty-seven intronless StAP genes were found, representing almost 44% of the total StAP genes (Figure 4). These were mostly grouped in clusters related to atypical APs, such as StAPV, StAPVI, StAPVII, and StAPVIII. The other four clusters contained StAP genes with numerous introns (up to 13), included mainly in StAPI (related with typical Aps), StAPII (atypical), and StAPIIIa and StAPIIIb (related with nucellin-like and atypical Aps, respectively). A similar intron–exon distribution among categories has been reported in rice [16], populus [53], and A. thaliana [77].
2.8. StAP Expression Analyses
In the last two decades, different studies have reported the participation of APs in plants exposed to adverse environmental conditions such as heat [83,84], osmotic and salt stresses [85], wounding [86], UV-B [13,87,88], and drought [13,21,42,45].
The expression analyses of 62 StAPs were retrieved from the PGSC DM1.3 RNASeq database to characterize StAPs and to identify candidate genes that could be involved in different abiotic stresses.
Twenty-six out of sixty-two genes presented a constitutive expression in different tissues or organs (Figure 4). Some StAP genes revealed preferential expression and tissue-specific patterns. High constitutive expression levels in most plant organs and tissues were associated with clusters StAPI (with StAP1.6, StAP7.5), StAPIV (with StAP8.7), StAPV (with StAP1.4, StAP7.5, StAP11.1), and StAPVII (StAP1.3 and StAP3.5), while the highest transcript level was related to cluster StAPVI (StAP1.4, StAP11.1) on petiole, to cluster StAPVIII (with StAP8.10) on fruits, and to cluster StAPVII (with StAP10.1 and StAP6.6) on mesocarp and endocarp tissues.
There were five genes without expression in any of the reported conditions or tissues (StAP7.1, StAP8.2, StAP8.3, StAP8.4, StAP8.6), which could mean that they might be induced in other conditions such as other plant physiological stages or stress sources. Alternatively, the lack of expression might suggest that they could in fact be pseudogenes. Interestingly, four of these genes were grouped in StAPVIII with the other five genes that were expressed in mesocarp, endocarp, and/or fruit, which points out a tight regulation.
An exception to the above-mentioned behavior of group StAPVIII was StAP5.6 (homolog to Nepenthesin I), which showed expression in several tissues including root, stem, callus, immature fruit, and upregulation after abiotic stress (including salt, mannitol, and temperature treatments) and under hormonal treatment mainly by Gibberellic acid (GA3) and 6-Benzylaminopurine (BAP). This broad expression profile could be related to a role as processing and/or degradative enzyme [8].
Approximately 21% of StAPs modified their gene expression profile under salt, osmotic, or temperature treatments. It was observed that those StAP genes with increased expression after NaCl and mannitol treatments also presented increased transcription profiles after GA3, BAP, and Abscisic Acid (ABA) treatments. The induced genes included StAP1.6, which showed the highest induction of transcript levels under NaCl and mannitol treatments, together with StAP11.1, StAP5.3, and StAP8.11. On the other hand, StAP1.4, StAP12.1, StAP1.3, StAP8.7, StAP4.4, and StAP10.1 were repressed.
Regarding heat treatments, 13 genes presented changes in their expression profile. Most of these genes were repressed except for StAP7.5, StAP2.6, StAP2.1, and StAP8.7, which exhibited different degrees of induction. Interestingly, some AP genes, orthologous to the ones described in this work, were expressed in similar environmental stress conditions. In a transcriptomic study, it was reported that an AtAP gene At1g62290 (a possible ortholog to StAp7.5) was induced in plants overexpressing the heat shock transcription factor HSfA2 [89]. Interestingly, HSfA2 overexpressing plants were tolerant to high temperatures, as well as to both salt and osmotic stress conditions. The authors did not analyze the role of these AP genes, but they concluded that genes that resulted regulated by HSfA2 could be implicated in multiple stress tolerance responses. Expression analysis of the AP gene superfamily in grape (V. vinifera.) identified a total of 19 VvAPs that responded to at least one abiotic stress condition. Specifically, the VvAP21 gene was upregulated in plants treated with NaCl 250 mM during 24 h, in accordance with its possible ortholog, StAP1.6 [9]. In another work, the same group overexpressed VvAP17, obtaining enhanced salt and drought stress-tolerant plants [85]. In detail, overexpressed VvAP17 lines exhibited an incremented level of the ABA hormone, a reduced stomatal aperture size, longer primary roots, as well as higher activities of several antioxidant enzymes such as catalase, superoxide dismutase, and peroxidase. A possible ortholog to VvAP17 is StAP4.3, which also showed upregulation after salt and mannitol treatment. These results suggest that StAPs hormonal expression profiles are related to abiotic stress.
The double monoploid reference genome DM1-3516-R44 is susceptible to P. infestans [50]. The analyses of the biotic stress expression profile revealed at least 16 genes that were affected by P. infestans treatment. Most of them were downregulated and included StAP1.6, StAP7.5, StAP2.2 (from StAPI), StAP8.11, StAP5.2 (from StAPVI), StAP2.6 (from StAPII), StAP11.1, StAP2.3 (from StAPV), StAP9.2 (from StAPIIIb), and StAP8.9 (from StAPIV). Typical aspartic proteases have been implicated in plant defense against pathogens. Thus, StAP3 protein (homolog to StAP1.6) accumulates in intercellular fluids after pathogen infection, and its levels are higher in a field-tolerant cultivar (cv. Pampeana) compared to a susceptible cultivar (cv. Bintje). It has also been shown that StAP3 inhibits the germination of P. infestans and F. solani spores in in vitro experiments [11,48]. The PSI domain StAP and procardosin A has been responsible for membrane interaction and permeabilization at low pH [90,91] and in alipid-composition-dependent manner [65,91].
Five genes grouped in cluster StAPIII had a slight expression increase after P. infestans infection (StAP11.2, StAP6.1, StAP6.2, StAP6.3, StAP3.5), and only StAP8.7 (StAPIV) showed strong upregulation. Some of these genes, StAP8.7 and StAP6.1, were also upregulated after BABA treatment, related to pathogen defense induction [92]. While StAP2.6, StAP9.2, and StAP7.5 presented an increase in their expression profile under BABA treatment, an opposite effect was observed with P. infestans infection and BTH treatment (a salicylic acid analog with SAR response signaling) [93].
Regarding BTH treatment, we observed upregulation in StAP1.6 and StAP2.2 (cluster StAPI), StAP10.1, StAP1.3 (cluster StAPVII), StAP11.2 (cluster StAPIII), and StAP8.9 (cluster StAPIV), while downregulation was observed in at least six genes (StAP8.11, StAP5.2, StAP2.6, StAP7.5, and StAP8.7) for this treatment (Figure 4).
2.9. Cis-Acting Regulatory Elements
Cis-acting elements were identified in silico in promoters of S. tuberosum aspartic proteases (Figure 5, Table S8), several of which correspond to abiotic and biotic stress responsiveness: anaerobic response element (ARE), enhancer element involved in specific anoxic inducibility (GC-motif), low-temperature response element (LTR), “cis-acting element involved in dehydration, low-temp, salt stresses” (DRE), water response elements (Myb, AT-rich element), drought response elements (DRE core, Dre1, as-1, MyC), MYB binding site involved in drought inducibility (MBS) and its drought stress-related MYB recognition site, “heat, osmotic stress, low pH, nutrient starvation stresses response element” (STRE), repetitions rich in TC involved in defense and response to stress (TC-rich repeats), and wound and pathogen response elements (WUN-motif, WRE3, W-box, Box S). In addition to stress-related elements, aspartic proteases presented motifs associated with phytohormone response: elements sensitive to abscisic acid (ABRE, ABRE2, ABRE3a, ABRE4), ABA and GA response element (CARE), cis-acting elements involved gibberellin response (GARE motif, P-box, and TATC box), auxin-responsive elements (AuxRR-core, TGA-box, TGA element), cis-acting regulatory element involved in the MeJA response (CGTCA-motif), ethylene-responsive element (ERE), and cis-acting elements involved in the response to salicylic acid (SARE, TCA, TCA element).
The majority of StAP coding genes presented more than one stress-related and phytohormonal-responsive cis-element in their promoter regions, some of them in very high copy numbers. In general, their presence correlates with the StAPs expression profiles observed in Figure 4.
StAPs that expressed differentially under NaCl and mannitol treatments (such as StAP1.6, StAP11.1, StAP5.3, and StAP8.11, among others) presented several cis-acting elements related to water, drought, salt, and osmotic stresses (AT-rich element, AT-rich sequence, DRE, DRE1, as-1, MBS, Myb, MYB recognition site, MyC, STRE). In addition, StAP1.6, StAP8.11, and StAP11.1 were also induced by ABA and presented cis-acting elements sensitive to abscisic acid, some of them in multiple copies (StAP1.6 has two copies of ABRE in the positive strand and three copies of ABRE, two of ABRE2, and one ABRE3a element in the negative strand; StAP8.11 has one ABRE4 in the positive strand and ABRE and ABRE3a in the negative strand, and StAP11.1 has one ABRE element in the negative strand). StAP1.6, which presented the highest copy number of abscisic acid-responsive elements, exhibited major induction. It has been described that more than one copy of ABRE is necessary for ABA-dependent gene expression [94].
Regarding heat stress response, promoters of StAPs induced by heat (StAP7.5, StAP2.6, StAP2.1, and StAP8.7) presented one or more copies of the heat, osmotic stress, low pH, and nutrient starvation stresses response element (STRE). P. infestans-responsive StAPs presented one, or in some cases two, cis-acting elements involved in defense, wound and pathogen responses (Box S, W box, WRE3, WUN-motif, TC-rich repeats). Salicylic acid-responsive elements (SARE, TCA, TCA-element) were present in three out of six upregulated and in five out of six downregulated genes under BTH treatments.
The presence of cis-acting regulatory elements found in StAP promoter regions is in agreement with previously reported roles of APs in water deficit, osmotic stress, and ABA signaling [45,85].
3. Materials and Methods
3.1. AP Identification and Characterization
A profile of Hidden Markov Model (HMMER 3.1b2, [95]) was created based on 59 aspartic protease amino acid sequences from A. thaliana [10]. These sequences were aligned with Multiple Alignment using Fast Fourier Transform (MAFFT) [96], and five of them were eliminated in order to improve the profile. The resulting alignment was trimmed on both ends, removing regions where the alignment was poor due to the variable nature. The HMMER profile was validated on the Arabidopsis TAIR database [97] (TAIR10_pep_20101214/) in accordance with sequences described [10,77].
Using this profile, an HMMER search was performed [98] against the Potato Genome Sequence Consortium reference genome Double monoploid DM1-3516-R44 [50] peptide database (Version 3.4). The HMMER logo was obtained with Skylign [99].
StAPs were named according to their chromosome localization with a first number that indicates the chromosome and a second number that accounts for its relative position within the chromosome (i.e., the aspartic protease localized on the most upstream position of chromosome 2 is named StAP2.1).
3.2. Phylogenetic Analyses
Phylogenetic analyses were carried out with StAP amino acid sequences identified from S. tuberosum (PGSC DM peptide database, Version 3.4), and 8 reference AP plant amino acid sequences: Nucellin (U87148) [57], Procardosin A (CAB40134) [56], StAP3 (AAT77954) [11], OsAP2 (AAK81699) [27], Nepenthesin I (AB114914) [59], Phytepsin (CAA39602) [55], NP_563851.1 [23], and NP_565219.1 [10]. A second phylogenetic analysis included these sequences and AtAP amino acids sequences (2005) [10].
Alignments and phylogenetic trees were performed with MAFFT [96]. A phylogenetic tree was generated using Simple Phylogeny with the neighbor-joining clustering method with 100 bootstrap replications and default parameters [100].
3.3. Domain Structure Analysis, MEME Motifs, and Subcellular Localization
AP protein domains and active sites were analyzed in Pfam database [101] and Scanprosite tool [102]. The Multiple EM for Motif Elicitation (MEME) online tool v5.2.0 [103] was employed to identify conserved motifs of StAP proteins. Parameter settings were as follows: optimum motif width was set to 6–50, the maximum number of motifs was set to 20, with expectation motif site distribution set to any number of repetitions and all other parameters set by default. SignalP V4.1 was used for identifying putative signal peptides [104]. PredoTar V1.04 [105], TargetP4.0 [106], ProtComp 9.0 [107] and YLOC [108,109] with default settings were used to predict signal sequences to organelles and other subcellular localizations [110].
3.4. Gene Structure Analysis
AP genes were checked for intron and exon structure based on primary transcripts from the PGSC [56].
3.5. Chromosome Localization and Duplication Analyses
A physical map was constructed with potato StAPs gene localization obtained from the potato sequenced genome PGSC database (
For duplication analyses, genome tandem duplications were defined for those StAP genes falling within 50 Kb of one another [53]. While duplicated blocks were downloaded from the Plant Genome Duplication Database (PGDD) [80], the number of Ka and Ks values of non-synonymous and synonymous nucleotide substitutions, respectively, were extracted from PGDD, and Ka/Ks ratio was calculated. The duplication time was estimated according to the formula: T = Ks/2λ, assuming a clock-like rate (λ) was 6.9 × 10 −9 substitutions/site/year for S. tuberosum [111].
3.6. Expression Profile
RNA-Seq expression data from double monoploid S. tuberosum group Phureja DM1-3 516 R44 clone [50,74] were integrated to local software (AgroBiotechTools from Agrobiotechnology Laboratory, INTA UIB-Balcarce) and represented as described in the graphical display. The data retrieved include libraries from a) abiotic stress conditions (salt, manitol and temperature); b) biotic stress (detached leaves inoculated with P. infestans (Pi), DL-β-aminobutyric acid (BABA), or benzoic-(1,2,3)-thiadiazol-7-carbotioic acid S-methyl ester (BTH) at 0 hs (biotic control) or a pool of samples at 24, 48, and 72 h post-inoculation (hpi) (biotic treatment)); c) hormonal treatments (indole acetic acid (IAA), gibberellic acid (GA3), 6-Benzylaminopurine (BAP), abscisic acid (ABA), or water (control)); and d) different tissues and plant organs (root, stem, leaf, stolon, primary and secondary tuber, flower, petiole, petal, carpel, sepal, callus, immature and mature fruits, and mesocarp and endocarp). Expression data were expressed in “Number of fragments per kilobase of exon per million fragments mapped” or FPKM [50,74].
3.7. Analysis of Cis-Acting Elements
Cis-acting elements were detected in the promoters of StAP by in silico analysis using the web software PlantCARE [112,113]. Upstream sequences of AP genes, -3000 bp to +103 bp from the start codon, were retrieved from the PGSC v4.03 database [74].
3.8. Graphical Displays
Protein domain structures, MEME motifs, expression analysis, cis-acting elements, and phylogenetic trees were plotted with iTOL program [81,82].
4. Conclusions
In the present study, we presented the first identification and description of 62 aspartic protease genes on the potato genome. The data obtained provide new knowledge about gene structure, chromosome localization, expression, and the evolution of the aspartic protease family in S. tuberosum, as well as S. tuberosum aspartic protease (StAP) protein organization, phylogenetic relations, and tissue and subcellular predicted localization.
Phylogenetic analysis revealed that StAPs are distributed in eight groups, named from StAPI to StAPVIII, that were differentiated into typical (StAPI), nucellin-like (StAPIIIa), and atypical aspartic proteases (StAPII, StAPIIIb to StAPVIII). Specific and common domains and Multiple EM for Motif Elicitation (MEME) motifs were defined among these groups.
Most StAPs were predicted to be multi-located, most localized on the vacuole (29 StAPs), plasma membrane (27 StAPs), extracellular space (32 StAPs), Golgi apparatus (21 StAPs), and cytoplasm (14 StAPs), while a few were in the nucleus (9 StAPs), chloroplast (3 StAPs), peroxisome (2 StAPs), and mitochondria (1 StAPs).
StAPs were distributed along all 12 potato chromosomes. Segmental and tandem duplications evidenced in the present work, might have contributed to the expansion of the potato AP gene family.
A total of 27 intronless StAP genes were found, representing almost 44% of the total StAP mainly on atypical AP.
The RNAseq analysis revealed that 26 out of 62 genes presented a constitutive expression in different tissues or organs. Some StAPs revealed preferential expression and tissue-specific patterns. Approximately 21% of StAPs modified their gene expression profile under salt, osmotic, or temperature treatments. Gene expression was consistent with the presence of cis-acting regulatory elements on StAP promoter regions.
This study leads to a better understanding of the StAP family and provides the baseline material for future studies tending to determine the functional assignment of the identified genes, as well as the redundancy or diversification of functions in duplicate and tandem genes. Additionally, the knowledge of the diversity of StAP genes and their potential functions can be a source of information for future researchers whose objective is potato breeding through new biotechnology techniques such as gene editing mediated by CRISPR.
Conceptualization, M.G.G. and N.S.N.; methodology, M.F.R.B., M.A.C. and N.S.N.; software, M.A.C., M.F.R.B., C.A.D.O. and N.S.N.; validation, M.F.R.B., N.S.N. and C.A.D.O.; formal analysis, M.F.R.B., N.S.N. and S.D.; investigation, M.F.R.B. and N.S.N.; resources, M.G.G. and S.E.F.; data curation, M.F.R.B. and N.S.N.; writing—original draft preparation, N.S.N. and M.G.G.; writing—review and editing, M.F.R.B., N.S.N. and C.A.D.O.; visualization, M.F.R.B., N.S.N. and G.A.M.; supervision, M.G.G.; project administration, M.G.G. and S.E.F.; funding acquisition, M.G.G. and S.E.F. All authors have read and agreed to the published version of the manuscript.
M.F.R.B. is a fellow of the National Agency for Scientific and Technological Promotion, Argentina (ANPCyT), grant 2348- 2017. M.G.G. and S.D. are established researchers of the National Scientific and Technical Research Council, Argentina (CONICET). C.A.D.O., M.G.G., S.E.F., and N.S.N. are established researchers of the National Institute of Agrotechnology, Argentina (INTA). This work was suuported by, Agencia Nacional de Promoción Científica y Tecnológica, Grant 2348-2017; National Scientific and Technical Research Council, Grant PIP 2020- 2023; INTA Grant PE:I115 and INTA Grant, PE:I116.
Not applicable.
Not applicable.
All data generated or analyzed during this study are included in this published article [and its
The authors declare no conflict of interest.
AP: aspartic protease; PGSC: potato genome sequence consortium, MEME: Multiple EM for Motif Elicitation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 1. Phylogenetic tree of representative aspartic proteases. The phylogenetic tree was built with a total of 62 and 51 AP aminoacid sequences from S. tuberosum and A. thaliana, respectively, and 6 reference sequences from O. sativa L. (OsAP2), S. tuberosum L. (StAP3), A. thaliana L., C. cardunculos L. (Cardosin A), N. gracilis (Nepenthesin I), and H. vulgare L. signed in bold. Cluster names include typical (a), atypical (b, c, e, f, g, and h), and nucellin-like (d) aspartic proteases.
Enlarge this image.Figure 2. Phylogenetic tree of 62 representative aspartic proteases from S. tuberosum. The tree includes 8 reference aspartic proteases from O. sativa L. (OsAP2), S. tuberosum L. (StAP3), A. thaliana L. (NP_565219.1 and NP_563851.1), C. cardunculos. (Procardosin A), N. gracilis (NepenthesinI), and H. vulgare (Phytepsin) indicated in red. Group names and branches are indicated with different colors. (A) Pfam Domains; (B) MEME motifs (specified in Table S6).
![View Image - Figure 3. Chromosome distribution of 62 StAP genes. Potato physical map (in bp) built by Agrobiotech tools based on PGSC DM3.4 sequence [50]. Blue lines connect genes in segmental duplications blocks previously reported [80], and green fine lines represent possible cases of segmental duplication detected.](https://media.proquest.com/media/hms/THUM/6/tjpDM?_a=ChgyMDI1MTAxODE3MTgzNzYwNzo3MDEwNDISBzEzODMzMjMaCk9ORV9TRUFSQ0giDTIxNi43My4yMTYuODcqBzIwMzIzNDcyCjI2MzMyNTkzNzI6C0lubGluZUltYWdlQgEwUgZPbmxpbmVaA0lMSWIEVEhVTWoKMjAyMi8wMS8wMXIKMjAyMi8xMi8zMXoAggEgUC0xMDAwMDg0LTIwMTM5NS1TVEFGRi1udWxsLW51bGySAQZPbmxpbmXKAWdNb3ppbGxhLzUuMCBBcHBsZVdlYktpdC81MzcuMzYgKEtIVE1MLCBsaWtlIEdlY2tvOyBjb21wYXRpYmxlOyBDbGF1ZGVCb3QvMS4wOyArY2xhdWRlYm90QGFudGhyb3BpYy5jb20p0gESU2Nob2xhcmx5IEpvdXJuYWxzmgIHUHJlUGFpZKoCG09TOkVNUy1GdWxsVGV4dC1nZXRGdWxsVGV4dMoCD0FydGljbGV8RmVhdHVyZeICAPICAPoCAVmCAwNXZWKKAxxDSUQ6MjAyNTEwMTgxNzE4Mzc2MDc6MjIyODkw&_s=gNdrgzsx9oY65MRAZYCE2A3coMU%3D "View Image - Figure 3. Chromosome distribution of 62 StAP genes. Potato physical map (in bp) built by Agrobiotech tools based on PGSC DM3.4 sequence [50]. Blue lines connect genes in segmental duplications blocks previously reported [80], and green fine lines represent possible cases of segmental duplication detected.")
Enlarge this image.Figure 3. Chromosome distribution of 62 StAP genes. Potato physical map (in bp) built by Agrobiotech tools based on PGSC DM3.4 sequence [50]. Blue lines connect genes in segmental duplications blocks previously reported [80], and green fine lines represent possible cases of segmental duplication detected.
![View Image - Figure 4. Aspartic proteases exonic and intronic structure and gene expression heatmap. The phylogenetic tree was built with 62 StAP amino acid sequences. (a) Exon and intron structures were built based on the predicted primary longest transcript structure, and (b) RNA-Seq data from different libraries from double monoploid S. tuberosum group Phureja DM1-3 516 R44 clone (DM) were indicated as “number of fragments per kilobase of exon per million fragments mapped” (FPKM) by color scale [50]. The phylogram was built with the iTol program [81,82]. Names of genes are colored based on the group membership in Figure 1.](https://media.proquest.com/media/hms/THUM/8/tjpDM?_a=ChgyMDI1MTAxODE3MTgzNzYwNzo3MDEwNDISBzEzODMzMjMaCk9ORV9TRUFSQ0giDTIxNi43My4yMTYuODcqBzIwMzIzNDcyCjI2MzMyNTkzNzI6C0lubGluZUltYWdlQgEwUgZPbmxpbmVaA0lMSWIEVEhVTWoKMjAyMi8wMS8wMXIKMjAyMi8xMi8zMXoAggEgUC0xMDAwMDg0LTIwMTM5NS1TVEFGRi1udWxsLW51bGySAQZPbmxpbmXKAWdNb3ppbGxhLzUuMCBBcHBsZVdlYktpdC81MzcuMzYgKEtIVE1MLCBsaWtlIEdlY2tvOyBjb21wYXRpYmxlOyBDbGF1ZGVCb3QvMS4wOyArY2xhdWRlYm90QGFudGhyb3BpYy5jb20p0gESU2Nob2xhcmx5IEpvdXJuYWxzmgIHUHJlUGFpZKoCG09TOkVNUy1GdWxsVGV4dC1nZXRGdWxsVGV4dMoCD0FydGljbGV8RmVhdHVyZeICAPICAPoCAVmCAwNXZWKKAxxDSUQ6MjAyNTEwMTgxNzE4Mzc2MDc6MjIyODkw&_s=WOx98sPVYNSoWsDZnM6%2BE%2FmtH%2FM%3D "View Image - Figure 4. Aspartic proteases exonic and intronic structure and gene expression heatmap. The phylogenetic tree was built with 62 StAP amino acid sequences. (a) Exon and intron structures were built based on the predicted primary longest transcript structure, and (b) RNA-Seq data from different libraries from double monoploid S. tuberosum group Phureja DM1-3 516 R44 clone (DM) were indicated as “number of fragments per kilobase of exon per million fragments mapped” (FPKM) by color scale [50]. The phylogram was built with the iTol program [81,82]. Names of genes are colored based on the group membership in Figure 1.")
Enlarge this image.Figure 4. Aspartic proteases exonic and intronic structure and gene expression heatmap. The phylogenetic tree was built with 62 StAP amino acid sequences. (a) Exon and intron structures were built based on the predicted primary longest transcript structure, and (b) RNA-Seq data from different libraries from double monoploid S. tuberosum group Phureja DM1-3 516 R44 clone (DM) were indicated as “number of fragments per kilobase of exon per million fragments mapped” (FPKM) by color scale [50]. The phylogram was built with the iTol program [81,82]. Names of genes are colored based on the group membership in Figure 1.
Enlarge this image.Figure 5. In silico identification of cis-acting elements in promoters of StAP. Colored boxes represent each element position in the coding strand.
Supplementary Materials
The following are available online at
References
1. Cruz, R.; Huesgen, P.; Riley, S.P.; Wlodawer, A.; Faro, C.; Overall, C.M.; Martinez, J.J.; Simões, I. RC1339/APRc from Rickettsia conorii Is a Novel Aspartic Protease with Properties of Retropepsin-Like Enzymes. PLoS Pathog.; 2014; 10, 13. [DOI: https://dx.doi.org/10.1371/journal.ppat.1004324]
2. Davies, D.R. The Structure and Function of the Aspartic Proteinases. Annu. Rev. Biophys. Biophys. Chem.; 1990; 19, pp. 189-215. [DOI: https://dx.doi.org/10.1146/annurev.bb.19.060190.001201]
3. Dunn, B.M. Structure and mechanism of the pepsin-like family of aspartic peptidases. Chem. Rev.; 2002; 102, pp. 4431-4458. [DOI: https://dx.doi.org/10.1021/cr010167q] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12475196]
4. Rawlings, N.D.; Barrett, A. Introduction: Peptidases of Unknown Catalytic Type. Handb. Proteolytic Enzym.; 2013; 3, pp. 3747-3749.
5. Domsalla, A.; Melzig, M.F. Occurrence and properties of proteases in plant latices. Planta Med.; 2008; 74, pp. 699-711. [DOI: https://dx.doi.org/10.1055/s-2008-1074530]
6. Siezen, R.J.; Leunissen, J.A.M. Subtilases: The superfamily of subtilisin-like serine proteases. Protein Sci.; 1997; 6, pp. 501-523. [DOI: https://dx.doi.org/10.1002/pro.5560060301]
7. Rawlings, N.D.; Barrett, A.J.; Thomas, P.D.; Huang, X.; Bateman, A.; Finn, R.D. The MEROPS database of proteolytic enzymes, their substrates and inhibitors in 2017 and a comparison with peptidases in the PANTHER database. Nucleic Acids Res.; 2018; 46, [DOI: https://dx.doi.org/10.1093/nar/gkx1134]
8. Simões, I.; Faro, C. Structure and function of plant aspartic proteinases. Eur. J. Biochem.; 2004; 271, pp. 2067-2075. [DOI: https://dx.doi.org/10.1111/j.1432-1033.2004.04136.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15153096]
9. Guo, R.; Xu, X.; Carole, B.; Li, X.; Gao, M.; Zheng, Y.; Wang, X. Genome-wide identification, evolutionary and expression analysis of the aspartic protease gene superfamily in grape. BMC Genom.; 2013; 14, 1. [DOI: https://dx.doi.org/10.1186/1471-2164-14-554] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23945092]
10. Faro, C.; Gal, S. Aspartic proteinase content of the Arabidopsis genome. Curr. Protein Pept. Sci.; 2005; 6, pp. 493-500. [DOI: https://dx.doi.org/10.2174/138920305774933268] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16381599]
11. Guevara, M.G.; Almeida, C.; Mendieta, J.R.; Faro, C.J.; Veríssimo, P.; Pires, E.V.; Daleo, G.R. Molecular cloning of a potato leaf cDNA encoding an aspartic protease (StAsp) and its expression after P. infestans infection. Plant Physiol. Biochem.; 2005; 43, pp. 882-889. [DOI: https://dx.doi.org/10.1016/j.plaphy.2005.07.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16198590]
12. Soares, A.; Niedermaier, S.; Faro, R.; Loos, A.; Manadas, B.; Faro, C.; Huesgen, P.F.; Cheung, A.Y.; Simões, I. An atypical aspartic protease modulates lateral root development in Arabidopsis thaliana. J. Exp. Bot.; 2019; 70, pp. 2157-2171. [DOI: https://dx.doi.org/10.1093/jxb/erz059] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30778561]
13. Timotijević, G.S.; Milisavljević, M.D.; Radović, S.R.; Konstantinović, M.M.; Maksimović, V.R. Ubiquitous aspartic proteinase as an actor in the stress response in buckwheat. J. Plant. Physiol.; 2010; 167, pp. 61-68. [DOI: https://dx.doi.org/10.1016/j.jplph.2009.06.017] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19643510]
14. Terauchi, K.; Asakura, T.; Ueda, H.; Tamura, T.; Tamura, K.; Matsumoto, I.; Misaka, T.; Hara-Nishimura, I.; Abe, K. Plant-specific insertions in the soybean aspartic proteinases, soyAP1 and soyAP2, perform different functions of vacuolar targeting. J. Plant Physiol.; 2006; 163, pp. 856-862. [DOI: https://dx.doi.org/10.1016/j.jplph.2005.08.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16777533]
15. Castanheira, P.; Samyn, B.; Sergeant, K.; Clemente, J.C.; Dunn, B.M.; Pires, E.; Van Beeumen, J.; Faro, C. Activation, proteolytic processing, and peptide specificity of recombinant cardosin A. J. Biol. Chem.; 2005; 280, pp. 13047-13054. [DOI: https://dx.doi.org/10.1074/jbc.M412076200]
16. Chen, J.; Ouyang, Y.; Wang, L.; Xie, W.; Zhang, Q. Aspartic proteases gene family in rice: Gene structure and expression, predicted protein features and phylogenetic relation. Gene; 2009; 442, pp. 108-118. [DOI: https://dx.doi.org/10.1016/j.gene.2009.04.021]
17. Figueiredo, J.; Sousa Silva, M.; Figueiredo, A. Subtilisin-like proteases in plant defence: The past, the present and beyond. Mol. Plant Pathol.; 2018; 19, pp. 1017-1028. [DOI: https://dx.doi.org/10.1111/mpp.12567]
18. Xia, Y.; Suzuki, H.; Borevitz, J.; Blount, J.; Guo, Z.; Patel, K.; Dixon, R.A.; Lamb, C. An extracellular aspartic protease functions in Arabidopsis disease resistance signaling. EMBO J.; 2004; 23, pp. 980-988. [DOI: https://dx.doi.org/10.1038/sj.emboj.7600086]
19. Gao, L.J.; Feng, J.Y.; Jin, B.; Zhang, Q.N.; Liu, T.Q.; Lun, Y.Q.; Wu, Z.J. Carbazole and hydroxy groups-tagged poly(aspartic acid) scale inhibitor for cooling water systems. Chem. Lett.; 2011; 40, pp. 1392-1394. [DOI: https://dx.doi.org/10.1246/cl.2011.1392]
20. Phan, H.A.; Iacuone, S.; Li, S.F.; Parish, R.W. The MYB80 transcription factor is required for pollen development and the regulation of tapetal programmed cell death in Arabidopsis thaliana. Plant Cell; 2011; 23, pp. 2209-2224. [DOI: https://dx.doi.org/10.1105/tpc.110.082651]
21. Yao, X.; Xiong, W.; Ye, T.; Wu, Y. Overexpression of the aspartic protease ASPG1 gene confers drought avoidance in Arabidopsis. J. Exp. Bot.; 2012; 63, pp. 2579-2593. [DOI: https://dx.doi.org/10.1093/jxb/err433] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22268147]
22. Paparelli, E.; Gonzali, S.; Parlanti, S.; Novi, G.; Giorgi, F.M.; Licausi, F.; Kosmacz, M.; Feil, R.; Lunn, J.E.; Brust, H. et al. Misexpression of a chloroplast aspartyl protease leads to severe growth defects and alters carbohydrate metabolism in Arabidopsis. Plant Physiol.; 2012; 160, pp. 1237-1250. [DOI: https://dx.doi.org/10.1104/pp.112.204016] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22987884]
23. Breitenbach, H.H.; Wenig, M.; Wittek, F.; Jordá, L.; Maldonado-Alconada, A.M.; Sarioglu, H.; Colby, T.; Knappe, C.; Bichlmeier, M.; Pabst, E. et al. Contrasting roles of the apoplastic aspartyl protease APOPLASTIC, ENHANCED DISEASE SUSCEPTIBILITY1-DEPENDENT1 and LEGUME LECTIN-LIKE PROTEIN1 in Arabidopsis systemic acquired resistance. Plant Physiol.; 2014; 165, pp. 791-809. [DOI: https://dx.doi.org/10.1104/pp.114.239665] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24755512]
24. Gao, H.; Zhang, Y.; Wang, W.; Zhao, K.; Liu, C.; Bai, L.; Li, R.; Guo, Y. Two membrane-anchored aspartic proteases contribute to pollen and ovule development [open]. Plant Physiol.; 2017; 173, pp. 219-239. [DOI: https://dx.doi.org/10.1104/pp.16.01719] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27872247]
25. Li, Y.; Kabbage, M.; Liu, W.; Dickman, M.B. Aspartyl protease-mediated cleavage of BAG6 is necessary for autophagy and fungal resistance in plants. Plant Cell; 2016; 28, pp. 233-247. [DOI: https://dx.doi.org/10.1105/tpc.15.00626] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26739014]
26. Prasad, B.D.; Creissen, G.; Lamb, C.; Chattoo, B.B. Overexpression of rice (Oryza sativa L.) OsCDRI leads to constitutive activation of defense responses in rice and Arabidopsis. Mol. Plant-Microbe Interact.; 2009; 22, pp. 1635-1644. [DOI: https://dx.doi.org/10.1094/MPMI-22-12-1635] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19888828]
27. Bi, X.; Khush, G.S.; Bennett, J. The Rice Nucellin Gene Ortholog OsAsp1 Encodes an Active Aspartic Protease without a Plant-specific Insert and is Strongly Expressed in Early Embryo. Plant Cell Physiol.; 2005; 46, pp. 87-98. [DOI: https://dx.doi.org/10.1093/pcp/pci002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15659452]
28. Niu, N.; Liang, W.; Yang, X.; Jin, W.; Wilson, Z.A.; Hu, J.; Zhang, D. EAT1 promotes tapetal cell death by regulating aspartic proteases during male reproductive development in rice. Nat. Commun.; 2013; 4, 11. [DOI: https://dx.doi.org/10.1038/ncomms2396]
29. Huang, J.; Zhao, X.; Cheng, K.; Jiang, Y.; Ouyang, Y.; Xu, C.; Li, X.; Xiao, J.; Zhang, Q. OsAP65, a rice aspartic protease, is essential for male fertility and plays a role in pollen germination and pollen tube growth. J. Exp. Bot.; 2013; 64, pp. 3351-3360. [DOI: https://dx.doi.org/10.1093/jxb/ert173]
30. Yamaya-Ito, H.; Shimoda, Y.; Hakoyama, T.; Sato, S.; Kaneko, T.; Hossain, M.S.; Shibata, S.; Kawaguchi, M.; Hayashi, M.; Kouchi, H. et al. Loss-of-function of ASPARTIC PEPTIDASE NODULE-INDUCED 1 (APN1) in Lotus japonicus restricts efficient nitrogen-fixing symbiosis with specific Mesorhizobium loti strains. Plant J.; 2018; 93, pp. 5-16. [DOI: https://dx.doi.org/10.1111/tpj.13759]
31. Milisavljevic, M.D.; Timotijevic, G.S.; Radovic, S.R.; Konstantinovic, M.M.; Maksimovic, V.R. Two types of aspartic proteinases from buckwheat seed-Gene structure and expression analysis. J. Plant Physiol.; 2008; 165, pp. 983-990. [DOI: https://dx.doi.org/10.1016/j.jplph.2007.03.016]
32. Milisavljević, M.D.; Timotijević, G.S.; Nikolić, D.B.; Samardžić, J.T.; Maksimović, V.R. Cell wall localization of the aspartic proteinase from buckwheat (FeAPL1) over-expressed in tobacco BY-2 cells. J. Serbian Chem. Soc.; 2011; 76, pp. 1229-1236. [DOI: https://dx.doi.org/10.2298/JSC110120106M]
33. Gong, L.; Zhang, H.; Gan, X.; Zhang, L.; Chen, Y.; Nie, F.; Shi, L.; Li, M.; Guo, Z.; Zhang, G. et al. Transcriptome profiling of the potato (Solanum tuberosum L.) plant under drought stress and water-stimulus conditions. PLoS ONE; 2015; 10, pp. 1-20. [DOI: https://dx.doi.org/10.1371/journal.pone.0128041]
34. Nakano, T.; Murakami, S.; Shoji, T.; Yoshida, S.; Yamada, Y.; Sato, F. A novel protein with DNA binding activity from tobacco chloroplast nucleoids. Plant Cell; 1997; 9, pp. 1673-1682. [DOI: https://dx.doi.org/10.1105/tpc.9.9.1673]
35. Athauda, S.B.P.; Matsumoto, K.; Rajapakshe, S.; Kuribayashi, M.; Kojima, M.; Kubomura-Yoshida, N.; Iwamatsu, A.; Shibata, C.; Inoue, H.; Takahashi, K. Enzymic and structural characterization of nepenthesin, a unique member of a novel subfamily of aspartic proteinases. Biochem. J.; 2004; 381, pp. 295-306. [DOI: https://dx.doi.org/10.1042/BJ20031575]
36. Butts, C.T.; Bierma, J.C.; Martin, R.W. Novel proteases from the genome of the carnivorous plant Drosera capensis: Structural prediction and comparative analysis. Proteins Struct. Funct. Bioinforma.; 2016; 84, pp. 1517-1533. [DOI: https://dx.doi.org/10.1002/prot.25095]
37. Schulze, W.X.; Sanggaard, K.W.; Kreuzer, I.; Knudsen, A.D.; Bemm, F.; Thøgersen, I.B.; Bräutigam, A.; Thomsen, L.R.; Schliesky, S.; Dyrlund, T.F. et al. The protein composition of the digestive fluid from the venus flytrap sheds light on prey digestion mechanisms. Mol. Cell. Proteom.; 2012; 11, pp. 1306-1319. [DOI: https://dx.doi.org/10.1074/mcp.M112.021006]
38. Olivares, J.E.; Díaz-Camino, C.; Estrada-Navarrete, G.; Alvarado-Affantranger, X.; Rodríguez-Kessler, M.; Zamudio, F.Z.; Olamendi-Portugal, T.; Márquez, Y.; Servín, L.E.; Sánchez, F. Nodulin 41, a novel late nodulin of common bean with peptidase activity. BMC Plant Biol.; 2011; 11, 13. [DOI: https://dx.doi.org/10.1186/1471-2229-11-134] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21985276]
39. Guo, R.; Tu, M.; Wang, X.; Zhao, J.; Wan, R.; Li, Z.; Wang, Y.; Wang, X. Ectopic expression of a grape aspartic protease gene, AP13, in Arabidopsis thaliana improves resistance to powdery mildew but increases susceptibility to Botrytis cinerea. Plant Sci.; 2016; 248, pp. 17-27. [DOI: https://dx.doi.org/10.1016/j.plantsci.2016.04.006]
40. D’Ippólito, S.; Rey-Burusco, M.F.; Feingold, S.E.; Guevara, M.G. Role of proteases in the response of plants to drought. Plant Physiol. Biochem.; 2021; 168, pp. 1-9. [DOI: https://dx.doi.org/10.1016/j.plaphy.2021.09.038] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34607206]
41. Contour-Ansel, D.; Torres-Franklin, M.L.; Zuily-Fodil, Y.; de Carvalho, M.H.C. An aspartic acid protease from common bean is expressed “on call” during water stress and early recovery. J. Plant Physiol.; 2010; 167, pp. 1606-1612. [DOI: https://dx.doi.org/10.1016/j.jplph.2010.06.018]
42. De Carvalho, M.H.C.; D’Arcy-Lameta, A.; Roy-Macauley, H.; Gareil, M.; El Maarouf, H.; Pham-Thi, A.T.; Zuily-Fodil, Y. Aspartic protease in leaves of common bean (Phaseolus vulgaris L.) and cowpea (Vigna unguiculata L. Walp): Enzymatic activity, gene expression and relation to drought susceptibility. FEBS Lett.; 2001; 492, pp. 242-246. [DOI: https://dx.doi.org/10.1016/S0014-5793(01)02259-1]
43. Sharma, S.K.; Bolser, D.; de Boer, J.; Sønderkær, M.; Amoros, W.; Carboni, M.F.; D’Ambrosio, J.M.; de la Cruz, G.; Di Genova, A.; Douches, D.S. et al. Construction of reference chromosome-scale pseudomolecules for potato: Integrating the potato genome with genetic and physical maps. G3 Genes Genomes Genet.; 2013; 3, pp. 2031-2047. [DOI: https://dx.doi.org/10.1534/g3.113.007153]
44. Nallamothu, J.; Motukuri, S.R.K.; Uppluri, L.S.; Ponnapati, S. Comparative expression profiling reveals the role of the hot pepper apoplast under drought conditions. Biologia; 2021; 76, pp. 375-380. [DOI: https://dx.doi.org/10.2478/s11756-020-00579-8]
45. D’´Ippolito, S.; Fiol, D.F.; Dale, G.R.; Guevara, M.G. Overexpression of Arabidopsis aspartic protease APA1 gene confers drought tolerance. Plant Sci.; 2020; 292, 110406. [DOI: https://dx.doi.org/10.1016/j.plantsci.2020.110406]
46. Guevara, M.G.; Oliva, C.R.; Machinandiarena, M.; Daleo, G.R. Purification and properties of an aspartic protease from potato tuber that is inhibited by a basic chitinase. Physiol. Plant.; 1999; 106, pp. 164-169. [DOI: https://dx.doi.org/10.1034/j.1399-3054.1999.106203.x]
47. Guevara, M.G.; Daleo, G.R.; Oliva, C.R. Purification and characterization of an aspartic protease from potato leaves. Physiol. Plant.; 2001; 112, pp. 321-326. [DOI: https://dx.doi.org/10.1034/j.1399-3054.2001.1120304.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11473688]
48. Guevara, M.G.; Oliva, C.R.; Huarte, M.; Daleo, G.R. An aspartic protease with antimicrobial activity is induced after infection and wounding in intercellular fluids of potato tubers. Eur. J. Plant Pathol.; 2002; 108, pp. 131-137. [DOI: https://dx.doi.org/10.1023/A:1015049629736]
49. Guevara, M.G.; Veríssimo, P.; Pires, E.; Faro, C.; Daleo, G.R. Potato aspartic proteases: Induction, antimicrobial activity and substrate specificity. J. Plant Pathol.; 2004; 86, pp. 233-238.
50. Xu, X.; Pan, S.; Cheng, S.; Zhang, B.; Mu, D.; Ni, P.; Zhang, G.; Yang, S.; Li, R.; Wang, J. et al. Genome sequence and analysis of the tuber crop potato. Nature; 2011; 475, pp. 189-195. [DOI: https://dx.doi.org/10.1038/NATURE10158] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21743474]
51. Hawkes, J. The Potato: Evolution, Biodiversity and Genetic Resources; Belhaven Press: London, UK, 1990; 259p. ISBN 1-85293-054-4
52. Takahashi, K.; Niwa, H.; Yokota, N.; Kubota, K.; Inoue, H. Widespread tissue expression of nepenthesin-like aspartic protease genes in Arabidopsis thaliana. Plant Physiol. Biochem.; 2008; 46, pp. 724-729. [DOI: https://dx.doi.org/10.1016/j.plaphy.2008.04.007]
53. Cao, S.; Guo, M.; Wang, C.; Xu, W.; Shi, T.; Tong, G.; Zhen, C.; Cheng, H.; Yang, C.; Elsheery, N.I. et al. Genome-wide characterization of aspartic protease (AP) gene family in Populus trichocarpa and identification of the potential PtAPs involved in wood formation. BMC Plant Biol.; 2019; 19, pp. 1-17. [DOI: https://dx.doi.org/10.1186/s12870-019-1865-0]
54. Gruis, D.; Selinger, D.A.; Curran, J.M.; Jung, R. Redundant proteolytic mechanisms process seed storage proteins in the absence of seed-type members of the vacuolar processing enzyme family of cysteine proteases. Plant Cell; 2002; 14, pp. 2863-2882. [DOI: https://dx.doi.org/10.1105/tpc.005009]
55. Runeberg-Roos, P.; Tormakangas, K.; Ostman, A. Primary structure of a barley-grain aspartic proteinase. A plant aspartic proteinase resembling mammalian cathepsin D. J. Biol. Inorg. Chem.; 1991; 202, pp. 1021-1027. [DOI: https://dx.doi.org/10.1111/j.1432-1033.1991.tb16465.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/1722454]
56. Faro, C.; Ramalho-Santos, M.; Vieira, M.; Mendes, A.; Simões, I.; Andrade, R.; Veríssimo, P.; Lin, X.L.; Tang, J.; Pires, E. Cloning and characterization of cDNA encoding cardosin A, an RGD-containing plant aspartic proteinase. J. Biol. Chem.; 1999; 274, pp. 28724-28729. [DOI: https://dx.doi.org/10.1074/jbc.274.40.28724] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10497243]
57. Chen, F.; Foolad, M.R. Molecular organization of a gene in barley which encodes a protein similar to aspartic protease and its specific expression in nucellar cells during degeneration. Plant Mol. Biol.; 1997; 35, pp. 821-831. [DOI: https://dx.doi.org/10.1023/A:1005833207707] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9426602]
58. Li, S.; Wang, N.; Ji, D.; Xue, Z.; Yu, Y.; Jiang, Y.; Liu, J.; Liu, Z.; Xiang, F. Evolutionary and functional analysis of membrane-bound NAC transcription factor genes in soybean. Plant Physiol.; 2016; 172, pp. 1804-1820. [DOI: https://dx.doi.org/10.1104/pp.16.01132]
59. Takahashi, K.; Athauda, S.; Matsumoto, K.; Rajapakshe, S.; Kuribayashi, M.; Kojima, M.; Kubomura-Yoshida, N.; Iwamatsu, A.; Shibata, C.; Inoue, H. Nepenthesin, a Unique Member of a Novel Subfamily of Aspartic Proteinases: Enzymatic and Structural Characteristics. Curr. Protein Pept. Sci.; 2005; 6, pp. 513-525. [DOI: https://dx.doi.org/10.2174/138920305774933259]
60. Wang, J.F.; Ho, F.I.; Chen, Y.Y.; Lin, Y.M.; Cheng, C.P. A tobacco rattle virus-induced gene silencing system for a soil-borne vascular pathogen Ralstonia solanacearum. Bot. Stud.; 2009; 50, pp. 413-424.
61. Fierens, K.; Brijs, K.; Courtin, C.M.; Gebruers, K.; Goesaert, H.; Raedschelders, G.; Robben, J.; Van Campenhout, S.; Volckaert, G.; Delcour, J.A. Molecular identification of wheat endoxylanase inhibitor TAXI-I 1, member of a new class of plant proteins. FEBS Lett.; 2003; 540, pp. 259-263. [DOI: https://dx.doi.org/10.1016/S0014-5793(03)00276-X]
62. Pollet, A.; Sansen, S.; Raedschelders, G.; Gebruers, K.; Rabijns, A.; Delcour, J.A.; Courtin, C.M. Identification of structural determinants for inhibition strength and specificity of wheat xylanase inhibitors TAXI-IA and TAXI-IIA. FEBS J.; 2009; 276, pp. 3916-3927. [DOI: https://dx.doi.org/10.1111/j.1742-4658.2009.07105.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19769747]
63. Sansen, S.; De Ranter, C.J.; Gebruers, K.; Brijs, K.; Courtin, C.M.; Delcour, J.A.; Rabijns, A. Structural basis for inhibition of Aspergillus niger xylanase by Triticum aestivum xylanase inhibitor-I. J. Biol. Chem.; 2004; 279, pp. 36022-36028. [DOI: https://dx.doi.org/10.1074/jbc.M404212200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15166216]
64. Bekalu, Z.E.; Dionisio, G.; Brinch-Pedersen, H. Molecular properties and new potentials of plant nepenthesins. Plants; 2020; 9, 13. [DOI: https://dx.doi.org/10.3390/plants9050570] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32365700]
65. Bryksa, B.C.; Bhaumik, P.; Magracheva, E.; De Moura, D.C.; Kurylowicz, M.; Zdanov, A.; Dutcher, J.R.; Wlodawer, A.; Yada, R.Y. Structure and mechanism of the saposin-like domain of a plant aspartic protease. J. Biol. Chem.; 2011; 286, pp. 28265-28275. [DOI: https://dx.doi.org/10.1074/jbc.M111.252619]
66. Almeida, C.M.; Manso, J.A.; Figueiredo, A.C.; Antunes, L.; Cruz, R.; Manadas, B.; Bur, D.; Pereira, P.J.B.; Faro, C.; Simões, I. Functional and structural characterization of synthetic cardosin B-derived rennet. Appl. Microbiol. Biotechnol.; 2017; 101, [DOI: https://dx.doi.org/10.1007/s00253-017-8445-8]
67. Pagano, M.R.; Mendieta, J.R.; Muñoz, F.F.; Daleo, G.R.; Guevara, M.G. Roles of glycosylation on the antifungal activity and apoplast accumulation of StAPs (Solanum tuberosum aspartic proteases). Int. J. Biol. Macromol.; 2007; 41, pp. 512-520. [DOI: https://dx.doi.org/10.1016/j.ijbiomac.2007.07.009]
68. Simões, I.; Faro, R.; Bur, D.; Faro, C. Characterization of recombinant CDR1, an Arabidopsis aspartic proteinase involved in disease resistance. J. Biol. Chem.; 2007; 282, pp. 31358-31365. [DOI: https://dx.doi.org/10.1074/jbc.M702477200]
69. Runeberg-Roos, P.; Kervinen, J.; Kovaleva, V.; Raikhel, N.V.; Gal, S. The Aspartic Proteinase of Barley Is a Vacuolar Enzyme That Processes Probarley Lectin In vitro. Plant Physiol.; 1991; 105, pp. 321-329. [DOI: https://dx.doi.org/10.1104/pp.105.1.321]
70. Hiraiwa, N.; Kondo, M.; Nishimura, M.; Hara-Nishimura, I. An aspartic endopeptidase is involved in the breakdown of propeptides of storage proteins in protein-storage vacuoles of plants. Eur. J. Biochem.; 1997; 246, pp. 133-141. [DOI: https://dx.doi.org/10.1111/j.1432-1033.1997.00133.x]
71. Mutlu, A.; Chen, X.; Reddy, S.M.; Gal, S. The aspartic proteinase is expressed in Arabidopsis thaliana seeds and localized in the protein bodies. Seed Sci. Res.; 1999; 9, pp. 75-84. [DOI: https://dx.doi.org/10.1017/S0960258599000082]
72. Rodrigo, I.; Vera, P.; Conejero, V. Degradation of tomato pathogenesis-related proteins by an endogenous 37-kDa aspartyl endoproteinase. Eur. J. Biochem.; 1989; 184, pp. 663-669. [DOI: https://dx.doi.org/10.1111/j.1432-1033.1989.tb15064.x]
73. Rodrigo, I.; Vera, P.; Van Loon, L.C.; Conejero, V. Degradation of tobacco pathogenesis-related proteins: Evidence for conserved mechanisms of degradation of pathogenesis-related proteins in plants. Plant Physiol.; 1991; 95, pp. 616-622. [DOI: https://dx.doi.org/10.1104/pp.95.2.616] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16668027]
74. Spud Data Base. Available online: http://spuddb.uga.edu/ (accessed on 25 October 2020).
75. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol.; 2004; 4, pp. 1-21. [DOI: https://dx.doi.org/10.1186/1471-2229-4-10] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15171794]
76. Wang, P.; Moore, B.M.; Panchy, N.L.; Meng, F.; Lehti-Shiu, M.D.; Shiu, S.H. Factors influencing gene family size variation among related species in a plant family, solanaceae. Genome Biol. Evol.; 2018; 10, pp. 2596-2613. [DOI: https://dx.doi.org/10.1093/gbe/evy193] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30239695]
77. Beers, E.P.; Jones, A.M.; Dickerman, A.W. The S8 serine, C1A cysteine and A1 aspartic protease families in Arabidopsis. Phytochemistry; 2004; 65, pp. 43-58. [DOI: https://dx.doi.org/10.1016/j.phytochem.2003.09.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14697270]
78. Gu, X.; Huang, W. Testing the Parsimony Test of Genome Duplications: A Counterexample. Genome Res.; 2002; 12, pp. 1-2. [DOI: https://dx.doi.org/10.1101/gr.214402] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11779825]
79. Zhang, Z.; Li, J.; Zhao, X.Q.; Wang, J. KaKs_Calculator: Calculating Ka and Ks through model selection and model averaging. Genom. Proteom. Bioinform.; 2006; 4, pp. 259-263. [DOI: https://dx.doi.org/10.1016/S1672-0229(07)60007-2]
80. Lee, T.H.; Tang, H.; Wang, X.; Paterson, A.H. PGDD: A database of gene and genome duplication in plants. Nucleic Acids Res.; 2013; 41, pp. 1152-1158. [DOI: https://dx.doi.org/10.1093/nar/gks1104]
81. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL): An online tool for phylogenetic tree display and annotation. Bioinformatics; 2007; 23, pp. 127-128. [DOI: https://dx.doi.org/10.1093/bioinformatics/btl529]
82. Letunic, I.; Bork, P. Interactive Tree of Life v2: Online annotation and display of phylogenetic trees made easy. Nucleic Acids Res.; 2011; 39, pp. 475-478. [DOI: https://dx.doi.org/10.1093/nar/gkr201]
83. Szucs, A.; Jäger, K.; Jurca, M.E.; Fábián, A.; Bottka, S.; Zvara, Á.; Barnabás, B.; Fehér, A. Histological and microarray analysis of the direct effect of water shortage alone or combined with heat on early grain development in wheat (Triticum aestivum). Physiol. Plant.; 2010; 140, pp. 174-188. [DOI: https://dx.doi.org/10.1111/j.1399-3054.2010.01394.x]
84. Echevarría-Zomeño, S.; Fernández-Calvino, L.; Castro-Sanz, A.B.; López, J.A.; Vázquez, J.; Castellano, M.M. Dissecting the proteome dynamics of the early heat stress response leading to plant survival or death in Arabidopsis. Plant Cell Environ.; 2016; 39, pp. 1264-1278. [DOI: https://dx.doi.org/10.1111/pce.12664] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26580143]
85. Guo, R.; Zhao, J.; Wang, X.; Guo, C.; Li, Z.; Wang, Y.; Wang, X. Constitutive expression of a grape aspartic protease gene in transgenic Arabidopsis confers osmotic stress tolerance. Plant Cell Tissue Organ Cult.; 2014; 121, pp. 275-287. [DOI: https://dx.doi.org/10.1007/s11240-014-0699-6]
86. Schaller, A.; Ryan, C.A. Molecular cloning of a tomato leaf cDNA encoding an aspartic protease, a systemic wound response protein. Plant Mol. Biol.; 1996; 31, pp. 1073-1077. [DOI: https://dx.doi.org/10.1007/BF00040725] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8843949]
87. Trentin, A.R.; Pivato, M.; Mehdi, S.M.M.; Barnabas, L.E.; Giaretta, S.; Fabrega-Prats, M.; Prasad, D.; Arrigoni, G.; Masi, A. Proteome readjustments in the apoplastic space of Arabidopsis thaliana ggt1 mutant leaves exposed to UV-B radiation. Front. Plant Sci.; 2015; 6, pp. 1-12. [DOI: https://dx.doi.org/10.3389/fpls.2015.00128]
88. Masi, A.; Trentin, A.R.; Arrigoni, G. Leaf apoplastic proteome composition in UV-B treated Arabidopsis thaliana mutants impaired in extracellular glutathione degradation. Data Brief; 2016; 6, pp. 368-377. [DOI: https://dx.doi.org/10.1016/j.dib.2015.12.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26862584]
89. Ogawa, D.; Yamaguchi, K.; Nishiuchi, T. High-level overexpression of the Arabidopsis HsfA2 gene confers not only increased themotolerance but also salt/osmotic stress tolerance and enhanced callus growth. J. Exp. Bot.; 2007; 58, pp. 3373-3383. [DOI: https://dx.doi.org/10.1093/jxb/erm184] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17890230]
90. Muñoz, F.F.; Mendieta, J.R.; Pagano, M.R.; Paggi, R.A.; Daleo, G.R.; Guevara, M.G. The swaposin-like domain of potato aspartic protease (StAsp-PSI) exerts antimicrobial activity on plant and human pathogens. Peptides; 2010; 31, pp. 777-785. [DOI: https://dx.doi.org/10.1016/j.peptides.2010.02.001]
91. Egas, C.; Lavoura, N.; Resende, R.; Brito, R.M.M.; Pires, E.; Pedroso De Lima, M.C.; Faro, C. The saposin-like domain of the plant aspartic proteinase precursor is a potent inducer of vesicle leakage. J. Biol. Chem.; 2000; 275, pp. 38190-38196. [DOI: https://dx.doi.org/10.1074/jbc.M006093200]
92. Olivieri, F.P.; Lobato, M.C.; Altamiranda, E.G.; Daleo, G.R.; Huarte, M.; Guevara, M.G.; Andreu, A.B. BABA effects on the behaviour of potato cultivars infected by Phytophthora infestans and Fusarium solani. Eur. J. Plant Pathol.; 2009; 123, pp. 47-56. [DOI: https://dx.doi.org/10.1007/s10658-008-9340-z]
93. Tripathi, D.; Raikhy, G.; Kumar, D. Chemical elicitors of systemic acquired resistance—Salicylic acid and its functional analogs. Curr. Plant Biol.; 2019; 17, pp. 48-59. [DOI: https://dx.doi.org/10.1016/j.cpb.2019.03.002]
94. Yamaguchi-Shinozaki, K.; Shinozaki, K. Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters. Trends Plant Sci.; 2005; 10, pp. 88-94. [DOI: https://dx.doi.org/10.1016/j.tplants.2004.12.012]
95. HMMER. Available online: http://hmmer.org/ (accessed on 15 April 2020).
96. MAFFT <Multiple Sequence Alignment <EMBL-EBI. Available online: https://www.ebi.ac.uk/Tools/msa/mafft/ (accessed on 20 March 2021).
97. TAIR Home Page. Available online: https://www.arabidopsis.org/ (accessed on 29 October 2021).
98. Botella, M.A.; Parker, J.E.; Frost, L.N.; Bittner-eddy, P.D.; Beynon, J.L. Three Genes of the Arabidopsis. Society; 1998; 10, pp. 1847-1860.
99. Wheeler, T.J.; Clements, J.; Finn, R.D. Skylign: A tool for creating informative, interactive logos representing sequence alignments and profile hidden Markov models. BMC Bioinformatics; 2014; 15, 7. [DOI: https://dx.doi.org/10.1186/1471-2105-15-7] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24410852]
100. Simple Phylogenetic Tree <Phylogeny <EMBL-EBI. Available online: https://www.ebi.ac.uk/Tools/phylogeny/simple_phylogeny/ (accessed on 29 October 2021).
101. Punta, M.; Coggill, P.C.; Eberhardt, R.; Mistry, J.; Tate, J.; Boursnell, C.; Pang, N.; Forslund, S.K.; Ceric, G.; Clements, J. et al. The Pfam protein families database. Nucleic Acids Res.; 2011; 40, pp. 290-301. [DOI: https://dx.doi.org/10.1093/nar/gkr1065] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22127870]
102. ScanProsite. Available online: https://prosite.expasy.org/scanprosite/ (accessed on 29 October 2021).
103. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME Suite: Tools for motif discovery and searching. Nucleic Acids Res.; 2009; 37, pp. 202-208. [DOI: https://dx.doi.org/10.1093/nar/gkp335]
104. Nielsen, H. Predicting secretory proteins with signaIP. Methods Mol. Biol.; 2017; 1611, pp. 59-73. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28451972]
105. Small, I.; Peeters, N.; Legeai, F.; Lurin, C. Predotar: A tool for rapidly screening proteomes for N-terminal targeting sequences. Proteomics; 2004; 4, pp. 1581-1590. [DOI: https://dx.doi.org/10.1002/pmic.200300776]
106. Petersen, T.N.; Brunak, S.; Von Heijne, G.; Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. Nat. Methods; 2011; 8, pp. 785-786. [DOI: https://dx.doi.org/10.1038/nmeth.1701]
107. Sharma, P.; Lin, T.; Grandellis, C.; Yu, M.; Hannapel, D.J. The BEL1-like family of transcription factors in potato. J. Exp. Bot.; 2014; 65, pp. 709-723. [DOI: https://dx.doi.org/10.1093/jxb/ert432]
108. Briesemeister, S.; Rahnenführer, J.; Kohlbacher, O. Going from where to why-interpretable prediction of protein subcellular localization. Bioinformatics; 2010; 26, pp. 1232-1238. [DOI: https://dx.doi.org/10.1093/bioinformatics/btq115]
109. Briesemeister, S.; Rahnenführer, J.; Kohlbacher, O. YLoc-an interpretable web server for predicting subcellular localization. Nucleic Acids Res.; 2010; 38, pp. 497-502. [DOI: https://dx.doi.org/10.1093/nar/gkq477]
110. Emanuelsson, O.; Brunak, S.; von Heijne, G.; Nielsen, H. Locating proteins in the cell using TargetP, SignalP and related tools. Nat. Protoc.; 2007; 2, pp. 953-971. [DOI: https://dx.doi.org/10.1038/nprot.2007.131]
111. De Sá, M.M.; Drouin, G. Phylogeny and substitution rates of angiosperm actin genes. Mol. Biol. Evol.; 1996; 13, pp. 1198-1212. [DOI: https://dx.doi.org/10.1093/oxfordjournals.molbev.a025685]
112. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van De Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res.; 2002; 30, pp. 325-327. [DOI: https://dx.doi.org/10.1093/nar/30.1.325]
113. Search for Cis-Acting Regulatory Elements. Available online: http://bioinformatics.psb.ugent.be/webtools/plantcare/html/search_CARE.html (accessed on 29 October 2021).
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Sebastián D’Ippólito 3 ; Cecilia Andrea Décima Oneto 1 ; Massa, Gabriela Alejandra 2
; Castellote, Martín Alfredo 1 ; Feingold, Sergio Enrique 1 ; Guevara, María Gabriela 3
1 Laboratory of Agrobiotechnology IPADS (INTA—CONICET), Balcarce B7620, Argentina;
2 Laboratory of Agrobiotechnology IPADS (INTA—CONICET), Balcarce B7620, Argentina;
3 Institute of Biological Research, University of Mar del Plata (IIB-UNMdP), Mar del Plata B7600, Argentina;