Next Article in Journal
Newcastle Disease Virus Induced Pathologies Severely Affect the Exocrine and Endocrine Functions of the Pancreas in Chickens
Next Article in Special Issue
Genomic and Phenotypic Diversity of Cultivated and Wild Tomatoes with Varying Levels of Heat Tolerance
Previous Article in Journal
Genomics of the Parasitic Nematode Ascaris and Its Relatives
Previous Article in Special Issue
Comparative Transcriptome Analysis Reveals Key Genes and Pathways Involved in Prickle Development in Eggplant
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unravelling Differences in Candidate Genes for Drought Tolerance in Potato (Solanum tuberosum L.) by Use of New Functional Microsatellite Markers

1
Institute of Biological Sciences, Plant Genetics, University of Rostock, Albert-Einstein-Str. 3, 18059 Rostock, Germany
2
MPI für Molekulare Pflanzenphysiologie, Am Mühlenberg 1, 14476 Potsdam, Germany
3
Department of Food Safety, German Federal Institute for Risk Assessment, Max-Dohrn Straße 8-10, 10589 Berlin, Germany
4
Deceased, Formerly MPI für Molekulare Pflanzenphysiologie, Am Mühlenberg 1, 14476 Potsdam, Germany
5
Julius Kühn-Institut, Rudolf-Schick-Platz 3, 18190 Sanitz, Germany
6
Landwirtschaftskammer Niedersachsen, Dethlingen 14, 29633 Munster, Germany
7
PotatoConsult UG, Hiddinger Straße 33, 27374 Visselhövede, Germany
8
Ministry of Agriculture, Environment and Climate Protection, Henning-Von-Tresckow-Straße 2-13, 14467 Potsdam, Germany
*
Author to whom correspondence should be addressed.
Genes 2021, 12(4), 494; https://doi.org/10.3390/genes12040494
Submission received: 23 February 2021 / Revised: 23 March 2021 / Accepted: 24 March 2021 / Published: 28 March 2021
(This article belongs to the Special Issue Genetics and Genomics of Solanaceae)

Abstract

:
Potato is regarded as drought sensitive and most vulnerable to climate changes. Its cultivation in drought prone regions or under conditions of more frequent drought periods, especially in subtropical areas, requires intensive research to improve drought tolerance in order to guarantee high yields under limited water supplies. A candidate gene approach was used to develop functional simple sequence repeat (SSR) markers for association studies in potato with the aim to enhance breeding for drought tolerance. SSR primer combinations, mostly surrounding interrupted complex and compound repeats, were derived from 103 candidate genes for drought tolerance. Validation of the SSRs was performed in an association panel representing 34 mainly starch potato cultivars. Seventy-five out of 154 SSR primer combinations (49%) resulted in polymorphic, highly reproducible banding patterns with polymorphic information content (PIC) values between 0.11 and 0.90. Five SSR markers identified allelic differences between the potato cultivars that showed significant associations with drought sensitivity. In all cases, the group of drought-sensitive cultivars showed predominantly an additional allele, indicating that selection against these alleles by marker-assisted breeding might confer drought tolerance. Further studies of these differences in the candidate genes will elucidate their role for an improved performance of potatoes under water-limited conditions.

1. Introduction

Drought has become the predominant abiotic stress and major yield limiting factor in crops in the past years [1,2,3]. Severe implications for food production have been observed [4,5,6] confronting plant breeding with the challenge of securing food supplies at a global level under these changing growth conditions. As an effect, research into drought tolerance mechanisms has been intensified in most important crops, such as maize [7,8,9], wheat [2,10], and rice [11,12,13]. Drought tolerance involves coordination of heavily interwoven networks of biosynthesis and signaling pathways of phytohormones such as abscisic acid (ABA), ethylene, jasmonate and the production of osmoprotective compounds such as proline, glycine betaine and trehalose [14]. Numerous transcription factors like ABA-responsive element (ABRE) binding proteins, dehydration responsive element binding (DREB) proteins, and ethylene response factors (ERFs) were identified as major players in controlling the expression of genes involved in drought response in plants [15,16,17].
As a temperate-zone crop with a shallow active root zone, potato is sensitive to drought [18,19] and heavy yield losses are predicted based on climate change prognosis [20,21,22]. Morphological and physiological traits such as shoot height, leaf size, number of leaves, photosynthetic rate, biomass as well as tuber yield are affected by water shortage [23]. Damage caused by water deficit is dependent on the developmental stage, the duration and the intensity of the drought period [18,24]. Potato tuber bulking stage (early stress) is the most critical phase with regard to water demands and an inadequate water supply can severely affect tuber yield and tuber quality [18,25]. Evaluation of marketable tuber yield in a panel of 103 European potato cultivars with different maturity identified cultivars with better performance under non-irrigated conditions and pointed at a single nucleotide polymorphism (SNP) PotVar0030768 from the 14 Infinium SNP marker array to be significantly associated with this trait [26]. Potato cultivars with a higher root to shoot ratio due to deeper roots show a better performance under drought conditions [27]. Change of plant morphology towards small, open stem-type canopies combined with shallow but dense root systems may be an interesting approach for improving drought tolerance in potato [28].
Different systems have been used to study drought stress in potato: field trials (irrigated and non-irrigated), rain-out shelters and/or trials under fully controlled conditions [29,30,31] as well as in vitro systems using proxies, such as sorbitol or polyethylene glycol [32,33,34]. Drought stress indices such as the stress susceptibility index (SSI) [35], the stress tolerance index (STI) or the geometric mean productivity (GMP) [36] have been successfully applied to define stress levels in experimental designs. Using an artificial data set, a new drought stress index calculated as deviation of the relative starch yield under drought and control conditions from the experimental median (DRYM) outperformed the other three stress indices in the ability to differentiate between drought sensitive and drought tolerant potato cultivars detached from the yield potential [24].
Transcriptome analyses identified differentially upregulated drought response genes involving ABA and ethylene signaling transduction like protein phosphatase 2C (PP2C), homolog of abscisic acid receptor PYL4 (PYRABACTIN RESISTANCE-LIKE 4), 1-aminocyclopropane-1-carboxylate oxidase (ACO) and ethylene responsive transcription factors in potato [19,37,38,39]. Moreover, expansin and genes involved in osmotic adjustment and protection such as delta-1-pyrroline-5-carboxylate synthase (P5CS), galactinol synthase, inositol-3-phosphate synthase, raffinose synthase (RFS), osmotin-like (OSML) and late embryogenesis abundant proteins (LEA) showed increased expression levels under drought stress. Differentially expressed transcription factors including 10 gene families such as bHLH, bZIP, ERF, HD-ZIP, MYB, NAC and WRKY represented 4.7% of the drought responsive genes in potato [37]. Investigating drought response in stolon tips significant differences in the expression of heat shock 70 kDa protein, aquaporins, bidirectional sugar transporter, peroxidase and pectinesterase/pectinesterase inhibitor as well as lipid-transfer proteins (LTP) were observed [38]. Among genes that showed the highest induction or repression comparing drought stress response with re-watering and control treatment were those coding for TAS14, WRKY, ERF and bHLHL transcription factors, auxin response protein, and abscisic acid 8′-hydroxylases [39]. In addition, protein kinases, mostly receptor-like kinases, but also kinases like mitogen-activated protein kinase kinase (MAPKK), histidine kinases, cyclin-dependent kinases and Sucrose non-fermenting 1-related protein kinase 2 (SnRK2) were differentially expressed under drought stress in potato [39,40]. Comparing Gwaizda/Oberon and Taijun/Owacja as two pairs of drought tolerant and drought sensitive potato cultivars 22 novel drought response genes were identified [41]. For six genes (armadillo/β-catenin-like repeat-containing protein, carbohydrate transporter, MAPKKK 15, nitrate transporter 2.7, nonspecific LTP and serine carboxypeptidase-like 19 protein), corresponding Arabidopsis mutants showed altered gene expressions and improved drought tolerance [41]. Constitutive differences in gene expression under both control and drought stress conditions have also been observed between drought tolerant and sensitive potato cultivars [40]. These included among others S-adenosyl-L-methionine-dependent methyltransferase superfamily, cytochrome P450, UDP-glycosyl transferase family proteins, disease-related genes and receptor kinases [40]. A random forest model allowed prediction of drought tolerance for different potato cultivars based on eight metabolite (e.g., 2-oxo-glutaric acid, ribonic acid, galactonic acid) and 19 transcript markers (e.g., glucosyltransferase, O-methyltransferase, extensin, poly(ADP-ribose) glycohydrolase, lipoxygenase) with an accuracy of 82.6% [31].
Transcriptome studies in potato identified drought response genes by differential gene expression, but have so far not distinguished between possible different alleles of the individual genes. Potato cultivars are highly heterozygous tetraploids allowing the presence of up to four different alleles per gene. A rapid and easy conversion of identified drought response candidate genes into molecular markers detecting variation in these genes will be essential for successful drought tolerance breeding. Simple sequence repeats (SSRs), also known as microsatellites, are the most frequently used markers for genotyping potato landraces [42]. Apart from pure forms, SSRs can also occur in complex, compound, interrupted pure, interrupted compound and interrupted complex forms [43,44]. In potato, SSR markers based on pure, compound and complex repeats have been successfully developed [45,46,47,48,49,50]. A new potato genetic identity kit for fingerprinting consisting of 51 SSR markers was finally assembled out of previously published SSR markers and new ones [51]. In addition, a genetic reference map based on SSR markers was developed [51]. Another set of 277 SSR markers based on pure repeats was derived from the whole potato genome sequence [43]. SSRs have been successfully applied in genetic diversity studies in potato cultivars, landraces and germplasms [51,52,53,54], in population structure analyses of association panels [55,56] and for the development of linkage maps for QTL analyses [47,57,58].
In this study, we present the identification of new functional SSR markers derived from potential candidate genes for drought tolerance in potato as well as association studies performed with these SSRs in a panel of 34 mainly starch potato cultivars, which were ranked for drought tolerance using DRYM [24] as drought tolerance index. The approach allowed the development of 75 new SSR markers that represent a valuable tool for mapping and further genetic studies in potato. Five of the SSR markers detected allelic variations in the potato cultivars that were significantly associated with drought sensitivity. Future functional studies of these differences in the candidate genes will reveal the role in drought tolerance. Combined application of three of these SSR markers allowed the identification of drought-sensitive cultivars in potato.

2. Methods

2.1. Extraction of Genomic DNA from Potato

Fully developed leaves of 34 tetraploid potato cultivars (Supplementary Table S1) were harvested from a polytunnel trial at the Max Planck Institute of Molecular Plant Physiology in Potsdam. This association panel consisting of 34 mainly starch potato cultivars was provided by German breeding companies (Böhm-Nordkartoffel-Agrarproduktion, Strehlow; Norika Nordring-Kartoffelzucht- u. Vermehrungs-GmbH, Groß-Lüsewitz; SaKa Pflanzenzucht, Hamburg). Genomic DNA was isolated by the protocol of Doyle and Doyle [59].

2.2. Drought Tolerance Assessment

