Next Article in Journal
Lifestyle and BMI Changes after the Release of COVID-19 Restrictions: Do Humans Go ‘Back to Normal’?
Previous Article in Journal
Genetic Adaptations of Soybean to Cold Stress Reveal Key Insights Through Transcriptomic Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Expression Analysis for Drought Tolerance in Early Stage of Potato Plant Development

1
Institute of Plant Biology and Biotechnology, Almaty 050040, Kazakhstan
2
Tanir Research Laboratory, Almaty 050060, Kazakhstan
*
Author to whom correspondence should be addressed.
Biology 2024, 13(11), 857; https://doi.org/10.3390/biology13110857
Submission received: 27 September 2024 / Revised: 15 October 2024 / Accepted: 22 October 2024 / Published: 23 October 2024

Simple Summary

Every year drought causes losses in yield of many crops including potato and causes significant economic problems. One of the ways to overcome drought issues is developing varieties that can withstand water scarcity better. In this study, we simulated osmotic stress using PEG-6000 and compared morphological characteristics of different varieties to identify contrasting varieties in drought tolerance in the early plant development stage. As a result, two contrasting varieties were analysed for gene expression, where clear differences in gene expression profiles were shown. These experiments allow to identify drought tolerance properties of variety in the early development stage of potato. Furthermore, comparative analysis showed that expression profiles of subset of drought-tolerance genes used in this study were also conserved across different crops.

Abstract

Drought has increasingly affected the yield of Solanum tuberosum L. (potato) every year over the last decade, posing serious economic problems for the global agricultural industry. Therefore, it is important to research drought tolerance in plants and obtain more robust varieties of crops. The aim of the present work was to study the expression of drought-upregulated genes in drought-tolerant and drought-sensitive varieties of potato. Bioreactors were used to identify whether each variety was drought-tolerant or drought-sensitive; then, expression analysis was performed according to the morphological characteristics of the plantlets in two different media: Murashige and Skoog (MS) medium and MS medium with 20% PEG-6000 to simulate osmotic stress. Based on the quantitative parameters of six initial varieties, two varieties were selected (Gala and Aksor) for further gene expression analysis. The expression of genes commonly upregulated in drought (ER24, TAS14, DREB147315, PP2C, 102605413 and NF-YC4) was higher in the drought-tolerant variety than in the sensitive one. Therefore, the expression of these genes can be used to determine the drought tolerance of a potato variety in vitro in the early plant development stage. Moreover, comparative analysis showed that some of the targeted genes used to identify drought tolerance in this study are conserved across different plant species.

1. Introduction

Drought years have become more frequent in the last decade [1], causing significant losses in crop yield, some estimating up to 35% for moderately drought-tolerant and up to 16% for drought-tolerant varieties [2]. Drought significantly increases the costs of agricultural production. Research on plants’ tolerance to drought is becoming increasingly important to minimize the impact of global climate change on agriculture [3]. Most crops are grown in water-deficient regions of the world [4]. In particular, in Kazakhstan, potatoes are grown in drought conditions in the southern foothill regions and over significant flat areas in the north of the country. Breeding potato plants for drought tolerance is important because climate change will increase the frequency of drought.
Conventional methods of breeding potatoes for drought tolerance are very labor-intensive, season-dependent, time-consuming and costly [5]. Vegetation indices, such as the leaf area index and aboveground biomass, have been found to be the main determinants of yield in tested varieties under drought conditions [6]. It should also be considered that early drought reduces the number of tubers per plant across all genotypes, and later drought starting from the moment of tuber initiation differs from variety to variety and is more pronounced in early varieties [7]. In addition, the ability of plants to recover after the elimination of soil drought should be considered, as this is a good indicator of their sensitivity to soil drought and allows the prediction of potato tuber yield [8]. For breeding programs, it is useful to identify physiological traits in potato varieties that are closely related but that differ in their tolerance to drought. However, evaluation based on morphological traits alone is often not effective. Therefore, various modern biotechnological methods are often used in breeding practices worldwide, with success.
In particular, in vitro cultivation methods are widely used to breed crops for tolerance to abiotic and biotic stresses. As a rule, drought stress conditions are modeled using polyethylene glycol (PEG) [9,10,11], including for potato [12,13]. Using PEG in in vitro culture to breed for tolerance to stress factors in the growing environment is efficient, as this approach can be used all year round, regardless of the season [14]. In the breeding process for drought tolerance, a high correlation between the results of in vitro experiments and field trials on soybean [15], rice [16], as well as for the plant growth and tuber yield of potato, has been shown [17].
In addition, the use of molecular markers has greatly facilitated the breeding process of almost all crops. At the same time, the number of molecular markers associated with certain breeding-valuable traits of agricultural crops is rapidly increasing. In particular, allelic differences between potato varieties were identified for five SSR markers, which showed a significant association with drought sensitivity. In all cases, an additional allele was predominant in the group of drought-sensitive varieties, indicating that selection against these alleles by marker-assisted crossing may confer drought tolerance to the hybrid [18]. Putative drought tolerance genes were differentially expressed in leaves of genotypes under water stress. The induction of protein phosphatase 2C gene was positively associated with yield maintenance under drought conditions. In addition, tolerant varieties expressed the transcription factor DREB to a greater extent than more sensitive varieties [6].
Drought stress causes changes in the expression of a number of genes that control morphological, physiological and biochemical traits at the level of transcription [19]. Transcriptomic changes in plants during drought stress provide insight into the mechanisms by which plants stabilize their metabolic processes to cope with drought conditions. In potatoes, understanding such drought-induced transcriptomic changes is crucial because prolonged field drought interferes with tuber formation and the developmental period of potatoes, which ultimately affects yield [20]. At the same time, it has been determined that drought stress turns on genes that are involved in the biosynthesis of endogenous hormones, which in turn, affect the expression of downstream genes associated with heat shock proteins [21]. Next-generation sequencing and single-molecule real-time sequencing technologies were used to study potato transcription profiles in response to drought stress modeling (20% PEG). This resulted in a significant improvement in the annotation of drought-induced genes in potatoes and provided a transcriptomic insight into the molecular basis of potato drought tolerance [22]. It was determined that in tomato, overexpression of the TAS14 gene under drought conditions leads to tolerance to osmotic stress by reducing osmotic potential and accumulation of solute substances such as sugars and K+. Moreover, overexpression of the TAS14 gene increased Na+ accumulation only in adult leaves, whereas in young leaves the accumulated solutes were K+ and sugars, suggesting that plants with overexpression of the TAS14 gene are able to distribute Na+ accumulation between young and adult leaves over a long period under stress conditions [23]. Simultaneously, transcriptome analysis of genes associated with photosynthesis determined that different genotypes can respond differently to drought at the leaf level, but lose tuber productivity to the same extent [24].
It is advisable to conduct the most extensive studies of enzymatic activity in order to have a complete picture of the functional consequences observed in modeling drought in potatoes. The regulation of drought tolerance related genes (DRO, ERECTA, ERF, DREB and StMYB) has been shown to differ between potato varieties with high and low-stress tolerance index [17]. Studies also showed that EF1 α and sec3 were the most stably expressed genes in potatoes under drought and osmotic stress [25]. It was determined that genes related to proline degradation were inhibited by PEG-induced stress in in vitro seedling culture, which leads to proline accumulation. In this case, drought-tolerant potato genotypes accumulate less proline than drought-sensitive genotypes, and proline accumulation is inversely correlated with root length [26]. It was shown that drought stress was reflected not only in drought-responsive genes associated with altered cell structure and components, but also in enhanced gene expression or “memory” of drought-responsive genes [27].
During drought stress plants employ adaptive measures to withstand osmotic and oxidative stress by initiating signals that modify present metabolism to restore homeostasis. Ethylene is a phytohormone known for inducing plant response to various stresses including drought, heat shock, osmotic stress, oxidative stress, and wounding, and activates plant defense [28,29,30,31]. It is a gaseous plant hormone controlled by many factors including ER24 with the molecular function of promoting the growth of root systems and sustaining metabolism, helping plants sustain the drought environment [32]. Another transcription factor expressed during drought is Nuclear factor Y (NF-Y), well regulated in plants during microbe attack, root development and different stresses [33,34,35].
Previously, researchers studied potato transcriptional responses in drought and osmotic stress conditions [36,37,38], some studies focused on comparing different cultivar gene expression profiles between them using aeroponic systems [39] and mannitol [40] to simulate osmotic stress conditions. Studies were conducted to study expression profiles of particular gene families of potatoes during drought and osmotic stress such as CDPK [41], CIPK [42], WD40 [43], DUF668 [44] comparing homologous genes between each other. Similar studies comparing gene expression profiles of contrasting varieties were conducted on other crops like Saccharum officinarum (sugarcane) [45]. The aim of this work was to determine whether the expression patterns of commonly upregulated genes differed between drought-tolerant and drought-sensitive varieties and whether there were common, conserved genes whose overexpression improved drought tolerance across plant species.

2. Materials and Methods

2.1. Plant Materials

The potato variety Gala and local Kazakh varieties Yagodnyi-19, Aksor, Tyanshanskyi and Shagalaly from the in vitro potato collections of the Institute of Plant Biology and Biotechnology were used in the experiment. The origins of these cultivars are provided in Table 1. A brief description of each variety can be found in Supplementary Materials S3.

