Next Article in Journal
Effect of Sowing Date on the Growth Characteristics and Yield of Growth-Constrained Direct-Seeding Rice
Next Article in Special Issue
Pathogenesis-Related Proteins (PRs) with Enzyme Activity Activating Plant Defense Responses
Previous Article in Journal
Macro and Micro-Nutrient Accumulation and Partitioning in Soybean Affected by Water and Nitrogen Supply
Previous Article in Special Issue
Characterization of microRNAs and Target Genes in Musa acuminata subsp. burmannicoides, var. Calcutta 4 during Interaction with Pseudocercospora musae
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proteomic Insights of Cowpea Response to Combined Biotic and Abiotic Stresses

by
Daiane Gonzaga Ribeiro
1,
Ana Carolina Mendes Bezerra
2,
Ivonaldo Reis Santos
3,
Priscila Grynberg
2,
Wagner Fontes
4,
Mariana de Souza Castro
4,
Marcelo Valle de Sousa
4,
Maria Eugênia Lisei-de-Sá
2,
Maria Fatima Grossi-de-Sá
1,2,5,
Octávio Luiz Franco
1,6 and
Angela Mehta
2,*
1
Centro de Análises Proteômicas e Bioquímicas Programa de Pós-Graduação em Ciências Genômicas e Biotecnologia, Universidade Católica de Brasília (UCB), Brasília CEP 71966-700, DF, Brazil
2
Embrapa Recursos Genéticos e Biotecnologia, PBI, Av. W/5 Norte Final, Brasília CEP 70770-917, DF, Brazil
3
Programa de Pós-Graduação em Ciências Biológicas (Biologia Molecular), Instituto de Ciências Biológicas, Campus Universitário Darcy Ribeiro-UnB, Universidade de Brasília, Brasília CEP 70910-900, DF, Brazil
4
Laboratório de Bioquímica e Química de Proteínas, Departamento de Biologia Celular, Universidade de Brasília, Brasília CEP 70910-900, DF, Brazil
5
National Institute of Science and Technology, INCT PlantStress Biotech, Embrapa, Brasilia CEP 70770-917, DF, Brazil
6
S-Inova Biotech, Universidade Católica Dom Bosco (UCDB), Campo Grande CEP 79117-900, MS, Brazil
*
Author to whom correspondence should be addressed.
Plants 2023, 12(9), 1900; https://doi.org/10.3390/plants12091900
Submission received: 21 March 2023 / Revised: 19 April 2023 / Accepted: 25 April 2023 / Published: 6 May 2023
(This article belongs to the Special Issue Defense-Related Proteins of Higher Plants)

Abstract

:
The co-occurrence of biotic and abiotic stresses in agricultural areas severely affects crop performance and productivity. Drought is one of the most adverse environmental stresses, and its association with root-knot nematodes further limits the development of several economically important crops, such as cowpea. Plant responses to combined stresses are complex and require novel adaptive mechanisms through the induction of specific biotic and abiotic signaling pathways. Therefore, the present work aimed to identify proteins involved in the resistance of cowpea to nematode and drought stresses individually and combined. We used the genotype CE 31, which is resistant to the root-knot nematode Meloidogyne spp. And tolerant to drought. Three biological replicates of roots and shoots were submitted to protein extraction, and the peptides were evaluated by LC-MS/MS. Shotgun proteomics revealed 2345 proteins, of which 1040 were differentially abundant. Proteins involved in essential biological processes, such as transcriptional regulation, cell signaling, oxidative processes, and photosynthesis, were identified. However, the main defense strategies in cowpea against cross-stress are focused on the regulation of hormonal signaling, the intense production of pathogenesis-related proteins, and the downregulation of photosynthetic activity. These are key processes that can culminate in the adaptation of cowpea challenged by multiple stresses. Furthermore, the candidate proteins identified in this study will strongly contribute to cowpea genetic improvement programs.

1. Introduction

In agricultural areas, cultivated plants can be simultaneously exposed to a wide range of biotic and abiotic stress, such as pathogen infection and drought. Climate change scenarios foresee increases in the earth’s temperature and a decrease in the rainfall regime, which can drastically affect the yield of important crops, including cowpea [1]. Cowpea (Vigna unguiculata L. Walp.) is a legume widely consumed all over the world. Vigna grains have high amounts of protein and vitamins, representing an important food source for the populations in these regions. Cowpea also plays an important role in soil fertilization due to the plant’s ability to fix atmospheric nitrogen in association with rhizobial bacteria [2]. However, nematode infestation and water deficits are among the main factors limiting cowpea productivity. The disease caused by root-knot nematodes (RKN) can impair water and nutrient absorption, causing leaf chlorosis and damaging the root system [3]. A water deficit restricts several aspects of vegetative growth, causing a series of biochemical, physiological, and morphological responses [4]. The effects of cross-stress, therefore, can severely affect cowpea development and productivity.
Plants have developed several mechanisms to cope with stress, and the study of the responses to single stresses is quite evolved. However, in the case of simultaneous biotic and abiotic stresses, plant responses are more complex, as they include the interaction of two living organisms combined with stress conditions [5]. It has been reported that the molecular responses of plants subjected to combined stresses may be completely different from their responses to individual stresses [6,7]. Studies have shown that plant responses to multiple stresses involve the activation of kinases that regulate specific genes, transcription factors, reactive oxygen species, heat shock proteins, and the crosstalk of several hormonal signaling pathways [8,9,10,11].
V. unguiculata is, in general, resistant to M. incognita and tolerant to a water deficit [12,13]. Cowpea resistance to root-knot nematodes seems to be controlled by the dominant gene Rk, which was proven effective against some RKN isolates but did not show any effectiveness against all RKN populations [14,15]. Regarding drought tolerance, studies on cowpea have associated several genes and multiple quantitative trait loci (QTL) with drought tolerance [16,17,18,19]. However, these studies represent individual assessments of cowpea-nematode interactions and cowpea response to drought. Although key mechanisms involved in RKN resistance and drought tolerance, such as activation of R genes and oxidative stress/ubiquitination, respectively, were reported [20,21,22,23], to the best of our knowledge, there are still few studies investigating plant responses to both stresses simultaneously.
Prospective studies using proteomic approaches can bring important contributions to understanding the molecular mechanisms involved in the defense response of plants to co-applied stresses. Furthermore, the proteomics approach is a promising tool for crop improvement programs as it allows for the global evaluation of total proteins, the final product of gene expression, from samples under different biological conditions. Therefore, the objective of this work was to identify proteins potentially involved in V. unguiculata resistance to M. incognita and drought tolerance, individually and combined (cross-stress). To evaluate the plants, we performed physiological (gas exchange) and shotgun proteomic analyses. Discovering the complex modulation of plant responses to simultaneous biotic and abiotic stresses may contribute to the development of cultivars resistant to multiple stresses, which often occur in the field.

2. Results and Discussion

2.1. Phenotypic Evaluation and Gas Exchange of Vigna unguiculata in Response to M. incognita and a Water Deficit

Plants subjected to a water deficit showed leaf wilting and retardation of vegetative growth during progressive water loss irrespective of nematode inoculation as compared to the control (Figure 1a,c,d). Inoculation with nematodes, on the other hand, caused leaf yellowing regardless of a water deficit (Figure 1a,b,d), as expected [24,25].
One of the main physiological strategies to cope with drought and prevent water loss in plants is stomatal control. The low water availability in the soil induces the closure of stomata through the induction of root signals, consequently causing a reduction in gas exchange and carbon fixation rates [26]. According to our results, we can observe a reduction in transpiration (E), stomatal conductance (gs), and net photosynthetic rates (An) in nematode and drought-exposed plants (ND) and in plants only under drought (D) from six days after treatment (DAT) (Figure 2B–D). In addition, plants submitted exclusively to D showed even lower rates of An when compared to the combined stress treatment, showing reductions of ~30% at 10 and 72% at 12 DAT. It was possible to observe an increase in intercellular CO2 concentrations (Ci) in response to water limitation in the last days of measurement (Figure 2A). Increases in Ci as a response to drought were previously reported in severely water-stressed cowpea plants [27]. As stated before, plant responses to a single stress cannot be extended to combined stress in most cases. In the literature, studies report biotic and abiotic stresses either enhancing or reducing plants’ tolerance when co-applied [28,29]. Plant responses will depend on many factors, such as plant species, stage of development, and intensity of the stressor factor [7]. Similar to our results, other authors also reported the reduction of E, gs, and An in plants subjected to single and combined stresses [30,31]. However, contrary to our results, the effects of the cross-stress (drought + inoculation with cowpea severe mosaic virus) in cowpea at an early stage of infection led to a progressive reduction of Ci, whereas An was reduced by the stresses but remained similar when comparing the values of single and combined stresses [12].
In our study, nematode inoculation induced higher An rates when compared to control plants, presenting increases of ~1.4, 0.5, and 0.3-fold changes at 6, 8, and 10 DAT, respectively (Figure 2). It is known that M. incognita is an obligate parasite and, therefore, needs the cellular content of plants to feed. The presence of the parasite may stimulate the cellular machinery of the plant increasing energy production. On the other hand, control plants showed regular rates of the aforementioned physiological processes, producing energy for their development.

2.2. Protein Profile of Vigna unguiculata in Response to Combined Stress