Drought tolerance of 34 European potato cultivars was assessed based on performance data from drought tolerance trials on three managed field sites in Germany (Max Planck Institute of Molecular Plant Physiology Potsdam-Golm [F1, F3] and Potato Research Station Dethlingen [F2, F5] in 2011 and 2012; Julius Kühn-Institute Groß-Lüsewitz [F4], only in 2012) as described in Sprenger et al. [24], but the field trials of 2013 were omitted. All five investigated field trials represent early stress and corresponding watered controls, whereas from the three omitted field trials in 2013 two field trials represented late stress scenarios. Plants were grown from seed tubers supplied by the respective breeding companies (Supplementary Table S1, [24]) and cultivated under agricultural conditions with additional irrigation for the optimally watered control plants and without or reduced irrigation for the drought stress treated plants (Supplementary Table S2). All treatments received additional water from precipitation, with the exception of the drought treatment in experiment F4, which was grown under a rain-out shelter. Plant performance was assessed based on the tuber starch yield. Tubers were harvested and weighted at the end of the experiment and tuber starch content was measured gravimetrically. Starch yield is the product of tuber mass and starch content for each replicate.
Data evaluation was performed in SAS (SAS 9.3 Level 1MO, Cary, NC, USA). The effect of cultivar and treatment on starch yield was assessed by analysis of variance with procedure GLM. Based on the results, data from the treatment “30% Field Capacity (FC)” at the field site Dethlingen 2011 were excluded from the analysis as they were not significantly different from the control treatment “50% FC”. In 2012, data from treatment “30% FC” and “no irrigation” were taken together as stress treatment. Drought tolerance was calculated as deviation of the relative starch yield of a genotype from the median relative starch yield of the respective experiment (DRYM) (see Equation (1)). Relative starch yield (RelSYGxEi) is the ratio of the starch yield under drought stress divided by the average starch yield of the respective cultivar (Gx) under control conditions in the same experiment (Ei).
D R Y M G x , E i = R e l S Y G x , E i m e d i a n ( R e l S Y E i )
This calculation centers DRYM around zero, with drought tolerant cultivars showing positive values and sensitive genotypes negatives values.
Analysis of covariance (ANCOVA) was performed with procedure GLM, testing the effect of experiment and cultivar (random factors) and normalized starch yield as covariable on DRYM. Means comparison for cultivars was performed with Ryan–Einot–Gabriel–Welch-Test (REGWQ) that provides multiple testing corrections.

2.3. Mining for Simple Sequence Repeats (SSR) and Characterization