2.2. Plant Propagation and Osmotic Stress Application

The 30 internodes were planted using Murashige and Skoog (MS) medium (4.43 g/L of an MS nutrient medium with vitamins (PhytoTechnology Lab, Lenexa, KS, USA), containing 0.8% agar (Santa Cruz Biotechnology, Dallas, TX, USA), 3% sucrose (Sigma-Aldrich Chemie GmbH, Steinheim, Germany)) with cefotaxime antibiotics (100 µg/L) in a bioreactor (RITA® systems (Cirad, Vitropic, France)) one week before forming roots with an immersion frequency of 3 h for 10 min, and air pressure of 0.5 bars as optimized by another study [46]. The rooted plantlets were transferred to a liquid MS medium: control medium or medium with 20% PEG-6000. Each variety was grown for 14 days at 24 °C with a 16/8 photoperiod. During 14 days, the morphological traits of plantlets (plant height, number of internodes and number of leaves), roots’ attributes (length and number) and maximum quantum efficiency of PSII were recorded every day. The maximum quantum efficiency of photosystem (PS) II (Fv/Fm) of leaves was determined using a portable chlorophyll fluorometer (Hansatech Instruments Pocket PEA, Serial Number 1947, Norfolk, UK) after 30 min of dark adaptation in the second and third leaves (from the top) of each plantlet [47].

2.3. Selection of Drought-Tolerant and Drought-Sensitive Varieties

The mean values of the potato traits were related using principal component analysis (PCA). Pearson’s correlation coefficients were used to examine the relationships between variables. The values for each trait were compared between groups using analysis of variance (ANOVA). Student’s criterion test was performed to detect significantly different mean values (p ≤ 0.05) to select contrast varieties. Statistical analysis was performed using JMP 17 PRO (JMP Statistical Discovery LLC, Cary, NC, USA) and R (version 4.4.1). Data visualization was performed using ggplot2 (version 3.5.1) [48] from tidyverse (version 2.0.0) [49] package in R. Pairs plots were created using GGally (version 2.2.1) package and Duncan’s multiple range test for expression data was performed using agricolae (version 1.3-7).

2.4. Plant Material for Gene Expression Analysis

The drought-tolerant “Gala” and its contrasting variety “Aksor” were chosen for gene expression analysis using bioreactor (RITA® systems). Two-week-old plantlets were transferred to bioreactors containing MS medium and MS medium with 20% PEG-6000 for osmotic stress. Three biological replicates were used for the control and treatment. The roots and stems were collected at time points of 0, 1, 2, 4, 6, 8, 10 and 15 days, immediately frozen in liquid nitrogen and kept at −70 °C in a freezer for further RNA extraction and analysis [50].

2.5. Gene Expression Analysis Using qRT-PCR

Total RNA was extracted from cells or snap-frozen tissues using TRIzol Reagent, according to the manufacturer’s protocol. Extracted total RNA was treated with DNAse I (Thermo Fischer Scientific, Waltham, MA, USA) to remove any leftover contaminating DNA. PCR amplifications were performed using a RevertAid First Strand cDNA Synthesis Kit (Thermo ScientificTM, K1621). Prior to amplification, the RNA quality and concentration were determined using a NanoDrop 2000 spectrophotometer (Thermo Fischer Scientific, USA) [51]. Gene expression analysis was performed with the PowerTrackTM SYBR Green Master Mix (Thermo Fischer Scientific, Waltham, MA, USA) using oligo(dT)18 in 15 µL reaction tubes. The protocols outlined in the manufacturer’s guidelines were followed. Expression data were normalized to Perox, and Actin mRNA expression and fold changes were calculated using the 2 Δ Δ C t method. The primers are presented in Table 2.
The expression of the ER24, TAS14, DREB147315, PP2C, 102605413 and NF-YC4 genes was measured using quantitative real-time PCR with the CFX ConnectTM Real-Time System (Bio-Rad, Hercules, CA, USA) and normalized to the expression of Actin and Perox. The cycles for PCR amplification were set in accordance with the PowerTrackTM SYBRTM Green Master Mix manufacturer’s protocol. Among genes that were regulated by drought in potatoes, we focused on those belonging to key categories playing an important role in drought response modulation such as hormones, transcription factors and abiotic defense responses. The primers used for the experiment are provided in Table 2.

2.6. Searching Criteria for Comparative Analysis

To compare the obtained results with those from other studies on potato and other crops, namely, tomato and wheat, studies were identified using the PubMed search tool with the following keywords: [Solanum tuberosum L., Solanum lycopersicum L., Triticum aestivum L.], drought, RNA-Seq. Studies providing raw data, consisting of multiple organs and involving drought tolerance, were selected for comparative analysis. Raw data from selected studies were downloaded and extracted using the NCBI-provided SRA toolkit.
Studies researching crops under drought in different phases and in different organs were downloaded using the SRA toolkit. These studies were identified through a PubMed search with the following keywords: Solanum tuberosum L., potato, drought, RNA-Seq. The selection criteria were as follows: RNA-Seq analysis, consisting of multiple organs, provided raw data or gene count, and preferably consisting of drought-tolerant and intolerant varieties for comparison.
Raw data from studies have been downloaded using the NCBI-provided sra-toolkit (version 3.1.0) according to their accession number given by the original authors. Raw data were evaluated for quality using FastQC (version 0.12.1) [55] and adapter sequence and ends were trimmed using fastp (version 0.23.4) [56]. Clean data were obtained by removing low-quality reads (reads containing more than 10% unknown nucleotides (N) and reads containing greater than 50% low-quality (Q-value ≤ 20) bases). Cleaned reads were aligned to the reference genome using Hisat2 [57], transcript compilation was conducted using StringTie [58], differential expression was performed using Ballgown [59] if gene expression matrix was not provided. The following reference genomes were used for read alignment: potato (SolTub_3.0), and tomato (SL3.1) from the NCBI genome database using NCBI’s datasets tool.

3. Results and Discussion

3.1. Morphological Analysis

PEG-induced osmotic stress on plant morphological traits in vitro is often used for the preliminary screening of crop genotypes for identifying drought tolerance. In particular, it is noted that the length of the seedling and root are the most important traits by which one can find out the drought tolerance of the genotypes. Moreover, these traits are important for most crops, such as soybean (Glycine max) [15], sorghum (Sorghum bicolor) [60], rice (Oryza sativa) [61], rapeseed (Brassica napus L.) [62] and tomatoes (Solanum lycopersicum L.) [9]. In all these experiments, stress led to a decrease in the length of seedlings and roots in vitro. At the same time, it is known that early drought reduces the number of tubers per plant in all studied potato genotypes [7]. Therefore, it would be logical to assume that modeling drought at the early stages of plant development will allow us to distinguish drought-tolerant and drought-sensitive potato varieties. In this study, the growth of in vitro plantlets was determined based on plant height, root length, number of internodes, number of leaves, number of roots, and maximum quantum efficiency of PSII of potato varieties Yagodnyi-19, Gala, Aksor, Tyanshanskyi, and Shagalaly during 15 DAT (days after treatment). All parameters were decreased (Supplementary File) in all evaluated varieties, except the maximum quantum efficiency of PSII of Yagodnyi-19 (0.685) and Tyanshanskyi (0.741) compared to the control (0.656 and 0.690, respectively). The plant height (PH) decreased especially in varieties Shagalaly (3.8 cm) and Aksor (6.4 cm) in comparison to its controls. The length of roots (LR) and number of internodes (NI) were reduced in Aksor (2.6 cm and 6.9) and Tyanshanskyi (4.5 cm and 6.6). The number of leaves decreased in Yagodnyi-19 (7.9) and Aksor (8.0), while the number of roots decreased in Yagodnyi-19 (10.2) and Shagalaly (9.8 cm). The most affected genotype was Aksor according to decreased number of morphological characteristics; in contrast, Gala showed tolerance responses to PEG-induced water stress (Supplementary Tables S6 and S7).

3.1.1. Principal Component Analysis

In this study, principal component analysis (PCA) was used as a multivariate statistical analysis for possible differences in potatoes using the dependent variables, morphological parameters, which served as the main component (Supplementary Table S2). As shown in the Figure 1, PCA explained 83.1% of the total variation; PC1 explained 55.6% and PC2 explained 27.5%. This figure also shows that the group with no treatment and PEG-6000-treated groups formed separate clusters, and each group could be distinguished. The results show that all the morphological parameters were formed in the control condition, which shows that stressors influenced PEG-6000 treatment. The range of association between the first component (PC1) and all the evaluated variables ranged from 11% to 52%, and all of them were in the same direction. The second component (PC2) had a positive direction for maximum quantum efficiency of PSII, plant height and root length (68%, 42% and 36%, respectively) and a negative direction for other variables (Supplementary Figures S5 and S6). PCA revealed trait influences between the control and PEG-6000-treated groups, which is in agreement with previous studies [63,64].
Analyzing the reaction of potato plantlets treated with PEG-6000 in the nutrient medium, it was determined that there was variability in the studied traits in plants. At the same time, the number of leaves and root length changed but not dramatically (Figure 2), while the plant height, number of internodes and number of roots changed significantly (Figure 2).
Thus, all the studied traits can be used to characterize potato varieties for resistance to simulated osmotic stress at the level of in vitro plantlets. Many authors have noted the relationship between abiotic stresses and the work of photosynthetic apparatus [65,66,67]. Thus, the species Triticum dicoccum Schuebl. and Triticum aethiopicum Jakubz. are noted as having the most stable indicators of root system development and relatively high photosynthetic activity under conditions of modeling drought [68]. In addition, it was shown that stress-tolerant plants effectively improved growth and physiological and photosynthetic parameters in pistachio seedlings subjected to salinity and drought stress, which is explained by local adaptation under in vitro conditions [69]. Drought in potatoes is believed to significantly reduce the photosynthetic rate, transpiration rate, stomatal conductance and intercellular CO2 concentration in leaves, while stress reduces the net photosynthetic rate of potato plants [70].

