Next Article in Journal
Neuronal Ablation of CoA Synthase Causes Motor Deficits, Iron Dyshomeostasis, and Mitochondrial Dysfunctions in a CoPAN Mouse Model
Next Article in Special Issue
Salinity Effects on Guard Cell Proteome in Chenopodium quinoa
Previous Article in Journal
Brain Site-Specific Inhibitory Effects of the GLP-1 Analogue Exendin-4 on Alcohol Intake and Operant Responding for Palatable Food
Previous Article in Special Issue
Genome-Wide Gene Expression Profiles Analysis Reveal Novel Insights into Drought Stress in Foxtail Millet (Setaria italica L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distinct Preflowering Drought Tolerance Strategies of Sorghum bicolor Genotype RTx430 Revealed by Subcellular Protein Profiling

1
Earth and Biological Sciences Directorate, Pacific Northwest National Laboratories, Richland, WA 99354, USA
2
Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, WA 99354, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(24), 9706; https://doi.org/10.3390/ijms21249706
Submission received: 26 October 2020 / Revised: 12 December 2020 / Accepted: 13 December 2020 / Published: 19 December 2020
(This article belongs to the Special Issue Mechanisms of Drought, Temperature and Salinity Tolerance in Plants)

Abstract

:
Drought is the largest stress affecting agricultural crops, resulting in substantial reductions in yield. Plant adaptation to water stress is a complex trait involving changes in hormone signaling, physiology, and morphology. Sorghum (Sorghum bicolor (L.) Moench) is a C4 cereal grass; it is an agricultural staple, and it is particularly drought-tolerant. To better understand drought adaptation strategies, we compared the cytosolic- and organelle-enriched protein profiles of leaves from two Sorghum bicolor genotypes, RTx430 and BTx642, with differing preflowering drought tolerances after 8 weeks of growth under water limitation in the field. In agreement with previous findings, we observed significant drought-induced changes in the abundance of multiple heat shock proteins and dehydrins in both genotypes. Interestingly, our data suggest a larger genotype-specific drought response in protein profiles of organelles, while cytosolic responses are largely similar between genotypes. Organelle-enriched proteins whose abundance significantly changed exclusively in the preflowering drought-tolerant genotype RTx430 upon drought stress suggest multiple mechanisms of drought tolerance. These include an RTx430-specific change in proteins associated with ABA metabolism and signal transduction, Rubisco activation, reactive oxygen species scavenging, flowering time regulation, and epicuticular wax production. We discuss the current understanding of these processes in relation to drought tolerance and their potential implications.

1. Introduction

Drought is the primary abiotic factors negatively affecting agriculture [1], causing tremendous losses in yield worldwide [2,3]. Current climate models predict rising seasonal temperatures and incidence of drought [2], and human population growth is expected to reach 9.7 billion by 2050 [4]. To prevent global food shortages, we must develop a better understanding of drought adaptation strategies in plants. In response to drought, general adaptation strategies include accelerating the transition from vegetative to reproductive (escape), increasing internal water content (avoidance), and to continue growth with decreased internal water content (tolerance) [5]. Each of these adaptation strategies require complex molecular responses to drought, including hormone signaling [6] and production of osmoprotectants [7], as well as developmental responses including stomatal closure, changes to leaf morphology, and altered root system architecture [8]. A common strategy to better understand these adaptation mechanisms is to compare plants with differing drought tolerances; this has been successful in multiple crop plants including wheat, rice, maize, and sorghum [9].
Sorghum (Sorghum bicolor (L.) Moench) (referred to as sorghum hereafter) is among the top five major cereal crops; it is used for food, biofuel production, and livestock feed [10]. Sorghum is a C4 grass, has a small and sequenced diploid genome, and is particularly drought-tolerant, making it an ideal candidate to identify novel drought adaptation strategies. Recent studies have identified multiple quantitative trait loci (QTLs) that control drought tolerance in sorghum at multiple life stages [11,12]. Multiple reports have identified molecular differences between sorghum genotypes differing in water use efficiency and drought tolerance. For example, increased expression of dehydrins and heat shock proteins as well as accumulation of specific sugars and sugar alcohols were identified as unique drought tolerance mechanisms in drought-tolerant sorghum Samsorg 17 compared to drought-susceptible sorghum Samsorg 40 [13]. Similarly, a comparison of transcript-level drought responses of sorghum genotypes IS22330 and IS20351 (drought-tolerant and -susceptible, respectively) suggests proteins involved in ROS scavenging, biosynthesis of secondary metabolites (e.g., phenylpropanoids, flavonoids), and a constitutively higher expression of cuticular wax biosynthesis genes are partly responsible for IS22330′s drought tolerance [9].
Comparisons of protein-level drought response using sorghum genotypes with differing drought tolerance identified the putative Rubisco chaperone HSP60 as changing uniquely to drought-tolerant cultivars [14]. Similarly, protein-level drought responses specific to the drought-tolerant sorghum line EL9 and not the drought-sensitive line Tabat included Rubisco large chain, a heat-shock protein, a glycine-rich RNA binding protein, and chloroplastic malate dehydrogenase [15]. Two sorghum genotypes with different drought tolerance, RTx430 and BTx642, are characterized as preflowering drought-tolerant, and postflowering drought-tolerant, respectively [16,17]. A large-scale RNA-Seq study found that preflowering drought stress caused similar transcriptional responses, but that the magnitude of changes in gene expression were different between genotypes [18]. Varoquaux et al. also found RTx430 genotype-specific transcriptional changes during preflowering drought stress, suggesting specific adaptation strategies, including upregulation of genes involved in shikimate pathway (relative to BTx642), increased glutathione-S transferase (GST) activity, and increased proline content [18]. Interestingly, preflowering drought induced inhibitors of ABA signaling in both sorghum genotypes, suggesting sorghum adapts to drought in part by preventing excessive ABA responses [18,19]. Recently, a comparison of histone modifications from drought-stressed RTx430 and BTx642 also revealed distinct histone posttranslational modification profiles between the two genotypes [20]. While most of the literature uses transcriptomics, experiments at both the protein and transcript levels suggest that photosynthesis and other organelle-localized processes are important for drought adaptation. Despite this trend, little is known about the protein-level response to drought in organelles.
The abundance of cytosolic proteins in whole-cell data-dependent proteomics experiments likely precludes identification of lower-abundance organelle proteins and their changes during water deficit. This study aimed to better understand RTx430′s unique drought tolerance by performing protein profiling of cytosolic- and organelle-enriched cellular compartments (CCs) in response to drought and comparing them with the preflowering drought-sensitive genotype BTx642. We extracted proteins from sorghum leaf tissues collected from a field study with well-controlled drought stress [18]. Our successful enrichment of the different CCs was confirmed by the far higher abundances of chloroplastic and histone proteins in organelle-enriched samples, and higher levels of ribosomal proteins in cytosolic-enriched samples. In response to drought, we observed many protein-level responses in both cellular compartments and genotypes that are consistent with previous literature findings. We show that the drought responses of organelle-enriched samples differs between RTx430 and BTx642, while the response of cytosolic protein profiles was similar between genotypes. Organelle-enriched proteins whose abundance significantly changed exclusively in RTx430 upon drought stress suggest multiple drought adaptation strategies. These include an RTx430-specific change in proteins associated with ABA metabolism and signal transduction, Rubisco activation, and flowering time regulation. We discuss the current understanding of these processes in relation to drought tolerance and their potential implications.

2. Results

2.1. Protein Profiles Are Enriched for Organelle and Cytoplasmic Compartments

The S. bicolor genotype RTx430 exhibits a unique preflowering drought tolerance compared to the preflowering drought sensitive genotype BTx642. To identify the potential mechanisms by which RTx430 uniquely adapts to preflowering water deficit, we compared the RTx430 and BTx642 organelle and cytoplasmic protein profiles after 8 weeks of drought stress preanthesis. To evaluate whether our procedure effectively enriched for organelle and cytoplasmic protein profiles, we first compared the abundance of histone, chloroplastic, and cytosolic ribosomal proteins independent of drought stress (Figure 1). In total, our analysis resulted in confident identification of 3683 and 4130 proteins in cytoplasmic- and organelle-enriched samples, respectively. A total of 3103 proteins were identified in both sample types (Supplemental Table S1). We observed cytosolic ribosomal proteins, including ribosomal L13, S16, and L20 were more abundant in cytosolic protein samples (Figure 1A). Conversely, chloroplastic and histone proteins were predominately more abundant in organelle-enriched samples (Figure 1B,C). These data suggest a successful enrichment for protein profiles representative of the different cellular compartments.

2.2. Organelle-Enriched Protein Profiles Comprise Genotype-Specific Drought Signatures

To evaluate the extent to which drought impacted the proteome of organelle- and cytoplasm- enriched samples, we first performed principal component analysis (PCA). For organelle-enriched samples, principal component (PC) 1 comprises 36.8% of data variation caused by drought (Figure 2A). Similarly, PC2 comprises 21.9% of data variation and is caused by differences between genotypes. Similarly, for cytosolic-enriched samples, 25.9% of data variation is explained in PC1 and appears to be caused by drought (Figure 2B). Unlike organelle-enriched samples, cytosolic-enriched samples do not appear to separate by genotype on PC2. These data suggest that drought impacts both organelle- and cytoplasm-enriched samples, and that the drought response of organelles may be more genotype-specific compared with the cytosol. A statistical comparison of organelle-enriched samples identified 100 and 99 drought-responsive DEPs in RTx430 and BTx642, respectively, with an overlap of 38 (Figure 2D). Similarly, a statistical test of cytosolic-enriched samples using limma [21,22] identified 45 and 42 drought-responsive differentially expressed proteins (DEPs, adj. p-value < 0.05) in RTx430 and BTx642, respectively, with an overlap of 19 (Figure 2F). Previous reports suggest that the drought stress responses of these genotypes are largely similar, but that the magnitude of RTx430′s response is larger than BTx642′s [18]. Consistent with this observation, the cytosolic- and organelle-enriched drought responses showed positive correlation (r2 = 0.5, slope = 0.5, Figure 2C,E). However, the DEPs exclusive to RTx430′s stress response showed a much larger change during drought (Figure 2C,E).

2.3. Protein Profiles in Response to Drought Support Previous Findings

To evaluate whether our protein profiles accurately captured a drought response, we first performed Gene Ontology (GO) enrichment [23] using DEPs that increase significantly >2-fold in response to drought (limma adj. p-value < 0.05). Drought-relevant GO terms, including response to heat (GO:0009408), response to osmotic stress (GO:0006970), and response to abiotic stimulus (GO:0009628) were significantly enriched among upregulated DEPs (Fisher’s Exact, Bonferroni adj. p-value < 0.05) (Supplemental Table S1). Proteins belonging to these GO terms include multiple small heat shock proteins (SHSPs) (e.g., Sb03G003530 and Sb01G040030) [24], as well as the water stress-induced dehydrin C5YX70 all of which were significantly >2-fold more abundant in response to drought in both genotypes (adj. p-value < 0.05). We also found that our data support previous findings. For example, the SHSPs Sb06G017850, Sb03G006870, Sb07G028370, Sb03G003530, and Sb04G035130 were significantly >2-fold more abundant in response to drought in cytosolic-enriched RTx430 and BTx642 samples, and were similarly upregulated in the study described in Johnson et al. 2014 [25]. Similarly, the dehydrin protein Sb09G018420, thioredoxin protein Sb03G036980, and glucose-1-phosphate adenylyltransferase Sb09G029610 were significantly >2-fold more abundant in our experiments and elsewhere [25]. We also observed a significant >2-fold increase in the AAA domain-containing protein Sb03G001130, and ribosomal L16 protein Sb01G036330 in cytosolic-enriched RTx430 and BTx642 during drought stress, consistent with observations in Abdel-Ghany et al., 2020 [26]. A cross reference of our observations with literature findings can be found in Supplemental Table S1. These findings suggest our experiment corroborates previous research findings and accurately captured protein profiles representative of sorghum’s drought response.

2.4. Genotype-Specific Drought Responses Involved in Flowering Time Control, Starch Biosynthesis, and Rubisco Activation May Explain RTx430′s Drought Tolerance

A comparison of significant drought-induced protein abundance changes in RTx430 with BTx642 identified multiple proteins that changed exclusively in RTx430 (Table 1). For example, the likely CDC73/PHP protein Sb01G012610 was significantly >100-fold more abundant in RTx430 organelle-enriched samples in response to drought and did not change significantly in BTx642 (Figure 3). The nearest Arabidopsis ortholog of Sb01G012610 is AtCDC73/PHP (At3G22590), and loss of AtCDC73 function caused reduced H3K4me3 histone methylation patterns at the FLOWERING LOCUS C (FLC), misregulation of FLC expression, and an accelerated transition from vegetative to reproductive growth [27,28]. These data suggest the RTx430-specific increase in CDC73-like Sb01G012610 may function similarly to suppress FLC-mediated transition from vegetative to reproductive growth during drought.
We also observed a significant 65-fold drought-responsive increase in the ADP-Glucose pyrophosphorylase (Sb09G029610) exclusively in organelle-enriched RTx430 samples (Figure 3). The closest Arabidopsis ortholog of Sb09G029610 is the glucose-APL2 (At1G27680), which catalyzes the first rate-limiting step in starch synthesis [29]. Starch biosynthesis plays an active role in drought tolerance in Z. mays, and loss of function in similar ADP-glucose pyrophosphorylases (e.g., SH2) leads to soluble sugar starvation and reduced leaf growth rates [30]. We observed a significant 61-fold drought-responsive increase in the rubisco activase (RCA) protein Sb05G027880 exclusively in RTx430. In C3 and C4 grasses, RCA plays a critical role in maintaining photosynthesis during heat and drought stress [31]. Taken together, these data suggest that altered starch metabolism, Rubisco activation, and suppression of the transition from vegetative to reproductive growth may be preflowering drought adaptation strategies specific to S. bicolor RTx430.

2.5. Genotype-Specific Drought Responses Involved in Abscisic Acid Signaling May Explain RTx430′s Drought Tolerance

We observed a significant 34-fold increase in the likely abscisic acid receptor PYL9 Sb04G009280 exclusively in organelle-enriched RTx430 samples (Figure 3). Interestingly, drought was shown to decrease PYL9 Sorghum bicolor L. Moench [32] (annotated as PYL8 in Dalal et al. 2014), and slightly downregulated in S. bicolor BTx642 in response to saline–alkali stress [33] (annotated as PYL6 in Ma et al. 2019). Nearest orthologs of Sb04G009280 are PYL9 genes (Z. mays Zm00001d016105 and Arabidopsis At1G01360 [34]), and Z. mays PYL9 overexpression lines were more resistant to drought [35,36]. We also observed a significant 11-fold decrease in the zeaxanthin epoxidase (ZEP, Sb06G018220) exclusively in organelle-enriched RTx430 during drought. Zeaxanthin epoxidase catalyzes the first step in ABA biosynthesis and plays a critical role in xanthophyll production [37] (Figure 3). In model grasses, decreased ABA production correlated with drought tolerance during vegetative [57] and reproductive growth [58]. We observed a significant 6.7-fold increase in the likely abscisic acid stress-ripening (ASR) protein Sb06G016540 exclusively in RTx430 organelle-enriched samples during drought. The Sb06G016540 protein is orthologous to Z. mays ASR3 and confers drought tolerance when expressed in Arabidopsis [38] and tobacco [39]. These data suggest that part of RTx430′s preflowering drought tolerance includes prevention of an excessive stress-induced ABA response by decreasing ABA production, while increasing PYL9 abundance may function to maintain sensitivity to ABA and stomatal control.

2.6. Other RTx430-Specific Preflowering Drought Responses Include Proteins Involved in ROS Scavenging, HSPs, Epicuticular Wax Production, and Phospholipid Metabolism

We observed a significant 65-fold increase in the superoxide dismutase (SOD) protein Sb01G035350 in organelle-enriched RTx430 in response to preflowering drought (Figure 3). Interestingly, the same protein was significantly 18-fold decreased in organelle-enriched BTx642 during drought. Orthologs of SOD in Arabidopsis and rice include AtCSD1 (At1g08830) and SODCC1 (Os03g22810), and have been shown to play a role in adaptation to high light and drought [40]. Another ROS scavenging enzyme, Ferredoxin (Sb06G015570), was significantly >100-fold more abundant exclusively in RTx430 organelles during drought. Ferredoxin was also upregulated in the drought-tolerant sorghum line SC56 compared with the drought-sensitive line Tx7000 [41]; it confers heat stress tolerance in Arabidopsis [42] and salt tolerance in rice [43]. We also observed a significant 44-fold increase in the thioredoxin protein Sb03G036980 (NRX1) exclusively in RTx430 organelles during drought. The NRX1 protein was previously shown to be drought-inducible in sorghum [14] and orthologous proteins (e.g., Arabidopsis AtNRX1 (At1G60420) and Z. mays Zm00001d012591) have been characterized. The Arabidopsis AtNRX1 was shown to guard sensitive enzymes from oxidative stress, and knockouts had reduced catalase activity and increased sensitivity to oxidative stress [44].
Another oxidative stress relevant protein, the sHSP Sb04G035130 was significantly 45-fold increased exclusively in organelle-enriched RTx430 samples (Figure 3). Sequence homology suggests Sb04G035130 is orthologous to HSP17.4B (Arabidopsis ortholog At1G54050 and Z. mays ortholog Zm00001d018298), which was shown to be involved in oxidative stress adaptation [45]. Similarly, we observed a significant 21-fold increase in the likely HSP90 protein C5YBL4 exclusively in drought stressed organelle-enriched RTx430 samples. The orthologs of Sb06G000660 are HSP90 proteins in Arabidopsis and Z. mays (At5g52640 and Zm00001d024903, respectively) and it was shown that AtHSP90 partly controls stomatal ontogenesis and heat stress adaptation by interacting with YODA to transduce the heat stress response, including phosphorylation of MPK6 and SPCH [46].
We observed a significant RTx430-specific 34-fold increase in 3-ketoacyl-CoA synthase (UniProt identifier A0A1B6PJK5), which catalyzes the first step in very-long-chain fatty acids (VLCFAs) necessary for epicuticular wax production [47]. Epicuticular wax is important for drought tolerance in multiple species [48], and its abundance in grass crops correlated with grain yield under drought stress [49], including sorghum [50,51]. Multiple drought-relevant proteins became significantly less abundant exclusively in organelle-enriched RTx430 during water deficit. For example, we observed a significant 43-fold reduction in the likely phospholipase A2 protein Sb01G040430 (PLA2) in organelle-enriched RTx430 samples during drought (Figure 3). Phospholipases hydrolyze membrane phospholipids, and they play a role in maintaining membrane integrity as well as production of fatty acid precursors of drought-relevant signaling molecules [52]. Mutant Arabidopsis plants lacking PLA2 show delayed stomata opening, and are increased during drought stress in the wheat [53]. We also observed a significant 93-fold decrease in the HMA domain-containing protein Sb04G021590 exclusively in RTx430 organelle-enriched samples during drought. Orthologous proteins in Z. mays and Arabidopsis (At3G56240 and Zm00001d016691, respectively) encode a copper chaperone CCH [54,55] important in maintaining electron transport and photosynthesis [56], and are downregulated during drought at the transcript level in Sorghum bicolor [9].

3. Discussion

Drought is the primary stress causing crop yield loss. Comparing closely related genotypes with differing drought tolerance has proved to be a successful way to identify plant adaptation strategies to water deficit stress. Two Sorghum bicolor genotypes, RTx430 and BTx642, are characterized as preflowering drought-tolerant and -sensitive, respectively. Multiple reports have examined the differences between each genotype in response to water limitation but have relied on whole-organ and transcript-level measurements. To further understand RTx430′s unique drought tolerance, we report here a close examination of the protein-level drought response in organelle- and cytosolic-enriched cellular compartments (CCs). Our analysis identified >3000 proteins in each CC that are enriched for proteins characteristic of each compartment, indicating a successful capture of organelle- and cytosolic-protein profiles. Our data set supports previous findings, including upregulation of multiple drought-responsive proteins previously characterized in the literature. Further examination of our data suggests that organelle-enriched protein profiles respond to drought in a genotype-specific manner more than cytosolic-enriched proteins and may shed light on water stress adaptation strategies unique to RTx430.
For example, drought caused an RTx430-specific increase in Sb01G012610 in organelle-enriched samples. Sequence homology suggests Sb01G012610 is a nuclear-localized PLANT HOMOLOGOUS TO PARAFIBROMIN (PHP/CDC73) protein, which has been characterized best in Arabidopsis as a component of the Paf1c complex [27]. The Paf1c complex is involved in regulation of transcription by altering histone methylation status and RNA processing [59,60,61]. Evidence suggests that Paf1c participates in controlling the transition from vegetative to reproductive growth via histone methylation at the FLOWERING LOCUS C (FLC) locus. The FLC is a well-characterized repressor of flowering, and it was recently shown that CDC73 associates with and modifies the histone methylation status of chromatin at the FLC locus, resulting in increased FLC expression [27,28]. Mutants lacking CDC73 exhibited reduced FLC expression and accelerated flowering times [27]. Differences between RTx430 and BTx642 histone methylation profiles during drought stress were recently reported [20], and the RTx430-specific increase in the CDC73 ortholog Sb01G012610 during drought observed here suggests a possible mechanism by which RTx430 suppresses flowering during water limitation and remains in a vegetative growth phase.
Drought tolerance requires maintaining photosynthetic rates despite decreasing cellular water availability, a process carried out in part by the Rubisco chaperone Rubisco Activase (RCA). RCA activates Rubisco by catalyzing the ATP-dependent removal of inhibitory sugar phosphates bound to RuBisCo [62], and RCA’s activity is sensitive to heat, chloroplast redox status, and ATP/ADP ratios [31]. Arabidopsis mutants with decreased RCA activity exhibited decreased growth and photosynthetic rates [63], and RCA expression is upregulated in drought-tolerant rice genotypes compared with drought-sensitive lines [64]. Furthermore, RCA protein levels are increased by heat stress in wheat [65]. Reviewed by Feller 2016, increased thermal stability of RCA at higher temperatures correlated with increased temperature optimum of photosynthesis, and consequently, RCA is often a target for breeding drought-tolerant crops [66]. In sorghum, rubisco activase gene expression is downregulated under heat stress, particularly for drought-sensitive lines [67]. While Rubisco and RCA activity was not measured here, the observation that RCA increased exclusively in S. bicolor RTx430 suggests that Rubisco maintenance contributes to RTx430′s preflowering drought adaptation strategy.
Starch biosynthesis occurs during photoperiods and is catabolized into soluble sugars that support respiration and growth at night. Drought typically increases starch content in leaves, and soluble sugars act as a carbon source, signal molecule, and osmolyte [68]. We also observed an RTx430-specific increase in the starch biosynthesis protein ADP-glucose pyrophosphorylase Sb09G029610 (AtAPL2 (At1g27680) and ZmAPL2 (Zm00001d039131) ortholog) during drought. While drought-induced changes in APL transcript or protein abundance is plant-specific (reviewed by Thalman et al., 2017 [68]), mutants impaired in starch metabolism often show impaired growth and drought tolerance [69]. In other grass crops like maize, APL2 mutants showed decreased drought tolerance, drought recovery [30], and growth. Our findings suggest that RTx430 adapts to preflowering drought stress in part by increasing starch biosynthesis.
We also observed RTx430-specific changes in proteins involved in ABA metabolism. ABA is a well-characterized stress hormone typically synthesized de novo in response to drought [19,70]. The first step in ABA biosynthesis is the epoxidation of zeaxanthin to antheraxanthin and eventually violaxanthin by zeaxanthin epoxidase (ZEP) [71]. Multiple studies have shown that drought-tolerant grasses accumulate less ABA during drought at vegetative and reproductive stages [57,58], and lower ABA in foliar tissues is often used as selection criteria for breeding drought-tolerant crops like maize [72]. It was shown in Arabidopsis that ZEP protein is predominantly chloroplastically localized and undergoes degradation upon drought stress [73,74], but Schwarz et al. report no significant change in ABA content, and posit that the decrease in ZEP functions to prevent zeaxanthin consumption, leaving additional zeaxanthin for production of photoprotective xanthophylls. Because ABA signaling and crosstalk is complex, and no measurements of ABA was made in this study, it remains to be tested whether the observed decrease in ZEP during preflowering drought reduced ABA levels in S. bicolor RTx430 leaves. However, our data suggests that decreased ZEP abundance in RTx430 but not BTx642 may be a possible adaptation strategy wherein RTx430 reduces ABA production and prevents an excessive stress response.
The primary receptors responsible for ABA perception and signal transduction are the PYROBACTIN RESISTANCE 1-Like (PYL) family proteins [75,76]. In guard cells, ABA binding to PYL results in inhibition of PP2C phosphatases and subsequent activation of SnRK2 kinases that leads to stomatal closure and decreased transpiration [77]. The PYL receptors are transcriptionally downregulated in drought-sensitive S. bicolor lines [32,33] and rice [78] but drought-inducible in Z. mays [36] and ABA-inducible in Arabidopsis [75,76]. Overexpression of PYL9 in Arabidopsis [35], rice [78], and maize [36] resulted in increased drought tolerance or ABA sensitivity. Our observation that S. bicolor PYL9 (Sb04G009280) significantly increased in genotype RTx430 and not in BTx642 suggests that ABA receptors are involved in RTx430′s preflowering drought tolerance. The only likely SnRK2 [32] identified here, SB06G033990, did not change significantly in either genotype or CC, and no likely PP2C’s were detected. It may be that the increased PYL observed here functions to maintain sensitivity to decreased ABA caused by lower ZEP protein abundance, allowing proper stomatal closure.
We also observed an RTx430-specific increase in the likely ASR (abscisic acid-, stress-, and ripening-induced) protein during preflowering drought. The Z. mays ASR (ZmASR3) was found to confer additional drought tolerance and longer roots when overexpressed in Arabidopsis [38]. The improved drought tolerance of Arabidopsis ZmASR3 overexpressing lines was attributed in part to a decreased ROS by increasing superoxide dismutase and catalase activity in leaves [38]. Similarly, other ASR orthologs in wheat (TaASR1) [39], tomato (ASR1) [79], lily (LLA23) [80], rice (OsAsr1) [81], and plantain (MpAsr) [82] conferred additional drought tolerance when overexpressed in tobacco, Arabidopsis, rice, and Arabidopsis, respectively. The S. bicolor RTx430-specific increase in the likely ASR (Sb06G016540) observed here suggests additional ASR production contributes to RTx430′s preflowering drought tolerance.
Epicuticular waxes are a complex mixture of lipids composed mostly of VLCFAs, and their production is critical for drought tolerance and photoprotection under high light in multiple species including maize, wheat, and soybean [48]. The first step in VLCFA biosynthesis is catalyzed by β-ketoacyl CoA synthase (KCS), and mutants lacking KCS function exhibit decreased drought tolerance in multiple species including rice [47,83] and sorghum [50,51,84]. The observation here that RTx430 exclusively increased KCS during drought suggests epicuticular wax production contributes to RTx430′s unique drought tolerance. Other proteins identified here that are known to increase drought tolerance include the ferredoxins. Ferredoxins are important in ROS scavenging, and their overexpression confers increased drought tolerance in multiple systems [42,43]. Our observations suggest RTx430′s unique drought tolerance may involve producing additional ROS scavenging ferredoxins.
To our knowledge, this study constitutes the first comprehensive protein profiling of sorghum subcellular compartments in response to drought. Consistent with previous observations, our findings suggest that while RTx430 and BTx642 respond similarly to drought, the magnitude of RTx430′s response is greater than BTx642 and these differences may be responsible for RTx430′s preflowering drought tolerance. Our data also suggest that genotype-specific changes in organelle protein profiles may be greater compared with the cytoplasm. The proteins identified as changing only in organelle-enriched RTx430 samples during drought implicate multiple biological processes in RTx430′s unique preflowering drought tolerance and provides multiple avenues for breeding crops with improved drought tolerance.

4. Materials and Methods

4.1. Plant Growth and Drought Stress

Plants used in this study were grown and collected by Xu et al., 2018 [85]. Briefly, two sorghum (S. bicolor (L.) Moench) genotypes, RTx430 and BTx642, were grown in Parlier, CA in a randomized block design. The average monthly rainfall at the field site during the growth season is approximately zero, and soils are characterized as sandy loam with silky substratum. All watering was accurately and uniformly controlled via drip lines at the furrow of each plant. Watered control plants received consistent watering at 80% of calculated evapotranspiration, as measured by an onsite weather station ~1 km from the field site. Preflowering drought was imposed by halting irrigation upon seedling emergence until sample collection 8 weeks later. Drought stressed plants therefore received no additional water for 8 weeks upon seedling emergence from the soil, while stress-free watered plants received consistent drip line irrigation based on local evapotranspiration measurements. Drip line emitters were spaced 0.3 m apart with an output of 2 L/h, and watering was conducted once a week. Less than 1 mm of rainfall was recorded during the growth period at the experiment site. Control and drought stressed sorghum leaves were manually excised at 8 weeks postemergence and frozen in liquid nitrogen prior to protein isolation. Samples used in this study included 3 biological replicates of S. bicolor RTx430 drought and control, as well as 2 control and 3 drought biological replicates of S. bicolor BTx642. Each biological replicate is a pool of 10 individuals from each treatment group, collected from the same row of plants, at the same date and time. The percent of plant-available soil moisture was estimated using soil matric potential sensors (Watermark sensors, Irrometer Corp., Riverside, CA, USA) placed 15 cm beneath the soil and centered within plant rows. The degree of water stress was measured using the crop water stress index (CWSI) for both drought and watered treatments [85]. Additional detailed metrics used to evaluate the extent of drought as well as the resulting plant phenotype can be found in Xu et al. 2018 [85].

4.2. Organelle and Cytosolic Enrichment

Organelle and cytosolic protein fractions were obtained using a modified protocol adapted from O’Green et al. 2011 [86] and our histone purification protocol [20]. Briefly, leaves were homogenized by bead beating at −80 °C, and resuspended in buffer EB1 (0.44 M Sucrose, 10 mM Tris-HCl pH 8.0, 5 mM DTT) containing protease inhibitor cocktail (Sigma 5892791001, St. Louis, MO, USA). Plant debris was filtered through mesh 100 cloth (Sigma F6801, St. Louis, MO, USA) and organelles were pelleted by centrifugation at 3000× g for 10 min at 4 °C. The supernatant containing the cytosolic fraction was collected and frozen at −80 °C. The organelle pellet was resuspended in cold buffer EB2 (0.25 M sucrose, 10 mM MgCl2, 10 mM Tris-HCl pH 8.0, 1% Triton X100, 5 mM DTT) containing protease inhibitor cocktail and incubated on ice for 10 min with gentle mixing. Organelles were pelleted by centrifugation at 2100× g for 15 min at 4 °C, and supernatant was decanted. Pellet was then washed twice by resuspension and centrifugation in buffer EB2B (0.25 M sucrose, 10 mM MgCl2, 10 mM Tris-HCl pH 8.0, 5 mM DTT) at 2100× g for 15 min at 4 °C. Nuclei were lysed by resuspension in 200 µL of high salt nuclear lysis buffer (400 mM NaCl, 10 mM Tris-HCl pH 8, 1% Triton x100, 20 mM EDTA). Protein was frozen at −80 °C until quantification and digestion.

4.3. Protein Digestion

To remove plant metabolites that may interfere with downstream processes, protein from cytosolic-enriched samples were precipitated by MPLEx [87] and resuspended in 50 mM ammonium bicarbonate containing 8 M urea. The protein yield of organelle-enriched samples was relatively low and contained fewer contaminating metabolites; therefore, it was immediately resuspended in 50 mM ammonium bicarbonate with 8 M urea to prevent protein loss during MPLEx precipitation. Protein from all samples was then quantified by bicinchoninic acid assay (BCA), and disulfide bonds were reduced by addition of dithiothreitol to the final concentration of 10 mM and incubation at 60 °C for 30 min with shaking. Cysteine residues were alkylated by addition of iodoacetamide to a final concentration of 40 mM and incubation with shaking at 37 °C for 3 h. Protein was trypsin digested at a 50:1 (w/w) protein to trypsin ratio for 3 h at 37 °C with shaking at 800 rpm. Peptides were then desalted using C18 SPE Phenomenex (Torrance, CA, USA), dried down, and resuspended in water and stored at −80 °C.

4.4. LC-MS/MS Analysis

Peptides were separated using a Waters nanoAcquity LC (Waters Corporation, Milford, MA, USA) with binary solvent mobile phase A (0.1% (v/v) formic acid in water) and mobile phase B (0.1% FA (v/v) in water). 0.5 µg peptides were loaded onto a C18 reversed phase trap column (5 cm length, 100 μm i.d, 3 μm particle, 300 Å pore size, Phenomenex, Torrance, CA, USA) and washed with mobile phase A for 10 min at 3 µL/min. Then, the peptides are loaded onto a C18 analytical column (70 cm, 75 μm i.d, same material as the trap column). The analytical gradient was ramped from 1% to 45% mobile phase B over 100 min at 300 nL/min. Eluting peptides were injected in-line into a Q Exactive Orbitrap mass spectrometer (Thermo Fisher Scientific, Bremen, Germany) via electrospray ionization with data-dependent acquisition. The scan range was limited to 300–1800 m/z in positive mode, with resolution of 70,000. MS2 spectra were acquired with resolution of 17,500, scan range 200–2000 m/z, and 30 s dynamic exclusion. Raw spectra and associated data were submitted to the Mass Spectrometry Interactive Virtual Environment (MassIVE; [https://massive.ucsd.edu]) accession number MSV000086332.

4.5. Protein Search and Data Analysis

Peptide sequence identification was carried out using MaxQuant [88] Andromeda [89]. Search parameters included variable methionine oxidation, N-terminal acetylation, and fixed carbamidomethylation. Proteins were identified with 1% false discovery rate (FDR) using a comprehensive S. bicolor protein database from Uniprot [90] (downloaded September 2020) and iBAQ normalized. Data were then further processed in R using the methods described by Zhang et al., 2018 [22]. Data were prefiltered to remove contaminants and proteins identified by only a single unique peptide. Only proteins identified in all biological replicates of at least one treatment group (e.g., all RTx430 control, or all RTx430 drought, etc.) were included in subsequent analyses. Data were subjected to variance stabilized normalization (VSN, Supplementary Figure S1), and missingness was evaluated. Proteins with missing values were nonrandom and predominately occurred for proteins near the limit of detection (Supplementary Figures S2 and S3). We therefore chose to impute missing values by random sampling from a left-shifted Gaussian distribution (shift = 1.8, scale = 0.3) (Supplementary Figures S2 and S3). Differentially expressed proteins (DEPs) were identified by applying an empirical Bayes linear model using limma [21,91]. Statistical analysis also included principal component analysis (PCA) of cytoplasmic- and organelle-enriched samples. For comparison of cytoplasmic- and organelle-enriched samples for Figure 1, data were Z-score transformed in Perseus [92] (χ-µ/σ, where χ = observed value, µ = mean of all samples, σ = standard deviation of all samples). A comprehensive list of proteins identified here, their relative abundances, statistics, orthologous proteins, ontology enrichment, and alternative identifiers from UniProt, STRING, and Ensembl Plants can be found in Supplemental Table S1.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/24/9706/s1. Figure S1: Data normalization and goodness of fit; Figure S2: Missingness caused by low abundance proteins in organelle-enriched samples; Figure S3: Missingness caused by low abundance proteins in cytosolic-enriched samples; Table S1: Main data analysis table containing the normalized protein abundance values, relevant statistics, fold-changes, GO enrichment.

Author Contributions

Conceptualization, A.J.O., M.Z., and P.P.H.; data curation, A.J.O.; formal analysis, A.J.O.; funding acquisition, A.J.O. and P.P.H.; investigation, A.J.O., S.A., K.M.E., M.Z., and P.P.H.; methodology, A.J.O., S.A., K.M.E., M.Z., and P.P.H.; project administration, A.J.O., M.Z., and P.P.H.; resources, A.J.O., M.Z., and P.P.H.; supervision, A.J.O., M.Z., and P.P.H.; validation, A.J.O.; visualization, A.J.O.; writing—original draft, A.J.O.; writing—review and editing, A.J.O., K.M.E., M.Z., and P.P.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research is part of the Mission SEED Initiative at the Pacific Northwest National Laboratory (PNNL). It was conducted under the Laboratory Directed Research and Development Program at PNNL, a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830. This work was partially performed in the W. R. Wiley Environmental Molecular Sciences Laboratory (EMSL) (grid.436923.9), a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research and located at the Pacific Northwest National Laboratory (PNNL).

Acknowledgments

We would like to thank the EPICON project (award DE-SC0014081) for allowing us to use samples previously collected. We thank Thomas Fillmore and Ronald Moore for assistance with proteomics data collection, and Amir H. Ahkami for valuable discussions. We specially thank Peggy G. Lemaux for her generosity assisting with EPICON project guidance and her inspiration for collaborative research.

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.

Abbreviations

DEPDifferentially Expressed Proteins
ABAAbscisic Acid
GOGene Ontology
FLCFLOWERING LOCUS C
RCARubisco Activase
ASRAbscisic Acid Stress-Ripening
sHSPSmall Heat Shock Protein
ZEPZeaxanthin Epoxidase
SODSuperoxide Dismutase

References

  1. Gupta, A.; Rico-Medina, A.; Caño-Delgado, A.I. The physiology of plant responses to drought. Science 2020, 368, 266–269. [Google Scholar] [CrossRef] [PubMed]
  2. Dbattisti, D.S.; Naylor, R.L. Historical Warnings of Future Food Insecurity with Unprecedented Seasonal Heat. Science 2009, 323, 240–244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kuwayama, Y.; Thompson, A.; Bernknopf, R.; Zaitchik, B.; Vail, P. Estimating the Impact of Drought on Agriculture Using the U.S. Drought Monitor. Am. J. Agric. Econ. 2018, 101, 193–210. [Google Scholar] [CrossRef]
  4. United Nations, Department of Economic and Social Affairs. World Population Prospects Highlights; 2019; ISBN 978-92-1-148316-1. Available online: https://www.un.org/development/desa/publications/world-population-prospects-2019-highlights.html (accessed on 15 October 2020).
  5. Basu, S.; Ramegowda, V.; Kumar, A.; Pereira, A. Plant adaptation to drought stress. F1000Research 2016, 5, 1554. [Google Scholar] [CrossRef]
  6. Okamoto, M.; Peterson, F.C.; DeFries, A.; Park, S.-Y.; Endo, A.; Nambara, E.; Volkman, B.F.; Cutler, S.R. Activation of dimeric ABA receptors elicits guard cell closure, ABA-regulated gene expression, and drought tolerance. Proc. Natl. Acad. Sci. USA 2013, 110, 12132–12137. [Google Scholar] [CrossRef] [Green Version]
  7. Fàbregas, N.; Lozano-Elena, F.; Blasco-Escámez, D.; Tohge, T.; Martínez-Andújar, C.; Albacete, A.; Osorio, S.; Bustamante, M.; Riechmann, J.L.; Nomura, T.; et al. Overexpression of the vascular brassinosteroid receptor BRL3 confers drought resistance without penalizing plant growth. Nat. Commun. 2018, 9, 1–13. [Google Scholar] [CrossRef] [Green Version]
  8. Rodrigues, J.; Inzé, D.; Nelissen, H.; Saibo, N.J.M. Source–Sink Regulation in Crops under Water Deficit. Trends Plant Sci. 2019, 24, 652–663. [Google Scholar] [CrossRef]
  9. Fracasso, A.; Trindade, L.M.; Amaducci, S. Drought stress tolerance strategies revealed by RNA-Seq in two sorghum genotypes with contrasting WUE. BMC Plant Biol. 2016, 16, 1–18. [Google Scholar] [CrossRef]
  10. Ananda, G.K.S.; Myrans, H.; Norton, S.L.; Gleadow, R.; Furtado, A.; Henry, R.J. Wild Sorghum as a Promising Resource for Crop Improvement. Front. Plant Sci. 2020, 11, 1108. [Google Scholar] [CrossRef]
  11. Kebede, H.; Subudhi, P.K.; Rosenow, D.T.; Nguyen, H.T. Quantitative trait loci influencing drought tolerance in grain sorghum (Sorghum bicolor L. Moench). Theor. Appl. Genet. 2001, 103, 266–276. [Google Scholar] [CrossRef]
  12. Sabadin, P.K.; Malosetti, M.; Boer, M.P.; Tardin, F.D.; Santos, F.G.; Guimaraes, C.T.; Gomide, R.L.; Andrade, C.L.T.; Albuquerque, P.E.P.; Caniato, F.F.; et al. Studying the genetic basis of drought tolerance in sorghum by managed stress trials and adjustments for phenological and plant height differences. Theor. Appl. Genet. 2012, 124, 1389–1402. [Google Scholar] [CrossRef] [PubMed]
  13. Ogbaga, C.C.; Stępień, P.; Dyson, B.C.; Rattray, N.J.W.; Ellis, D.I.; Goodacre, R.; Johnson, G.N. Biochemical Analyses of Sorghum Varieties Reveal Differential Responses to Drought. PLoS ONE 2016, 11, e0154423. [Google Scholar] [CrossRef] [PubMed]
  14. Jedmowski, C.; Ashoub, A.; Beckhaus, T.; Berberich, T.; Karas, M.; Brüggemann, W. Comparative Analysis of Sorghum bicolor Proteome in Response to Drought Stress and following Recovery. Available online: https://www.hindawi.com/journals/ijpro/2014/395905/ (accessed on 28 September 2020).
  15. Fadoul, H.E.; El Siddig, M.A.; Abdalla, A.W.H.; Hussien, A.A. Physiological and proteomic analysis of two contrasting Sorghum bicolor genotypes in response to drought stress. Aust. J. Crop. Sci. 2018, 12, 1543–1551. [Google Scholar] [CrossRef]
  16. Thomas, H.; Howarth, C.J. Five ways to stay green. J. Exp. Bot. 2000, 51, 329–337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Smith, R.H.; Bhaskaran, S.; Miller, F.R. Screening for drought tolerance in Sorghum using cell culture. Vitr. Cell. Dev. Biol. Anim. 1985, 21, 541–545. [Google Scholar] [CrossRef]
  18. Varoquaux, N.; Cole, B.; Gao, C.; Pierroz, G.; Baker, C.R.; Patel, D.; Madera, M.; Jeffers, T.; Hollingsworth, J.; Sievert, J.; et al. Transcriptomic analysis of field-droughted sorghum from seedling to maturity reveals biotic and metabolic responses. Proc. Natl. Acad. Sci. USA 2019, 116, 27124–27132. [Google Scholar] [CrossRef] [Green Version]
  19. Sah, S.K.; Reddy, K.R.; Li, J. Abscisic Acid and Abiotic Stress Tolerance in Crop Plants. Front. Plant Sci. 2016, 7, 571. [Google Scholar] [CrossRef] [Green Version]
  20. Zhou, M.; Malhan, N.; Ahkami, A.H.; Engbrecht, K.; Myers, G.; Dahlberg, J.; Hollingsworth, J.; Sievert, J.A.; Hutmacher, R.; Madera, M.; et al. Top-down mass spectrometry of histone modifications in sorghum reveals potential epigenetic markers for drought acclimation. Methods 2019. [Google Scholar] [CrossRef]
  21. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  22. Zhang, X.; Smits, A.H.; van Tilburg, G.B.; Ovaa, H.; Huber, W.; Vermeulen, M. Proteome-wide identification of ubiquitin interactions using UbIA-MS. Nat. Protoc. 2018, 13, 530. [Google Scholar] [CrossRef]
  23. Thomas, P.D.; Muruganujan, A.; Ebert, D.; Huang, X.; Thomas, P.D. PANTHER version 14: More genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 2019, 47, D419–D426. [Google Scholar] [CrossRef]
  24. Nagaraju, M.; Reddy, P.S.; Kumar, S.A.; Kumar, A.; Rajasheker, G.; Rao, D.M.; Kishor, P.B.K. Genome-wide identification and transcriptional profiling of small heat shock protein gene family under diverse abiotic stress conditions in Sorghum bicolor (L.). Int. J. Biol. Macromol. 2020, 142, 822–834. [Google Scholar] [CrossRef] [PubMed]
  25. Johnson, S.M.; Lim, F.-L.; Finkler, A.; Fromm, H.; Slabas, A.R.; Knight, M.R. Transcriptomic analysis of Sorghum bicolor responding to combined heat and drought stress. BMC Genom. 2014, 15, 1–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Abdel-Ghany, S.E.; Ullah, F.; Ben-Hur, A.; Reddy, A.S. Transcriptome Analysis of Drought-Resistant and Drought-Sensitive Sorghum (Sorghum bicolor) Genotypes in Response to PEG-Induced Drought Stress. Int. J. Mol. Sci. 2020, 21, 772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Yu, X.; Michaels, S.D. The Arabidopsis Paf1c Complex Component CDC73 Participates in the Modification of FLOWERING LOCUS C Chromatin. Plant Physiol. 2010, 153, 1074–1084. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Park, S.; Oh, S.; Ek-Ramos, J.; Van Nocker, S. Plant Homologous to Parafibromin is a Component of the PAF1 Complex and Assists in Regulating Expression of Genes within H3K27ME3-Enriched Chromatin. Plant Physiol. 2010, 153, 821–831. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Ventriglia, T.; Kuhn, M.L.; Ruiz, M.T.; Ribeiro-Pedro, M.; Valverde, F.; Ballicora, M.A.; Preiss, J.; Romero, J.M. Two Arabidopsis ADP-Glucose Pyrophosphorylase Large Subunits (APL1 and APL2) Are Catalytic. Plant Physiol. 2008, 148, 65–76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. AbdelGawad, H.; Avramova, V.; Baggerman, G.; Van Raemdonck, G.; Valkenborg, D.; Van Ostade, X.; Guisez, Y.; Prinsen, E.; Asard, H.; Ende, W.V.D.; et al. Starch biosynthesis contributes to the maintenance of photosynthesis and leaf growth under drought stress in maize. Plant Cell Environ. 2020, 43, 2254–2271. [Google Scholar] [CrossRef]
  31. Perdomo, J.A.; Capó-Bauçà, S.; Carmo-Silva, A.E.; Galmés, J. Rubisco and Rubisco Activase Play an Important Role in the Biochemical Limitations of Photosynthesis in Rice, Wheat, and Maize under High Temperature and Water Deficit. Front. Plant Sci. 2017, 8, 490. [Google Scholar] [CrossRef] [Green Version]
  32. Dalal, M.; Inupakutika, M. Transcriptional regulation of ABA core signaling component genes in sorghum (Sorghum bicolor L. Moench). Mol. Breed. 2014, 34, 1517–1525. [Google Scholar] [CrossRef]
  33. Ma, S.; Lv, L.; Meng, C.; Zhou, C.; Fu, J.; Shen, X.; Zhang, S.C.; Li, Y. Genome-Wide Analysis of Abscisic Acid Biosynthesis, Catabolism, and Signaling in Sorghum Bicolor under Saline-Alkali Stress. Biomolecules 2019, 9, 823. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zhao, Y.; Zhang, Z.; Gao, J.; Wang, P.; Hu, T.; Wang, Z.; Hou, Y.-J.; Wan, Y.; Liu, W.; Xie, S.; et al. Arabidopsis Duodecuple Mutant of PYL ABA Receptors Reveals PYL Repression of ABA-Independent SnRK2 Activity. Cell Rep. 2018, 23, 3340–3351.e5. [Google Scholar] [CrossRef] [PubMed]
  35. Zhao, Y.; Chan, Z.; Gao, J.; Xing, L.; Cao, M.; Yu, C.; Hu, Y.; You, J.; Shi, H.; Zhu, Y.; et al. ABA receptor PYL9 promotes drought resistance and leaf senescence. Proc. Natl. Acad. Sci. USA 2016, 113, 1949–1954. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. 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]
  37. Shaw, A.K.; Bhardwaj, P.K.; Ghosh, S.; Roy, S.; Saha, S.; Sherpa, A.R.; Saha, S.K.; Hossain, Z. β-aminobutyric acid mediated drought stress alleviation in maize (Zea mays L.). Environ. Sci. Pollut. Res. 2015, 23, 2437–2453. [Google Scholar] [CrossRef]
  38. Liang, Y.; Jiang, Y.; Du, M.; Li, B.; Chen, L.; Chen, M.; Jin, D.; Wu, J. ZmASR3 from the Maize ASR Gene Family Positively Regulates Drought Tolerance in Transgenic Arabidopsis. Int. J. Mol. Sci. 2019, 20, 2278. [Google Scholar] [CrossRef] [Green Version]
  39. Hu, W.; Huang, C.; Deng, X.; Zhou, S.; Chen, L.; Li, Y.; Wang, C.; Ma, Z.; Yuan, Q.; Wang, Y.; et al. TaASR1, a transcription factor gene in wheat, confers drought stress tolerance in transgenic tobacco. Plant Cell Environ. 2013, 36, 1449–1464. [Google Scholar] [CrossRef]
  40. Xing, Y.; Cao, Q.; Zhang, Q.; Qin, L.; Jia, W.; Zhang, J. MKK5 Regulates High Light-Induced Gene Expression of Cu/Zn Superoxide Dismutase 1 and 2 in Arabidopsis. Plant Cell Physiol. 2013, 54, 1217–1227. [Google Scholar] [CrossRef] [Green Version]
  41. Azzouz-Olden, F.; Hunt, A.G.; Dinkins, R.D. Transcriptome analysis of drought-tolerant sorghum genotype SC56 in response to water stress reveals an oxidative stress defense strategy. Mol. Biol. Rep. 2020, 47, 3291–3303. [Google Scholar] [CrossRef] [Green Version]
  42. Lin, Y.-H.; Huang, L.-F.; Hase, T.; Huang, H.-E.; Feng, T.-Y. Expression of plant ferredoxin-like protein (PFLP) enhances tolerance to heat stress in Arabidopsis thaliana. New Biotechnol. 2015, 32, 235–242. [Google Scholar] [CrossRef]
  43. Huang, H.-E.; Ho, M.-H.; Chang, H.; Chao, H.-Y.; Ger, M.-J. Overexpression of plant ferredoxin-like protein promotes salinity tolerance in rice (Oryza sativa). Plant Physiol. Biochem. 2020, 155, 136–146. [Google Scholar] [CrossRef] [PubMed]
  44. Kneeshaw, S.; Keyani, R.; Delorme-Hinoux, V.; Imrie, L.H.; Loake, G.J.; Le Bihan, T.; Reichheld, J.-P.; Spoel, S.H. Nucleoredoxin guards against oxidative stress by protecting antioxidant enzymes. Proc. Natl. Acad. Sci. USA 2017, 114, 8414–8419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Sewelam, N.A.; Kazan, K.; Hüdig, M.; Maurino, V.G.; Schenk, P.M. The AtHSP17.4C1 Gene Expression Is Mediated by Diverse Signals that Link Biotic and Abiotic Stress Factors with ROS and Can Be a Useful Molecular Marker for Oxidative Stress. Int. J. Mol. Sci. 2019, 20, 3201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Samakovli, D.; Tichá, T.; Vavrdová, T.; Ovečka, M.; Luptovčiak, I.; Zapletalová, V.; Kuchařová, A.; Křenek, P.; Krasylenko, Y.; Margaritopoulou, T.; et al. YODA-HSP90 Module Regulates Phosphorylation-Dependent Inactivation of SPEECHLESS to Control Stomatal Development under Acute Heat Stress in Arabidopsis. Mol. Plant 2020, 13, 612–633. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, X.; Guan, Y.; Zhang, D.; Dong, X.; Tian, L.; Qu, L.Q. A β-Ketoacyl-CoA Synthase Is Involved in Rice Leaf Cuticular Wax Synthesis and Requires a CER2-LIKE Protein as a Cofactor. Plant Physiol. 2017, 173, 944–955. [Google Scholar] [CrossRef] [Green Version]
  48. Xue, D.; Zhang, X.; Lu, X.; Chen, G.; Chen, Z. Molecular and Evolutionary Mechanisms of Cuticular Wax for Plant Drought Tolerance. Front. Plant Sci. 2017, 8, 621. [Google Scholar] [CrossRef]
  49. Dhanyalakshmi, K.H.; Soolanayakanahally, R.Y.; Rahman, T.; Tanino, K.K.; Nataraja, K.N. Leaf Cuticular Wax, a Trait for Multiple Stress Resistance in Crop Plants. Abiotic Biotic Stress Plants 2019. [Google Scholar] [CrossRef] [Green Version]
  50. Jordan, W.R.; Shouse, P.J.; Blum, A.; Miller, F.R.; Monk, R.L. Environmental Physiology of Sorghum. II. Epicuticular Wax Load and Cuticular Transpiration. Crop. Sci. 1984, 24, 1168–1173. [Google Scholar] [CrossRef]
  51. Awika, H.O.; Hays, D.B.; Mullet, J.E.; Rooney, W.L.; Weers, B.D. QTL mapping and loci dissection for leaf epicuticular wax load and canopy temperature depression and their association with QTL for staygreen in Sorghum bicolor under stress. Euphytica 2017, 213, 207. [Google Scholar] [CrossRef] [Green Version]
  52. Matos, A.R.; Pham-Thi, A.-T. Lipid deacylating enzymes in plants: Old activities, new genes. Plant Physiol. Biochem. 2009, 47, 491–503. [Google Scholar] [CrossRef]
  53. Verlotta, A.; Liberatore, M.T.; Cattivelli, L.; Trono, D. Secretory Phospholipases A2 in Durum Wheat (Triticum durum Desf.): Gene Expression, Enzymatic Activity, and Relation to Drought Stress Adaptation. Int. J. Mol. Sci. 2013, 14, 5146–5169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Chai, L.-X.; Dong, K.; Liu, S.-Y.; Zhang, Z.; Zhang, X.-P.; Tong, X.; Zhu, F.-F.; Zou, J.-Z.; Wang, X.-B. A putative nuclear copper chaperone promotes plant immunity in Arabidopsis. J. Exp. Bot. 2020, 71, 6684–6696. [Google Scholar] [CrossRef] [PubMed]
  55. Shin, L.-J.; Lo, J.-C.; Yeh, K. Copper Chaperone Antioxidant Protein1 Is Essential for Copper Homeostasis. Plant Physiol. 2012, 159, 1099–1110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Burkhead, J.L.; Reynolds, K.A.G.; Abdel-Ghany, S.E.; Cohu, C.M.; Pilon, M. Copper homeostasis. New Phytol. 2009, 182, 799–816. [Google Scholar] [CrossRef] [PubMed]
  57. Wang, Z.; Huang, B.; Bonos, S.A.; Meyer, W.A. Abscisic Acid Accumulation in Relation to Drought Tolerance in Kentucky Bluegrass. HortScience 2004, 39, 1133–1137. [Google Scholar] [CrossRef]
  58. Ji, X.; Dong, B.; Shiran, B.; Talbot, M.J.; Edlington, J.E.; Hughes, T.; White, R.G.; Gubler, F.; Dolferus, R. Control of Abscisic Acid Catabolism and Abscisic Acid Homeostasis Is Important for Reproductive Stage Stress Tolerance in Cereals. Plant Physiol. 2011, 156, 647–662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. He, Y.; Doyle, M.R.; Amasino, R.M. PAF1-complex-mediated histone methylation of FLOWERING LOCUS C chromatin is required for the vernalization-responsive, winter-annual habit in Arabidopsis. Genes Dev. 2004, 18, 2774–2784. [Google Scholar] [CrossRef] [Green Version]
  60. Xu, L.; Zhao, Z.; Dong, A.; Soubigou-Taconnat, L.; Renou, J.-P.; Steinmetz, A.; Shen, W.-H. Di- and Tri- but Not Monomethylation on Histone H3 Lysine 36 Marks Active Transcription of Genes Involved in Flowering Time Regulation and Other Processes in Arabidopsis thaliana. Mol. Cell. Biol. 2007, 28, 1348–1360. [Google Scholar] [CrossRef] [Green Version]
  61. Oh, S.; Park, S.; Van Nocker, S. Genic and Global Functions for Paf1C in Chromatin Modification and Gene Expression in Arabidopsis. PLoS Genet. 2008, 4, e1000077. [Google Scholar] [CrossRef] [Green Version]
  62. Archie, R.P. Rubisco activase–Rubisco’s catalytic chaperone. Photosynth. Res. 2003, 75, 11–27. [Google Scholar] [CrossRef]
  63. Eckardt, N.A.; Snyder, C.W.; Jr, A.R.P.; Ogren, W.L. Growth and Photosynthesis under High and Low Irradiance of Arabidopsis thaliana Antisense Mutants with Reduced Ribulose-1,5-Bisphosphate Carboxylase/Oxygenase Activase Content. Plant Physiol. 1997, 113, 575–586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Fukayama, H.; Ueguchi, C.; Nishikawa, K.; Katoh, N.; Ishikawa, C.; Masumoto, C.; Hatanaka, T.; Misoo, S. Overexpression of Rubisco Activase Decreases the Photosynthetic CO2 Assimilation Rate by Reducing Rubisco Content in Rice Leaves. Plant Cell Physiol. 2012, 53, 976–986. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Law, R.; Crafts-Brandner, S.J. High Temperature Stress Increases the Expression of Wheat Leaf Ribulose-1,5-Bisphosphate Carboxylase/oxygenase Activase Protein. Arch. Biochem. Biophys. 2001, 386, 261–267. [Google Scholar] [CrossRef] [PubMed]
  66. Feller, U. Drought stress and carbon assimilation in a warming climate: Reversible and irreversible impacts. J. Plant Physiol. 2016, 203, 84–94. [Google Scholar] [CrossRef] [Green Version]
  67. Djanaguiraman, M.; Prasad, P.V.V.; Murugan, M.; Perumal, R.; Reddy, U.K. Physiological differences among sorghum (Sorghum bicolor L. Moench) genotypes under high temperature stress. Environ. Exp. Bot. 2014, 100, 43–54. [Google Scholar] [CrossRef]
  68. Thalmann, M.; Santelia, D. Starch as a determinant of plant fitness under abiotic stress. New Phytol. 2017, 214, 943–951. [Google Scholar] [CrossRef] [Green Version]
  69. Stitt, M.; Zeeman, S.C. Starch turnover: Pathways, regulation and role in growth. Curr. Opin. Plant Biol. 2012, 15, 282–292. [Google Scholar] [CrossRef]
  70. Roychoudhury, A.; Paul, S.; Basu, S. Cross-talk between abscisic acid-dependent and abscisic acid-independent pathways during abiotic stress. Plant Cell Rep. 2013, 32, 985–1006. [Google Scholar] [CrossRef]
  71. Yamamoto, H.; Nakayama, T.; Chichester, C. Studies on the light and dark interconversions of leaf xanthophylls. Arch. Biochem. Biophys. 1962, 97, 168–173. [Google Scholar] [CrossRef]
  72. Landi, P.; Sanguineti, M.; Conti, S.; Tuberosa, R. Direct and Correlated Responses to Divergent Selection for Leaf Abscisic Acid Concentration in Two Maize Populations. Crop. Sci. 2001, 41, 335–344. [Google Scholar] [CrossRef]
  73. Schwarz, N.; Armbruster, U.; Iven, T.; Brückle, L.; Melzer, M.; Petit-Houdenot, Y.; Jahns, P. Tissue-Specific Accumulation and Regulation of Zeaxanthin Epoxidase in Arabidopsis Reflect the Multiple Functions of the Enzyme in Plastids. Plant Cell Physiol. 2014, 56, 346–357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. North, H.M.; Frey, A.; Boutin, J.-P.; Sotta, B.; Marion-Poll, A. Analysis of xanthophyll cycle gene expression during the adaptation of Arabidopsis to excess light and drought stress: Changes in RNA steady-state levels do not contribute to short-term responses. Plant Sci. 2005, 169, 115–124. [Google Scholar] [CrossRef]
  75. Park, S.-Y.; Fung, P.; Nishimura, N.; Jensen, D.R.; Fujii, H.; Zhao, Y.; Lumba, S.; Santiago, J.; Rodrigues, A.; Chow, T.-F.F.; et al. Abscisic Acid Inhibits Type 2C Protein Phosphatases via the PYR/PYL Family of START Proteins. Science 2009, 324, 1068–1071. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Ma, Y.; Szostkiewicz, I.; Korte, A.; Moes, D.; Yang, Y.; Christmann, A.; Grill, E. Regulators of PP2C Phosphatase Activity Function as Abscisic Acid Sensors. Science 2009, 324, 1064–1068. [Google Scholar] [CrossRef]
  77. Lee, S.C.; Lim, C.W.; Lan, W.; He, K.; Luan, S. ABA Signaling in Guard Cells Entails a Dynamic Protein–Protein Interaction Relay from the PYL-RCAR Family Receptors to Ion Channels. Mol. Plant 2013, 6, 528–538. [Google Scholar] [CrossRef] [Green Version]
  78. Tian, X.; Wang, Z.Y.; Li, X.; Lv, T.; Liu, H.; Wang, L.; Niu, H.; Bu, Q. Characterization and Functional Analysis of Pyrabactin Resistance-Like Abscisic Acid Receptor Family in Rice. Rice 2015, 8, 1–13. [Google Scholar] [CrossRef] [Green Version]
  79. Kalifa, Y.; Perlson, E.; Gilad, A.; Konrad, Z.; Scolnik, P.A.; Bar-Zvi, D. Over-expression of the water and salt stress-regulated Asr1 gene confers an increased salt tolerance. Plant Cell Environ. 2004, 27, 1459–1468. [Google Scholar] [CrossRef]
  80. Yang, C.-Y.; Chen, Y.-C.; Jauh, G.Y.; Wang, C. A Lily ASR Protein Involves Abscisic Acid Signaling and Confers Drought and Salt Resistance in Arabidopsis. Plant Physiol. 2005, 139, 836–846. [Google Scholar] [CrossRef] [Green Version]
  81. Kim, S.-J.; Lee, S.-C.; Hong, S.K.; An, K.; An, G.; Kim, S.-R. Ectopic expression of a cold-responsive OsAsr1 cDNA gives enhanced cold tolerance in transgenic rice plants. Mol. Cells 2009, 27, 449–458. [Google Scholar] [CrossRef]
  82. Liu, H.-Y.; Dai, J.-R.; Feng, D.-R.; Liu, B.; Wang, H.; Wang, J.-F. Characterization of a Novel Plantain Asr Gene, MpAsr, that is Regulated in Response to Infection of Fusarium oxysporum f. sp. cubense and Abiotic Stresses. J. Integr. Plant Biol. 2010, 52, 315–323. [Google Scholar] [CrossRef]
  83. Yu, D.; Ranathunge, K.; Huang, H.; Pei, Z.; Franke, R.; Schreiber, L.; He, C. Wax Crystal-Sparse Leaf1 encodes a β–ketoacyl CoA synthase involved in biosynthesis of cuticular waxes on rice leaf. Planta 2008, 228, 675–685. [Google Scholar] [CrossRef] [PubMed]
  84. Burow, G.B.; Franks, C.D.; Xin, Z. Genetic and Physiological Analysis of an Irradiated Bloomless Mutant (Epicuticular Wax Mutant) of Sorghum. Crop. Sci. 2008, 48, 41–48. [Google Scholar] [CrossRef] [Green Version]
  85. Xu, L.; Naylor, D.; Dong, Z.; Simmons, T.; Pierroz, G.; Hixson, K.K.; Kim, Y.-M.; Zink, E.M.; Engbrecht, K.M.; Wang, Y.; et al. Drought delays development of the sorghum root microbiome and enriches for monoderm bacteria. Proc. Natl. Acad. Sci. USA 2018, 115, E4284–E4293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. O’Geen, H.; Echipare, L.; Farnham, P.J. Using ChIP-Seq Technology to Generate High-Resolution Profiles of Histone Modifications. In Epigenetics Protocols; Tollefsbol, T.O., Ed.; Methods in Molecular Biology; Humana Press: Totowa, NJ, USA, 2011; pp. 265–286. ISBN 978-1-61779-316-5. [Google Scholar]
  87. Nakayasu, E.S.; Nicora, C.D.; Sims, A.C.; Burnum-Johnson, K.E.; Kim, Y.-M.; Kyle, J.E.; Matzke, M.M.; Shukla, A.K.; Chu, R.K.; Schepmoes, A.A.; et al. MPLEx: A Robust and Universal Protocol for Single-Sample Integrative Proteomic, Metabolomic, and Lipidomic Analyses. mSystems 2016, 1, e00043-16. [Google Scholar] [CrossRef] [Green Version]
  88. Tyanova, S.; Temu, T.; Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protoc. 2016, 11, 2301–2319. [Google Scholar] [CrossRef]
  89. Cox, J.; Neuhauser, N.; Michalski, A.; Scheltema, R.A.; Olsen, J.V.; Mann, M. Andromeda: A Peptide Search Engine Integrated into the MaxQuant Environment. J. Proteome Res. 2011, 10, 1794–1805. [Google Scholar] [CrossRef]
  90. The UniProt Consortium UniProt: A hub for protein information. Nucleic Acids Res. 2015, 43, D204–D212. [CrossRef]
  91. Smyth, G.K. Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments. Stat. Appl. Genet. Mol. Biol. 2004, 3, 1–25. [Google Scholar] [CrossRef]
  92. Tyanova, S.; Temu, T.; Sinitcyn, P.; Carlson, A.; Hein, M.Y.; Geiger, T.; Mann, M.; Cox, J. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods 2016, 13, 731–740. [Google Scholar] [CrossRef]
Figure 1. Evaluation of cytosolic and organelle enrichment. Ribosomal proteins abundant in the cytoplasm were more abundant in cytosolic-enriched samples (A), while photosynthesis-related chloroplastic (B) and nuclear histone proteins (C) were more abundant in organelle-enriched samples. Color and scale bars represent average Z-score transformed protein abundances of all biological replicates for each sample. Protein identifiers correspond to STRING or UniProt accessions. Wat., Watered; Drt., Drought.
Figure 1. Evaluation of cytosolic and organelle enrichment. Ribosomal proteins abundant in the cytoplasm were more abundant in cytosolic-enriched samples (A), while photosynthesis-related chloroplastic (B) and nuclear histone proteins (C) were more abundant in organelle-enriched samples. Color and scale bars represent average Z-score transformed protein abundances of all biological replicates for each sample. Protein identifiers correspond to STRING or UniProt accessions. Wat., Watered; Drt., Drought.
Ijms 21 09706 g001
Figure 2. Principal component analysis of organelle- (A) and cytosolic-enriched protein profiles (B) in response to drought. Correlation between RTx430 and BTx642 drought responses in (C) organelle- and (E) cytoplasm-enriched samples. Venn diagram summary of proteins with adj. p-value < 0.05 in each genotype for (D) organelle- and (F) cytosolic-enriched samples. Purple data points are significantly changed in either genotype, while red points correspond to RTx430-specific changes (adj. p-value < 0.05). D, Drought; W, Watered; PC, Principal Component.
Figure 2. Principal component analysis of organelle- (A) and cytosolic-enriched protein profiles (B) in response to drought. Correlation between RTx430 and BTx642 drought responses in (C) organelle- and (E) cytoplasm-enriched samples. Venn diagram summary of proteins with adj. p-value < 0.05 in each genotype for (D) organelle- and (F) cytosolic-enriched samples. Purple data points are significantly changed in either genotype, while red points correspond to RTx430-specific changes (adj. p-value < 0.05). D, Drought; W, Watered; PC, Principal Component.
Ijms 21 09706 g002
Figure 3. Comparison of genotype-specific drought responses in organelle-enriched samples of RTx430 (A) and BTx642 (B). Red datapoints passed a limma adjusted p-value threshold of 0.05. Pink circles indicate proteins that changed significantly only in RTx430 (limma adj. p-value < 0.05). Protein names correspond to STRING identifiers, or when unavailable, UniProt accessions.
Figure 3. Comparison of genotype-specific drought responses in organelle-enriched samples of RTx430 (A) and BTx642 (B). Red datapoints passed a limma adjusted p-value threshold of 0.05. Pink circles indicate proteins that changed significantly only in RTx430 (limma adj. p-value < 0.05). Protein names correspond to STRING identifiers, or when unavailable, UniProt accessions.
Ijms 21 09706 g003
Table 1. RTx430-specific drought responsive proteins.
Table 1. RTx430-specific drought responsive proteins.
IdentifierProteinAbbr.Ens. OrthologFunctionRef.
Sb01g012610Plant Homologous to ParafribominCDC73/PHPAt3g22590, Zm00001d013597, Os03g0664700Histone methylation, FLC-regulation[27,28]
Sb09g029610ADP-Glc pyrophosphorylaseAPL2At1G27680, Zm00001d039131, Os05g0580000Starch biosynthesis[29,30]
Sb05g027880Rubisco activase alpha isoformRCAAt2G39730, Zm00001d048592, Photosynthesis[31]
Sb04g009280Pyrobactin Resistance 1-likePYL9At1G01360, Zm00001d016105, Os02g0255500ABA Signal Transduction[32,33,34,35,36]
Sb06g018220Zeaxanthin epoxidaseZEPAt5G67030, Os04g0448900ABA metabolism[37]
Sb06g016540Abscisic Acid Stress Ripening ASR2Zm00001d003712, Os04g0423400ABA metabolism[38,39]
Sb01g035350Superoxide dismutaseSODAt1G08830, Zm00001d047479, Os03g0351500ROS Scavenging[40]
Sb06g015570FerredoxinFdAt2G27510, Zm00001d003797, Os04g0412200ROS Scavenging[41,42,43]
Sb03g036980NRX1-like Thioredoxin domain-containingNRXZm00001d012591, Os01g0794400Oxidative Stress[14,44]
Sb04g035130Small heatshock protein 17sHSP17At1G54050, Zm00001d018298, Os02g0782500Protein Folding[45]
Sb06g000660Heatshock protein 90HSP90At5G52640, Zm00001d024903, Os04g0107900Protein Folding[46]
A0A1B6PJK53-ketoacyl-CoA synthaseKCSAt1G25450, Os04g0116800Epicuticular Wax Production[47,48,49,50,51]
Sb01g040430Phospholipase A2PLA2Zm00001d028505Fatty acid and membrane integrity[52,53]
Sb04G021590HMA domain-containing proteinHMAAt3G56240, Zm00001d016691Copper Ion Homeostasis[9,54,55,56]
Sb06g033990Sucrose non-fermenting 1-related prot. Kinase 2SnRK2At1G60940, Zm00001d026690, Os04g0691100ABA Signal Transduction[32]
Identifiers consist of STRING and UniProt accessions. Ensemble Orthologs, when available, are reported for Arabidopsis thaliana, Zea mays, and Oryza sativa. Ens., Ensemble; Prot., Protein; Ref., Reference; Abbr., Abbreviated Name; ADP, Adenosine diphosphate; Glc, Glucose.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ogden, A.J.; Abdali, S.; Engbrecht, K.M.; Zhou, M.; Handakumbura, P.P. Distinct Preflowering Drought Tolerance Strategies of Sorghum bicolor Genotype RTx430 Revealed by Subcellular Protein Profiling. Int. J. Mol. Sci. 2020, 21, 9706. https://doi.org/10.3390/ijms21249706

AMA Style

Ogden AJ, Abdali S, Engbrecht KM, Zhou M, Handakumbura PP. Distinct Preflowering Drought Tolerance Strategies of Sorghum bicolor Genotype RTx430 Revealed by Subcellular Protein Profiling. International Journal of Molecular Sciences. 2020; 21(24):9706. https://doi.org/10.3390/ijms21249706

Chicago/Turabian Style

Ogden, Aaron J., Shadan Abdali, Kristin M. Engbrecht, Mowei Zhou, and Pubudu P. Handakumbura. 2020. "Distinct Preflowering Drought Tolerance Strategies of Sorghum bicolor Genotype RTx430 Revealed by Subcellular Protein Profiling" International Journal of Molecular Sciences 21, no. 24: 9706. https://doi.org/10.3390/ijms21249706

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