Plants have developed a series of morphological, physiological, and molecular adaptations to survive in adverse conditions. The higher complexity of plants’ responses to combined biotic and abiotic stresses is due to the involvement of two living organisms, the plant and the pathogen, accompanied by the abiotic factor. Therefore, to better understand the mechanisms of resistance of cowpea to combined stress, the protein profiles of the roots and shoots were investigated. Proteomic analysis of cowpea roots and shoots revealed, respectively, 570 and 470 proteins differentially abundant in the three conditions evaluated (N, D, and ND). Principal component analysis and heatmaps were carried out for stress and control conditions, showing greater variability in the biological replicas under N (roots) and ND (shoots) (Figure 3). Gene ontology analysis showed several important biological processes in which differentially abundant proteins are involved, such as homeostatic processes, cell wall organization, transport, photosynthesis, and response to stimuli, among others (Figure 4 and Figure 5). In this study, important proteins from cowpea roots reported to be involved in common defense responses to biotic and abiotic stresses, such as the activation of the antioxidant mechanisms, production of NB-LRR proteins, and accumulation of secondary metabolites, were identified. Interestingly, some specific defense responses of cowpea subjected to cross-stress were also observed, such as the activation of the jasmonic acid hormonal signaling pathway and the intense production of PR proteins. These and other proteins appear to be important in activating efficient mechanisms of adaptation to cross-stress and will be discussed in more detail below.

2.2.1. Defense Proteins Orchestrate Cowpea Adaptation Mechanisms to Cross-Stress

The NBS-LRR proteins (nucleotide-binding site leucine-rich repeat) are encoded by one of the largest families of genes involved in plant disease resistance and recognize pathogen-derived effector proteins [32,33]. In addition to NBS-LRR, pathogenesis-related proteins (PR) are also essential for plant defense against pathogens. PRs tend to accumulate in the infected tissue, protecting plants from further infections. In this study, increased LRR and PR proteins were identified in ND plants (Table 1 and Table 2), both in the shoots and roots (LRR receptor-like serine/threonine-protein kinase FLS2, Log2FC (FC): 12; chitinase, FC: 6; pathogenesis-related protein, FC: 5; thaumatin, FC: 3). Previous studies have identified NBS-LRR and PRs proteins involved in resistance to both biotic and abiotic stresses when evaluated individually [34,35,36], but little is known about the presence of these proteins under cross-stresses [10,37]. A recent study reported that the overexpression of an A. stenosperma endochitinase (AsECHI) induces resistance to RKN and a water deficit, isolated and combined. The higher drought tolerance of these plants was achieved concomitantly with a reduction of ~30% in infection by RKN [38]. In this context, the increased NB-LRR proteins modulated in ND may activate cascades of effective defenses against both stresses. ND plants also seem to increase the production of PRs to minimize the oxidative effects caused by a water deficit and nematode infection. Interestingly, we identified a chitinase protein only under ND. This protein can be a key candidate potentially involved in RKN resistance and drought tolerance. This protein was not identified in cowpea inoculated only with nematodes. The proteins mentioned above are included in the stimulus response category (Figure 4 and Figure 5).
Other proteins related to cowpea resistance were increased in the roots (Table 1) of N, D, and ND plants (CPRD86, CPRD14, and CPRD2) and in the shoots (Table 2) of ND plants (CPRD2 and dehydrin). Cowpea clones responsive to dehydration (CPRD) proteins have already been identified in other cowpea varieties, such as EPACE [39,40]. Dehydrin proteins play an important role in dehydration tolerance, as they help to increase water uptake and water holding capacity, favoring the maintenance of photosynthetic rates and consequently reducing the production of reactive oxygen species (ROS) [41,42]. Even though the responsiveness to dehydration proteins is primarily related to single drought resistance, our results showed its increase under both cross and single stresses, appearing to be a common cowpea resistance response. Overall, our results suggest that some resistance mechanisms of cowpea subjected to drought and cross-stress seem to be shared and that other mechanisms seem to be unique to ND, such as the intense activation of PRs proteins identified only in ND. The imposed water deficit may increase cowpea resistance to RKN. The proteins involved in dehydration are included in the protein binding category (Figure 4 and Figure 5).

2.2.2. Antioxidant System Drives Plant Adaptation to Combined Stress

Reactive oxygen species are important signaling molecules of the signal transduction pathway that triggers stress defense responses, but the exacerbated production of these molecules can cause cellular damage [43]. Biotic and abiotic stresses cause dysregulation of the plant metabolism, which results in increased ROS generation affecting cellular redox homeostasis. ROS react with different cellular structures, such as the nucleus, proteins, and membranes, impairing their integrity, so its excessive production must be controlled by antioxidant systems to prevent damage and cell death. The activation of antioxidant enzymes crucial in the efficient elimination of ROS has become an essential marker of plant adaptation [44]. In this study, many antioxidant proteins were identified both in the shoots and in the roots of cowpea, especially under conditions of D and ND (Figure 4 and Figure 5—oxidoreductase activity category). Among the increased antioxidant proteins, we can highlight L-ascorbate peroxidase, glutathione transferase, and superoxide dismutase with fold-changes raging between two and nine (Table 1 and Table 2). Ascorbate peroxidase and glutathione transferase are antioxidant enzymes responsible for protecting the plasma membrane of cells [45,46]. In addition, it has been reported that one of the main antioxidant compounds strongly linked to plant adaptation to drought is ascorbate [47]. Other studies also reported the identification of several antioxidant proteins in plants subjected to both drought [48,49] and combined abiotic stresses [50,51] but there are few reports of antioxidant proteins in combined biotic and abiotic stresses [52,53]. Given the importance of antioxidant defense for stress resistance, it is notable that the activation of the ROS detoxification pathway in cowpea plants is key to preventing cellular damage caused by drought and root-knot nematodes.

2.2.3. Defense Responses Mediated by Hormonal Signaling in Cowpea

Plants do not have an adaptive immune system, but they do have a signaling network to ensure growth and development as well as adaptation to adverse conditions. Growth regulators play a crucial role in the reprogramming of complex mechanisms of adaptation to stresses, such as abscisic acid (ABA), jasmonic acid (JA), salicylic acid (SA), and ethylene (ET) [37,54,55]. ABA and JA are strongly related to stress tolerance since ABA-mediated signaling allows the regulation of stomatal closure, and JA-mediated signaling triggers plant defense responses through its accumulation in leaves [56,57]. Although JA is known to be related to biotic stress, many other studies reinforce the role of this phytohormone in response to drought, showing a significant and rapid increase of JA in plants under water stress conditions [58,59]. Our study identified the increase of some proteins involved in signaling via JA (Lipoxygenases, FC: 2–3) and ABA pathways (protein containing Bet_v_1 domain, FC: 2–6 and annexin, FC: 1.8) in plants subjected to D and ND in both the roots and shoots (Table 1 and Table 2). Studies on the evaluation of cross and individual stresses in Arabidopsis and tomato reported that the hormonal signaling pathways and the expression of defense genes are different for each stress [5,60]. Another report involving the same treatments in peanuts showed that JA and ABA were activated in response to RKN and drought stress, respectively, while ET was activated in Arachis plants subjected to cross-stress (nematode and drought) [38]. Here, we can observe that through the crosstalk network, JA and ABA phytohormones may be working together to regulate cowpea tolerance to drought and combined stress.
Interestingly, our study identified several proteins related to the abscisic acid pathway that were highlighted in cowpea tolerance to drought and cross-stress. These data corroborate the physiological analysis of stomatal conductance (gs), which showed reduced rates on all evaluated days, with a significant reduction on the last day.

2.2.4. Other Proteins’ Important in the Process of Adaptation to Drought and Multiple Stresses

During PTI and ETI induction, some signaling cascades are rapidly initiated, such as the activation of MAP kinases, production of reactive oxygen species, and transcriptional reprogramming, among others [61]. In plants, stimuli caused by adverse conditions are detected by specific receptors on the plasma membrane, initiating the signaling of the lipid cascade. Phospholipid-based signaling cascades are extremely important in signal transduction as phospholipids activate communication with the immune system preventing invasion by pathogens [62,63]. The main components of this signaling cascade are phospholipase C (PLC) and phospholipase D (PLD), which are also responsible for the turnover of phospholipids along with diacylglycerol kinase (DAK) [64]. Some studies have identified phospholipases involved in plant resistance against fungi and bacteria (mainly in effector recognition) [65,66], as well as in ABA-induced water stress tolerance [67,68,69]. In this study, some phospholipase C and phospholipase D (FC: 1.1-5, respectively) proteins were identified in cowpea plants subjected to D and ND in both roots and shoots (Table 1 and Table 2). Thus, it is likely that activation of signaling via phospholipids in cowpea plants is directly related to stomatal closure and reduced transpiration, as well as cell wall-based defenses against pathogenic infections. Protein phospholipases are included in the catalytic activity category (Figure 4 and Figure 5).
Secondary metabolites also play an important role in plant defense responses. Flavonoids are a large class of secondary metabolites present in plants and have several functions, such as plant development through the control of auxin, antioxidants, chemoattractants, and defense compounds, among others. Studies report that flavonoids can induce quiescence by decreasing the movement of nematodes, modifying their migration towards the root, repelling, and even killing them [70,71]. Flavonoids are also important in response to a water deficit, and their accumulation is often associated with drought resistance [72,73]. In our study, proteins involved in the flavonoid biosynthesis pathway were increased only in cowpea roots subjected to D (chalcone synthase, isoflavone reductase) and ND (flavonol 3-O-methyltransferase). Some authors associated the increase in plant flavonoids with abiotic stresses tolerance [74,75], but reports of alterations in metabolites involved in response to combined stresses are rare [76]. Cowpea plants subjected to D and ND seemed to accumulate secondary metabolites, such as flavonoids, important to eliminate cells infected by the parasite and damaged by drought. The proteins mentioned above were not identified in the cowpea-nematode interaction (N) in our study.

2.2.5. Reduction of Photosynthetic Proteins in Response to Combined Stress in Cowpea

Photosynthesis is one of the primary processes altered by exposure to stress. Due to its direct relation with energy production, reductions in photosynthesis, depending on the severity, may lead to reduced growth, yield, and even plant death [77]. In the shoots of cowpea, D differentially induced the increase of proteins associated with leaf photochemistry (Table 2), such as chlorophyll a-b binding protein (FC: 4.47), involved with the capture and delivery of excitation energy to the photosystems and a protein identified as photosystem II stability/assembly factor (FC: 1.3), essential for photosystem II (PSII) biogenesis. It is known that PSII is highly susceptible to excessive excitation energy and elevated levels of ROS and the rates of its damage and repair of PSII control photoinhibition [78]. Photoinhibition is reversible in its early stages, and the repair of PSII is mainly associated with the de novo biosynthesis of PSII photosynthetic proteins [79]. In contrast, two other leaf photosynthesis-related proteins were decreased due to the exposure of plants to drought: the photosystem I (PSI) subunit PsaN (FC: 1.8), involved in the docking of the mobile electron carrier plastocyanin [80,81], and the subunit 1 of chlorophyll a-b binding protein (FC: 1.3), also involved in the light-harvesting complex [82].
The simultaneous application of the stresses (ND), however, induced a completely different response. Among the differentially abundant proteins in the shoots of V. unguiculata under combined stress, 172 were decreased (Supplementary Table S4), with the coordinated decrease of several major photosynthetic proteins involved in the light reactions and carbon assimilation with at least 2-fold change (Table 2), including several proteins involved in the biosynthesis of chlorophyll, such as: magnesium-protoporphyrin IX monomethyl ester (oxidative) cyclase (FC 1.8), magnesium protoporphyrin IX methyltransferase activity (FC: 2.72), NADPH-protochlorophyllide oxidoreductase (FC: 2.2), and the enzyme geranylgeranyl reductase (FC: 2.4); two proteins involved in energy capture and transfer chlorophyll a/b binding proteins in the light-harvesting complex I (FC: 3.3 and 3.2, respectively); PSII oxygen-evolving enhancer protein units 1 (FC: 1.2), 2 (FC: 3.3), and 3 (FC: 2.2); photosystem I chlorophyll a apoproteins A1, PSI-subunit VI, and the PSI iron-sulfur center were down-regulated with FC of 2.4, 2.0, and 1.5, respectively; two components of the electron carrier cytochrome b6-f complex, cytochrome f (petA; FC: 1.5) and cytochrome b6 (petB, FC: 1.8); the chloroplastic enzyme ferredoxin-NADP reductase (FC: 1.79); and the subunits alpha and beta of ATP synthase (FC: 1.1 and 0.9); among other. A fragment of the molecular chaperone of rubisco, rubisco activase, and ribulose bisphosphate carboxylase large chain were also reduced by FC 1.4 and 0.8, respectively. A progressive reduction of rubisco activity was reported under high temperatures (>35 °C) due to a decrease in rubisco activase activity [83]. The proteins involved in photosynthesis are represented in the cellular metabolic process category (Figure 5).
Thus, it was possible to identify that under single and cross-stress, cowpea plants adopted different strategies to cope. While under single stress, plants seemed to increase sink capacity and invest in the repair of PSII photodamaged units to try and maintain carbon fixation rates, while under combined stress, we observed a coordinated decrease of light and assimilatory reactions. The coordinated decrease of photosynthetic transcripts involved in leaf photochemistry, carbon reduction, and pigment synthesis was reported as a defense mechanism of plants against biotic stresses. The authors hypothesized that the turnover of photosynthetic proteins allows the redirection of plants’ energy and reallocation of nitrogen to the induction of defense proteins against the pathogen while maintaining the losses in carbon assimilation to only a moderate level [84]. In our study, the downregulation of photosynthetic proteins was observed exclusively under the ND condition (Table 2). Even though the reduced activity of photosynthesis lowered carbon fixation rates and subsequently negatively affected plant productivity, in terms of plant biology, it may also represent an important strategy to cope with the stresses, indicating that the damage caused by the production and accumulation of toxic products of photorespiration and/or reactive intermediates threatening to cause permanent photodamage could outweigh the benefit of sustaining the high carbon fixation rates [85]. Both hypotheses are consistent with our gas exchange results. On the one hand, we can observe a reduction in photosynthetic assimilatory rates (An) in all plants submitted to water restriction (D and ND) when compared to control plants from 10 and 12 DAT concomitant with changes in the protein profile, characterized by the increased expression of pathogen defense-related proteins under ND. On the other hand, levels of An were severely reduced at 12 DAT under D when compared to ND (showing reductions of 85%), indicating that in the long term, the coordinated reduction of photosynthetic proteins might be a more efficient strategy to cope with drought and avoid photodamage.

3. Material and Methods

3.1. Plant Material and Growing Conditions

The cultivar of V. unguiculata CE 31, which is resistant to M. incognita race 3 and tolerant to water deficits [12,13], was used in this study. These authors reported that cowpea plants were challenged with M incognita (race 3) and collected after sixty days of cultivation. The egg mass index analysis was scored according to [86], and the degree of resistance was determined according to [87]. Furthermore, cowpea CE 31 has been reported as tolerant to drought in other works and used in the evaluation of other combined stresses [12,22].
Seeds were disinfected with 3% (v/v) sodium hypochlorite for 5 min and soaked in distilled water for 20 min. Germination was carried out using Germitest® paper (neutral pH, 28 × 38 cm), and the seeds were kept in the dark for the first three days until the seedlings’ emergence. Seedlings were transplanted into 500 mL plastic cups containing a mixture of autoclaved substrate and sand (1:1) and kept in a growth chamber under a photoperiod of 16 h/8 h (light/dark) and a temperature of 25 °C. During the first 20 days after transplantation, the plants were irrigated with tap water, and from the 21st day on, plants were watered only with Hoagland’s nutrient solution (1:10) [88]. We evaluated plants inoculated with nematodes (N), plants submitted to a water deficit (D), plants inoculated with nematodes and a water deficit simultaneously (ND), and control plants (C). Control plants were not inoculated with nematodes and received a water regime to maintain 70% of FC (field capacity). Treatments were initiated 25 days after transplantation, and three biological replicas were used for each treatment. In the first three days after nematode inoculation and interruption of irrigation, no physiological measurements were performed, as this time was estimated for parasite penetration into plant roots. Gas exchange measurements were performed during the last nine days of the experiment, and all treatments were evaluated and compared to the control. Whole plants were collected for proteomic analysis. Plants were separated into shoots (all leaves and stems) and roots (whole root) and collected 12 days after treatment when they presented an average of 2 to 3 trifoliate leaves.

3.2. Plant-Nematode Interaction Assay (N)

The Meloidogyne incognita race 3 population was previously cultivated in tomato roots (cv. Santa Cruz) in a greenhouse. To extract the nematodes, the roots were triturated in a blender with 0.5% (v/v) sodium hypochlorite, as previously described [89]. Second-stage juvenile nematodes (J2) were collected using modified Baermann funnels, and nematode counting was performed with a light microscope using Peter’s counting slides. Cowpea plants were inoculated with 1000 J2 25 days after transplantation to evaluate the plant-pathogen interaction. Shoots and roots of plants inoculated with nematodes were collected 12 days after inoculation and stored at −80 °C for further evaluation.

3.3. Water Deficit Assay (D)

Field capacity (FC) was determined considering the weight difference between the wet soil after saturation and free drainage and the dry soil. FC maintenance was measured daily in all plants by weighing the cups and replacing the volume of water lost through transpiration using a scale with a capacity of 20 kg. Twenty-five days after transplantation, irrigation was interrupted for 12 consecutive days. During this period, soil water gradually reduced to 25% of FC. Shoots and roots of the plants subjected to drought were collected 12 days after irrigation interruption and stored at −80 °C for further evaluation.

3.4. Cross-Stress Assay (Nematode and Drought; ND)

Cowpea plants at 25 days of age were inoculated with 1000 J2 and subjected to drought simultaneously. Irrigation interruption was performed for 12 consecutive days until the plants reached 25% of FC. The shoots and roots of the plants were collected 12 days after nematode inoculation and irrigation interruption and stored at −80 °C for further evaluation.

3.5. Evaluation of Gas Exchange

Gas exchange measurements were performed in the second fully expanded trefoil of V. unguiculata in the mornings (from 8 am to 12 pm.) using a portable infrared gas analyzer system (IRGA LCpro-SD—ADC BioScientific, Hoddesdon, England). Intercellular CO2 partial pressure (Ci), transpiration rate (E), stomatal conductance (gs), and net photosynthesis rate (An) were recorded. The environmental conditions inside the IRGA chamber during the evaluation were: photosynthetic photon flux density (PPFD), 1000 µmol m−2 s−1; temperature, 28 °C; and ambient CO2 concentration. Gas exchange measurements were performed from the sixth day after nematode inoculation, a period sufficient for the penetration of the parasite into the roots of the plants. We presented data on gas exchange parameters at 6, 8, 10, and 12 days after treatment (DAT).

3.6. Statistical Analysis

The homogeneity of variances was previously analyzed using the Cochran homogeneity test. Whenever the data distribution was not homogeneous, box-cox transformation was used. After that, data were submitted to one-way analysis of variance (ANOVA), and in the presence of significant differences, means were compared using Tukey’s test with a 5% confidence level (p ≤ 0.05).

3.7. Protein Extraction, Digestion, and Desalting

Protein extraction was performed following Ribeiro et al., 2022 [21]. For the extraction of total proteins, 300 mg of plant tissue (root and shoots) were solubilized in an extraction buffer, and later phenol was added. The samples were agitated for 15 min and centrifuged at 10,000× g for 3 min. Proteins were precipitated in 0.1 M ammonium acetate in methanol and washed with 80% acetone. For protein digestion, the extracted proteins were solubilized with 60 µL of 50 mM ammonium bicarbonate (NH4HCO3 pH 8.5) and 25 µL of RapiGestTM SF—Waters (0.2% v/v). The samples were heated at 80 °C for 15 min, under agitation, and 2.5 µL of 100 mM dithiothreitol was added. For protein alkylation, 2.5 µL of 300 mM iodoacetamide was added, and the samples were kept at room temperature for 30 min in the dark. After reduction and alkylation, the proteins (approximately 80 µg) were hydrolyzed with trypsin (1 µg) and incubated in an oven at 37 °C for 19 h. Later, 10 µL of 5% trifluoroacetic acid (TFA) was added, and the samples were incubated again at 37 °C for 90 min. Samples were centrifuged at 10 °C for 30 min, and the supernatant was recovered and dried. Protein quantification was performed with Qubit® Fluorometer (Invitrogen), according to the manufacturer’s manual. Subsequently, the samples were desalted following the protocol of Rappsilber et al., 2007 [90] with modifications [91]. First, the tips were prepared (P200 µL) with Empore® C18 disks, with one disk in each tip. Second, the tips containing the disks were washed with 20 µL of methanol and centrifuged at 2000× g for 30 s. Third, 20 µL of POROS® R2 resin in solvent B (0.1% formic acid/98% acetonitrile) was added, and the tips were centrifuged at 2000× g for 30 s. Each tip was washed two times with 20 µL of solvent A (0.1% formic acid/2% acetonitrile) followed by centrifugation. After tip preparation, the samples were solubilized in solvent A, inserted in the tips, and centrifuged at 2000× g for 2 min. The peptides were eluted with 20 µL of solvent B and dried in a vacuum centrifuge for 45 min. Then, the peptides were solubilized with 0.1% formic acid and injected into the ESI LC-MS/MS.

3.8. Chromatography and Mass Spectrometry Analysis

The peptides obtained were injected into a chromatographic system (Dionex Ultimate 3000 RSLCnano UPLC, Thermo, Waltham, MA, USA) configured with a trap column (3 cm × 100 µm) containing C18 particles with 5 µm, 120 Å (ReprosilPur, Dr. Maich GmbH, Ammerbuch, Germany). The samples were injected to obtain 2 µg in the column and submitted to a linear gradient of elution between solvents A (0.1% formic acid in 2% acetonitrile/water) and B (0.1% formic acid in 80% acetonitrile/water) from 2% B to 35% B for 155 min. The fractions separated in the chromatographic system were eluted directly at the ionization source of an Orbitrap Elite mass spectrometer (Thermo, USA) configured for data-dependent acquisition mode. MS1 spectra were acquired on the orbitrap analyzer with a resolution of 120,000 and a range between 300 and 1650 m/z. The 15 most intense ions were fragmented, generating MS2 spectra with a resolution of 15,000 [92,93]. The reanalysis of already fragmented ions was inhibited by dynamic exclusion, favoring the identification of less abundant peptides.

3.9. Quantitative and Qualitative Analysis of Proteins

Chromatogram alignment and peptide quantification were performed using the Progenesis QI for Proteomics software. Protein identification was performed using Peaks® 7.0 software (Bioinformatics Solutions Inc., Waterloo, ON, Canada) with the following parameters: database obtained from the UniProt repository (Universal protein), filtered for Vigna unguiculata (Taxon ID 3917), and submitted to the removal of redundant sequences using the FASTAtools software (http://lbqp.unb.br/LBQPtools/; accessed on 15 July 2021). The search was performed based on de novo and PSM sequencing with the following parameters: tolerance for precursor mass of 10 ppm, fragments of 0.05 Da, tolerance of up to 2 missed cleavages, carbamidomethylation of cysteines as fixed modification, and oxidation of methionine as a variable modification. Protein identification was considered significant at a false discovery rate of less than 1% (FDR < 1%). Differentially regulated proteins were detected using the ANOVA test (p < 0.05) and fold change ≥1.1. Abundant proteins were submitted to multivariate statistical analysis—PCA and grouped according to relative abundance profiles. The gene ontology was obtained using pfam2go software (https://rdrr.io/github/missuse/ragp/man/pfam2go.html, accessed on 10 June 2022). Proteomic data is available via ProteomeXchange (PXD031824 and PXD031808).

4. Conclusions

We performed gas exchange and proteomic analysis of cowpea plants inoculated with nematodes and submitted to a water deficit, individually and simultaneously. According to the photosynthesis analysis, plants submitted to a water deficit and combined stress showed reductions in several gas exchange parameters such as E, gs, and An when compared to irrigated plants. Meanwhile, plants only inoculated with nematode presented higher values of An compared to control plants until 10 DAT. We also investigated the proteins and genes expressed in resistant cowpea (CE31) during M. incognita infection and a water deficit, and several differentially abundant candidates were identified. To the best of our knowledge, this is the first study that evaluates the protein profile of resistant cowpea subjected to biotic and abiotic stresses simultaneously, aiming to identify the main processes involved in plant resistance and elucidate the main defense strategies against the nematode and drought tolerance.
In the proteomic analysis of cowpea subjected to isolated and combined stresses, it was possible to observe shared mechanisms of resistance between the stresses, such as the activation of NB-LRR proteins and antioxidant activity. We also identified specific strategies induced by the cross-stress, such as the activation of the signaling mediated by jasmonic acid, the intense production of PR proteins, and a coordinated down-regulation of photosynthetic proteins. The amount and abundance of PRs proteins identified in cowpea subjected to cross-stress may indicate that the imposition of a water deficit can increase the resistance of cowpea to RKN. The regulation of biotic and abiotic responses in plants is complex; however, the differentially abundant proteins here identified (shared and/or specific) as a response to the stresses can help draw a clearer picture of the commonalities and differences in cowpea responses to single and cross stresses. Furthermore, our results suggest that the abiotic stress (drought) had a predominant effect over the biotic stress (nematode). It is noteworthy that the target proteins identified in this study could contribute to the genetic improvement of cowpea in order to develop cultivars resistant to multiple stresses. In addition, farmers could have a lower cost in grain production due to the reductions in irrigation demands and the use of pesticides.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12091900/s1. Differentially abundant proteins from cowpea roots inoculated with nematodes (Table S1), subjected to drought (Table S2) and combined stress (Table S3) when compared to control. Differentially abundant proteins from shoots of cowpea inoculated with nematodes (Table S4), subjected to drought (Table S5) and combined stress (Table S6) when compared to control.

Author Contributions

D.G.R. performed the bioassay, protein extraction, proteomic analysis, and wrote the paper. A.C.M.B. assisted in the drought stress analyses and wrote the paper. I.R.S. performed the proteomic analysis. P.G. assisted in the functional annotation of proteins. W.F. and M.V.d.S. assisted in the acquisition of proteomic data from the mass spectrometer. M.d.S.C. assisted in the analysis of quantification and identification of proteins. M.E.L.-d.-S. assisted in the experimental design and inoculation of nematodes. M.F.G.-d.-S. and O.L.F. supervised the work and corrected the manuscript. A.M. conceived the study and was in charge of overall direction and planning. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by INCT PlantStress Biotech, Conselho Nacional de Desenvolvimento Científico e Tecnológico—CNPq, Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—CAPES (Financial Code 001), Empresa Brasileira de Pesquisa Agropecuária—EMBRAPA, Fundação de Apoio à Pesquisa do Distrito Federal—FAPDF, and Fundação de Apoio a Pesquisa do Estado de Mato Grosso do Sul (FUNDECT).

Data Availability Statement

Proteomic data is available via ProteomeXchange (PXD031824 and PXD031808).

Conflicts of Interest

The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.

References

  1. Raza, A.; Razzaq, A.; Mehmood, S.S.; Zou, X.; Zhang, X.; Lv, Y.; Xu, J. Impact of climate change on crops adaptation and strategies to tackle its outcome: A review. Plants 2019, 8, 34. [Google Scholar] [CrossRef] [PubMed]
  2. Boukar, O.; Massawe, F.; Muranaka, S.; Franco, J.; Maziya-Dixon, B.; Singh, B.; Fatokun, C. Evaluation of cowpea germplasm lines for protein and mineral concentrations in grains. Plant Genet. Resour. 2011, 9, 515–522. [Google Scholar] [CrossRef]
  3. Moens, M.; Perry, R.N.; Starr, J.L. Meloidogyne species—A diverse group of novel and important plant parasites. Root-Knot Nematodes 2009, 1, 483. [Google Scholar]
  4. Shao, H.-B.; Chu, L.-Y.; Jaleel, C.A.; Manivannan, P.; Panneerselvam, R.; Shao, M.-A. Understanding water deficit stress-induced changes in the basic metabolism of higher plants–biotechnologically and sustainably improving agriculture and the ecoenvironment in arid regions of the globe. Crit. Rev. Biotechnol. 2009, 29, 131–151. [Google Scholar] [CrossRef] [PubMed]
  5. Kissoudis, C.; Sunarti, S.; Van De Wiel, C.; Visser, R.G.F.; van der Linden, C.G.; Bai, Y. Responses to combined abiotic and biotic stress in tomato are governed by stress intensity and resistance mechanism. J. Exp. Bot. 2016, 67, 5119–5132. [Google Scholar] [CrossRef] [PubMed]
  6. Ramegowda, V.; Senthil-Kumar, M. The interactive effects of simultaneous biotic and abiotic stresses on plants: Mechanistic understanding from drought and pathogen combination. J. Plant Physiol. 2015, 176, 47–54. [Google Scholar] [CrossRef]
  7. Zhang, H.; Sonnewald, U. Differences and commonalities of plant responses to single and combined stresses. Plant J. 2017, 90, 839–855. [Google Scholar] [CrossRef] [PubMed]
  8. Ku, Y.-S.; Sintaha, M.; Cheung, M.-Y.; Lam, H.-M. Plant hormone signaling crosstalks between biotic and abiotic stress responses. Int. J. Mol. Sci. 2018, 19, 3206. [Google Scholar] [CrossRef]
  9. Bahuguna, R.N.; Gupta, P.; Bagri, J.; Singh, D.; Dewi, A.K.; Tao, L.; Islam, M.; Sarsu, F.; Singla-Pareek, S.L.; Pareek, A. Forward and reverse genetics approaches for combined stress tolerance in rice. Indian J. Plant Physiol. 2018, 23, 630–646. [Google Scholar] [CrossRef]
  10. Bai, Y.; Kissoudis, C.; Yan, Z.; Visser, R.G.F.; van der Linden, G. Plant behaviour under combined stress: Tomato responses to combined salinity and pathogen stress. Plant J. 2018, 93, 781–793. [Google Scholar] [CrossRef]
  11. Muthusamy, S.K.; Lenka, S.K.; Katiyar, A.; Chinnusamy, V.; Singh, A.K.; Bansal, K.C. Genome-wide identification and analysis of biotic and abiotic stress regulation of C4 photosynthetic pathway genes in rice. Appl. Biochem. Biotechnol. 2019, 187, 221–238. [Google Scholar] [CrossRef]
  12. Silva, R.G.G.; Vasconcelos, I.M.; Martins, T.F.; Varela, A.L.N.; Souza, P.F.N.; Lobo, A.K.M.; Silva, F.D.A.; Silveira, J.A.G.; Oliveira, J.T.A. Drought increases cowpea (Vigna unguiculata [L.] Walp.) susceptibility to cowpea severe mosaic virus (CPSMV) at early stage of infection. Plant Physiol. Biochem. 2016, 109, 91–102. [Google Scholar] [CrossRef]
  13. Oliveira, J.T.A.; Andrade, N.C.; Martins-Miranda, A.S.; Soares, A.A.; Gondim, D.M.F.; Araújo-Filho, J.H.; Freire-Filho, F.R.; Vasconcelos, I.M. Differential expression of antioxidant enzymes and PR-proteins in compatible and incompatible interactions of cowpea (Vigna unguiculata) and the root-knot nematode Meloidogyne incognita. Plant Physiol. Biochem. 2012, 51, 145–152. [Google Scholar] [CrossRef] [PubMed]
  14. Huynh, B.-L.; Matthews, W.C.; Ehlers, J.D.; Lucas, M.R.; Santos, J.R.P.; Ndeve, A.; Close, T.J.; Roberts, P.A. A major QTL corresponding to the Rk locus for resistance to root-knot nematodes in cowpea (Vigna unguiculata L. Walp.). Theor. Appl. Genet. 2016, 129, 87–95. [Google Scholar] [CrossRef]
  15. Roberts, P.A.; Matthews, W.C.; Ehlers, J.D. New resistance to virulent root-knot nematodes linked to the Rk locus of cowpea. Crop Sci. 1996, 36, 889–894. [Google Scholar] [CrossRef]
  16. D’Arcy-Lameta, A.; Ferrari-Iliou, R.; Contour-Ansel, D.; Pham-Thi, A.-T.; Zuily-Fodil, Y. Isolation and characterization of four ascorbate peroxidase cDNAs responsive to water deficit in cowpea leaves. Ann. Bot. 2006, 97, 133–140. [Google Scholar] [CrossRef]
  17. Contour-Ansel, D.; Torres-Franklin, M.L.; Cruz De Carvalho, M.H.; D’Arcy-Lameta, A.; Zuily-Fodil, Y. Glutathione reductase in leaves of cowpea: Cloning of two cDNAs, expression and enzymatic activity under progressive drought stress, desiccation and abscisic acid treatment. Ann. Bot. 2006, 98, 1279–1287. [Google Scholar] [CrossRef]
  18. El-Maarouf, H.; d’Arcy-Lameta, A.; Gareil, M.; Zuily-Fodil, Y.; Pham-Thi, A.-T. Cloning and expression under drought of cDNAscoding for two PI-PLCs in cowpea leaves. Plant Physiol. Biochem. 2001, 39, 167–172. [Google Scholar] [CrossRef]
  19. Muchero, W.; Ehlers, J.D.; Close, T.J.; Roberts, P.A. Mapping QTL for drought stress-induced premature senescence and maturity in cowpea [Vigna unguiculata (L.) Walp.]. Theor. Appl. Genet. 2009, 118, 849–863. [Google Scholar] [CrossRef] [PubMed]
  20. Villeth, G.R.C.; Carmo, L.S.T.; Silva, L.P.; Fontes, W.; Grynberg, P.; Saraiva, M.; Brasileiro, A.; Carneiro, R.; Oliveira, J.T.A.; Grossi-de-Sá, M.F. Cowpea–Meloidogyne incognita interaction: Root proteomic analysis during early stages of nematode infection. Proteomics 2015, 15, 1746–1759. [Google Scholar] [CrossRef]
  21. Ribeiro, D.G.; Mota, A.P.Z.; Santos, I.R.; Arraes, F.B.M.; Grynberg, P.; Fontes, W.; de Souza Castro, M.; de Sousa, M.V.; Lisei-de-Sá, M.E.; Grossi-de-Sá, M.F. NBS-LRR-WRKY genes and protease inhibitors (PIs) seem essential for cowpea resistance to root-knot nematode. J. Proteom. 2022, 261, 104575. [Google Scholar] [CrossRef]
  22. Carvalho, M.; Castro, I.; Moutinho-Pereira, J.; Correia, C.; Egea-Cortines, M.; Matos, M.; Rosa, E.; Carnide, V.; Lino-Neto, T. Evaluating stress responses in cowpea under drought stress. J. Plant Physiol. 2019, 241, 153001. [Google Scholar] [CrossRef] [PubMed]
  23. Mofokeng, M.A.; Mashingaidze, K. Breeding and genetic management of drought in cowpea: Progress and technologies. Aust. J. Crop Sci. 2019, 13, 1920–1926. [Google Scholar] [CrossRef]
  24. Thuy, T.T.T.; Yen, N.T.; Tuyet, N.T.A.; Te, L.L.; De Waele, D. Plant-parasitic nematodes and yellowing of leaves associated with black pepper plants in Vietnam. Arch. Phytopathol. Plant Prot. 2012, 45, 1183–1200. [Google Scholar] [CrossRef]
  25. Yigezu Wendimu, G. Biology, Taxonomy, and Management of the Root-Knot Nematode (Meloidogyne incognita) in Sweet Potato. Adv. Agric. 2021, 2021, 8820211. [Google Scholar] [CrossRef]
  26. Else, M.A.; Janowiak, F.; Atkinson, C.J.; Jackson, M.B. Root signals and stomatal closure in relation to photosynthesis, chlorophyll a fluorescence and adventitious rooting of flooded tomato plants. Ann. Bot. 2009, 103, 313–323. [Google Scholar] [CrossRef]
  27. Singh, S.K.; Raja Reddy, K. Regulation of photosynthesis, fluorescence, stomatal conductance and water-use efficiency of cowpea (Vigna unguiculata [L.] Walp.) under drought. J. Photochem. Photobiol. B 2011, 105, 40–50. [Google Scholar] [CrossRef]
  28. Mayek-PÉrez, N.; GarcÍa-Espinosa, R.; LÓpez-CastaÑeda, C.; Acosta-Gallegos, J.A.; Simpson, J. Water relations, histopathology and growth of common bean (Phaseolus vulgaris L.) during pathogenesis of Macrophomina phaseolina under drought stress. Physiol. Mol. Plant Pathol. 2002, 60, 185–195. [Google Scholar] [CrossRef]
  29. Xu, P.; Chen, F.; Mannas, J.P.; Feldman, T.; Sumner, L.W.; Roossinck, M.J. Virus infection improves drought tolerance. New Phytol. 2008, 180, 911–921. [Google Scholar] [CrossRef]
  30. Hosseini, S.A.; Reza, Z.G.; Mostafa, Y.Z.; Mehdi, K. Effects of Cucumber Mosaic Virus infection and drought tolerance of tomato plants under greenhouse conditions: Preliminary results. J. Berry Res. 2018, 8, 129–136. [Google Scholar] [CrossRef]
  31. Vieira, E.A.; das Graças Silva, M.; Moro, C.F.; Laura, V.A. Physiological and biochemical changes attenuate the effects of drought on the Cerrado species Vatairea macrocarpa (Benth.) Ducke. Plant Physiol. Biochem. 2017, 115, 472–483. [Google Scholar] [CrossRef] [PubMed]
  32. Dubey, N.; Singh, K. Role of NBS-LRR proteins in plant defense. In Molecular Aspects of Plant-Pathogen Interaction; Springer: Berlin/Heidelberg, Germany, 2018; pp. 115–138. [Google Scholar]
  33. Elmore, J.M.; Lin, Z.-J.D.; Coaker, G. Plant NB-LRR signaling: Upstreams and downstreams. Curr. Opin. Plant Biol. 2011, 14, 365–371. [Google Scholar] [CrossRef] [PubMed]
  34. Akbudak, M.A.; Yildiz, S.; Filiz, E. Pathogenesis related protein-1 (PR-1) genes in tomato (Solanum lycopersicum L.): Bioinformatics analyses and expression profiles in response to drought stress. Genomics 2020, 112, 4089–4099. [Google Scholar] [CrossRef]
  35. Sahoo, D.K.; Hegde, C.; Bhattacharyya, M.K. Cold Tolerance is Governed by Diverse Genetic Mechanisms Including Those Regulated by NB-LRR-type Receptor Proteins in Arabidopsis. bioRxiv 2022. [Google Scholar] [CrossRef]
  36. Sinha, R.K.; Verma, S.S.; Rastogi, A. Role of Pathogen-Related Protein 10 (PR 10) under abiotic and biotic stresses in plants. Phyton (B Aires) 2020, 89, 167. [Google Scholar] [CrossRef]
  37. Kissoudis, C.; Seifi, A.; Yan, Z.; Islam, A.T.M.; Van der Schoot, H.; van de Wiel, C.; Visser, R.G.F.; Van Der Linden, C.G.; Bai, Y. Ethylene and abscisic acid signaling pathways differentially influence tomato resistance to combined powdery mildew and salt stress. Front. Plant Sci. 2017, 7, 2009. [Google Scholar] [CrossRef]
  38. Mota, A.P.Z.; Brasileiro, A.C.M.; Vidigal, B.; Oliveira, T.N.; da Cunha Quintana Martins, A.; Saraiva, M.A.d.P.; de Araújo, A.C.G.; Togawa, R.C.; Grossi-de-Sá, M.F.; Guimaraes, P.M. Defining the combined stress response in wild Arachis. Sci. Rep. 2021, 11, 11097. [Google Scholar] [CrossRef]
  39. Agbicodo, E.M.; Fatokun, C.A.; Muranaka, S.; Visser, R.G.F. Breeding drought tolerant cowpea: Constraints, accomplishments, and future prospects. Euphytica 2009, 167, 353–370. [Google Scholar] [CrossRef]
  40. Diouf, D. Recent advances in cowpea [Vigna unguiculata (L.) Walp.] “omics” research for genetic improvement. Afr. J. Biotechnol. 2011, 10, 2803–2810. [Google Scholar]
  41. Riyazuddin, R.; Nisha, N.; Singh, K.; Verma, R.; Gupta, R. Involvement of dehydrin proteins in mitigating the negative effects of drought stress in plants. Plant Cell Rep. 2021, 41, 519–533. [Google Scholar] [CrossRef]
  42. Tiwari, P.; Indoliya, Y.; Singh, P.K.; Singh, P.C.; Chauhan, P.S.; Pande, V.; Chakrabarty, D. Role of dehydrin-FK506-binding protein complex in enhancing drought tolerance through the ABA-mediated signaling pathway. Environ. Exp. Bot. 2019, 158, 136–149. [Google Scholar] [CrossRef]
  43. Gill, S.S.; Tuteja, N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 2010, 48, 909–930. [Google Scholar] [CrossRef] [PubMed]
  44. Choudhury, F.K.; Rivero, R.M.; Blumwald, E.; Mittler, R. Reactive oxygen species, abiotic stress and stress combination. Plant J. 2017, 90, 856–867. [Google Scholar] [CrossRef] [PubMed]
  45. Pandey, P.; Irulappan, V.; Bagavathiannan, M.V.; Senthil-Kumar, M. Impact of combined abiotic and biotic stresses on plant growth and avenues for crop improvement by exploiting physio-morphological traits. Front. Plant Sci. 2017, 8, 537. [Google Scholar] [CrossRef]
  46. Anjum, N.A.; Panda, B.B.; Pareek, A.; Talukdar, D.; Pandey, P.; Singh, J.; Mohan, V.; Achary, M.; Reddy, M.K. Redox homeostasis via gene families of ascorbate-glutathione pathway. Front. Environ. Sci. 2015, 3, 25. [Google Scholar]
  47. Laxa, M.; Liebthal, M.; Telman, W.; Chibani, K.; Dietz, K.J. The role of the plant antioxidant system in drought tolerance. Antioxidants 2019, 8, 94. [Google Scholar] [CrossRef] [PubMed]
  48. Easwar Rao, D.; Viswanatha Chaitanya, K. Changes in the antioxidant intensities of seven different soybean (Glycine max (L.) Merr.) cultivars during drought. J. Food Biochem. 2020, 44, e13118. [Google Scholar] [CrossRef]
  49. Zhang, L.; Peng, J.; Chen, T.T.; Zhao, X.H.; Zhang, S.P.; Liu, S.D.; Dong, H.L.; Feng, L.; Yu, S.X. Effect of drought stress on lipid peroxidation and proline content in cotton roots. JAPS J. Anim. Plant Sci. 2014, 24, 1729–1736. [Google Scholar]
  50. Demirevska, K.; Simova-Stoilova, L.; Fedina, I.; Georgieva, K.; Kunert, K. Response of oryzacystatin I transformed tobacco plants to drought, heat and light stress. J. Agron. Crop Sci. 2010, 196, 90–99. [Google Scholar] [CrossRef]
  51. Sekmen, A.H.; Ozgur, R.; Uzilday, B.; Turkan, I. Reactive oxygen species scavenging capacities of cotton (Gossypium hirsutum) cultivars under combined drought and heat induced oxidative stress. Environ. Exp. Bot. 2014, 99, 141–149. [Google Scholar] [CrossRef]
  52. Kissoudis, C.; van de Wiel, C.; Visser, R.G.F.; van der Linden, G. Enhancing crop resilience to combined abiotic and biotic stress through the dissection of physiological and molecular crosstalk. Front. Plant Sci. 2014, 5, 207. [Google Scholar] [CrossRef] [PubMed]
  53. Nostar, O.; Ozdemir, F.; Bor, M.; Turkan, I.; Tosun, N. Combined effects of salt stress and cucurbit downy mildew (Pseudoperospora cubensis Berk. and Curt. Rostov.) infection on growth, physiological traits and antioxidant activity in cucumber (Cucumis sativus L.) seedlings. Physiol. Mol. Plant Pathol. 2013, 83, 84–92. [Google Scholar] [CrossRef]
  54. Ali, S.; Hayat, K.; Iqbal, A.; Xie, L. Implications of abscisic acid in the drought stress tolerance of plants. Agronomy 2020, 10, 1323. [Google Scholar] [CrossRef]
  55. Dong, T.; Park, Y.; Hwang, I. Abscisic acid: Biosynthesis, inactivation, homoeostasis and signalling. Essays Biochem. 2015, 58, 29–48. [Google Scholar]
  56. Ruan, J.; Zhou, Y.; Zhou, M.; Yan, J.; Khurshid, M.; Weng, W.; Cheng, J.; Zhang, K. Jasmonic Acid Signaling Pathway in Plants. Int. J. Mol. Sci. 2019, 20, 2479. [Google Scholar] [CrossRef]
  57. Danquah, A.; de Zelicourt, A.; Colcombet, J.; Hirt, H. The role of ABA and MAPK signaling pathways in plant abiotic stress responses. Biotechnol. Adv. 2014, 32, 40–52. [Google Scholar] [CrossRef]
  58. de Ollas, C.; Hernando, B.; Arbona, V.; Gómez-Cadenas, A. Jasmonic acid transient accumulation is needed for abscisic acid increase in citrus roots under drought stress conditions. Physiol. Plantarum. 2013, 147, 296–306. [Google Scholar] [CrossRef]
  59. Balbi, V.; Devoto, A. Review Summary Tansley review Jasmonate signalling network in Arabidopsis thaliana: Crucial regulatory nodes and new physiological scenarios. New Phytol. 2008, 177, 301–318. [Google Scholar] [CrossRef]
  60. Davila Olivas, N.H.; Kruijer, W.; Gort, G.; Wijnen, C.L.; van Loon, J.J.A.; Dicke, M. Genome-wide association analysis reveals distinct genetic architectures for single and combined stress responses in Arabidopsis thaliana. New Phytol. 2017, 213, 838–851. [Google Scholar] [CrossRef]
  61. Kiba, A.; Nakano, M.; Hosokawa, M.; Galis, I.; Nakatani, H.; Shinya, T.; Ohnishi, K.; Hikichi, Y. Phosphatidylinositol-phospholipase C2 regulates pattern-triggered immunity in Nicotiana benthamiana. J. Exp. Bot. 2020, 71, 5027–5038. [Google Scholar] [CrossRef]
  62. Rodas-Junco, B.A.; Nic-Can, G.I.; Muñoz-Sánchez, A.; Teresa Hernández-Sotomayor, S.M. Phospholipid signaling is a component of the salicylic acid response in plant cell suspension cultures. Int. J. Mol. Sci. 2020, 21, 5285. [Google Scholar] [CrossRef] [PubMed]
  63. Singh, S.; Singh, B.; Singh, A.P. Nematodes: A threat to sustainability of agriculture. Procedia Environ. Sci. 2015, 29, 215–216. [Google Scholar] [CrossRef]
  64. González-Mendoza, V.M.; Sánchez-Sandoval, M.E.; Castro-Concha, L.A.; Teresa Hernández-Sotomayor, S.M. Phospholipases c and d and their role in biotic and abiotic stresses. Plants 2021, 10, 921. [Google Scholar] [CrossRef]
  65. Johansson, O.N.; Fahlberg, P.; Karimi, E.; Nilsson, A.K.; Ellerström, M.; Andersson, M.X. Redundancy among phospholipase D isoforms in resistance triggered by recognition of the Pseudomonas syringae effector AvrRpm1 in Arabidopsis thaliana. Front. Plant Sci. 2014, 5, 639. [Google Scholar] [CrossRef]
  66. Zhao, J.; Devaiah, S.P.; Wang, C.; Li, M.; Welti, R.; Wang, X. Arabidopsis phospholipase Dβ1 modulates defense responses to bacterial and fungal pathogens. New Phytol. 2013, 199, 228–240. [Google Scholar] [CrossRef]
  67. Zhang, Q.; Van Wijk, R.; Shahbaz, M.; Roels, W.; Van Schooten, B.; Vermeer, J.E.M.; Zarza, X.; Guardia, A.; Scuffi, D.; García-Mata, C.; et al. Arabidopsis phospholipase C3 is involved in lateral root initiation and ABA responses in seed germination and stomatal closure. Plant Cell Physiol. 2018, 59, 469–486. [Google Scholar] [CrossRef] [PubMed]
  68. Georges, F.; Das, S.; Ray, H.; Bock, C.; Nokhrina, K.; Kolla, V.A.; Keller, W. Over-expression of Brassica napus phosphatidylinositol-phospholipase C2 in canola induces significant changes in gene expression and phytohormone distribution patterns, enhances drought tolerance and promotes early flowering and maturation. Plant Cell Environ. 2009, 32, 1664–1681. [Google Scholar] [CrossRef]
  69. Liu, Q.; Zhou, Y.; Li, H.; Liu, R.; Wang, W.; Wu, W.; Yang, N.; Wang, S. Osmotic stress-triggered stomatal closure requires Phospholipase Dδ and hydrogen sulfide in Arabidopsis thaliana. Biochem. Biophys. Res. Commun. 2021, 534, 914–920. [Google Scholar] [CrossRef] [PubMed]
  70. Faizi, S.; Fayyaz, S.; Bano, S.; Yawar Iqbal, E.; Siddiqi, H.; Naz, A. Isolation of nematicidal compounds from Tagetes patula L. yellow flowers: Structure–activity relationship studies against cyst nematode Heterodera zeae infective stage larvae. J. Agric. Food Chem. 2011, 59, 9080–9093. [Google Scholar] [CrossRef]
  71. Wuyts, N.; Swennen, R.; De Waele, D. Effects of plant phenylpropanoid pathway products and selected terpenoids and alkaloids on the behaviour of the plant-parasitic nematodes Radopholus similis, Pratylenchus penetrans and Meloidogyne incognita. Nematology 2006, 8, 89–101. [Google Scholar] [CrossRef]
  72. Ma, D.; Sun, D.; Wang, C.; Li, Y.; Guo, T. Expression of flavonoid biosynthesis genes and accumulation of flavonoid in wheat leaves in response to drought stress. Plant Physiol. Biochem. 2014, 80, 60–66. [Google Scholar] [CrossRef]
  73. Yuan, Y.; Liu, Y.; Wu, C.; Chen, S.; Wang, Z.; Yang, Z.; Qin, S.; Huang, L. Water deficit affected flavonoid accumulation by regulating hormone metabolism in Scutellaria baicalensis Georgi roots. PLoS ONE 2012, 7, e42946. [Google Scholar] [CrossRef]
  74. Ibrahim, W.; Zhu, Y.; Chen, Y.; Qiu, C.; Zhu, S.; Wu, F. Genotypic differences in leaf secondary metabolism, plant hormones and yield under alone and combined stress of drought and salinity in cotton genotypes. Physiol. Plant 2019, 165, 343–355. [Google Scholar] [CrossRef] [PubMed]
  75. Schulz, E.; Tohge, T.; Winkler, J.B.; Albert, A.; Schäffner, A.R.; Fernie, A.R.; Zuther, E.; Hincha, D.K. Natural variation among Arabidopsis accessions in the regulation of flavonoid metabolism and stress gene expression by combined UV radiation and cold. Plant Cell Physiol. 2021, 62, 502–514. [Google Scholar] [CrossRef] [PubMed]
  76. Tarafdar, A.; Rani, T.S.; Chandran, U.S.; Ghosh, R.; Chobe, D.R.; Sharma, M. Exploring combined effect of abiotic (soil moisture) and biotic (Sclerotium rolfsii Sacc.) stress on collar rot development in chickpea. Front. Plant Sci. 2018, 9, 1154. [Google Scholar] [CrossRef]
  77. Morales, F.; Ancín, M.; Fakhet, D.; González-Torralba, J.; Gámez, A.L.; Seminario, A.; Soba, D.; Ben Mariem, S.; Garriga, M.; Aranjuelo, I. Photosynthetic metabolism under stressful growth conditions as a bases for crop breeding and yield improvement. Plants 2020, 9, 88. [Google Scholar] [CrossRef]
  78. Vass, I. Molecular mechanisms of photodamage in the Photosystem II complex. Biochim. Biophys. Acta Bioenerg. 2012, 1817, 209–217. [Google Scholar] [CrossRef] [PubMed]
  79. Chaves, M.M.; Flexas, J.; Pinheiro, C. Photosynthesis under drought and salt stress: Regulation mechanisms from whole plant to cell. Ann. Bot. 2009, 103, 551–560. [Google Scholar] [CrossRef]
  80. Haldrup, A.; Simpson, D.J.; Scheller, H.V. Down-regulation of the PSI-F subunit of photosystem I (PSI) in Arabidopsis thaliana. The PSI-F subunit is essential for photoautotrophic growth and contributes to antenna function. J. Biol. Chem. 2000, 275, 31211–31218. [Google Scholar] [CrossRef] [PubMed]
  81. Schöttler, M.A.; Albus, C.A.; Bock, R. Photosystem I: Its biogenesis and function in higher plants. J. Plant Physiol. 2011, 168, 1452–1461. [Google Scholar] [CrossRef]
  82. Jansson, S. The light-harvesting chlorophyll ab-binding proteins. Biochim. Biophys. Acta (BBA)-Bioenerg. 1994, 1184, 1–19. [Google Scholar] [CrossRef]
  83. Crafts-Brandner, S.J.; Salvucci, M.E. Sensitivity of Photosynthesis in a C4 Plant, Maize, to Heat Stress. Plant Physiol. 2002, 129, 1773–1780. [Google Scholar] [CrossRef] [PubMed]
  84. Bilgin, D.D.; Zavala, J.A.; Zhu, J.; Clough, S.J.; Ort, D.R.; Delucia, E.H. Biotic stress globally downregulates photosynthesis genes. Plant Cell Environ. 2010, 33, 1597–1613. [Google Scholar] [CrossRef] [PubMed]
  85. Osei-Bonsu, I.; McClain, A.M.; Walker, B.J.; Sharkey, T.D.; Kramer, D.M. The roles of photorespiration and alternative electron acceptors in the responses of photosynthesis to elevated temperatures in cowpea. In Plant Cell and Environment; Blackwell Publishing Ltd.: Hoboken, NJ, USA, 2021; Volume 44, pp. 2290–2307. [Google Scholar]
  86. Taylor, A.L.; Sasser, J.N. Biology, Identification and Control of Root-Knot Nematodes; North Carolina State University Graphics: Raleigh, NC, USA, 1978; 111p. [Google Scholar]
  87. Sasser, J.N.; Carter, C.C.; Hartman, K.M. Standardization of Host Suitability Studies and Reporting of Resistance to Root-Knot Nematodes; North Carolina State University Graphics: Raleigh, NC, USA, 1984. [Google Scholar]
  88. Hoagland, D.R.; Arnon, D.I. The Water-Culture Method for Growing Plants without Soil; Circular 347; California Agricultural Experiment Station: Berkeley, CA, USA, 1950. [Google Scholar]
  89. Hussey, R.S. A comparison of methods of collecting inocula of Meloidogyne spp., including a new technique. Plant Dis. Rep. 1973, 57, 1025–1028. [Google Scholar]
  90. Rappsilber, J.; Mann, M.; Ishihama, Y. Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. Nat. Protoc. 2007, 2, 1896–1906. [Google Scholar] [CrossRef] [PubMed]
  91. Ribeiro, D.G.; de Almeida, R.F.; Fontes, W.; de Souza Castro, M.; de Sousa, M.V.; Ricart, C.A.O.; da Cunha, R.N.V.; Lopes, R.; Scherwinski-Pereira, J.E.; Mehta, A. Stress and cell cycle regulation during somatic embryogenesis plays a key role in oil palm callus development. J. Proteom. 2019, 192, 137–146. [Google Scholar] [CrossRef] [PubMed]
  92. Kalli, A.; Smith, G.T.; Sweredoski, M.J.; Hess, S. Evaluation and Optimization of Mass Spectrometric Settings during Data-dependent Acquisition Mode: Focus on LTQ-Orbitrap Mass Analyzers. J. Proteome Res. 2013, 12, 3071–3086. [Google Scholar] [CrossRef]
  93. Hawkridge, A.M. Practical considerations and current limitations in quantitative mass spectrometry-based proteomics. In Quantitative Proteomics; Eyers, C.E., Gaskell, S., Eds.; Royal Society of Chemistry: Cambridge, UK, 2014; pp. 3–25. [Google Scholar]
Figure 1. Morphological aspects of cowpea plants at 12 days after inoculation with nematodes (1000 J2) and interruption of irrigation (25% CF). Control CE 31 plants (a), inoculated with nematodes (b), submitted to water deficit (c), and inoculated with nematode and submitted to drought (d).
Figure 1. Morphological aspects of cowpea plants at 12 days after inoculation with nematodes (1000 J2) and interruption of irrigation (25% CF). Control CE 31 plants (a), inoculated with nematodes (b), submitted to water deficit (c), and inoculated with nematode and submitted to drought (d).
Plants 12 01900 g001
Figure 2. Intercellular CO2 concentrations, Ci (A), carbon transpiration rate, E (B), stomatal conductance, gs (C), and net photosynthesis rate, An (D) in leaves of V. unguiculata CE 31 inoculated with M. incognita (1000 J2) and/or submitted to water deficit (25% of field capacity) evaluated at 6, 8, 10, and 12 DAT. Each bar represents the mean of three biological replicates with standard deviation. Small caption letters represent comparisons among treatments using Tukey’s test with a 5% confidence level (p ≤ 0.05).
Figure 2. Intercellular CO2 concentrations, Ci (A), carbon transpiration rate, E (B), stomatal conductance, gs (C), and net photosynthesis rate, An (D) in leaves of V. unguiculata CE 31 inoculated with M. incognita (1000 J2) and/or submitted to water deficit (25% of field capacity) evaluated at 6, 8, 10, and 12 DAT. Each bar represents the mean of three biological replicates with standard deviation. Small caption letters represent comparisons among treatments using Tukey’s test with a 5% confidence level (p ≤ 0.05).
Plants 12 01900 g002
Figure 3. Principal component analysis (PCA) and heat maps of differentially abundant proteins in roots (A and C, respectively) and shoots (B and D, respectively) of cowpea.
Figure 3. Principal component analysis (PCA) and heat maps of differentially abundant proteins in roots (A and C, respectively) and shoots (B and D, respectively) of cowpea.
Plants 12 01900 g003
Figure 4. Bubble plot of the biological processes, molecular functions, and cellular components in which differentially abundant proteins of roots are involved. Increased proteins are represented by red bubbles and decreased by blue bubbles.
Figure 4. Bubble plot of the biological processes, molecular functions, and cellular components in which differentially abundant proteins of roots are involved. Increased proteins are represented by red bubbles and decreased by blue bubbles.
Plants 12 01900 g004
Figure 5. Bubble plot of biological processes, molecular functions, and cellular components in which differentially abundant proteins from shoots are involved. Increased proteins are represented by red bubbles and decreased proteins by blue bubbles.
Figure 5. Bubble plot of biological processes, molecular functions, and cellular components in which differentially abundant proteins from shoots are involved. Increased proteins are represented by red bubbles and decreased proteins by blue bubbles.
Plants 12 01900 g005
Table 1. Differentially abundant proteins of cowpea roots inoculated with nematodes and submitted to drought, individually and simultaneously, when compared to the control.
Table 1. Differentially abundant proteins of cowpea roots inoculated with nematodes and submitted to drought, individually and simultaneously, when compared to the control.
Biological ProcessProtein ID (UniProt)DescriptionScoreAnova (p ≤ 0.05)Log2FCCross/ControlDrought/ControlNematode/Control
Defense responseA0A4D6LDE6LRR receptor-like serine/threonine-protein kinase FLS2188.440.001-Exclusive
A0A4D6LDX5LRR receptor-like serine/threonine-protein kinase FLS236.160.023-Exclusive
A0A4D6M0M0Chitinase68.290.0086.29Increased
A0A4D6KNF4Pathogenesis-related protein 1350.079.42 × 10−55.58Increased
Dehydration responseQ9AYM8CPRD2 protein119.930.0013.22Increased
P93700CPRD14 protein1000.730.0021.7 Increased
Q9FS23CPRD86 protein (Fragment)121.70.0256.11 Increased
Oxidative stressA0A4D6LCC8Glutathione S-transferase50.390.0179.43Increased
Q9M7R2Superoxide dismutase510.360.0013.74 Increased
A0A4D6MLJ0Glutathione peroxidase221.880.0032.15 Increased
Jasmonic acid biosynthetic processA0A4D6MNE8Lipoxygenase132.891.36 × 10−63.72Increased
A0A4D6MNU1Lipoxygenase84.890.0033.53Increased
Response to abscisic acidA0A4D6LN47Bet_v_1 domain-containing protein284.530.0193.61Increased
A0A4D6LBT9Bet_v_1 domain-containing protein29.690.0196.73 Increased
Lipid metabolic processA0A4D6N7S0Phospholipase D1345.210.0311.11Increased
O04865Phospholipase D alpha 1424.290.0131.06Increased
A0A4D6MIW6Phosphoinositide phospholipase C51.840.0095.77 Increased
Flavonoid biosynthetic processA0A4D6NPD2Chalcone synthase90.180.0002.31 Increased
A0A4D6NRT4Isoflavone reductase840.960.0060.77 Increased
A0A4D6NPC3Flavonol 3-O-methyltransferase407.450.0023.86Increased
Table 2. Differentially abundant proteins of cowpea shoots inoculated with nematodes and submitted to drought, individually and simultaneously, when compared to the control.
Table 2. Differentially abundant proteins of cowpea shoots inoculated with nematodes and submitted to drought, individually and simultaneously, when compared to the control.
Biological ProcessProtein ID (UniProt)DescriptionScoreAnova (p ≤ 0.05)Log2FCCross/ControlDrought/ControlNematode/Control
Defense responseA0A4D6M0I2LRR receptor-like serine/threonine-protein kinase FLS226.997.07 × 10−512.52Increased
A0A4D6LBY3Chitinase247.880.0052.14Increased
A0A4D6LT51Thaumatin120.727.80 × 10−43.62Increased
Dehydration responseP93700CPRD14 protein566.231.98 × 10−41.64Increased
A0A4D6NTQ2Dehydrin327.310.031.95Increased
Oxidative stressA0A4D6MLG9Glutathione transferase320.252.174 × 10−52.94Increased
A0A4D6LKZ9Glutathione S-transferase41.730.0012.25 Increased
Q41712L-ascorbate peroxidase623.750.021.24Increased
Q5QIA9L-ascorbate peroxidase221.680.0012.01 Increased
A0A4D6N707Superoxide dismutase329.210.0411.20 Increased
Jasmonic acid biosynthetic processA0A4D6MNU1Lipoxygenase86.990.0392.72Increased
A0A4D6MNL5Lipoxygenase406.251.88 × 10−43.81 Increased
Response to abscisic acidA0A4D6LBT9Bet_v_1 domain-containing protein29.090.0492.55 Increased
A0A4D6LB15Annexin1459.331.87 × 10−51.83Increased
Lipid metabolic processA0A4D6N7S0Phospholipase D1182.195.98 × 10−61.68Increased
A0A4D6MCA1Phospholipase A167.940.0172.87 Increased
PhotosynthesisA0A4D6L6M3Chlorophyll a-b binding protein, chloroplastic288.640.0081.32 Increased
A0A4D6LJ91Photosystem II stability/assembly factor730.690.0141.86 Increased
A0A4D6L8F7Photosystem I subunit PsaN82.180.021.80 Decreased
A0A4D6M6X8Magnesium-protoporphyrin IX monomethyl ester (oxidative) cyclase50.371.49 × 10−22.72Decreased
A0A4D6L976Magnesium-protoporphyrin O-methyltransferase252.455.47 × 10−42.26Decreased
A0A4D6N7U7NADPH-protochlorophyllide oxidoreductase281.881.81 × 10−52.45Decreased
A0A4D6KW86Geranylgeranyl reductase58.534.78 × 10−43.34Decreased
A0A4D6MNX0Chlorophyll a-b binding protein, chloroplastic703.791.40 × 10−43.21Decreased
A0A4D6NHQ4Chlorophyll a-b binding protein, chloroplastic95.340.043.32Decreased
A0A4D6N658Photosystem II oxygen-evolving enhancer protein 2118.024.66 × 10−32.26Decreased
A0A4D6NAD0Photosystem II oxygen-evolving enhancer protein 3232.750.031.27Decreased
A0A4D6L0K2Photosystem II oxygen-evolving enhancer protein 154.410.052.47Decreased
A0A4D6M3U4Photosystem I P700 chlorophyll a apoprotein A11060.971.40 × 10−21.51Decreased
I2E2T7Photosystem I iron-sulfur center322.985.73 × 10−41.87Decreased
J7EX90Cytochrome b6806.963.33 × 10−31.54Decreased
I2E2Q5Cytochrome f341.35.15 × 10−51.79Decreased
A0A4D6LQE0Ferredoxin--NADP reductase, chloroplastic2375.032.84 × 10−51.13Decreased
I2E2P9ATP synthase subunit alpha3858.042.73 × 10−30.97Decreased
I2E2M7ATP synthase subunit beta388.272.78 × 10−62.04Decreased
A0A4D6MKZ9Photosystem I subunit VI637.475.69 × 10−31.42Decreased
A6H596Putative rubisco activase (Fragment)38.235.76 × 10−54.47Decreased
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

Ribeiro, D.G.; Bezerra, A.C.M.; Santos, I.R.; Grynberg, P.; Fontes, W.; de Souza Castro, M.; de Sousa, M.V.; Lisei-de-Sá, M.E.; Grossi-de-Sá, M.F.; Franco, O.L.; et al. Proteomic Insights of Cowpea Response to Combined Biotic and Abiotic Stresses. Plants 2023, 12, 1900. https://doi.org/10.3390/plants12091900

AMA Style

Ribeiro DG, Bezerra ACM, Santos IR, Grynberg P, Fontes W, de Souza Castro M, de Sousa MV, Lisei-de-Sá ME, Grossi-de-Sá MF, Franco OL, et al. Proteomic Insights of Cowpea Response to Combined Biotic and Abiotic Stresses. Plants. 2023; 12(9):1900. https://doi.org/10.3390/plants12091900

Chicago/Turabian Style

Ribeiro, Daiane Gonzaga, Ana Carolina Mendes Bezerra, Ivonaldo Reis Santos, Priscila Grynberg, Wagner Fontes, Mariana de Souza Castro, Marcelo Valle de Sousa, Maria Eugênia Lisei-de-Sá, Maria Fatima Grossi-de-Sá, Octávio Luiz Franco, and et al. 2023. "Proteomic Insights of Cowpea Response to Combined Biotic and Abiotic Stresses" Plants 12, no. 9: 1900. https://doi.org/10.3390/plants12091900

APA Style

Ribeiro, D. G., Bezerra, A. C. M., Santos, I. R., Grynberg, P., Fontes, W., de Souza Castro, M., de Sousa, M. V., Lisei-de-Sá, M. E., Grossi-de-Sá, M. F., Franco, O. L., & Mehta, A. (2023). Proteomic Insights of Cowpea Response to Combined Biotic and Abiotic Stresses. Plants, 12(9), 1900. https://doi.org/10.3390/plants12091900

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