3.1.2. Correlation Analysis

Figure 3 shows the Pearson correlation coefficients between the quantitative traits in the control group. These correlations showed significantly positive correlations between plant height and root length (r = 0.7725), and the number of internodes and number of leaves (r = 0.8085). Figure 4 shows a correlation between quantitative variables combing both control and PEG-6000 treated groups with each variety separated (Supplementary Table S1). Although root length is related to drought tolerance, the authors of the study examining the relationship between root length and drought tolerance suggest that correlations are too weak and effect sizes too small to explain the differences in drought tolerance observed among potato varieties [71]. The maximum quantum efficiency of PSII had significant negative correlations with the number of leaves (r = −0.9732). Root length had a significant negative correlation with the number of internodes (r = −0.9013). The maximum quantum efficiency of PSII had no positive correlation between all parameters except root length (r = 0.3516) (Supplementary Figure S1). It is known that various stresses reduce leaf area, stem area and plant height [72], respectively, reducing the photosynthetic activity of potato [73].
Pearson correlation coefficients between the main parameters in PEG-6000-treated plantlets showed significantly positive correlations between the plant height and maximum quantum efficiency of PSII (r = 0.8558) and the number of internodes and number of leaves (r = 0.8054), as shown in the table (Supplementary Material). The maximum quantum efficiency of PSII showed significant negative correlations with the number of leaves (r = −0.9732). In correlation analysis, no significant negative correlations were found for the PEG-6000-treated group, but the number of leaves showed negative correlations with maximum quantum efficiency of PSII (r = −0.7704) and plant height (r = −0.6589). The results also showed that the maximum quantum efficiency of PSII had no positive correlation with any parameter except plant height (r = 0.8558). According to the correlation data obtained, it should be noted that the maximum quantum efficiency of PSII was generally negatively correlated with all the parameters in both the control and PEG-6000-treated groups, except for root length in the control and plant height in the PEG-6000-treated group, while in the control group, the maximum quantum efficiency of PSII, plant height and root length were negatively correlated with the number of internodes, number of leaves and number of roots. In the PEG-6000-treated group, only root length showed a weak positive correlation with the maximum quantum efficiency of PSII.The PEG-6000-treated group clearly showed higher correlations for some varieties (Yagodnyi-19: 0.708), while a very small change in correlation between the control and PEG-6000-treated groups was observed in some cases (Gala: −0.149; Shagalaly: −0.061) (Supplementary Figures S2–S4).
In our case, in vitro drought modeling also affected the maximum quantum efficiency of PSII, and the different genotypes studied responded differently to the stress in this regard. All the varieties showed greater upward and downward variability in the maximum quantum efficiency of PSII under stress than the control variant (Figure 2). Since all the investigated traits showed reliable variability in in vitro drought modeling with PEG-6000, we had an opportunity to distinguish the studied genotypes in terms of their sensitivity to drought conditions at the early development stage of the plantlets.

3.2. Selection of Drought-Tolerant and Drought-Sensitive Varieties

The presence of PEG-6000 in the nutrient medium affected almost all the potato varieties involved in the experiment (Table 3). Firstly, all the measured quantitative traits decreased. At the same time, the maximum quantum efficiency of PSII changed to varying degrees in different varieties. Thus, in the varieties Yagodnyi-19 and Tyanshanskyi, it increased in the presence of PEG-6000, and in other varieties, it decreased. The results showed that in the Gala variety, the number of internodes, length of roots and their quantity did not change significantly in PEG-6000-induced drought, while in the Aksor and Tyanshanskyi varieties, all the measured quantitative parameters decreased significantly in the presence of PEG-6000 compared to the control (Supplementary Table S3).
Plantlets of six potato varieties were grown in bioreactors with 20% PEG-6000 and bioreactors to compare the morphological properties between them and determine which varieties were tolerant and sensitive to osmotic stress. Prior plantlets were grown in vitro until root formation. Based on the results obtained (Table 3), Aksor was selected as the most sensitive to osmotic stress and Gala as a tolerant to osmotic stress variety to further study the expression of major genes expressed during drought to assess the difference in expression of these genes.Aksor showed a highly noticeable difference in the control and PEG-6000-treated groups in morphological properties compared to Gala, for which no significant difference was observed in the length of roots, number of internodes or number of roots. In conclusion, the variety Gala is the most tolerant among the studied varieties in terms of osmotic stress, while the Aksor variety can be considered the least tolerant.

3.3. Gene Expression Analysis

The quantitative real-time PCR was performed under two conditions (control and PEG-6000) to determine the drought tolerance difference at the molecular level between the two varieties (Gala and Aksor). Figure 5 shows the expression of the genes DREB147315, ER24, NF-YC4, 102605413, PP2C and TAS14 in the roots and stems of control and PEG-6000-treated plantlets of the two varieties with contrasting osmotic stress tolerance properties. Results for stem samples at 0, 1, 2, 4, 6, 8 and 15 days after treatment (DAT) and for roots at 2, 4, 6, 8, 10 and 15 DAT are shown. Overall, all six selected upregulated genes showed a clear difference between the two varieties, the osmotic stress-tolerant Gala and osmotic stress-sensitive Aksor, in both the roots and stems across all days.The expression of DREB147315 in stem samples reached a high level in Gala at 1 DAT (with a relative expression level of 3.46) and showed little deviation across all days, while it increased and decreased a small amount in Aksor. This pattern was also noticeable in root samples. The overall expression of DREB147315 was higher in root samples compared to stem samples. ER24 in both stems and roots showed high expression right after PEG-6000 treatment, declined slowly after 2 DAT in Gala (1 DAT: 5.4; 2 DAT: 4.83, 4 DAT: 4.56, 5 DAT: 3.58, etc.) and had low expression across all days in Aksor. NF-YC4 gene expression was high one day after treatment with PEG-6000 in Gala, and started slowly decreasing in stems, slowly increased until 8 DAT and showed a further decline in root samples, while both the roots and stems of Aksor showed inconsistent and low NF expression. The PP2C gene showed overall high and similar expression across all days in Gala, while Aksor showed inconsistent and low expression. The expression of the 102605413 gene in the stems of Gala increased up to day 6 and, on 15 DAT, had a similar expression level to that on the first day. TAS14’s expression in stem samples of the Gala variety was the highest on 1 DAT, started slowly decreasing until 4 DAT and afterward showed consistent expression. No genes exhibited a noticeable change in the Aksor variety control and PEG-6000-treated samples (Supplementary Tables S4 and S5).
The sixth DAT showed the most significant change in the transcription of major genes in Gala. It can be seen from the figure that all the major drought-responsive genes are much more strongly expressed in the more osmotic stress-tolerant variety (Gala) than in the sensitive variety (Aksor). Leaf and root samples were selected and examined using targeted qRT-PCR for major genes regulated during drought (Supplementary Figures S8 and S9).
It was shown before that overexpression of the DREB genes in potato [74], wheat [75] and sugarcane [76] improves the drought tolerance of crops affected with a low amount of water. It is well known that genes involved in ethylene response (ER24) and NF-YC4 genes were also upregulated in drought and other abiotic stresses in potato [77] and tomato [78,79]. The potato study also revealed the upregulation of the TAS14 gene which was also selected for quantitative gene expression analysis. PP2C is a protein phosphatase gene, also commonly known to be upregulated during drought stress in potato [27,80] and other crops such as wheat [81]. Generally, presented genes are up-regulated in both varieties, but clearly in the drought-tolerant variety it had higher expression levels. This validates the claim that overexpression of genes up-regulated during droughts such as DREB, TAS14 and others improves the drought tolerance of crops [23,82,83,84,85].
The major plant signal transduction pathway in abiotic stresses such as drought is the abscisic acid (ABA)-signaling pathway. Among many mechanisms controlling this pathway, phosphorylation and dephosphorylation, regulated by protein phosphatases are considered important signal transduction mechanisms [86,87]. Protein phosphatases type 2C (PP2Cs)counteract SnRK2 kinases by physical interaction, and thereby inhibit the activation of the transcription factors that mediate ABA-responsive gene expression. ABA modulation via negative feedback regulates the maintenance of plant homeostatis [88]. Stomatal closure control and respective reduction in the leaf area where water loss can occur are considered the main physiological mechanisms by which the ABA signaling pathway plays a role in plants during drought [89].

3.4. Comparison with Other Studies

Based on the criteria described in the Materials and Methods, three studies were selected involving potato [22], tomato [90] and wheat [91]. Using the criteria listed in the Materials and Methods, a number of studies were selected for comparative analysis. After aligning and annotating aligned reads, FPKM values were used for comparison. A Solanum tuberosum L. study involved potato seedlings as in the present study and applied drought stress. Using the raw data provided and methods described above, a table of expression counts was created for each gene that was aligned to the reference genome. Genes in the table (Supplementary Material) that were used in the current study were extracted, namely, DREB147315, ER24, NF-YC4, 102605413, PP2C and TAS14. All six genes involved in the study were upregulated in drought studies of potato and tomato, which may be because they belong to the same group. Among the genes that were regulated by drought in seedlings, we focused on those belonging to key categories playing an important role in drought response modulation such as hormones, transcription factors and abiotic defense responses. In addition, the DREB147315 gene, which is known as a major transcription factor involved in the responses to most abiotic stresses, was also found to be upregulated in all of the studies involving seedlings [92,93,94]. The most upregulated genes in this study were those responsible for the modulation of hormones, transcription factors and abiotic defense responses. In particular, DREB147315, ER24, TAS14 and PP2C were upregulated in all the studies, including those on wheat.
Meta-analyses identifying key genes involved in drought across different species [95] and individual species such as wheat [96] and rice [97] have previously been published. This study used a similar approach, expanding the analysis across different crop species. The aim of this study was to investigate molecular markers of drought tolerance in potato varieties while also investigating markers conserved across crops instead of identifying specific responses in each crop. This study demonstrated that the drought tolerance of a variety can be identified using qRT-PCR on a small subset of genes responsible for drought stress. The study also shows that some genes upregulated during drought stress in potatoes can also be used for screening in other crop species such as wheat and tomato.
In the study [95], authors conducted a meta-analysis of transcriptome comparison across seven plant species, in particular, Arabidopsis thaliana, Solanum lycopersicum L., Zea mays, Vitis vinifera, Malus domestica, Solanum tuberosum L., Triticum aestivum L. including multiple organs and multiple developmental stages of plants.As a result of the comparative analysis, authors found 934 genes commonly regulated among three crops (Solanum lycopersicum L., Malus domestica, Arabidopsis thaliana) in seedlings and 132 genes commonly regulated among four species (Solanum tuberosum L., Zea mays, Triticum aestivum L., Vitis vinifera) in mature leaves.Also, the authors found that 22 genes are commonly regulated among both leaves and seedlings. This shows that there are molecular mechanisms against drought conserved across species and there is the possibility of choosing a small subset of genes to improve the time and efficiency of selection using molecular markers. For example, genes involved in ethylene signaling, DREB family genes and other abiotic stress response genes are commonly regulated in multiple crops.
The present work showed that, at least in the Solanaceae family, the chosen genes (DREB147315, ER24, NF-YC4, 102605413, PP2C and TAS14) can be used to determine the drought tolerance of varieties and comparatively analyze multiple varieties using gene expression analysis. Expanding this gene set and selecting other sets with better performance are future possibilities. Also, combining gene expression analysis with other morphological analyses could improve the selection results. A previous study in wheat drought tolerance also showed that a higher stem dry weight and content of water-soluble carbohydrates correlated with drought tolerance [98].

4. Conclusions

The treatment with PEG-6000 in the nutrient medium affected almost all the potato varieties involved in the experiment. Most of the quantitative traits decreased, while the maximum quantum efficiency of PSII changed to varying degrees in different varieties. Gala was selected as the most osmotic stress-tolerant variety, showing no significant change in internode and root length between the control and PEG-6000-treated plantlets, while the Aksor variety showed significant decreases in all the quantitative parameters. Therefore, Aksor was selected as the sensitive variety and Gala as the tolerant to osmotic stress variety. Quantitative real-time PCR showed clear differences in drought tolerance gene expression regulation between Gala and Aksor in both root and stem samples, and almost all the genes upregulated during drought showed higher expression in Gala than in Aksor. The results obtained from the study show that the drought tolerance of potato varieties can be determined in the early stage of plant development through the gene expression analysis of genes commonly upregulated during drought (DREB147315, TAS14, NF-YC4, ER24, PP2C and 102605413) as simulated with PEG-6000 medium. This study demonstrated that varieties with contrasting osmotic stress tolerance properties show vastly different expression patterns of the main drought-tolerance genes, which can be used for variety selection in the early stage. Comparative analyses were also carried out, showing similar levels of expression of drought-tolerance genes in different plant species, thereby validating claims that the overexpression of genes upregulated in drought improves the drought tolerance of varieties. In the future, morpho-physiological and molecular assessments of drought tolerance in contrasting (arid and non-arid) field conditions will be studied to confirm the hypotheses of the present studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology13110857/s1, Table S1: Pearson’s correlation coefficients between morphological characteristics; Table S2: Eigenvector and variance explained associated with principal component; Table S3: Effect of the PEG-6000 in the nutrient medium; Table S4: relative expression level expression data of each sample; Table S5: Duncan’s test for expression data; Tables S6 and S7: Raw data of morphological characteristics of each group; Figure S1: Correlation analysis plot of morphological characteristics between all samples; Figure S2: Pairs plot between all samples; Figure S3: Pairs plot between control samples; Figure S4: Pairs plot between PEG-6000 treated samples; Figure S5: PCA plot containing replicate level averaged samples; Figure S6: PCA components weight distribution; Figure S7: Boxplot of maximum qunatum efficiency of PSII of Aksor and Gala varieties; Figures S8 and S9: Relative expression levels (fold) of drought-response genes faceted by days.

Author Contributions

Conceptualization, K.Z. and M.S.; data curation, R.K., Z.A. and M.T.; formal analysis, M.T.; funding acquisition, K.Z.; investigation, A.D., Z.S. and N.R.; methodology, D.D., N.R. and Z.S.; software, M.T. and R.K.; project administration, K.Z.; resources, M.S.; supervision, K.Z.; visualization, D.D. and M.T.; writing—original draft, K.Z. and R.K.; writing—review and editing, K.Z. and Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research and APC were funded by the Committee of Science of the Ministry of Science and Higher Education of the Republic of Kazakhstan (grant No. AP14871451).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article are provided in Supplementary Materials.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mukherjee, S.; Mishra, A.; Trenberth, K.E. Climate Change and Drought: A Perspective on Drought Indices. Curr. Clim. Chang. Rep. 2018, 4, 145–163. [Google Scholar] [CrossRef]
  2. Dietz, K.J.; Zörb, C.; Geilfus, C.M. Drought and Crop Yield. Plant Biol. 2021, 23, 881–893. [Google Scholar] [CrossRef] [PubMed]
  3. Sapakhova, Z.; Raissova, N.; Daurov, D.; Zhapar, K.; Daurova, A.; Zhigailov, A.; Zhambakin, K.; Shamekova, M. Sweet Potato as a Key Crop for Food Security under the Conditions of Global Climate Change: A Review. Plants 2023, 12, 2516. [Google Scholar] [CrossRef] [PubMed]
  4. Jain, S.K.; von Wettberg, E.J.; Punia, S.S.; Parihar, A.K.; Lamichaney, A.; Kumar, J.; Gupta, D.S.; Ahmad, S.; Pant, N.C.; Dixit, G.P.; et al. Genomic-Mediated Breeding Strategies for Global Warming in Chickpeas (Cicer arietinum L.). Agriculture 2023, 13, 1721. [Google Scholar] [CrossRef]
  5. Boguszewska-Mańkowska, D.; Ruszczak, B.; Zarzyńska, K. Classification of Potato Varieties Drought Stress Tolerance Using Supervised Learning. Appl. Sci. 2022, 12, 1939. [Google Scholar] [CrossRef]
  6. Schafleitner, R.; Gutierrez, R.; Espino, R.; Gaudin, A.; Pérez, J.; Martínez, M.; Domínguez, A.; Tincopa, L.; Alvarado, C.; Numberto, G.; et al. Field Screening for Variation of Drought Tolerance in Solanum tuberosum L. by Agronomical, Physiological and Genetic Analysis. Potato Res. 2007, 50, 71–85. [Google Scholar] [CrossRef]
  7. Lahlou, O.; Ouattar, S.; Ledent, J.F. The Effect of Drought and Cultivar on Growth Parameters, Yield and Yield Components of Potato. Agronomie 2003, 23, 257–268. [Google Scholar] [CrossRef]
  8. Boguszewska-Mańkowska, D.; Pieczyński, M.; Wyrzykowska, A.; Kalaji, H.M.; Sieczko, L.; Szweykowska-Kulińska, Z.; Zagdańska, 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]
  9. Kumar, P.A.; Reddy, N.; Lakshmi, N. PEG Induced Screening for Drought Tolerance in Tomato Genotypes. Int. J. Curr. Microbiol. Appl. Sci. 2017, 6, 168–181. [Google Scholar] [CrossRef]
  10. Reyes, J.A.O.; Casas, D.E.; Gandia, J.L.; Parducho, M.J.L.; Renovalles, E.M.; Quilloy, E.P.; Delfin, E.F. Polyethylene Glycol-Induced Drought Stress Screening of Selected Philippine High-Yielding Sugarcane Varieties. J. Agric. Food Res. 2023, 14, 100676. [Google Scholar] [CrossRef]
  11. Khan, M.; Khattak, A.M.; Ahmad, S.; Khalil, S.A.; Ahmad, N.; Shah, S.A.A.; Ali, S.; Moon, Y.S.; Hamayun, M.; Alrefaei, A.F. Polyethylene Glycol-Stimulated Drought Stress Enhanced the Biosynthesis of Steviol Glycosides in Stevia Rebaudiana. Pak. J. Bot. 2024, 56, 1645–1652. [Google Scholar] [CrossRef] [PubMed]
  12. Wilmer, L.; Tränkner, M.; Pawelzik, E.; Naumann, M. Sufficient Potassium Supply Enhances Tolerance of Potato Plants to PEG-induced Osmotic Stress. Plant Stress 2022, 5, 100102. [Google Scholar] [CrossRef]
  13. Wang, D.; Song, J.; Lin, T.; Yin, Y.; Mu, J.; Liu, S.; Wang, Y.; Kong, D.; Zhang, Z. Identification of Potato Lipid Transfer Protein Gene Family and Expression Verification of Drought Genes StLTP1 and StLTP7. Plant Direct 2023, 7, e491. [Google Scholar] [CrossRef] [PubMed]
  14. Dobránszki, J.; Magyar-Tábori, K.; Takács-Hudák, Á. Growth and Developmental Responses of Potato to Osmotic Stress Under in Vitro Conditions. Acta Biol. Hung. 2003, 54, 365–372. [Google Scholar] [CrossRef] [PubMed]
  15. Yan, C.; Song, S.; Wang, W.; Wang, C.; Li, H.; Wang, F.; Li, S.; Sun, X. Screening Diverse Soybean Genotypes for Drought Tolerance by Membership Function Value Based on Multiple Traits and Drought-Tolerant Coefficient of Yield. BMC Plant Biol. 2020, 20, 321. [Google Scholar] [CrossRef]
  16. Fukai, S.; Pantuwan, G.; Jongdee, B.; Cooper, M. Screening for Drought Resistance in Rainfed Lowland Rice. Field Crops Res. 1999, 64, 61–74. [Google Scholar] [CrossRef]
  17. Zaki, H.E.M.; Radwan, K.S.A. Response of Potato (Solanum tuberosum L.) Cultivars to Drought Stress under in Vitro and Field Conditions. Chem. Biol. Technol. Agric. 2022, 9, 1. [Google Scholar] [CrossRef]
  18. Shiri, M. Identification of Informative Simple Sequence Repeat (SSR) Markers for Drought Tolerance in Maize. Afr. J. Biotechnol. 2011, 10, 16414–16420. [Google Scholar] [CrossRef]
  19. Guo, P.; Baum, M.; Grando, S.; Ceccarelli, S.; Bai, G.; Li, R.; von Korff, M.; Varshney, R.K.; Graner, A.; Valkoun, J. Differentially Expressed Genes between Drought-Tolerant and Drought-Sensitive Barley Genotypes in Response to Drought Stress during the Reproductive Stage. J. Exp. Bot. 2009, 60, 3531–3544. [Google Scholar] [CrossRef]
  20. Aliche, E.B.; Gengler, T.; Hoendervangers, I.; Oortwijn, M.; Bachem, C.W.B.; Borm, T.; Visser, R.G.F.; van der Linden, C.G. Transcriptomic Responses of Potato to Drought Stress. Potato Res. 2022, 65, 289–305. [Google Scholar] [CrossRef]
  21. 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, e0128041. [Google Scholar] [CrossRef] [PubMed]
  22. Jian, H.; Sun, H.; Liu, R.; Zhang, W.; Shang, L.; Wang, J.; Khassanov, V.; Lyu, D. Construction of Drought Stress Regulation Networks in Potato Based on SMRT and RNA Sequencing Data. BMC Plant Biol. 2022, 22, 381. [Google Scholar] [CrossRef] [PubMed]
  23. Muñoz-Mayor, A.; Pineda, B.; Garcia-Abellán, J.O.; Antón, T.; Garcia-Sogo, B.; Sanchez-Bel, P.; Flores, F.B.; Atarés, A.; Angosto, T.; Pintor-Toro, J.A.; et al. Overexpression of Dehydrin Tas14 Gene Improv. Osmotic Stress Impos. Drought Salin. Tomato. J. Plant Physiol. 2012, 169, 459–468. [Google Scholar] [CrossRef] [PubMed]
  24. Evers, D.; Lefèvre, I.; Legay, S.; Lamoureux, D.; Hausman, J.F.; Rosales, R.O.G.; Marca, L.R.T.; Hoffmann, L.; Bonierbale, M.; Schafleitner, R. Identification of Drought-Responsive Compounds in Potato through a Combined Transcriptomic and Targeted Metabolite Approach. J. Exp. Bot. 2010, 61, 2327–2343. [Google Scholar] [CrossRef] [PubMed]
  25. Tang, X.; Zhang, N.; Si, H.; Calderón-Urrea, A. Selection and Validation of Reference Genes for RT-qPCR Analysis in Potato under Abiotic Stress. Plant Methods 2017, 13, 85. [Google Scholar] [CrossRef]
  26. Liu, Y.; Wang, L.; Li, Y.; Li, X.; Zhang, J. Proline Metabolism-Related Gene Expression in Four Potato Genotypes in Response to Drought Stress. Biol. Plant. 2019, 63, 757–764. [Google Scholar] [CrossRef]
  27. 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]
  28. Arraes, F.B.M.; Beneventi, M.A.; Lisei de Sa, M.E.; Paixao, J.F.R.; Albuquerque, E.V.S.; Marin, S.R.R.; Purgatto, E.; Nepomuceno, A.L.; Grossi-de-Sa, M.F. Implications of Ethylene Biosynthesis and Signaling in Soybean Drought Stress Tolerance. BMC Plant Biol. 2015, 15, 213. [Google Scholar] [CrossRef]
  29. Zhang, Z.; Li, F.; Li, D.; Zhang, H.; Huang, R. Expression of Ethylene Response Factor JERF1 in Rice Improves Tolerance to Drought. Planta 2010, 232, 765–774. [Google Scholar] [CrossRef]
  30. Wang, Y.; Yuan, M.; Li, Z.; Niu, Y.; Jin, Q.; Zhu, B.; Xu, Y. Effects of Ethylene Biosynthesis and Signaling on Oxidative Stress and Antioxidant Defense System in Nelumbo Nucifera G. under Cadmium Exposure. Environ. Sci. Pollut. Res. 2020, 27, 40156–40170. [Google Scholar] [CrossRef]
  31. Deikman, J. Molecular Mechanisms of Ethylene Regulation of Gene Transcription. Physiol. Plant. 1997, 100, 561–566. [Google Scholar] [CrossRef]
  32. Nazir, F.; Peter, P.; Gupta, R.; Kumari, S.; Nawaz, K.; Khan, M.I.R. Plant Hormone Ethylene: A Leading Edge in Conferring Drought Stress Tolerance. Physiol. Plant. 2024, 176, e14151. [Google Scholar] [CrossRef]
  33. Chen, M.; Zhao, Y.; Zhuo, C.; Lu, S.; Guo, Z. Overexpression of a—Transcription Factor from Bermudagrass Confers Tolerance to Drought and Salinity in Transgenic Rice. Plant Biotechnol. J. 2015, 13, 482–491. [Google Scholar] [CrossRef] [PubMed]
  34. Zanetti, M.E.; Rípodas, C.; Niebel, A. Plant NF-Y Transcription Factors: Key Players in Plant-Microbe Interactions, Root Development and Adaptation to Stress. Biochim. Biophys. Acta (BBA)-Gene Regul. Mech. 2017, 1860, 645–654. [Google Scholar] [CrossRef]
  35. Swain, S.; Myers, Z.A.; Siriwardana, C.L.; Holt, B.F. The Multifaceted Roles of NUCLEAR FACTOR-Y in Arabidopsis thaliana Development and Stress Responses. Biochim. Biophys. Acta (BBA)-Gene Regul. Mech. 2017, 1860, 636–644. [Google Scholar] [CrossRef]
  36. Li, Q.; Qin, Y.; Hu, X.; Jin, L.; Li, G.; Gong, Z.; Xiong, X.; Wang, W. Physiology and Gene Expression Analysis of Potato (Solanum tuberosum L.) in Salt Stress. Plants 2022, 11, 1565. [Google Scholar] [CrossRef]
  37. Li, Q.; Qin, Y.; Hu, X.; Li, G.; Ding, H.; Xiong, X.; Wang, W. Transcriptome Analysis Uncovers the Gene Expression Profile of Salt-Stressed Potato (Solanum tuberosum L.). Sci. Rep. 2020, 10, 5411. [Google Scholar] [CrossRef]
  38. Dekomah, S.D.; Wang, Y.; Qin, T.; Xu, D.; Sun, C.; Yao, P.; Liu, Y.; Bi, Z.; Bai, J. Identification and Expression Analysis of Calcium-Dependent Protein Kinases Gene Family in Potato Under Drought Stress. Front. Genet. 2022, 13, 874397. [Google Scholar] [CrossRef]
  39. Ponce, O.P.; Torres, Y.; Prashar, A.; Buell, R.; Lozano, R.; Orjeda, G.; Compton, L. Transcriptome Profiling Shows a Rapid Variety-Specific Response in Two Andigenum Potato Varieties under Drought Stress. Front. Plant Sci. 2022, 13, 1003907. [Google Scholar] [CrossRef]
  40. Qin, T.; Sun, C.; Kazim, A.; Cui, S.; Wang, Y.; Richard, D.; Yao, P.; Bi, Z.; Liu, Y.; Bai, J. Comparative Transcriptome Analysis of Deep-Rooting and Shallow-Rooting Potato (Solanum tuberosum L.) Genotypes under Drought Stress. Plants 2022, 11, 2024. [Google Scholar] [CrossRef]
  41. Bi, Z.; Wang, Y.; Li, P.; Sun, C.; Qin, T.; Bai, J. Evolution and Expression Analysis of CDPK Genes under Drought Stress in Two Varieties of Potato. Biotechnol. Lett. 2021, 43, 511–521. [Google Scholar] [CrossRef] [PubMed]
  42. Ma, R.; Liu, W.; Li, S.; Zhu, X.; Yang, J.; Zhang, N.; Si, H. Genome-Wide Identification, Characterization and Expression Analysis of the CIPK Gene Family in Potato (Solanum tuberosum L.) and the Role of StCIPK10 in Response to Drought and Osmotic Stress. Int. J. Mol. Sci. 2021, 22, 13535. [Google Scholar] [CrossRef] [PubMed]
  43. Liu, Z.; Liu, Y.; Coulter, J.A.; Shen, B.; Li, Y.; Li, C.; Cao, Z.; Zhang, J. The WD40 Gene Family in Potato (Solanum tuberosum L.): Genome-Wide Analysis and Identification of Anthocyanin and Drought-Related WD40s. Agronomy 2020, 10, 401. [Google Scholar] [CrossRef]
  44. Liu, E.; Li, Z.; Luo, Z.; Xu, L.; Jin, P.; Ji, S.; Zhou, G.; Wang, Z.; Zhou, Z.; Zhang, H. Genome-Wide Identification of DUF668 Gene Family and Expression Analysis under Drought and Salt Stresses in Sweet Potato [Ipomoea batatas (L.) Lam]. Genes 2023, 14, 217. [Google Scholar] [CrossRef]
  45. Contiliani, D.F.; de Oliveira Nebó, J.F.C.; Ribeiro, R.V.; Andrade, L.M.; Peixoto Júnior, R.F.; Lembke, C.G.; Machado, R.S.; Silva, D.N.; Belloti, M.; de Souza, G.M.; et al. Leaf Transcriptome Profiling of Contrasting Sugarcane Genotypes for Drought Tolerance under Field Conditions. Sci. Rep. 2022, 12, 9153. [Google Scholar] [CrossRef]
  46. Daurov, D.; Daurova, A.; Sapakhova, Z.; Kanat, R.; Akhmetzhanova, D.; Abilda, Z.; Toishimanov, M.; Raissova, N.; Otynshiyev, M.; Zhambakin, K.; et al. The Impact of the Growth Regulators and Cultivation Conditions of Temporary Immersion Systems (TISs) on the Morphological Characteristics of Potato Explants and Microtubers. Agronomy 2024, 14, 1782. [Google Scholar] [CrossRef]
  47. Daurov, D.; Daurova, A.; Karimov, A.; Tolegenova, D.; Volkov, D.; Raimbek, D.; Zhambakin, K.; Shamekova, M. Determining Effective Methods of Obtaining Virus-Free Potato for Cultivation in Kazakhstan. Am. J. Potato Res. 2020, 97, 367–375. [Google Scholar] [CrossRef]
  48. Villanueva, R.A.M.; Chen, Z.J. Ggplot2: Elegant Graphics for Data Analysis. Meas. Interdiscip. Res. Perspect. 2019, 17, 160–167. [Google Scholar] [CrossRef]
  49. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.D.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  50. Daurov, D.; Argynbayeva, A.; Daurova, A.; Zhapar, K.; Sapakhova, Z.; Zhambakin, K.; Shamekova, M. Monitoring the Spread of Potato Virus Diseases in Kazakhstan. Am. J. Potato Res. 2023, 100, 63–70. [Google Scholar] [CrossRef]
  51. Daurov, D.; Lim, Y.H.; Park, S.U.; Kim, Y.H.; Daurova, A.; Sapakhova, Z.; Zhapar, K.; Abilda, Z.; Toishimanov, M.; Shamekova, M.; et al. Selection and Characterization of Lead-Tolerant Sweetpotato Cultivars for Phytoremediation. Plant Biotechnol. Rep. 2024, 18, 327–339. [Google Scholar] [CrossRef]
  52. Sultana, N. Molecular and Genetic Analysis of Drought Tolerance in Potato; Wageningen University: Wageningen, The Netherlands, 2014. [Google Scholar]
  53. 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. [Google Scholar] [CrossRef]
  54. Zhang, N.; Liu, B.; Ma, C.; Zhang, G.; Chang, J.; Si, H.; Wang, D. Transcriptome Characterization and Sequencing-Based Identification of Drought-Responsive Genes in Potato. Mol. Biol. Rep. 2014, 41, 505–517. [Google Scholar] [CrossRef]
  55. Brown, J.; Pirrung, M.; McCue, L.A. FQC Dashboard: Integrates FastQC Results into a Web-Based, Interactive, and Extensible FASTQ Quality Control Tool. Bioinformatics 2017, 33, 3137–3139. [Google Scholar] [CrossRef]
  56. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. Fastp: An Ultra-Fast All-in-One FASTQ Preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  57. Zhang, Y.; Park, C.; Bennett, C.; Thornton, M.; Kim, D. Rapid and Accurate Alignment of Nucleotide Conversion Sequencing Reads with HISAT-3N. Genome Res. 2021, 31, 1290–1295. [Google Scholar] [CrossRef]
  58. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. StringTie Enables Improved Reconstruction of a Transcriptome from RNA-seq Reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef]
  59. Frazee, A.C.; Pertea, G.; Jaffe, A.E.; Langmead, B.; Salzberg, S.L.; Leek, J.T. Ballgown Bridges the Gap between Transcriptome Assembly and Expression Analysis. Nat. Biotechnol. 2015, 33, 243–246. [Google Scholar] [CrossRef]
  60. Smith, R.H.; Bhaskaran, S.; Miller, F.R. Screening for Drought Tolerance in Sorghum Using Cell Culture. In Vitro Cell. Dev. Biol. 1985, 21, 541–545. [Google Scholar] [CrossRef]
  61. de Datta, S.K.; Malabuyoc, J.A.; Aragon, E.L. A Field Screening Technique for Evaluating Rice Germplasm for Drought Tolerance during the Vegetative Stage. Field Crops Res. 1988, 19, 123–134. [Google Scholar] [CrossRef]
  62. Daurova, A.; Daurov, D.; Volkov, D.; Zhapar, K.; Raimbek, D.; Shamekova, M.; Zhambakin, K. Doubled Haploids of Interspecific Hybrids between Brassica Napus and Brassica Rapa for Canola Production with Valuable Breeding Traits. OCL 2020, 27, 45. [Google Scholar] [CrossRef]
  63. Bakhshandeh, E.; Gholamhossieni, M. Quantification of Soybean Seed Germination Response to Seed Deterioration under PEG-induced Water Stress Using Hydrotime Concept. Acta Physiol. Plant. 2018, 40, 126. [Google Scholar] [CrossRef]
  64. Tajaragh, R.P.; Rasouli, F.; Giglou, M.T.; Zahedi, S.M.; Hassanpouraghdam, M.B.; Aazami, M.A.; Adámková, A.; Mlček, J. Morphological and Physiological Responses of In Vitro-Grown Cucurbita sp. Landraces Seedlings under Osmotic Stress by Mannitol and PEG. Horticulturae 2022, 8, 1117. [Google Scholar] [CrossRef]
  65. Balla, K.; Bencze, S.; Bónis, P.; Árendás, T.; Veisz, O. Changes in the Photosynthetic Efficiency of Winter Wheat in Response to Abiotic Stress. Cent. Eur. J. Biol. 2014, 9, 519–530. [Google Scholar] [CrossRef]
  66. Sharma, A.; Kumar, V.; Shahzad, B.; Ramakrishnan, M.; Singh Sidhu, G.P.; Bali, A.S.; Handa, N.; Kapoor, D.; Yadav, P.; Khanna, K.; et al. Photosynthetic Response of Plants Under Different Abiotic Stresses: A Review. J. Plant Growth Regul. 2020, 39, 509–531. [Google Scholar] [CrossRef]
  67. Muhammad, I.; Shalmani, A.; Ali, M.; Yang, Q.H.; Ahmad, H.; Li, F.B. Mechanisms Regulating the Dynamics of Photosynthesis Under Abiotic Stresses. Front. Plant Sci. 2021, 11, 615942. [Google Scholar] [CrossRef]
  68. Terletskaya, N.; Zobova, N.; Stupko, V.; Shuyskaya, E. Growth and Photosynthetic Reactions of Different Species of Wheat Seedlings under Drought and Salt Stress. Period. Biol. 2017, 119, 37–45. [Google Scholar] [CrossRef]
  69. Khalilpour, M.; Mozafari, V.; Abbaszadeh-Dahaji, P. Tolerance to Salinity and Drought Stresses in Pistachio (Pistacia vera L.) Seedlings Inoculated with Indigenous Stress-Tolerant PGPR Isolates. Sci. Hortic. 2021, 289, 110440. [Google Scholar] [CrossRef]
  70. Wang, L.; Liu, Y.; Feng, S.; Yang, J.; Li, D.; Zhang, J. Roles of Plasmalemma Aquaporin Gene StPIP1 in Enhancing Drought Tolerance in Potato. Front. Plant Sci. 2017, 8, 616. [Google Scholar] [CrossRef]
  71. Hill, D.; Nelson, D.; Hammond, J.; Bell, L. Morphophysiology of Potato (Solanum Tuberosum) in Response to Drought Stress: Paving the Way Forward. Front. Plant Sci. 2021, 11, 597554. [Google Scholar] [CrossRef]
  72. Deblonde, P.M.K.; Ledent, J.F. Effects of Moderate Drought Conditions on Green Leaf Number, Stem Height, Leaf Length and Tuber Yield of Potato Cultivars. Eur. J. Agron. 2001, 14, 31–41. [Google Scholar] [CrossRef]
  73. Allen, E.J.; Scott, R.K. An Analysis of Growth of the Potato Crop. J. Agric. Sci. 1980, 94, 583–606. [Google Scholar] [CrossRef]
  74. Jia, X.X.; Li, Y.T.; Qi, E.F.; Ma, S.; Hu, X.Y.; Wen, G.H.; Wang, Y.H.; Li, J.W.; Zhang, X.H.; Wang, H.M.; et al. Overexpression of the Arabidopsis DREB1A Gene Enhances Potato Drought-Resistance. Russ. J. Plant Physiol. 2016, 63, 523–531. [Google Scholar] [CrossRef]
  75. Wang, J.W.; Yang, F.P.; Chen, X.Q.; Liang, R.Q.; Zhang, L.Q.; Geng, D.M.; Zhang, X.D.; Song, Y.Z.; Zhang, G.S. Induced Expression of DREB Transcriptional Factor and Study on Its Physiological Effects of Drought Tolerance in Transgenic Wheat. Acta Genet. Sin. 2006, 33, 468–476. [Google Scholar] [CrossRef]
  76. Reis, R.R.; Andrade Dias Brito da Cunha, B.; Martins, P.K.; Martins, M.T.B.; Alekcevetch, J.C.; Chalfun-Júnior, A.; Andrade, A.C.; Ribeiro, A.P.; Qin, F.; Mizoi, J.; et al. Induced Over-Expression of AtDREB2A CA Improves Drought Tolerance in Sugarcane. Plant Sci. 2014, 221–222, 59–68. [Google Scholar] [CrossRef]
  77. van Muijen, D.; Anithakumari, A.; Maliepaard, C.; Visser, R.G.F.; 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]
  78. Zegzouti, H.; Jones, B.; Frasse, P.; Marty, C.; Maitre, B.; Latché, A.; Pech, J.C.; Bouzayen, M. Ethylene-Regulated Gene Expression in Tomato Fruit: Characterization of Novel Ethylene-Responsive and Ripening-Related Genes Isolated by Differential Display. Plant J. 1999, 18, 589–600. [Google Scholar] [CrossRef]
  79. Hommel, M.; Khalil-Ahmad, Q.; Jaimes-Miranda, F.; Mila, I.; Pouzet, C.; Latché, A.; Pech, J.C.; Bouzayen, M.; Regad, F. Over-Expression of a Chimeric Gene of the Transcriptional Co-Activator MBF1 Fused to the EAR Repressor Motif Causes Developmental Alteration in Arabidopsis and Tomato. Plant Sci. 2008, 175, 168–177. [Google Scholar] [CrossRef]
  80. Wang, Y.f.; Liao, Y.q.; Wang, Y.p.; Yang, J.w.; Zhang, N.; Si, H.j. Genome-Wide Identification and Expression Analysis of StPP2C Gene Family in Response to Multiple Stresses in Potato (Solanum tuberosum L.). J. Integr. Agric. 2020, 19, 1609–1624. [Google Scholar] [CrossRef]
  81. Yu, X.; Han, J.; Li, L.; Zhang, Q.; Yang, G.; He, G. Wheat PP2C-a10 Regulates Seed Germination and Drought Tolerance in Transgenic Arabidopsis. Plant Cell Rep. 2020, 39, 635–651. [Google Scholar] [CrossRef]
  82. Chen, J.Q.; Meng, X.P.; Zhang, Y.; Xia, M.; Wang, X.P. Over-Expression of OsDREB Genes Lead to Enhanced Drought Tolerance in Rice. Biotechnol. Lett. 2008, 30, 2191–2198. [Google Scholar] [CrossRef] [PubMed]
  83. Cui, M.; Zhang, W.; Zhang, Q.; Xu, Z.; Zhu, Z.; Duan, F.; Wu, R. Induced Over-Expression of the Transcription Factor OsDREB2A Improves Drought Tolerance in Rice. Plant Physiol. Biochem. 2011, 49, 1384–1391. [Google Scholar] [CrossRef] [PubMed]
  84. Hu, Z.; Li, Y.; Yang, J.; Song, S.; Li, X.; Xiong, C.; Yi, P.; Liu, C.; Hu, R.; Huang, X. The Positive Impact of the NtTAS14-like1 Gene on Osmotic Stress Response in Nicotiana Tabacum. Plant Cell Rep. 2023, 43, 25. [Google Scholar] [CrossRef] [PubMed]
  85. Han, X.; Tang, S.; An, Y.; Zheng, D.C.; Xia, X.L.; Yin, W.L. Overexpression of the Poplar NF-YB7 Transcription Factor Confers Drought Tolerance and Improves Water-Use Efficiency in Arabidopsis. J. Exp. Bot. 2013, 64, 4589–4601. [Google Scholar] [CrossRef]
  86. Chuong, N.N.; Hoang, X.L.T.; Nghia, D.H.T.; Dai, T.N.T.; Thi, V.A.L.; Thao, N.P. Protein Phosphatase Type 2C Functions in Phytohormone-Dependent Pathways and in Plant Responses to Abiotic Stresses. Curr. Protein Pept. Sci. 2021, 22, 430–440. [Google Scholar] [CrossRef]
  87. Kim, T.H. Mechanism of ABA Signal Transduction: Agricultural Highlights for Improving Drought Tolerance. J. Plant Biol. 2014, 57, 1–8. [Google Scholar] [CrossRef]
  88. Jung, C.; Nguyen, N.H.; Cheong, J.J. Transcriptional Regulation of Protein Phosphatase 2C Genes to Modulate Abscisic Acid Signaling. Int. J. Mol. Sci. 2020, 21, 9517. [Google Scholar] [CrossRef]
  89. Davies, W.J.; Wilkinson, S.; Loveys, B. Stomatal Control by Chemical Signalling and the Exploitation of This Mechanism to Increase Water Use Efficiency in Agriculture. New Phytol. 2002, 153, 449–460. [Google Scholar] [CrossRef]
  90. Liu, M.; Yu, H.; Zhao, G.; Huang, Q.; Lu, Y.; Ouyang, B. Profiling of Drought-Responsive microRNA and mRNA in Tomato Using High-Throughput Sequencing. BMC Genom. 2017, 18, 481. [Google Scholar] [CrossRef]
  91. Chu, C.; Wang, S.; Paetzold, L.; Wang, Z.; Hui, K.; Rudd, J.C.; Xue, Q.; Ibrahim, A.M.H.; Metz, R.; Johnson, C.D.; et al. RNA-seq Analysis Reveals Different Drought Tolerance Mechanisms in Two Broadly Adapted Wheat Cultivars ‘TAM 111’ and ‘TAM 112’. Sci. Rep. 2021, 11, 4301. [Google Scholar] [CrossRef]
  92. Khan, M.S. The Role of Dreb Transcription Factors in Abiotic Stress Tolerance of Plants. Biotechnol. Biotechnol. Equip. 2011, 25, 2433–2442. [Google Scholar] [CrossRef]
  93. Singh, K.; Chandra, A. DREBs-potential Transcription Factors Involve in Combating Abiotic Stress Tolerance in Plants. Biologia 2021, 76, 3043–3055. [Google Scholar] [CrossRef]
  94. Lata, C.; Prasad, M. Role of DREBs in Regulation of Abiotic Stress Responses in Plants. J. Exp. Bot. 2011, 62, 4731–4748. [Google Scholar] [CrossRef]
  95. Benny, J.; Pisciotta, A.; Caruso, T.; Martinelli, F. Identification of Key Genes and Its Chromosome Regions Linked to Drought Responses in Leaves across Different Crops through Meta-Analysis of RNA-Seq Data. BMC Plant Biol. 2019, 19, 194. [Google Scholar] [CrossRef] [PubMed]
  96. Saidi, M.N.; Mahjoubi, H.; Yacoubi, I. Transcriptome Meta-Analysis of Abiotic Stresses-Responsive Genes and Identification of Candidate Transcription Factors for Broad Stress Tolerance in Wheat. Protoplasma 2023, 260, 707–721. [Google Scholar] [CrossRef]
  97. Buti, M.; Baldoni, E.; Formentin, E.; Milc, J.; Frugis, G.; Lo Schiavo, F.; Genga, A.; Francia, E. A Meta-Analysis of Comparative Transcriptomic Data Reveals a Set of Key Genes Involved in the Tolerance to Abiotic Stresses in Rice. Int. J. Mol. Sci. 2019, 20, 5662. [Google Scholar] [CrossRef]
  98. Hou, J.; Huang, X.; Sun, W.; Du, C.; Wang, C.; Xie, Y.; Ma, Y.; Ma, D. Accumulation of Water-Soluble Carbohydrates and Gene Expression in Wheat Stems Correlates with Drought Resistance. J. Plant Physiol. 2018, 231, 182–191. [Google Scholar] [CrossRef]
Figure 1. Principal component analysis (PCA) between main indicators of control and PEG-6000-treated plantlets. PEG-6000-treated and control plantlets form distinct clusters in two components. (A) PCA plot with two components (PC1: 52.8%, PC2: 25.9%), (B) Biplot maximum quantum efficiency of photosystem II (Fv/Fm), plant height in cm (PH), length of roots in cm (LR), number of internodes (NI), number of leaves (NL), number of roots (NR).
Figure 1. Principal component analysis (PCA) between main indicators of control and PEG-6000-treated plantlets. PEG-6000-treated and control plantlets form distinct clusters in two components. (A) PCA plot with two components (PC1: 52.8%, PC2: 25.9%), (B) Biplot maximum quantum efficiency of photosystem II (Fv/Fm), plant height in cm (PH), length of roots in cm (LR), number of internodes (NI), number of leaves (NL), number of roots (NR).
Biology 13 00857 g001
Figure 2. Boxplot showing difference between morphological characteristics of plantlets treated with polyethylene glycol (PEG)-6000 and control. Outliers shown as red dots.
Figure 2. Boxplot showing difference between morphological characteristics of plantlets treated with polyethylene glycol (PEG)-6000 and control. Outliers shown as red dots.
Biology 13 00857 g002
Figure 3. Pearson’s correlation analysis between morphological characteristics of control and PEG-6000-treated potato plantlets. Maximum quantum efficiency of photosystem II (Fv/Fm), plant height in cm (PH), length of roots in cm (LR), number of internodes (NI), number of leaves (NL), and number of roots (NR). (A) shows tcorrelation plot for control group, (B) shows correlation plot for PEG-6000 treated group.
Figure 3. Pearson’s correlation analysis between morphological characteristics of control and PEG-6000-treated potato plantlets. Maximum quantum efficiency of photosystem II (Fv/Fm), plant height in cm (PH), length of roots in cm (LR), number of internodes (NI), number of leaves (NL), and number of roots (NR). (A) shows tcorrelation plot for control group, (B) shows correlation plot for PEG-6000 treated group.
Biology 13 00857 g003
Figure 4. Pairs plot showing Pearson’s correlation coefficients between the morphological characteristics of the different potato varieties used in this study. A pairs plot showing separate correlations for the control and PEG-6000-treated groups can be found in the Supplementary Material. Maximum quantum efficiency of photosystem II (Fv/Fm), plant height in cm (PH), length of roots in cm (LR), number of internodes (NI), number of leaves (NL), and number of roots (NR). p-values of correlation tests are represented as follows: “***”—p-value is < 0.001; “**”—p-value is < 0.01; “*”—p-value is < 0.05; “ ”—otherwise.
Figure 4. Pairs plot showing Pearson’s correlation coefficients between the morphological characteristics of the different potato varieties used in this study. A pairs plot showing separate correlations for the control and PEG-6000-treated groups can be found in the Supplementary Material. Maximum quantum efficiency of photosystem II (Fv/Fm), plant height in cm (PH), length of roots in cm (LR), number of internodes (NI), number of leaves (NL), and number of roots (NR). p-values of correlation tests are represented as follows: “***”—p-value is < 0.001; “**”—p-value is < 0.01; “*”—p-value is < 0.05; “ ”—otherwise.
Biology 13 00857 g004
Figure 5. Relative expression levels of DREB147315, ER24, NF-YC4, 102605413, PP2C and TAS14 genes in Aksor and Gala varieties in two different MS media, with the other medium containing PEG-6000 to simulate osmotic stress. (A) Stem samples (B) Root samples.
Figure 5. Relative expression levels of DREB147315, ER24, NF-YC4, 102605413, PP2C and TAS14 genes in Aksor and Gala varieties in two different MS media, with the other medium containing PEG-6000 to simulate osmotic stress. (A) Stem samples (B) Root samples.
Biology 13 00857 g005
Table 1. Solanum tuberosum L. varieties used in this study.
Table 1. Solanum tuberosum L. varieties used in this study.
VarietiesOrigin
GalaNORIKA GMBH, Germany
Yagodnyi-19Northwest Agricultural Research Center
Aksor“Kazakh Research Institute of Potato and Vegetables Growing”, LLP
Tyanshanskyi“Kazakh Research Institute of Potato and Vegetables Growing”, LLP
Shagalaly“Kazakh Research Institute of Potato and Vegetables Growing”, LLP × A.I. Barayev Research and Production Center for Grain Farming
Table 2. Primers used for the expression analysis of the six chosen genes.
Table 2. Primers used for the expression analysis of the six chosen genes.
GeneSequenceReference
ActinAGGAGCATCCTGTCCTCCTAA[25]
CACCATCACCAGAGTCCAACA
PeroxAGCACTGATCCATCACATCCC[52]
TGGTAGTTGGAAAGATTGAGAAGC
ER24GAATCAGGCATTGCGAGCTGAATCAGGCATT[52]
GCCGCCTTCTTGTTCAATCGCCGCCTTCTTG
TAS14CAACAGCAGCTTCGTCGATCAACAGCAGCTT[52]
CATGTCCTCCTCCTGGCATCATGTCCTCCTC
DREB147315TGTTCATGGATGAGGAAGCG[6]
AACATTGGGGAGGAGGTAGCAT
PP2CTCACCGATTGCTCGAGACA[53]
GTCCCAATTCCTCTGTCCA
102605413AGTGAAGCATTTCGTAGAGCCA[54]
ACGATGAGTCATGGTTCTGCTT
NF-YC4CGGAGATACCCACCAACTCCGGAGATACCCA[52]
AAAGCTCGGTGGAACTAGCAAAGCTCGGTGG
Table 3. Effect of the PEG-6000 in the nutrient medium during potato seedling growth. NS: not significant; the p-values are for comparisons between the control and PEG-6000-treated groups of each genotype. Maximum quantum efficiency of photosystem II (Fv/Fm), plant height in cm (PH), length of roots in cm (LR), number of internodes (NI), number of leaves (NL), and number of roots (NR).
Table 3. Effect of the PEG-6000 in the nutrient medium during potato seedling growth. NS: not significant; the p-values are for comparisons between the control and PEG-6000-treated groups of each genotype. Maximum quantum efficiency of photosystem II (Fv/Fm), plant height in cm (PH), length of roots in cm (LR), number of internodes (NI), number of leaves (NL), and number of roots (NR).
Varieties(Fv/Fm)PHLRNINLNR
concp-Valueconcp-Valueconcp-Valueconcp-Valueconcp-Valueconcp-Value
Yagodnyi-19-C0.6560.015611.30.00013.80.00158.6NS13.40.017614.40.0005
Yagodnyi-19-PEG0.685 8.2 2.9 7.5 7.9 10.2
1-13 Gala-C0.6850.00059.70.00024.5NS9.5NS11.20.04612.9NS
Gala-PEG0.635 7.0 4.3 8.5 10.1 12.3
Aksor-C0.6860.00029.10.00014.10.00019.60.000111.10.000114.60.0001
Aksor-PEG0.635 6.4 2.6 6.9 8.0 10.6
Tyanshanskyi-C0.6900.000112.30.00017.80.00018.30.00039.40.000112.60.0275
Tyanshanskyi-PEG0.741 9.1 4.5 6.6 7.3 10.3
Shagalaly-C0.6900.000112.00.00015.2NS8.90.000110.6NS15.70.0001
Shagalaly-PEG0.619 3.8 4.6 7.3 9.6 9.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kanat, R.; Shamekova, M.; Sapakhova, Z.; Toishimanov, M.; Daurov, D.; Raissova, N.; Abilda, Z.; Daurova, A.; Zhambakin, K. Gene Expression Analysis for Drought Tolerance in Early Stage of Potato Plant Development. Biology 2024, 13, 857. https://doi.org/10.3390/biology13110857

AMA Style

Kanat R, Shamekova M, Sapakhova Z, Toishimanov M, Daurov D, Raissova N, Abilda Z, Daurova A, Zhambakin K. Gene Expression Analysis for Drought Tolerance in Early Stage of Potato Plant Development. Biology. 2024; 13(11):857. https://doi.org/10.3390/biology13110857

Chicago/Turabian Style

Kanat, Rakhim, Malika Shamekova, Zagipa Sapakhova, Maxat Toishimanov, Dias Daurov, Nurgul Raissova, Zhanar Abilda, Ainash Daurova, and Kabyl Zhambakin. 2024. "Gene Expression Analysis for Drought Tolerance in Early Stage of Potato Plant Development" Biology 13, no. 11: 857. https://doi.org/10.3390/biology13110857

APA Style

Kanat, R., Shamekova, M., Sapakhova, Z., Toishimanov, M., Daurov, D., Raissova, N., Abilda, Z., Daurova, A., & Zhambakin, K. (2024). Gene Expression Analysis for Drought Tolerance in Early Stage of Potato Plant Development. Biology, 13(11), 857. https://doi.org/10.3390/biology13110857

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