Genomic sequences of the candidate genes for drought tolerance were retrieved from different databases as Phytozome (https://phytozome.jgi.doe.gov/pz/portal.html (accessed on 27 March 2021)), EnsemblPlants ([60], http://plants.ensembl.org/index.html (accessed on 27 March 2021)), Spud DB Potato Genomics Resource ([43], http://potatogenomics.plantbiology.msu.edu (accessed on 27 March 2021)), Sol Genomics Network (https://solgenomics.net/ (accessed on 27 March 2021)) and NCBI (https://www.ncbi.nlm.nih.gov/ (accessed on 27 March 2021)). The retrieved sequences were then mined for repeats using the default parameters in the program Microsatellite repeats finder (http://insilico.ehu.es/mini_tools/microsatellites/ (accessed on 27 March 2021)). In case, no microsatellites could be detected in the gene sequence (exons and introns), sequences of the 5′- and 3′-UTR were studied. If still no repeats could be observed the flanking downstream and upstream sequences were investigated for microsatellites, which were then termed intergenic SSRs. The microsatellites are given as present in the coding strand as this is the custom for naming repeats in coding regions [44]. Primer combinations were then derived in sequence segments of about 240 bp around the repeat sequence using the program Primer3 (http://bioinfo.ut.ee/primer3-0.4.0/ (accessed on 27 March 2021)). All 154 primer combinations are given in Supplementary Table S3. The polymerase chain reaction (PCR) was performed as described in Sajer et al. [61]. However, the DNA concentration was increased from 10 ng/µL to 20 ng/µL considering the heterozygous tetraploid genetic constitution of cultivated potatoes. The PCR reaction volume was reduced to 15 µL. By using IRD-labelled M13 primers (IRD700 and IRD800) (Biomers, Ulm, Germany) PCR products were fluorescence labelled for the detection on the DNA Analyzer Model 4300 (LI-COR Biosciences, Lincoln, NE, USA). The newly derived SSR primer combinations were analyzed in the association panel of 34 potato cultivars. The estimation of the polymorphic information content (PIC) values and of the expected heterozygosity (Hexp) was performed according to Bali et al. [52].
In order to determine the statistical significance of the banding patterns concerning drought tolerance the Fisher’s exact test was calculated online in a case-control type study dividing the association panel into two groups of 17 drought tolerant and 17 drought sensitive cultivars (https://www.langsrud.com/stat/Fishertest.htm (accessed on 27 March 2021)). A significant association was reached at a (two-sided) p value (calculated probability) < 0.05.

2.4. AFLP Analyses

The potato DNA was digested with EcoRI and MseI and ligated to EcoRI and MseI adapters as described in Vos et al. [62]. AFLP analyses were performed based on a pre-amplification with E01 (5′-GACTGCGTACCAATTCA-3′) and M02 (5′-GATGAGTCCTGAGTAAC-3′) as primers. For the selective amplification, the EcoRI primers E32–E35 were combined with MseI primers M47–M49 or M60–M62 (Supplementary Table S4). In addition, E32M50 was used. All primer sequences and numbers were used according to http://wheat.pw.usda.gov/ggpages/keygeneAFLPs.html (accessed on 27 March 2021) (Keygene, N.V., Wageningen, NL, USA). IRD-labelled EcoRI primers (IRD700 and IRD800) from Biomers (Ulm, Germany) were used for fluorescence labeling of the selective amplification products, which were separated on a DNA Analyzer Model 4300 (LI-COR Biosciences, Lincoln, NE, USA).

2.5. Population Structure

In total, 54 linkage group specific SSR primer combinations (Supplementary Table S5) were used to determine the population structure of the association panel, resulting in at least four SSR primer combinations per linkage group. In addition, data from 19 AFLP primer combinations were included to have better genome coverage (Supplementary Table S4). The program NTSYSpc 2.2 (Numerical Taxonomy System for personal computer) was used to identify and reveal structures in multivariate data [63]. The similarity matrix was based on the Jaccard’s coefficient to reflect the genetic similarity of the 34 different potato cultivars. For clustering, UPGMA (Unweighted Pair Group Method with Arithmetic Mean) was applied to obtain a dendrogram. Mantel test was conducted, and a cophenetic correlation coefficient computed, which allowed a statement regarding the reliability of the UPGMA-dendrogram in representing the similarity matrix.

3. Results

3.1. Characterization of Functional SSR Markers Derived from Candidate Genes for Drought Tolerance in Potato

In total, 154 SSR primer combinations were derived from 103 candidate genes for drought tolerance (Figure 1; Supplementary Table S3). Eighty of the candidate genes were taken from previous publications (reviewed in Krannich et al. [14]) and 23 were identified in a drought-related transcriptome study as top transcript markers for drought tolerance prediction (Supplementary Table S6, [39]). Some candidate genes are involved in biosynthesis and signaling pathways of plant hormones like ethylene and ABA as well as the biosynthesis of osmolytes and starch (Figure 1).
Other potential candidate genes for drought tolerance were transcription factors, kinases/phosphatases and transporters. Another group covered the area of plant development, especially tuber formation. In addition, genes encoding proteins involved in cell wall formation, defense, protective cellular structures, protein modification and DNA repair as well as detoxification were included.
Most of the detected microsatellites showed an interrupted complex structure of the type (repeats)k X (repeats)n(repeats)m (42%) or in 18% a compound structure such as (repeats)n (repeats)m, where the repeats were di-, tri-, tetra-, penta- and hexanucleotide motifs, the number of repeats n + m or k + n+m added up to at least six and the number of nucleotides X between the repeats did not exceed four, with nine exceptions (Figure 2A, Supplementary Table S3). Most of the microsatellite repeats were located in intron regions (44%), followed by 28% in exon regions. A lower number of SSRs was present in 5′-UTRs (13%), in 3′-UTRs (8%), and in intergenic regions (7%) (Figure 2a). From 154 microsatellites tested and optimized in an association panel of 34 tetraploid starch potato cultivars, 75 SSR primer combinations (49%) gave reproducible polymorphic banding patterns (Supplementary Table S7). For these 75 SSR markers, a total of 278 alleles were observed resulting in an average of 3.7 alleles per SSR marker. The minority of these polymorphic SSR markers represented complex repeats (9%) and pure repeats (10%), followed by interrupted pure repeats (12%) (Figure 2b). Compound (17%), interrupted compound (20%) and interrupted complex (32%) repeats were the most prominent SSR types.
Only seven candidate genes for drought tolerance (PGSC0003DMG400021476, PGSC0003DMG400000284, PGSC0003DMG400021426, PGSC0003DMG400009719, PGSC0003DMG400021651, PGSC0003DMG402026767, PGSC0003DMG400014417) contained pure repeats. These repeats were (AT)n, (GA)n, (ACT)n, (TTA)n and (CGG)n. In three of the candidate genes, 1-aminocyclopropane-1-carboxylate synthase 4 (PGSC0003DMG400021651, ([ACT]n], delta-1-pyrroline-5-carboxylate synthetase (PGSC0003DMG402026767; [TTA]n) and ethylene response factor ERF12 (PGSC0003DMG400014417; [CGG]n), SSRs had been previously described [43], but had not been screened with regard to drought tolerance. The polymorphic information content (PIC) for the 75 SSR markers ranged between 0.11 and 0.90 and for the expected heterozygosity (He) between 0.11 and 0.91 (Supplementary Table S7). We also tested exceptions from X > 4 in nine cases by allowing X = 5 (HRO_BDGEH_1A, HRO_BDGEH_1B, HRO_BFNLRP_2C, HRO_GLUCT_7B, HRO_STPKPT_1A) or X = 6 (HRO_ETOP1_3, HRO_ETP2_2, HRO_ETR2_1, HRO_GLUCT_7), because there were no other possibilities for exploring microsatellites in these genes. In four cases, this resulted also in polymorphic banding patterns and might be considered as a possible option. Polymorphic SSR markers were mostly derived from intron sequences (45%), but also to a considerable percentage from 5′- and 3′-untranslated regions (29%), where changes in the repeat number might have an influence on gene expression. Only a low percentage of polymorphic SSR markers (15%) were obtained from exon regions as in case of PP2C (PGSC0003DMG400011321) and ethylene responsive transcription factor (PGSC0003DMG400002185). Here a change in the number of repeats will lead to changes in the encoded proteins.
A single locus was amplified by nearly all developed SSR primer combinations, resulting in 1-4 bands for each potato cultivar. The exception was the primer combination for aldehyde dehydrogenase (ALDH), which detected at least two loci. All 12 chromosomes carried candidate genes for drought tolerance that allowed development of functional SSR markers. A clustering of candidate genes was observed in the distal parts of linkage group (LG)1, LG2, LG5 and LG6 as well as the proximal parts of LG1 and LG12 (Figure 3).

3.2. Ranking of the 34 Starch Potato Cultivars According to Drought Tolerance

Drought tolerance of 34 potato cultivars (Supplementary Table S1) was assessed from the tuber starch yield obtained under optimal and reduced water supply in five field experiments. These experiments were performed for two years (2011 and 2012) on three sites (details see [24]). Analysis of covariance indicated that starch yield was significantly affected by the drought treatment in all experiments [24]. The drought tolerance index DRYM was calculated as the deviation of the relative starch yield from the experimental median of the starch yield (Figure 4). Drought tolerant potato cultivars show positive DRYM values, whereas drought sensitive genotypes have negative values. Analysis of covariance indicated a significant effect of cultivar on DRYM (p < 0.0001), suggesting a significant genetic variance of drought tolerance within the population. The means test (see bars at the bottom of Figure 4) indicates a more or less continuous variation in tolerance. Tolerance was significantly higher (α = 0.1) in the cultivars with the ranks 1 to 14 compared to those of rank 20 to 34. Subsequent comparison of means (REGWQ-Test, α = 0.1) yielded the ranking shown in Supplementary Table S1.

3.3. Population Structure of the Investigated 34 Potato Cultivars

The 34 potato cultivars were analyzed with 19 AFLP primer combinations (Supplementary Table S4) and 54 linkage group-specific SSR primers (Supplementary Table S5). The population structure shown in an UPGMA dendrogram was based on the Jaccard’s coefficient (Supplementary Figure S1). The Mantel test resulted in a cophenetic correlation coefficient of r = 0.62, indicating a moderate fit. No clustering of the potato cultivars was observed that would lead to a false association of markers regarding drought tolerance.

3.4. Identification of SSR Markers Associated with Drought Tolerance

Association studies applying a case-control type approach were performed within the association panel of 34 potato cultivars, which had been previously ranked according to their drought tolerance using DRYM as drought tolerance index (1t = most tolerant and 34t = most sensitive) (Figure 4, Supplementary Table S1). After splitting the panel into two groups, one of drought tolerant (1t–17t) and one of drought-sensitive (18t–34t) cultivars, all bands were tested for association within these two groups using the Fisher’s Exact test (Supplementary Table S8). This association study allowed the identification of differences in five SSR markers significantly associated with drought sensitivity (Figure 5).
The first primer combination HRO_ACS3 was derived from the 1-aminocyclopropane-1-carboxylate synthase 3 (ACS3, PGSC0003DMG400021426, located on LG2) representing a rate-limiting step of the ethylene biosynthesis. This primer combination amplified 1-4 alleles in the different cultivars (Figure 5a). The band of 173 bp (HRO_ACS3_D) was significantly associated with drought sensitivity (p = 0.0366). This band was present in 13 drought sensitive cultivars. In addition, most of the more drought tolerant cultivars showed only one to two bands, whereas the more sensitive cultivars tended to display two to three alleles (Table 1).
The second primer combination HRO_PP2C amplified two alleles, a monomorphic band of 205 bp and a smaller band of 196 bp that was predominantly present in the drought sensitive group (HRO_PP2C_1_B, p = 0.0366) (Figure 5b). This SSR detected differences in the gene encoding a protein phosphatase 2C (PP2C, PGSC0003DMG400011321, located on LG1) as part of the ABA signaling pathway.
The third primer combination HRO_ALDH was derived from the aldehyde dehydrogenase (ALDH, PGSC0003DMG400034597, located on LG9). This primer combination amplified more than one locus showing on average 6.4 bands per cultivar. One band of 184 bp (HRO_ALDH_H, p = 0.0366) was significantly associated with drought sensitivity, a second one of 172 bp (HRO_ALDH_N, p = 0.0854) was only weakly associated (Figure 5c), but proved to be important for further analyses.
With the fourth primer combination HRO_ETRTF_5a derived from an ethylene responsive transcription factor (ERF, PGSC0003DMG400002185, located on LG 11) two to four alleles were amplified (Figure 5d). Again, the more drought sensitive cultivars showed a more complex pattern. The band HRO_ETRTF_5a_D (217 bp) was associated with drought sensitivity (p = 0.0366). The band was present in 11 drought sensitive cultivars, but only in four tolerant cultivars.
The last SSR marker (HRO_PARG_1A) was derived from a gene encoding a poly (ADP-ribose) glycohydrolase (PARG, PGSC0003DMG400029361, located on LG 12), which is involved in DNA repair. The primers amplified one to two alleles per cultivar with the lower band of 147 bp showing association with drought sensitivity (p = 0.0324) (Figure 5e).
The banding pattern of all five candidate genes was characterized by an additional allele in the group of drought sensitive cultivars, which was significantly associated with drought sensitivity indicating that in all cases these bands might represent alleles inferring less drought tolerance (Table 2).

3.5. Selection Using Functional SSR Markers Associated with Drought Tolerance

The banding pattern obtained by 75 polymorphic SSR primer combinations derived from candidate genes for drought tolerance was used to estimate the genetic similarity among 34 potato cultivars based on the Jaccard’s coefficient (Figure 6). The clustering of the 34 potato cultivars as shown in the UPGMA dendrogram did not reveal any grouping corresponding to drought tolerance. However, when the SSR markers HRO_ACS3_D, HRO_PP2C_B, HRO_ALDH_H and HRO_ALDH_N associated with drought sensitivity were used for marker-assisted selection, 15 out of 17 drought-sensitive cultivars were correctly identified. Only 33t and 25t were not detected as drought-sensitive and only 16t was falsely identified as drought sensitive, even though the ranking put this cultivar in the drought tolerant group. The marker-assisted selection was performed in the way that a cultivar was considered drought-sensitive if both bands, HRO_ACS3_D as well as HRO_PP2C_B, were present or if the cultivar was characterized by the presence of both bands H and N produced by the primer combination HRO_ALDH (Figure 5). The band HRO_ALDH_N had only been weakly associated with drought tolerance using the Fisher’s exact test. However, it proved to be essential for the selection of drought sensitive cultivars using the three primer combinations HRO_ACS3, HRO_PP2C and HRO_ALDH.
The alleles for drought sensitivity detected by the three SSR primer combinations should be eliminated by breeding to obtain more drought tolerant cultivars. Including the two additional primer combinations HRO_ETRTF_5a and HRO_PARG_1A in the selection procedure did not improve the selection efficiency, even though most drought sensitive cultivars (14) contained at least one of these negative alleles, six even both negative alleles detectable by these two SSR primer combinations (Supplementary Table S8).

4. Discussion

4.1. Development of Drought Tolerance Associated SSR Markers for Potato

In the present study, 75 functional SSR markers were developed starting from 103 potential candidate genes for drought tolerance derived from own transcriptome analyses [39] and literature. SSR markers derived from five candidate genes (aldehyde dehydrogenase, PGSC0003DMG400034597; 1-aminocyclopropane-1-carboxylate synthase 3, PGSC0003DMG400021426; ethylene responsive transcription factor, PGSC0003DMG400002185; poly (ADP-ribose) glycohydrolase, PGSC0003DMG 400029361; protein phosphatase 2C, PGSC0003DMG400011321) were significantly associated with drought sensitivity. These five SSR-markers identified allelic differences between potato cultivars that are crucial for sensitivity towards drought. Selection for drought tolerance using three of the significantly associated SSR markers allowed in 91 % of the cases a correct grouping into drought tolerant or drought sensitive potato cultivars. Now, it would be interesting to test these SSR markers in a larger international potato association panel.
Our approach to derive SSR markers from candidate genes for drought tolerance made mostly use of the presence of interrupted pure, compound, interrupted compound, complex and interrupted complex repeats. This strategy differed from previous studies, which concentrated mostly on the use of pure as well as compound and complex SSRs [45,46,47,48,49,50,51]. Our results showed that the use of SSR markers derived from potential candidate genes worked well and can be used to detect genetic diversity in potato cultivars. In addition, the number of repeats was lowered from eight to six as compared to previous studies, which still resulted in highly polymorphic banding patterns. Comparing the locations for all derived SSR primer combinations versus only the polymorphic SSR markers the percentage of polymorphic SSR markers located in exon regions was considerably reduced to only 15% compared to 28% for all SSR primer combinations. On the other hand, the percentages of polymorphic SSR markers derived from 5′-UTR, 3′-UTR and intergenic regions and were slightly increased with 18%, 11% and 11%, respectively.

4.2. Role of PP2C as Part of the ABA Signaling Pathway

Our SSR analyses showed a significant association of drought sensitivity with allelic variations in PP2C (PGSC0003DMG400011321) located on LG1. The interrupted complex SSR (GAA)kA(AT)mCAT(GAT)n used for detection is located in exon 3 of the gene. Changes in the repeat numbers will have an effect on the encoded protein. PP2Cs are known drought responsive genes in potato. Induction of expression was shown for drought exposed moderately drought tolerant native Andean potato clones [19]. In addition, comparing transcriptome profiles of drought treated and re-watered plants of the potato strain Ningshu4 also revealed a contrasting expression pattern for PP2C 30 (PGSC0003DMG400027174, located on LG5) [38]. In diploid potato, a different PP2C (PGSC0003DMG400027196, also located on LG5) was identified as a central point in the gene regulatory network [64]. PP2Cs are part of the ABA signaling pathway, which consists of three major players: ABA receptors, PP2C and SnRK2 [65]. The ABA signaling pathway leads to phosphorylation of SnRK2 by inhibiting dephosphorylation by PP2C in the presence of ABA. Auto- and trans-activation of SnRK2 then results in phosphorylation and activation of downstream transcription factors [66]. SnRK2 modulates gene expression by targeting basic-zipper class transcription factors (ABI5) for arrest of seed germination, but also phosphorylates SLAC1 (slow anion channel-associated 1) to induce stomatal closure under limiting water conditions [67,68,69].
Arabidopsis mutants of PP2C like ABI1 (ABA insensitive 1) and its homolog ABI2 revealed PP2C role as negative regulator of ABA trigged responses such as inhibition of germination, root elongation, and induction of stomatal closure [70]. Furthermore, PP2C mutant sag113 (senescence associated gene 113) showed higher ABA sensitivity for stomatal movement and reduced leaf senescence [71]. In wheat, a genome-wide study of the TaPP2C gene family identified 257 homoeologs of 95 TaPP2C genes distributed over the A, B and D sub genomes [72]. Phylogenetic analyses separated the TaPP2C genes into 13 groups A-M. TaPP2Cs of group A were upregulated by ABA treatment and showed interaction with TaSnRK2.1 and TaSnRK2.2. Overexpression of TaPP2C135 in Arabidopsis increased tolerance of germinating seedlings towards ABA [72].
Enhanced ABA synthesis during drought stress induces signaling cascades that activate plant stress responses and regulate the water status of plants mainly through management of stomatal aperture and maintenance of root growth [73]. Therefore, allelic variations in PP2C (PGSC0003DMG400011321) as part of the ABA signaling pathway could play a major role in drought tolerance in potato and further functional studies of the different allelic compositions are required to better understand the causes behind drought tolerance in potato.

4.3. Damage Control of Drought Induced Reactive Oxygen

The production of reactive oxygen species (ROS) under drought stress leads to oxidative stress and an excessive accumulation of oxidized toxic compounds such as aldehydes [74]. ALDHs act as aldehyde scavengers by irreversibly oxidizing aldehydes into the respective carboxylic acids [75]. Our SSR analyses in potato showed that one allelic difference in ALDH (PGSC0003DMG400034597, LG9) was significantly associated with drought sensitivity and another one also relevant for selection. This ALDH demonstrates the highest similarity in Arabidopsis thaliana to AT1G54100 (84.1% similarity according to Spud DB), which is annotated as ALDH family 7 member B4. In A. thaliana, members of ALDH gene families 3 and 7 (ALDH3I1 and ALDH7B4), show the strongest transcriptional activation under osmotic stress [76]. Cytosolic overexpression of ALDH7B4 confers enhanced detoxification and turgor response, demonstrating its importance under stress conditions [77]. As ALDH7B4 plays a major role in drought tolerance in A. thaliana, association of variation in this enzyme with drought sensitivity in potato suggests a functional relevance, which needs to be further investigated. In addition, the upregulation of this ALDH (PGSC0003DMG400034597) under water-limited conditions in greenhouse and field trials supports its importance in drought stress response in potato [39].

4.4. Damage Repair of DNA under Drought

Drought stress can lead to DNA lesions including double strand breaks (DSBs) as a major form [78,79,80]. DNA damage surveillance and repair processes have to recruit a number of repair proteins including poly (ADP-ribose) polymerases (PARPs), which recognize DNA SSBs and DSBs, assemble DNA repair complexes and auto-PARylate [81,82]. As the PARylation process is reversible, ADP-ribose chains are eventually hydrolyzed by the enzyme poly (ADP-ribose) glycohydrolase (PARG) [83]. As DNA damage repair is required for recovery of plants from drought stress, our SSR results implicating a role of allelic variation in PARG (PGSC0003DMG400029361, LG12) with drought tolerance in potato seem relevant. PARG (PGSC0003DMG400029361) was also used as one of the 27 predictors in forecasting drought tolerant potato based on transcripts and metabolites [31].
The PARG gene (PGSC0003DMG400029361) in potato shows the highest similarity on protein level to PARG1 (AT2G31870) (73.2% similarity according to Spud DB), one of two PARG genes present in tandem array on chromosome 2 in Arabidopsis [84]. An Arabidopsis parg1 mutant exhibits loss of drought, osmotic, and oxidative stress tolerance, coupled with enhanced cell damage and loss of stomatal closure under drought stress [85].

4.5. Role of ACC Synthase and Ethylene Response Factors in Drought Tolerance

Ethylene represents a small two carbon gaseous hormone, which is involved in abiotic and biotic stress reactions, seed germination, fruit ripening, cell elongation, cell fate, leaf abscission and senescence [86]. Synthesis of ethylene involves two enzymes, 1-aminocyclopropane-1-carboxylic acid (ACC) synthase (ACS) and ACC-oxidase (ACO) [87,88].
For potato, we identified allelic differences in ACC synthase 3 (StACS3, PGSC0003DMG400021426, located on LG2) to be significantly associated with drought sensitivity. StACS3 shows the highest similarity to the Arabidopsis homolog AtACS8 (85% similarity according to Spud DB) encoded by AT4G37770. In roots, higher expression levels have been observed for this gene under drought conditions (TAIR Arabidopsis eFP Browser [89]). In Arabidopsis, the ACS gene family originally contained 12 annotated members ACS1-12 [90]. A large scale of single as well as multiple mutants (acs1, acs2, acs4, asc5, acs6 and acs9) demonstrated interactions between different ACS isoforms [91]. The ACS family is grouped into three types I-III based on the C-terminal end [92]. StACS3 in potato represents a type II ACS as ACC synthases encoded in Arabidopsis by AtACS4, AtACS5, AtACS8, AtACS9 and AtACS11 [93]. These ACS are characterized by a shorter conserved C-terminal domain (TOE domain) around a serine as one potential phosphorylation site. Phosphorylation of ACS of type II by casein kinase 1.8 (CK1.8) leads to interaction with the E3 ubiquitin ligase Ethylene Overproducer 1 (ETO1) and ETO1-like (EOL), which results in degradation of ACS via 26S proteasome [88]. Binding of ACS by 14-3-3 protein represses this interaction and destabilizes ETO1/EOL [88]. In maize, ZmACS6-deficient mutants showed an inhibition of drought-induced leaf senescence, an increased transpiration rate and higher CO2 assimilation under drought stress conditions [94]. Variations in StACS3 may be interesting for drought tolerance breeding in potato as ACS represents a potentially rate-limiting enzyme in ethylene biosynthesis. Ethylene invokes responses in both young and mature potato leaves [18].
Ethylene acts as central integrator linking and reprogramming complex stress-adaptive signaling cascades through activation of downstream ERFs [95]. ERFs show conserved binding preference for the dehydration responsive element (DRE) (CCGAC) or the GCC motif (GCCGCC) in promoter sequences that control expression of abiotic and biotic stress-responsive genes [96,97,98,99]. Some ERFs can also bind to both motives, DRE and GCC.
A genome-wide study of ERFs in potato identified 155 genes, which could be mainly distributed into two groups, the CBF/DREB and the ERF group [100]. Comparison of RNA-seq data revealed that most ERF genes are controlled by abiotic stress conditions (drought, salt, low temperature, ABA treatment). Overexpression of StDREB1/StERF186, as well as of StDREB2/StERF137, improves tolerance to salt stress in transgenic potatoes [101,102]. Overexpression of the tomato jasmonate and ethylene response factor 1 (JERF1) results in enhanced drought tolerance and ABA biosynthesis in tobacco and rice [103]. The ERF gene (PGSC0003DMG400002185) significantly associated with drought sensitivity in our SSR study corresponds to StERF8 [100]. As the SSR (CCA)n(ACA)k is located in exon 1, the observed allelic variation will have an effect on the amino acid sequence of the encoded transcription factor. Further studies will be necessary to reveal the specific effects of these allelic variations of StERF8 for drought tolerance in potato.

5. Conclusions

In potato, the candidate gene approach proved to be efficient to derive new functional SSR primer combinations and to identify SSR markers detecting allelic variations significantly associated with drought tolerance or sensitivity. Future testing of these SSR-markers in a larger panel of potato cultivars from all over the world will give a better idea about their versatility. Our association studies identified allelic differences in five of the candidate genes (aldehyde dehydrogenase, PGSC0003DMG400034597; 1-aminocyclopropane-1-carboxylate synthase 3 PGSC0003DMG400021426; ethylene responsive transcription factor, PGSC0003DMG400002185; poly (ADP-ribose) glycohydrolase, PGSC0003DMG400029361; protein phosphatase 2C, PGSC0003DMG400011321) as significantly associated with drought sensitivity. Further investigations will elucidate the role of these allelic variations in drought response of cultivated tetraploid potato and will give a better insight into mechanisms involved in drought tolerance in Solanaceae. The identified SSR markers are now available for genetic diversity studies and for marker-assisted selection in potato breeding programs. Interestingly, in our investigation breeding for drought tolerance in potato seems to require elimination of certain alleles from drought sensitive cultivars.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12040494/s1, Table S1. Characterization of potato cultivars. Table S2. Field trials for ranking the potato cultivars using DRYM as drought stress index. Table S3. Overview over all 154 derived SSR primer combinations. Table S4. Overview of 19 AFLP primer combinations used in the UPGMA dendrogram (Figure S1). Table S5. Linkage group specific SSR primer combinations. Table S6. Selected 23 transcripts identified in comparative transcriptome analysis [39]. Table S7. Characterization of 75 polymorphic SSR primer combinations. Table S8. Association studies of the SSR markers applying the Fisher Exact Test in the potato association panel (1t most drought tolerant cultivar, 34t most drought sensitive cultivar). Figure S1. Population structure in the investigated association panel of 34 starch potato cultivars based on 54 linkage group specific SSR primer combinations and 19 AFLP primer combinations. The UPGMA dendrogram based on the Jaccard’s coefficient is shown. The potato cultivars are shown with the drought tolerance rank from 1t = most tolerant to 34t = most sensitive.

Author Contributions

C.S. development and testing of SSRs derived from candidate genes, AFLP analysis and linkage group specific SSRs for populations structure analysis using NTSYS, statistical analysis (exact Fisher Test, PICs and Heterozygosity calculations); C.T.K., L.M. development and testing of SSRs; R.H. design and supervision of the project part in VALDIS Trost, writing of the manuscript, planning of the marker development, selection of candidate genes; K.K., S.S., R.P., H.S., and M.H. organized, coordinated and performed the field and greenhouse trials, ranked the potato varieties; S.H. minor participation in draft manuscript writing; J.K., D.K.H., E.Z., H.S. provided data from the transcript analysis. All authors read, corrected and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Work on drought tolerance in potato was funded in our laboratory by the Bundesmini-sterium für Ernährung und Landwirtschaft, Germany, (grant No. 22014713) as well as by the Gemeinschaft zur Förderung von Pflanzeninnovation e.V. (GFPi), Germany, (K 79/14 NR) and accompanied by the Fachagentur Nachwachsende Rohstoffe, Germany.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Genomic sequences used are available in EnsemblPlants (Potato Genome Sequencing Consortium 2011), NCBI, Phytozome, Sol Genomics Network and Spud DB Potato Genomics Resource.

Acknowledgments

We like to acknowledge Monja Sundt and Hilke Brandt at the Institute of Biological Sciences, Plant Genetics, University of Rostock, for the excellent technical support in the work with the SSR analyses. We are also grateful to Bernd Truberg, Katja Muders and Josef Strahwald for providing and propagating the potato varieties.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Fahad, S.; Bajwa, A.A.; Nazir, U.; Anjum, S.A.; Farooq, A.; Zohaib, A.; Sadia, S.; Nasim, W.; Adkins, S.; Saud, S.; et al. Crop production under drought and heat stress: Plant responses and management options. Front. Plant Sci. 2017, 8, 1147. [Google Scholar] [CrossRef] [Green Version]
  2. Fleury, D.; Jefferies, S.; Kuchel, H.; Langridge, P. Genetic and genomic tools to improve drought tolerance in wheat. J. Exp. Bot. 2010, 61, 3211–3222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Legay, S.; Lefèvre, I.; Lamoureux, D.; Barreda, C.; Luz, R.T.; Gutierrez, R.; Quiroz, R.; Hoffmann, L.; Hausman, J.-F.; Bonierbale, M.; et al. Carbohydrate metabolism and cell protection mechanisms differentiate drought tolerance and sensitivity in advanced potato clones (Solanum tuberosum L.). Funct. Integr. Genom. 2011, 11, 275–291. [Google Scholar] [CrossRef] [PubMed]
  4. Leisner, C.P. Review: Climate change impacts on food security—Focus on perennial cropping systems and nutritional value. Plant Sci. 2020, 293, 110412. [Google Scholar] [CrossRef] [PubMed]
  5. Raza, A.; Razzaq, A.; Mehmood, S.S.; Zou, Y.; Zhang, X.; Lv, Y.; Xu, J. Impact of climate change on crops adaptation and strategies to tackle its outcome: A review. Plants 2019, 8, 34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Thornton, P.K.; Ericksen, P.J.; Herrero, M.; Challinor, A.J. Climate variability and vulnerability to climate change: A review. Global Chang. Biol. 2014, 20, 3313–3328. [Google Scholar] [CrossRef] [PubMed]
  7. Adee, E.; Roozeboom, K.; Balboa, G.R.; Schlegel, A.; Ciampitti, I.A. Drought-tolerant corn hybrids yield more in drought-stressed environments with no penalty in non-stressed environments. Front. Plant Sci. 2016, 7, 1534. [Google Scholar] [CrossRef] [Green Version]
  8. He, Z.; Zhong, J.; Sun, X.; Wang, B.; Terzaghi, W.; Dai, M. The maize ABA receptors ZMPYL8, 9, and 12 facilitate plant drought resistance. Front. Plant Sci. 2018, 9, 422. [Google Scholar] [CrossRef] [PubMed]
  9. Tumova, L.; Tarkowska, D.; Rehorova, K.; Markova, H.; Kocova, M.; Rothova, O.; Cecetka, P.; Hola, D. Drought-tolerant and drought-sensitive genotypes of maize (Zea mays L.) differ in contents of endogenous brassinosteroids and their drought-induced changes. PLoS ONE 2018, 13, e0197870. [Google Scholar] [CrossRef] [Green Version]
  10. Mwadzingeni, L.; Shimmelis, H.; Dube, E.; Laing, M.D.; Tsilo, T.J. Breeding wheat for drought tolerance: Progress and technologies. J. Integr. Agricult. 2016, 15, 935–943. [Google Scholar] [CrossRef] [Green Version]
  11. Lawas, L.M.F.; Shi, W.; Yoshimoto, M.; Hasegawa, T.; Hincha, D.K.; Zuther, E.; Jagadish, S.V.K. Combined drought and heat stress impact during flowering and grain filling in contrasting rice cultivars grown under field conditions. Field Crop. Res. 2018, 229, 66–77. [Google Scholar] [CrossRef]
  12. Oladosu, Y.; Rafii, M.Y.; Samuel, C.; Fatai, A.; Magaji, U.; Kareem, I.; Kamarudin, Z.S.; Muhammad, I.; Kolapo, K. Drought resistance in rice from conventional to molecular breeding: A review. Int. J. Mol. Sci. 2019, 20, 3519. [Google Scholar] [CrossRef] [Green Version]
  13. Das, G.; Patra, J.K.; Baek, K.H. Insight into MAS: A molecular tool for development of stress resistant and quality of rice through gene stacking. Front. Plant Sci. 2017, 8, 985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Krannich, C.T.; Maletzki, L.; Kurowsky, C.; Horn, R. Network candidate genes in breeding for drought tolerant crops. Int. J. Mol. Sci. 2015, 16, 16378–16400. [Google Scholar] [CrossRef]
  15. Joshi, R.; Wani, S.H.; Singh, B.; Bohra, A.; Dar, Z.A.; Lone, A.A.; Pareek, A.; Singla-Pareek, S.L. Transcription factors and plants response to drought stress: Current understanding and future directions. Front. Plant Sci. 2016, 7, 1029. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Shinozaki, K.; Yamaguchi-Shinozaki, K. Gene networks involved in drought stress response and tolerance. J. Exp. Bot. 2007, 58, 221–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Singh, D.; Laxmi, A. Transcriptional regulation of drought response: A tortuous network of transcriptional factors. Front. Plant Sci. 2015, 6, 895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Obidiegwu, J.E.; Bryan, G.J.; Jones, H.G.; Prashar, A. Coping with drought: Stress and adaptive responses in potato and perspectives for improvement. Front. Plant Sci. 2015, 6, 542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Schafleitner, R.; Rosales, R.O.G.; Gaudin, A.; Aliaga, C.A.A.; Martinez, G.N.; Marca, L.R.T.; Bolivar, L.A.; Delgado, F.M.; Simon, R.; Bonierbale, M. Capturing candidate drought tolerance traits in two native Andean potato clones by transcription profiling of field grown plants under water stress. Plant Physiol. Biochem. 2007, 45, 673–690. [Google Scholar] [CrossRef]
  20. Dahal, K.; Li, X.Q.; Tai, H.; Creelman, A.; Bizimungu, B. Improving potato stress tolerance and tuber yield under a climate change scenario—A current overview. Front Plant Sci. 2019, 10, 563. [Google Scholar] [CrossRef] [PubMed]
  21. Hijmans, R.J. The effect of climate change on global potato production. Am. J. Potato Res. 2003, 80, 271–280. [Google Scholar] [CrossRef]
  22. Raymundo, R.; Asseng, S.; Robertson, R.; Petsakos, A.; Hoogenboom, G.; Quiroz, R.; Hareau, G.; Wolf, J. Climate change impact on global potato production. Eur. J. Agron. 2018, 100, 87–98. [Google Scholar] [CrossRef]
  23. Anithakumari, A.M.; Nataraja, K.N.; Visser, R.G.F.; van der Linden, C.G. Genetic dissection of drought tolerance and recovery potential by quantitative trait locus mapping of a diploid potato population. Mol. Breed. 2012, 30, 1413–1429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Sprenger, H.; Rudack, K.; Schudoma, C.; Neumann, A.; Seddig, S.; Peters, R.; Zuther, E.; Kopka, J.; Hincha, D.K.; Walther, D.; et al. Assessment of drought tolerance and its potential yield penalty in potato. Funct. Plant Biol. 2015, 42, 655–667. [Google Scholar] [CrossRef] [PubMed]
  25. Rykaczewska, K. Impact of heat and drought stresses on size and quality of the potato yield. Plant Soil Environ. 2017, 63, 40–46. [Google Scholar]
  26. Aliche, E.B.; Oortwijn, M.; Theeuwen, T.P.J.M.; Bachem, C.W.B.; Van Eck, H.J.; Visser, R.G.F.; van der Linden, C.G. Genetic mapping of tuber size distribution and marketable yield under drought stress in potatoes. Euphytica 2019, 251, 186. [Google Scholar] [CrossRef] [Green Version]
  27. Zarzynska, K.; Boguszewska-Mankowska, D.; Nosalewicz, A. Differences in size and architecture of the potato cultivars root system and their tolerance to drought stress. Plant Soil Environ. 2017, 63, 159–164. [Google Scholar]
  28. Hill, D.; Nelson, D.; Hammond, J.; Bell, L. Morphophysiology of potato (Solanum tuberosum) in repsonse to drought stress: Paving the way forward. Front. Plant Sci. 2021, 11, 597554. [Google Scholar] [CrossRef] [PubMed]
  29. Aliche, E.B.; Oortwijn, M.; Theeuwen, T.P.J.M.; Bachem, C.W.B.; Visser, R.G.F.; van der Linden, C.G. Drought response in field grown potatoes and the interactions between canopy growth and yield. Agric. Water Manag. 2018, 206, 20–30. [Google Scholar] [CrossRef]
  30. Boguszewska-Mankowska, D.; Pieczynski, M.; Wyrzykowska, A.; Kalaji, H.M.; Sieczko, L.; Szweykowska-Kulinska, Z.; Zagdanska, B. Divergent strategies displayed by potato (Solanum tuberosum L.) cultivars to cope with soil drought. J. Agron. Crop Sci. 2018, 204, 13–30. [Google Scholar] [CrossRef] [Green Version]
  31. Sprenger, H.; Erban, A.; Seddig, S.; Rudack, K.; Thalhammer, A.; Le, M.Q.; Walther, D.; Zuther, E.; Köhl, K.I.; Kopka, J.; et al. Metabolite and transcript markers for the prediction of potato drought tolerance. Plant Biotech. J. 2018, 16, 939–950. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ambrosone, A.; Batelli, G.; Bostan, H.; D’Agostino, N.; Chiusano, M.L.; Perrotta, G.; Leone, A.; Grillo, S.; Costa, A. Distinct gene networks drive differential response to abrupt or gradual water deficit in potato. Gene 2017, 597, 30–39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Bündig, C.; Vu, T.H.; Meise, P.; Seddig, S.; Schum, A.; Winkelmann, T. Variability in osmotic stress tolerance of starch potato genotypes (Solanum tuberosum L.) as revealed by an in vitro screening: Role of proline, osmotic adjustment and drought resposne in pot trials. J. Agron. Crop Sci. 2016, 203, 206–218. [Google Scholar] [CrossRef]
  34. Moon, K.B.; Ahn, D.J.; Park, J.S.; Jung, W.Y.; Cho, H.S.; Kim, H.R.; Jeon, J.H.; Park, Y.I.; Kim, H.S. Transcriptome profiling and characterization of drought-tolerant potato plant (Solanum tuberosum L.). Mol. Cells 2018, 41, 979–992. [Google Scholar] [PubMed]
  35. Fischer, R.A.; Maurer, R. Drought resistance in spring wheat cultivars. I Grain yield responses. Aust. J. Agric. Res. 1978, 29, 897–912. [Google Scholar] [CrossRef]
  36. Fernandez, G.C.J. Effective selection criteria for assessing stress tolerance. In Adaptation of Food Crops to Temperature and Water Stress; Kuo, C.G., Ed.; Asian Vegetable Research and Development Center: Shanhuan, China, 1992; pp. 257–270. [Google Scholar]
  37. Chen, Y.; Li, C.; Yi, J.; Yang, Y.; Lei, C.; Gong, M. Transcriptome response to drought, rehydration and re-dehydration in potato. Int. J. Mol. Sci. 2020, 21, 159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. 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, e012804. [Google Scholar] [CrossRef]
  39. Sprenger, H.; Kurowsky, C.; Horn, R.; Erban, A.; Seddig, S.; Rudack, K.; Fischer, A.; Walther, D.; Zuther, E.; Köhl, K.; et al. The drought response of potato reference cultivars with contrasting tolerance. Plant Cell Environ. 2016, 39, 2370–2389. [Google Scholar] [CrossRef] [PubMed]
  40. Yang, X.; Liu, J.; Xu, J.; Duan, S.; Wang, Q.; Li, G.; Jin, L. Transcriptome profiling reveals effects of drought stress on gene expression in diploid potato genotype P3-198. Int. J. Mol. Sci. 2019, 20, 852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Pieczynski, M.; Wyrzykowska, A.; Milanowska, K.; Boguszewska-Mankowska, D.; Zagdanska, B.; Karlowski, W.; Jarmolowski, A.; Szweykowska-Kulinska, Z. Genomewide identification of genes involved in the potato response to drought indicates functional evolutionary conservation with Arabidopsis plants. Plant Biotech. J. 2018, 16, 603–614. [Google Scholar] [CrossRef] [Green Version]
  42. Spooner, D.M.; Nunez, J.; Trujillo, G.; Herrera Mdel, R.; Guzman, F.; Ghislain, M. Extensive simple sequence repeat genotyping of potato landraces supports a major reevaluation of their gene pool structure and classification. Proc. Natl. Acad. Sci. USA 2007, 104, 19398–19403. [Google Scholar] [CrossRef] [Green Version]
  43. Sharma, S.K.; Bolser, D.; de Boer, J.; Sonderkaer, 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, 2031–2047. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Chambers, G.F.; MacAvoy, E.S. Microsatellites: Consensus and controversy. Comp. Biochem. Physiol. 2000, 126, 455–476. [Google Scholar] [CrossRef]
  45. Feingold, S.; Lloyd, J.; Norero, N.; Bonierbale, M.; Lorenzen, J. Mapping and characterization of new EST-derived microsatellites for potato (Solanum tuberosum L.). Theor. Appl. Genet. 2005, 111, 456–466. [Google Scholar] [CrossRef]
  46. Ghislain, M.; Spooner, D.M.; Rodríguez, F.; Villamón, F.; Núñez, J.; Vásquez, C.; Waugh, R.; Bonierbale, M. Selection of highly informative and user-friendly microsatellites (SSRs) for genotyping of cultivated potato. Theor. Appl. Genet. 2004, 108, 881–890. [Google Scholar] [CrossRef] [PubMed]
  47. Milbourne, D.; Meyer, R.C.; Collins, A.J.; Ramsay, L.D.; Gebhardt, C.; Waugh, R. Isolation, characterisation and mapping of simple sequence repeat loci in potato. Mol. Gen. Genet. 1998, 259, 233–245. [Google Scholar] [CrossRef] [PubMed]
  48. Provan, J.; Powell, W.; Waugh, R. Microsatellite analysis of relationships within cultivated potato (Solanum tuberosum). Theor. Appl. Genet. 1996, 92, 1078–1084. [Google Scholar] [CrossRef]
  49. Tang, J.; Baldwin, S.J.; Jacobs, J.M.; van der Linden, C.G.; Voorrips, R.E.; Am Leunissen, J.; van Eck, H.; Vosman, B. Large-scale identification of polymorphic microsatellites using an in silico approach. BMC Bioinform. 2008, 9, 374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Veilleux, R.E.; Shen, L.Y.; Paz, M.M. Analysis of the genetic composition of anther-derived potato by randomly amplified polymorphic DNA and simple sequence repeats. Genome 1995, 38, 1153–1162. [Google Scholar] [CrossRef]
  51. Ghislain, M.; Núñez, J.; del Rosario Herrera, M.; Pignataro, J.; Guzman, F.; Bonierbale, M.; Spooner, D.M. Robust and highly informative microsatellite-based genetic identity kit for potato. Mol. Breed. 2009, 23, 377–388. [Google Scholar] [CrossRef] [Green Version]
  52. Bali, S.; Sathuvalli, V.; Brown, C.; Novy, R.; Ewing, L.; Debons, J.; Douches, D.; Coombs, J.; Navarre, D.; Withworth, J.; et al. Genetic fingerprinting of potato varieties from the Northwest potato variety development program. Am. J. Potato Res. 2017, 94, 54–63. [Google Scholar] [CrossRef]
  53. Monte, M.N.; Burusco, M.F.R.; Carboni, M.F.; Castellote, M.A.; Sucar, S.; Norero, N.S.; Colman, S.L.; Massa, G.A.; Colavita, M.L.; Feingold, S.E. Genetic diversity in Argentine Andean potatoes by means of functional markers. Am. J. Potato Res. 2018, 95, 286–300. [Google Scholar] [CrossRef]
  54. Song, X.; Zhang, C.; Li, Y.; Feng, S.; Yang, Q.; Huang, S. SSR analysis of genetic diversity among 192 diploid potato cultivars. Hort. Plant J. 2016, 2, 163–171. [Google Scholar] [CrossRef] [Green Version]
  55. Schönhals, E.M.; Ortega, F.; Barandalla, L.; Aragones, A.; Ruiz de Galarreta, J.I.; Liao, J.-C.; Sanetomo, R.; Walkemeier, B.; Tacke, E.; Ritter, E.; et al. Identification and reproducibility of diagnostic DNA markers for tuber starch and yield optimization in a novel association mapping population. Theor. Appl. Genet. 2016, 129, 767–785. [Google Scholar] [CrossRef] [Green Version]
  56. Wang, Y.; Rahsid, M.A.R.; Li, X.; Yao, C.; Lu, L.; Bai, J.; Li, Y.; Xu, N.; Ynag, Q.; Zhang, L.; et al. Collection and evaluation of genetic diversity and population structure of potato landraces and varieties in China. Front. Plant Sci. 2019, 10, 139. [Google Scholar] [CrossRef] [Green Version]
  57. Danan, S.; Veyrieras, J.-B.; Lefebvre, V. Construction of a potato consensus map and QTL meta-analysis offer new insights into the genetic architecture of late blight resistance and plant maturity traits. BMC Plant Biol. 2011, 11, 16. [Google Scholar] [CrossRef] [Green Version]
  58. Watanabe, K. Potato genetics, genomics, and applications. Breed. Sci. 2015, 65, 53–68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Doyle, J.J.; Doyle, J.L. Isolation of plant DNA from fresh tissue. Focus 1990, 12, 13–15. [Google Scholar]
  60. Potato Genome Sequencing Consortium. Genome sequence and analysis of the tuber crop potato. Nature 2011, 475, 189–195. [Google Scholar] [CrossRef] [Green Version]
  61. Sajer, O.; Scorza, R.; Dardick, C.; Zhebentyayeva, T.; Abbott, A.G.; Horn, R. Development of sequence-tagged site markers linked to the pillar growth type in peach (Prunus persica). Plant Breed. 2012, 131, 186–192. [Google Scholar] [CrossRef]
  62. Vos, P.; Hogers, R.; Bleeker, M.; Reijans, M.; van de Lee, T.; Hornes, M.; Frijters, A.; Pot, J.; Peleman, J.; Kuiper, M.; et al. AFLP: A new technique for DNA fingerprinting. Nucleic Acids Res. 1995, 23, 4407–4414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Rohlf, F.J. NTSYSpc—Numerical Taxonomy and Multivariate Analysis System; Version 2.1, User Guide, Exeter Software; Applied Biostatistics Inc.: Setauket, NY, USA, 2004; pp. 1–43. [Google Scholar]
  64. Van Muijen, D.; Anithakumari, A.M.; Maliepaard, C.; Visser, R.G.; van der Linden, C.G. Systems genetics reveals key genetic elements of drought induced gene regulation in diploid potato. Plant Cell Environ. 2016, 39, 1895–1908. [Google Scholar] [CrossRef] [PubMed]
  65. Weiner, J.J.; Peterson, F.C.; Volkman, B.F.; Cutler, S.R. Structural and functional insights into core ABA signaling. Curr. Opin. Plant Biol. 2010, 13, 495–502. [Google Scholar] [CrossRef] [Green Version]
  66. Umezawa, T.; Sugiyama, N.; Mizoguchi, M.; Hayashi, S.; Myouga, F.; Yamaguchi-Shinozaki, K.; Ishihama, Y.; Hirayama, T.; Shinozaki, K. Type 2C protein phosphatases directly regulate abscisic acid-activated protein kinases in Arabidopsis. Proc. Natl. Acad. Sci. USA 2009, 106, 17588–17593. [Google Scholar] [CrossRef] [Green Version]
  67. Brocard, I.M.; Lynch, T.J.; Finkelstein, R.R. Regulation and role of the Arabidopsis abscisic acid-insensitive 5 gene in abscisic acid, sugar, and stress response. Plant Physiol. 2002, 129, 1533–1543. [Google Scholar] [CrossRef] [Green Version]
  68. Munemasa, S.; Hauser, F.; Park, J.; Waadt, R.; Brandt, B.; Schroeder, J.I. Mechanism of abscisic acid-mediated control of stomatal aperture. Curr. Opin. Plant Biol. 2015, 28, 154–162. [Google Scholar] [CrossRef] [Green Version]
  69. Zhang, A.; Ren, H.M.; Tan, Y.Q.; Qi, G.N.; Yao, F.Y.; Wu, G.L.; Yang, L.W.; Hussain, J.; Sun, S.J.; Wang, Y.F. S-Type anion channels SLAC1 and SLAH3 function as essential negative regulators of inward K+ channels and stomatal opening in Arabidopsis. Plant Cell 2016, 28, 949–965. [Google Scholar] [CrossRef] [PubMed]
  70. Merlot, S.; Gosti, F.; Guerrier, D.; Vavasseur, A.; Giraudat, J. The ABI1 and ABI2 protein phosphatases 2C act in a negative feedback regulatory loop of the abscisic acid signalling pathway. Plant J. 2001, 25, 295–303. [Google Scholar] [CrossRef]
  71. Zhang, K.W.; Xia, X.Y.; Zhang, Y.Y.; Gan, S.S. An ABA-regulated and Golgi-localized protein phosphatase controls water loss during leaf senescence in Arabidopsis. Plant J. 2012, 69, 667–678. [Google Scholar] [CrossRef] [PubMed]
  72. Yu, X.; Han, J.; Wang, E.; Xiao, J.; Hu, R.; Yang, G.; He, G. Genome-wide identification and homoeologous expression analysis of PP2C genes in wheat (Triticum aestivum L.). Front. Genet. 2019, 10, 561. [Google Scholar] [CrossRef] [Green Version]
  73. Umezawa, T.; Nakashima, K.; Miyakawa, T.; Kuromori, T.; Tanokura, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Molecular basis of the core regulatory network in ABA responses: Sensing, signaling and transport. Plant Cell Physiol. 2010, 51, 1821–1839. [Google Scholar] [CrossRef]
  74. Das, K.; Roychoudhury, A. Reactive oxygen sypecies (ROS) and response of antioxidants as ROS-scavengers during environmental stress in plants. Front. Plant Sci. 2014, 2, 53. [Google Scholar]
  75. Singh, S.; Brocker, C.; Koppaka, V.; Chen, Y.; Jackson, B.C.; Matsumoto, A.; Thompson, D.C.; Vasiliou, V. Aldehyde dehydrogenases in cellular responses to oxidative/ electrophilic stress. Free Radic. Biol. Med. 2013, 56, 89–101. [Google Scholar] [CrossRef] [Green Version]
  76. Kirch, H.H.; Schlingensiepen, S.; Kotchoni, S.; Sunkar, R.; Bartels, D. Detailed expression analysis of selected genes of the aldehyde dehydrogenase (ALDH) gene superfamily in Arabidopsis thaliana. Plant Mol. Biol. 2005, 57, 315–332. [Google Scholar] [CrossRef]
  77. Kotchoni, S.O.; Kuhns, C.; Ditzer, A.; Kirch, H.H.; Bartels, D. Over-expression of different aldehyde dehydrogenase genes in Arabidopsis thaliana confers tolerance to abiotic stress and protects plants against lipid peroxidation and oxidative stress. Plant Cell Environ. 2006, 29, 1033–1048. [Google Scholar] [CrossRef]
  78. Nisa, M.-U.; Huang, Y.; Benhamed, M.; Raynaud, C. The plant DNA damage response: Signaling pathway leading to growth inhibition and putative role in response to stress conditions. Front. Plant Sci. 2019, 10, 653. [Google Scholar] [CrossRef] [Green Version]
  79. Roy, S. Maintenance of genome stability in plants: Repairing DNA double strand breaks and chromatin structure stability. Front. Plant Sci. 2014, 5, 487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Yao, Y.; Bilichak, A.; Titov, V.; Golubov, A.; Kovalchuk, I. Genome stability of Arabidopsis atm, ku80 and rad51b mutants: Somatic and transgenerational responses to stress. Plant Cell Physiol. 2013, 54, 982–989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Lamb, R.S.; Citarelli, M.; Teotia, S. Functions of the poly(ADP-ribose) polymerase superfamily in plants. Cell Mol. Life Sci. 2012, 69, 175–189. [Google Scholar] [CrossRef] [PubMed]
  82. Rissl, D.; Peiter, E. Poly(ADP-ribose) polymerase in plants and their human counterparts: Parallels and pecularities. Int. J. Mol. Sci. 2019, 20, 1638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Hatakeyama, K.; Nemoto, Y.; Ueda, K.; Hayaishi, O. Purification and characterization of poly(ADP-ribose) glycohydrolase. Different modes of action on large and small poly(ADP-ribose). J. Biol. Chem. 1986, 261, 14902–14911. [Google Scholar] [CrossRef]
  84. Zhang, H.L.; Gu, Z.Y.; Wu, Q.; Yang, L.F.; Liu, C.F.; Ma, H.; Xia, Y.J.; Ge, X.C. Arabidopsis PARG1 is the key factor promoting cell survival among the enzymes regulating post-translational poly(ADP-ribosyl)ation. Sci. Rep. 2015, 5, 15892. [Google Scholar] [CrossRef]
  85. Adams-Phillips, L.; Briggs, A.G.; Bent, A.F. Disruption of poly(ADP-ribosyl)ation mechanisms alters responses of Arabidopsis to biotic stress. Plant Physiol. 2010, 152, 267–280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Iqbal, N.; Khan, N.A.; Ferrante, A.; Trivellini, A.; Francini, A.; Khan, M.I.R. Ethylene role in plant growth, development and senescence: Interaction with other phytohormones. Front. Plant Sci. 2017, 8, 475. [Google Scholar] [CrossRef] [Green Version]
  87. Houben, M.; Van de Poel, B. 1-Aminocyclopropane-1-carboxylic acid aoxidase (ACO): The enzyme that makes the plant hormone ethylene. Front. Plant Sci. 2019, 10, 695. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  88. Pattyn, J.; Vaughan-Hirsch, J.; Van de Poel, B. The regulation of ethylene biosynthesis: A complex multilevel control circuitry. New Phytol. 2021, 229, 770–782. [Google Scholar] [CrossRef]
  89. Winter, D.; Vinegar, B.; Nahal, H.; Ammar, R.; Wilson, G.V.; Provart, N.J. An “Electronic Fluorescent Pictograph” browser for exploring and analyzing large scale biological data sets. PLoS ONE 2007, 8, e718. [Google Scholar] [CrossRef] [PubMed]
  90. Yamagami, T.; Tsuchisaka, A.; Yamada, K.; Haddon, W.F.; Harden, L.A.; Theologis, A. Biochemical diversity among the 1-amino-cyclopropane-1-carboxylate synthase isozymes encoded by the Arabidopsis gene family. J. Biol. Chem. 2003, 278, 49102–49112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Tsuchisaka, A.; Yu, G.; Jin, H.; Alonso, J.M.; Ecker, J.R.; Zhang, X.; Gao, S.; Theologis, A. A combinatorial interplay among the 1-aminocyclopropane-1-carboxylate isoforms regulates ethylene biosynthesis in Arabidopsis thaliana. Genetics 2009, 183, 979–1003. [Google Scholar] [CrossRef] [Green Version]
  92. Chae, H.S.; Kieber, J.J. Eto Brute? Role of ACS turnover in regulating ethylene biosynthesis. Trends Plant Sci. 2005, 10, 291–296. [Google Scholar] [CrossRef]
  93. Eun, H.-D.; Ali, S.; Jung, H.; Kim, K.; Kim, W.-C. Profiling of ACC synthase gene (ACS11) expression in Arabidopsis induced by abiotic stress. Appl. Biol. Chem. 2019, 62, 42. [Google Scholar] [CrossRef] [Green Version]
  94. Young, T.E.; Meeley, R.B.; Gallie, D.R. ACC synthase expression regulates leaf performance and drought tolerance in maize. Plant J. 2004, 40, 813–825. [Google Scholar] [CrossRef] [PubMed]
  95. Moschen, S.; Di Rienzo, J.A.; Higgins, J.; Tohge, T.; Watanabe, M.; Gonzalez, S.; Rivarola, M.; García-García, F.; Dopazo, J.; Hopp, H.E.; et al. Integration of transcriptomic and metabolic data reveals hub transcription factors involved in drought stress response in sunflower (Helianthus annuus L.). Plant Mol. Biol. 2017, 94, 549–564. [Google Scholar] [CrossRef] [PubMed]
  96. Berrocal-Lobo, M.; Molina, A.; Solano, R. Constitutive expression of ETHYLENE-RESPONSE-FACTOR1 in Arabidopsis confers resistance to several necrotrophic fungi. Plant J. 2002, 29, 23–32. [Google Scholar] [CrossRef] [PubMed]
  97. Cheng, M.C.; Liao, P.M.; Kuo, W.W.; Lin, T.P. The Arabidopsis ETHYLENE RESPONSE FACTOR1 regulates abiotic stress-responsive gene expression by binding to different cis-acting elements in response to different stress signals. Plant Physiol. 2013, 162, 1566–1582. [Google Scholar] [CrossRef] [Green Version]
  98. Müller, M.; Munné-Bosch, S. Ethylene response factors: A key regulatory hub in hormone and stress signaling. Plant Physiol. 2015, 169, 32–41. [Google Scholar] [CrossRef] [Green Version]
  99. Solano, R.; Stepanova, A.; Chao, Q.M.; Ecker, J.R. Nuclear events in ethylene signaling: A transcriptional cascade mediated by ETHYLENE-INSENSITIVE3 and ETHYLENE-RESPONSE-FACTOR1. Gene Dev. 1998, 12, 3703–3714. [Google Scholar] [CrossRef] [Green Version]
  100. Charfeddine, M.; Saidi, M.N.; Charfeddine, S.; Hammami, A.; Bouzid, R.G. Genome-wide analysis and expression profiling of the ERF transcription factor family in potato (Solanum tuberosum L.). Mol. Biotechnol. 2015, 57, 348–358. [Google Scholar] [CrossRef]
  101. Bouaziz, D.; Pirrello, J.; Amor, H.B.; Hammami, A.; Charfeddine, M.; Dhieb, A.; Bouzayen, M.; Gargouri-Bouzid, R. Ectopic expression of dehydration responsive element binding proteins (StDREB2) confers higher tolerance to salt stress in potato. Plant Physiol. Biochem. 2012, 60, 98–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Bouaziz, D.; Pirrello, J.; Charfeddine, M.; Hammami, A.; Jbir, R.; Dhieb, A.; Bouzayen, M.; Gargouri-Bouzid, R. Overexpression of StDREB1 transcription factor increases tolerance to salt in transgenic potato plants. Mol. Biotechnol. 2013, 54, 803–817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Zhang, Z.J.; Li, F.; Li, D.J.; Zhang, H.W.; Huang, R.F. Expression of ethylene response factor JERF1 in rice improves tolerance to drought. Planta 2010, 232, 765–774. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Overview about processes involved in response to water deficiency showing the selected candidate genes for drought tolerance used in this study for the development of SSR markers. Abbreviations: ABA: Abscisic acid, ABA8H: ABA 8′-hydroxylase, ACO: 1-Aminocyclopropane-1-carboxylic acid oxidase, ACS: 1-Aminocyclopropane-1-carboxylic acid synthase, ALDH: Aldehyde dehydrogenase, APX: Ascorbate peroxidase, APY: Apyrase, ARF: Auxin response factor, BA: β-amylase, BDGEH: β-D-glucan exohydrolase, BADH: Betaine aldehyde dehydrogenase, BSDRP4: Bacterial spot disease resistance protein 4, BEAF: Boundary element associated factor, BED: BEAF and DREF, CHO: Carbohydrate, CTR1K: Constitutive Triple Response1 kinase, DBP: DNA binding protein, P5CS: Delta 1-pyrroline-5-carboxylate synthetase, DREB: Dehydration-Responsive Element-Binding, DREF: DNA replication-related element factor, EBF1: EIN3-binding F-box protein 1, EIN: Ethylene insensitive, EIL: Ethylene-Insensitive3-Like, EREBP: Ethylene-responsive element-binding protein, EREBPC2: Ethylene responsive element binding protein C2, ERF: Ethylene response factor, ERS 1: Ethylene response sensor 1, ET: Ethylene, ETO1: Ethylene-overproduction protein 1, ETP: EIN2 targeting protein, ETR1: Ethylene receptor 1, ETRHOM: Ethylene receptor homolog, EIX: Ethylene-inducing xylanase, FLS 2: Flagellin-sensing 2, GA3ox: Gibberellin 3-oxidase, GA: Gibberellin, GB: Glycine betaine, GR: Glutathione reductase, GST: Glutathione-S-transferase, GT: Glucosyltransferase, JA: Jasmonic acid, JA2: Jasmonic acid 2, JERF1: Jasmonate and Ethylene Response Factor 1, KAO: Ent-kaurenoic acid oxidase, KIN: Kinase, LEA5: Late embryogenesis abundant 5, LOX: Lipoxygenase, LRR: Leucine-rich repeat, MetRS: Methionyl-tRNA synthetase, MOCOS: Molybdenum Cofactor sulfurase, MRP-ABC: Multidrug resistance protein ATP binding cassette, NAM: No apical meristem, NBS: Nucleotide binding site, NI: Neutral invertase, NCED: 9-Cis-epoxycarotenoid dioxygenase, NDPK2: Nucleoside diphosphate protein kinase 2, P5CR: Pyrroline-5-carboxylate reductase, P5CS: Pyrroline-5-carboxylate synthase, PARG: Poly (ADP-ribose) glycohydrolase, PP2C: Protein phosphatase 2C, PRO: Proline, PRODH: Proline dehydrogenase, PYL4: Pyrabactin resistance-like 4, RPK: Receptor protein kinase, SBT4B: subtilase 4B, SCR: SCARECROW, SIPP: Soluble inorganic pyrophosphatase, SLSG6: S-locus specific glycoprotein S6, SN-1: Snaking-1, SOD: Superoxide dismutase, St: Solanum tuberosum, STPK: Serine/threonine protein kinase, TIR: Toll/interleukin-1 receptor-like protein, TNL-DR: TIR-NBS-LRR disease resistance, TPS1: Trehalose-6-phosphate synthase 1, TRE: Trehalose, UGGT: UDP-glucose glycosyltransferase, UGPase: UDP-glucose-1-phosphate uridylyltransferase, VPSP: Vacuolar protein sorting protein, WRKY: WRKY transcription factor, ZEP: Zeaxanthin epoxidase
Figure 1. Overview about processes involved in response to water deficiency showing the selected candidate genes for drought tolerance used in this study for the development of SSR markers. Abbreviations: ABA: Abscisic acid, ABA8H: ABA 8′-hydroxylase, ACO: 1-Aminocyclopropane-1-carboxylic acid oxidase, ACS: 1-Aminocyclopropane-1-carboxylic acid synthase, ALDH: Aldehyde dehydrogenase, APX: Ascorbate peroxidase, APY: Apyrase, ARF: Auxin response factor, BA: β-amylase, BDGEH: β-D-glucan exohydrolase, BADH: Betaine aldehyde dehydrogenase, BSDRP4: Bacterial spot disease resistance protein 4, BEAF: Boundary element associated factor, BED: BEAF and DREF, CHO: Carbohydrate, CTR1K: Constitutive Triple Response1 kinase, DBP: DNA binding protein, P5CS: Delta 1-pyrroline-5-carboxylate synthetase, DREB: Dehydration-Responsive Element-Binding, DREF: DNA replication-related element factor, EBF1: EIN3-binding F-box protein 1, EIN: Ethylene insensitive, EIL: Ethylene-Insensitive3-Like, EREBP: Ethylene-responsive element-binding protein, EREBPC2: Ethylene responsive element binding protein C2, ERF: Ethylene response factor, ERS 1: Ethylene response sensor 1, ET: Ethylene, ETO1: Ethylene-overproduction protein 1, ETP: EIN2 targeting protein, ETR1: Ethylene receptor 1, ETRHOM: Ethylene receptor homolog, EIX: Ethylene-inducing xylanase, FLS 2: Flagellin-sensing 2, GA3ox: Gibberellin 3-oxidase, GA: Gibberellin, GB: Glycine betaine, GR: Glutathione reductase, GST: Glutathione-S-transferase, GT: Glucosyltransferase, JA: Jasmonic acid, JA2: Jasmonic acid 2, JERF1: Jasmonate and Ethylene Response Factor 1, KAO: Ent-kaurenoic acid oxidase, KIN: Kinase, LEA5: Late embryogenesis abundant 5, LOX: Lipoxygenase, LRR: Leucine-rich repeat, MetRS: Methionyl-tRNA synthetase, MOCOS: Molybdenum Cofactor sulfurase, MRP-ABC: Multidrug resistance protein ATP binding cassette, NAM: No apical meristem, NBS: Nucleotide binding site, NI: Neutral invertase, NCED: 9-Cis-epoxycarotenoid dioxygenase, NDPK2: Nucleoside diphosphate protein kinase 2, P5CR: Pyrroline-5-carboxylate reductase, P5CS: Pyrroline-5-carboxylate synthase, PARG: Poly (ADP-ribose) glycohydrolase, PP2C: Protein phosphatase 2C, PRO: Proline, PRODH: Proline dehydrogenase, PYL4: Pyrabactin resistance-like 4, RPK: Receptor protein kinase, SBT4B: subtilase 4B, SCR: SCARECROW, SIPP: Soluble inorganic pyrophosphatase, SLSG6: S-locus specific glycoprotein S6, SN-1: Snaking-1, SOD: Superoxide dismutase, St: Solanum tuberosum, STPK: Serine/threonine protein kinase, TIR: Toll/interleukin-1 receptor-like protein, TNL-DR: TIR-NBS-LRR disease resistance, TPS1: Trehalose-6-phosphate synthase 1, TRE: Trehalose, UGGT: UDP-glucose glycosyltransferase, UGPase: UDP-glucose-1-phosphate uridylyltransferase, VPSP: Vacuolar protein sorting protein, WRKY: WRKY transcription factor, ZEP: Zeaxanthin epoxidase
Genes 12 00494 g001
Figure 2. Composition and distribution of microsatellites in 5′-UTR-regions, exons, introns and 3′-UTR-regions of the investigated 103 candidate genes. (a) SSR classes and distribution of SSR loci of all tested SSR markers; (b) SSR classes and distribution of SSR loci of polymorphic SSR markers in the potato association panel (34 tetraploid potato cultivars).
Figure 2. Composition and distribution of microsatellites in 5′-UTR-regions, exons, introns and 3′-UTR-regions of the investigated 103 candidate genes. (a) SSR classes and distribution of SSR loci of all tested SSR markers; (b) SSR classes and distribution of SSR loci of polymorphic SSR markers in the potato association panel (34 tetraploid potato cultivars).
Genes 12 00494 g002
Figure 3. Physical location of the candidate gene specific SSR markers according to the potato reference genome (The Potato Genome Consortium 2011).
Figure 3. Physical location of the candidate gene specific SSR markers according to the potato reference genome (The Potato Genome Consortium 2011).
Genes 12 00494 g003
Figure 4. Ranking of the 34 potato varieties of the association panel according to the drought tolerance. Drought tolerance index (DRYM) of 34 potato varieties against rank of cultivar according to drought tolerance (DRYM) (1 most drought tolerant, 34 most sensitive). The bars at the bottom of the figure indicate the grouping according to the means comparison (REGWQ, α = 0.1). Cultivars that are underlined by the same bar are not significantly different in DRYM, e.g., 1 (most tolerant group) is different from 2 to 4 that belong into the second most tolerant group.
Figure 4. Ranking of the 34 potato varieties of the association panel according to the drought tolerance. Drought tolerance index (DRYM) of 34 potato varieties against rank of cultivar according to drought tolerance (DRYM) (1 most drought tolerant, 34 most sensitive). The bars at the bottom of the figure indicate the grouping according to the means comparison (REGWQ, α = 0.1). Cultivars that are underlined by the same bar are not significantly different in DRYM, e.g., 1 (most tolerant group) is different from 2 to 4 that belong into the second most tolerant group.
Genes 12 00494 g004
Figure 5. SSR markers associated with drought sensitivity in an association panel of 34 mostly starch potato varieties ranked according to drought tolerance (1t most drought tolerant, 34t most sensitive). (a) SSR primer combination HRO_ACS3; (b) SSR primer combination HRO-PP2C_1; (c) SSR primer combination HRO_ALDH (only relevant bands are marked); (d) SSR primer combination HRO_ETRTF_5a; (e) SSR primer combination HRO_PARG_1A; M: 50–350 IRD-Marker; Marker bands are alphabetically labelled starting with the largest fragment. Bands associated with drought sensitivity are marked with asterisks.
Figure 5. SSR markers associated with drought sensitivity in an association panel of 34 mostly starch potato varieties ranked according to drought tolerance (1t most drought tolerant, 34t most sensitive). (a) SSR primer combination HRO_ACS3; (b) SSR primer combination HRO-PP2C_1; (c) SSR primer combination HRO_ALDH (only relevant bands are marked); (d) SSR primer combination HRO_ETRTF_5a; (e) SSR primer combination HRO_PARG_1A; M: 50–350 IRD-Marker; Marker bands are alphabetically labelled starting with the largest fragment. Bands associated with drought sensitivity are marked with asterisks.
Genes 12 00494 g005
Figure 6. Clustering of the 34 potato varieties based on the genetic similarity obtained by SSR analyses and selection with SSR markers. (a) UPGMA dendrogram of the 34 potato varieties (1t = most drought tolerant, 34t = most sensitive) based on all 75 SSR primer combinations derived from candidate genes for drought tolerance, Mantel coefficient r = 0.62; (b) Selection by applying four SSR markers associated with drought sensitivity. Black box= band is present, grey box= band is absent. The varieties are ordered according to the UPGMA dendrogram. Boxed varieties are selected by the presence of two marker bands either in the combination HRO_ACS3_D and HRO_PP2C_1_B or the combination HRO_ALDH_H and HRO-ALDH_N. Black Arrow head = falsely negative selected; Grey arrow head = falsely positive selected.
Figure 6. Clustering of the 34 potato varieties based on the genetic similarity obtained by SSR analyses and selection with SSR markers. (a) UPGMA dendrogram of the 34 potato varieties (1t = most drought tolerant, 34t = most sensitive) based on all 75 SSR primer combinations derived from candidate genes for drought tolerance, Mantel coefficient r = 0.62; (b) Selection by applying four SSR markers associated with drought sensitivity. Black box= band is present, grey box= band is absent. The varieties are ordered according to the UPGMA dendrogram. Boxed varieties are selected by the presence of two marker bands either in the combination HRO_ACS3_D and HRO_PP2C_1_B or the combination HRO_ALDH_H and HRO-ALDH_N. Black Arrow head = falsely negative selected; Grey arrow head = falsely positive selected.
Genes 12 00494 g006
Table 1. Scoring pattern of SSR markers significantly associated with drought tolerance in potato.
Table 1. Scoring pattern of SSR markers significantly associated with drought tolerance in potato.
Primer\Cultivar12345678910111213141516171819202122232425262728293031323334
HRO_ACS3_D0010111000000110011111100011111110
HRO_PP2C_1_B0100000101100000011101101011011100
HRO_ALDH_H1110000000000001001110010111110101
HRO_ETRTF_5a_D0110000100000010001101111101001110
HRO_PARGH_1A_B0000100100010000001111100111100001
Table 2. SSR Markers significantly associated with drought tolerance in potato.
Table 2. SSR Markers significantly associated with drought tolerance in potato.
SSR MarkerSizeLGGene ID PGSC0003Transcript ID PGSC0003Gene Annotationp-Value
HRO_ACS3_D173 bpLG 2DMG400021426DMT4000552031-aminocyclopropane-1-
carboxylate synthase 3
p = 0.0366
HRO_PP2C_1_B205 bpLG 1DMG400011321DMT400029441protein phosphatase 2Cp = 0.0366
HRO_ALDH_H184 bpLG 9DMG400034597DMT400085026aldehyde dehydrogenasep = 0.0366
HRO_ETRTF_5a_D217 bpLG 11DMG400002185DMT400005585ethylene responsive
transcription factor
p = 0.0366
HRO_PARG_1A_B147 bpLG 12DMG400029361DMT400075512poly (ADP-ribose)
glycohydrolase
p = 0.0324
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schumacher, C.; Krannich, C.T.; Maletzki, L.; Köhl, K.; Kopka, J.; Sprenger, H.; Hincha, D.K.; Seddig, S.; Peters, R.; Hamera, S.; et al. Unravelling Differences in Candidate Genes for Drought Tolerance in Potato (Solanum tuberosum L.) by Use of New Functional Microsatellite Markers. Genes 2021, 12, 494. https://doi.org/10.3390/genes12040494

AMA Style

Schumacher C, Krannich CT, Maletzki L, Köhl K, Kopka J, Sprenger H, Hincha DK, Seddig S, Peters R, Hamera S, et al. Unravelling Differences in Candidate Genes for Drought Tolerance in Potato (Solanum tuberosum L.) by Use of New Functional Microsatellite Markers. Genes. 2021; 12(4):494. https://doi.org/10.3390/genes12040494

Chicago/Turabian Style

Schumacher, Christina, Christoph Tim Krannich, Lisa Maletzki, Karin Köhl, Joachim Kopka, Heike Sprenger, Dirk Karl Hincha, Sylvia Seddig, Rolf Peters, Sadia Hamera, and et al. 2021. "Unravelling Differences in Candidate Genes for Drought Tolerance in Potato (Solanum tuberosum L.) by Use of New Functional Microsatellite Markers" Genes 12, no. 4: 494. https://doi.org/10.3390/genes12040494

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop