Cysteine protease and cystatin expression and activity during soybean nodule development and senescence

Defense and senescence interplay in legume nodules

Immunity and senescence play a crucial role in the functioning of the legume symbiotic nodules. The miss-regulation of one of these processes compromises the symbiosis leading to death of the endosymbiont and the arrest of the nodule functioning. The relationship between immunity and senescence has been extensively studied in plant organs where a synergistic response can be observed. However, the interplay between immunity and senescence in the symbiotic organ is poorly discussed in the literature and these phenomena are often mixed up. Recent studies revealed that the cooperation between immunity and senescence is not always observed in the nodule, suggesting complex interactions between these two processes within the symbiotic organ. Here, we discuss recent results on the interplay between immunity and senescence in the nodule and the specificities of this relationship during legume-rhizobium symbiosis.

Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis

The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

Effects of Elevated Temperature on Pisum sativum Nodule Development: I-Detailed Characteristic of Unusual Apical Senescence

Despite global warming, the influence of heat on symbiotic nodules is scarcely studied. In this study, the effects of heat stress on the functioning of nodules formed by Rhizobium leguminosarum bv. viciae strain 3841 on pea (Pisum sativum) line SGE were analyzed. The influence of elevated temperature was analyzed at histological, ultrastructural, and transcriptional levels. As a result, an unusual apical pattern of nodule senescence was revealed. After five days of exposure, a senescence zone with degraded symbiotic structures was formed in place of the distal nitrogen fixation zone. There was downregulation of various genes, including those associated with the assimilation of fixed nitrogen and leghemoglobin. After nine days, the complete destruction of the nodules was demonstrated. It was shown that nodule recovery was possible after exposure to elevated temperature for 3 days but not after 5 days (which coincides with heat wave duration). At the same time, the exposure of plants to optimal temperature during the night leveled the negative effects. Thus, the study of the effects of elevated temperature on symbiotic nodules using a well-studied pea genotype and Rhizobium strain led to the discovery of a novel positional response of the nodule to heat stress.

Nodule-specific cysteine-rich peptide 343 is required for symbiotic nitrogen fixation in Medicago truncatula

Symbiotic interactions between legumes and rhizobia lead to the development of root nodules and nitrogen fixation by differentiated bacteroids within nodules. Differentiation of the endosymbionts is reversible or terminal, determined by plant effectors. In inverted repeat lacking clade legumes, nodule-specific cysteine-rich (NCR) peptides control the terminal differentiation of bacteroids. Medicago truncatula contains ∼700 NCR-coding genes. However, the role of few NCR peptides has been demonstrated. Here, we report characterization of fast neutron 2106 (FN2106), a symbiotic nitrogen fixation defective (fix-) mutant of M. truncatula. Using a transcript-based approach, together with linkage and complementation tests, we showed that loss-of-function of NCR343 results in impaired bacteroid differentiation and/or maintenance and premature nodule senescence of the FN2106 mutant. NCR343 was specifically expressed in nodules. Subcellular localization studies showed that the functional NCR343-YFP fusion protein colocalizes with bacteroids in symbiosomes in infected nodule cells. Transcriptomic analyses identified senescence-, but not defense-related genes, as being significantly upregulated in ncr343 (FN2106) nodules. Taken together, results from our phenotypic and transcriptomic analyses of a loss-of-function ncr343 mutant demonstrate an essential role of NCR343 in bacteroid differentiation and/or maintenance required for symbiotic nitrogen fixation.

GmNAC039 and GmNAC018 activate the expression of cysteine protease genes to promote soybean nodule senescence

Root nodules are major sources of nitrogen for soybean (Glycine max (L.) Merr.) growth, development, production, and seed quality. Symbiotic nitrogen fixation is time-limited, as the root nodule senesces during the reproductive stage of plant development, specifically during seed development. Nodule senescence is characterized by the induction of senescence-related genes, such as papain-like cysteine proteases (CYPs), which ultimately leads to the degradation of both bacteroids and plant cells. However, how nodule senescence-related genes are activated in soybean is unknown. Here, we identified 2 paralogous NAC transcription factors, GmNAC039 and GmNAC018, as master regulators of nodule senescence. Overexpression of either gene induced soybean nodule senescence with increased cell death as detected using a TUNEL assay, whereas their knockout delayed senescence and increased nitrogenase activity. Transcriptome analysis and nCUT&Tag-qPCR assays revealed that GmNAC039 directly binds to the core motif CAC(A)A and activates the expression of 4 GmCYP genes (GmCYP35, GmCYP37, GmCYP39, and GmCYP45). Similar to GmNAC039 and GmNAC018, overexpression or knockout of GmCYP genes in nodules resulted in precocious or delayed senescence, respectively. These data provide essential insights into the regulatory mechanisms of nodule senescence, in which GmNAC039 and GmNAC018 directly activate the expression of GmCYP genes to promote nodule senescence.

CaCP15 Gene Negatively Regulates Salt and Osmotic Stress Responses in Capsicum annuum L

Salt and osmotic stress seriously restrict the growth, development, and productivity of horticultural crops in the greenhouse. The papain-like cysteine proteases (PLCPs) participate in multi-stress responses in plants. We previously demonstrated that salt and osmotic stress affect cysteine protease 15 of pepper (Capsicum annuum L.) (CaCP15); however, the role of CaCP15 in salt and osmotic stress responses is unknown. Here, the function of CaCP15 in regulating pepper salt and osmotic stress resistance was explored. Pepper plants were subjected to abiotic (sodium chloride, mannitol, salicylic acid, ethrel, methyl jasmonate, etc.) and biotic stress (Phytophthora capsici inoculation). The CaCP15 was silenced through the virus-induced gene silencing (VIGS) and transiently overexpressed in pepper plants. The full-length CaCP15 fragment is 1568 bp, with an open reading frame of 1032 bp, encoding a 343 amino acid protein. CaCP15 is a senescence-associated gene 12 (SAG12) subfamily member containing two highly conserved domains, Inhibitor 129 and Peptidase_C1. CaCP15 expression was the highest in the stems of pepper plants. The expression was induced by salicylic acid, ethrel, methyl jasmonate, and was infected by Phytophthora capsici inoculation. Furthermore, CaCP15 was upregulated under salt and osmotic stress, and CaCP15 silencing in pepper enhanced salt and mannitol stress resistance. Conversely, transient overexpression of CaCP15 increased the sensitivity to salt and osmotic stress by reducing the antioxidant enzyme activities and negatively regulating the stress-related genes. This study indicates that CaCP15 negatively regulates salt and osmotic stress resistance in pepper via the ROS-scavenging.

Soybean gene co-expression network analysis identifies two co-regulated gene modules associated with nodule formation and development

Gene co-expression network analysis is an efficient systems biology approach for the discovery of novel gene functions and trait-associated gene modules. To identify clusters of functionally related genes involved in soybean nodule formation and development, we performed a weighted gene co-expression network analysis. Two nodule-specific modules (NSM-1 and NSM-2, containing 304 and 203 genes, respectively) were identified. The NSM-1 gene promoters were significantly enriched in cis-binding elements for ERF, MYB, and C2H2-type zinc transcription factors, whereas NSM-2 gene promoters were enriched in cis-binding elements for TCP, bZIP, and bHLH transcription factors, suggesting a role of these regulatory factors in the transcriptional activation of nodule co-expressed genes. The co-expressed gene modules included genes with potential novel roles in nodulation, including those involved in xylem development, transmembrane transport, the ethylene signalling pathway, cytoskeleton organization, cytokinesis and regulation of the cell cycle, regulation of meristem initiation and growth, transcriptional regulation, DNA methylation, and histone modifications. Functional analysis of two co-expressed genes using TILLING mutants provided novel insight into the involvement of unsaturated fatty acid biosynthesis and folate metabolism in nodule formation and development. The identified gene co-expression modules provide valuable resources for further functional genomics studies to dissect the genetic basis of nodule formation and development in soybean.

Critical role for uricase and xanthine dehydrogenase in soybean nitrogen fixation and nodule development

De novo purine biosynthesis is required for the incorporation of fixed nitrogen in ureide exporting nodules, as formed on soybean [Glycine max (L.) Merr.] roots. However, in many cases, the enzymes involved in this pathway have been deduced strictly from genome annotations with little direct genetic evidence, such as mutant studies, to confirm their biochemical function or importance to nodule development. While efforts to develop large mutant collections of soybean are underway, research on this plant is still hampered by the inability to obtain mutations in any specific gene of interest. Using a forward genetic approach, as well as CRISPR/Cas9 gene editing via Agrobacterium rhizogenes-mediated hairy root transformation, we identified and characterized the role of GmUOX (Uricase) and GmXDH (Xanthine Dehydrogenase) in nitrogen fixation and nodule development in soybean. The gmuox knockout soybean mutants displayed nitrogen deficiency chlorosis and early nodule senescence, as exemplified by the reduced nitrogenase (acetylene reduction) activity in nodules, the internal greenish-white internal appearance of nodules, and diminished leghemoglobin production. In addition, gmuox1 nodules showed collapsed infected cells with degraded cytoplasm, aggregated bacteroids with no discernable symbiosome membranes, and increased formation of poly-β-hydroxybutyrate granules. Similarly, knockout gmxdh mutant nodules, generated with the CRISPR/Cas9 system, also exhibited early nodule senescence. These genetic studies confirm the critical role of the de novo purine metabolisms pathway not only in the incorporation of fixed nitrogen but also in the successful development of a functional, nitrogen-fixing nodule. Furthermore, these studies demonstrate the great utility of the CRISPR/Cas9 system for studying root-associated gene traits when coupled with hairy root transformation.

Insight into the control of nodule immunity and senescence during Medicago truncatula symbiosis

Medicago (Medicago truncatula) establishes a symbiosis with the rhizobia Sinorhizobium sp, resulting in the formation of nodules where the bacteria fix atmospheric nitrogen. The loss of immunity repression or early senescence activation compromises symbiont survival and leads to the formation of nonfunctional nodules (fix-). Despite many studies exploring an overlap between immunity and senescence responses outside the nodule context, the relationship between these processes in the nodule remains poorly understood. To investigate this phenomenon, we selected and characterized three Medicago mutants developing fix- nodules and showing senescence responses. Analysis of specific defense (PATHOGENESIS-RELATED PROTEIN) or senescence (CYSTEINE PROTEASE) marker expression demonstrated that senescence and immunity seem to be antagonistic in fix- nodules. The growth of senescence mutants on non-sterile (sand/perlite) substrate instead of sterile in vitro conditions decreased nodule senescence and enhanced defense, indicating that environment can affect the immunity/senescence balance. The application of wounding stress on wild-type (WT) fix+ nodules led to the death of intracellular rhizobia and associated with co-stimulation of defense and senescence markers, indicating that in fix+ nodules the relationship between the two processes switches from opposite to synergistic to control symbiont survival during response to the stress. Our data show that the immune response in stressed WT nodules is linked to the repression of DEFECTIVE IN NITROGEN FIXATION 2 (DNF2), Symbiotic CYSTEINE-RICH RECEPTOR-LIKE KINASE (SymCRK), and REGULATOR OF SYMBIOSOME DIFFERENTIATION (RSD), key genes involved in symbiotic immunity suppression. This study provides insight to understand the links between senescence and immunity in Medicago nodules.

Genome-wide association identifies candidate genes for drought tolerance in coast redwood and giant sequoia

Drought is a major limitation for survival and growth in plants. With more frequent and severe drought episodes occurring due to climate change, it is imperative to understand the genomic and physiological basis of drought tolerance to be able to predict how species will respond in the future. In this study, univariate and multitrait multivariate genome-wide association study methods were used to identify candidate genes in two iconic and ecosystem-dominating species of the western USA, coast redwood and giant sequoia, using 10 drought-related physiological and anatomical traits and genome-wide sequence-capture single nucleotide polymorphisms. Population-level phenotypic variation was found in carbon isotope discrimination, osmotic pressure at full turgor, xylem hydraulic diameter, and total area of transporting fibers in both species. Our study identified new 78 new marker × trait associations in coast redwood and six in giant sequoia, with genes involved in a range of metabolic, stress, and signaling pathways, among other functions. This study contributes to a better understanding of the genomic basis of drought tolerance in long-generation conifers and helps guide current and future conservation efforts in the species.

Plant protease as regulator and signaling molecule for enhancing environmental stress-tolerance

Proteases are ubiquitous in prokaryotes and eukaryotes. Plant proteases are key regulators of various physiological processes, including protein homeostasis, organelle development, senescence, seed germination, protein processing, environmental stress response, and programmed cell death. Proteases are involved in the breakdown of peptide bonds resulting in irreversible posttranslational modification of the protein. Proteases act as signaling molecules that specifically regulate cellular function by cleaving and triggering receptor molecules. Peptides derived from proteolysis regulate ROS signaling under oxidative stress in the plant. It degrades misfolded and abnormal proteins into amino acids to repair the cell damage and regulates the biological process in response to environmental stress. Proteases modulate the biogenesis of phytohormones which control plant growth, development, and environmental stresses. Protein homeostasis, the overall balance between protein synthesis and proteolysis, is required for plant growth and development. Abiotic and biotic stresses are major factors that negatively impact cellular survivability, biomass production, and reduced crop yield potentials. Therefore, the identification of various stress-responsive proteases and their molecular functions may elucidate valuable information for the development of stress-resilient crops with higher yield potentials. However, the understanding of molecular mechanisms of plant protease remains unexplored. This review provides an overview of proteases related to development, signaling, and growth regulation to acclimatize environmental stress in plants.

"The PLCP gene family of grapevine (Vitis vinifera L.): characterization and differential expression in response to Plasmopara Viticola"

BACKGROUND

Papain-like cysteine proteases (PLCPs), a large group of cysteine proteases, are structurally related to papain. The members belonging to PLCPs family contribute to plant immunity, senescence, and defense responses in plants. The PLCP gene family has been identified in Arabidopsis, rice, soybean, and cotton. However, no systematic analysis of PLCP genes has been undertaken in grapevine. Since Plasmopara viticola as a destructive pathogen could affect immunity of grapes in the field, we considered that the members belonged to PLCPs family could play a crucial role in defensive mechanisms or programmed cell death. We aimed to evaluate the role of PLCPs in 2 different varieties of grapevines and compared the changes of their expressions with the transcriptional data in response to P. viticola.

RESULTS

In this study, 23 grapevine PLCP (VvPLCP) genes were identified by comprehensive bioinformatics analysis. Subsequently, the chromosomal localizations, gene structure, conserved domains, phylogenetic relationship, gene duplication, and cis-acting elements were analyzed. Numerous cis-acting elements related to plant development, hormone, and stress responses were identified in the promoter of the VvPLCP genes. Phylogenetic analysis grouped the VvPLCP genes into nine subgroups. The transcription of VvPLCP in different inoculation time points and varieties indicated that VvPLCP may have vital functions in grapevine defense against Plasmopara viticola. According to transcriptome data and qPCR analysis, we observed the increasing expression levels of VvRD21-1 at 72 h after inoculation in resistant variety, inferring that it was related to grape downy mildew resistance. Meanwhile, 3 genes including VvXBCP1, VvSAG12-1, and VvALP1 showed higher expression at 24 h after pathogen inoculation in the susceptible variety and might be related to the downy mildew phenotype. We nominated these four genes to function during hypersensitive response (HR) process, inferring that these genes could be associated with downy mildew resistance in grapes.

CONCLUSIONS

Our results provide the reference for functional studies of PLCP gene family, and highlight its functions in grapevine defense against P. viticola. The results help us to better understand the complexity of the PLCP gene family in plant immunity and provide valuable information for future functional characterization of specific genes in grapevine.

Transcriptional control of local auxin distribution by the CsDFB1-CsPHB module regulates floral organogenesis in cucumber

Auxin is a key phytohormone influencing multiple aspects of plant development, including meristem maintenance, primordia initiation, floral organogenesis, and vascular differentiation. Local auxin biosynthesis and polar auxin transport are essential to establish and maintain auxin gradients that ensure proper plant development. Here, we demonstrate that CsDFB1, a member of the plant cystatin superfamily, which was previously implicated in defense responses, plays a critical role in regulating local auxin distribution and thus influences floral organogenesis in cucumber. Genetic and biochemical assays suggest that CsDFB1 affects local auxin distribution by acting as an attenuator that interacts with CsPHB and modulates CsPHB-mediated transcriptional control of CsYUC2 and CsPIN1. Our results shed light on the fine tuning of local auxin distribution in plants.

Plant cystatins are cysteine proteinase inhibitors that play key roles in defense responses. In this work, we describe an unexpected role for the cystatin-like protein DEFORMED FLORAL BUD1 (CsDFB1) as a transcriptional regulator of local auxin distribution in cucumber (Cucumis sativus L.). CsDFB1 was strongly expressed in the floral meristems, floral primordia, and vasculature. RNA interference (RNAi)-mediated silencing of CsDFB1 led to a significantly increased number of floral organs and vascular bundles, together with a pronounced accumulation of auxin. Conversely, accompanied by a decrease of auxin, overexpression of CsDFB1 resulted in a dramatic reduction in floral organ number and an obvious defect in vascular patterning, as well as organ fusion. CsDFB1 physically interacted with the cucumber ortholog of PHABULOSA (CsPHB), an HD-ZIP III transcription factor whose transcripts exhibit the same pattern as CsDFB1. Overexpression of CsPHB increased auxin accumulation in shoot tips and induced a floral phenotype similar to that of CsDFB1-RNAi lines. Furthermore, genetic and biochemical analyses revealed that CsDFB1 impairs CsPHB-mediated transcriptional regulation of the auxin biosynthetic gene YUCCA2 and the auxin efflux carrier PIN-FORMED1, and thus plays a pivotal role in auxin distribution. In summary, we propose that the CsDFB1-CsPHB module represents a regulatory pathway for local auxin distribution that governs floral organogenesis and vascular differentiation in cucumber.

Characteristics and Research Progress of Legume Nodule Senescence

Delaying the nodule senescence of legume crops can prolong the time of nitrogen fixation and attenuate the lack of fertilizer in the later stage of legume crop cultivation, resulting in improved crop yield and reduced usage of nitrogen fertilizer. However, effective measures to delay the nodule senescence of legume crops in agriculture are relatively lacking. In the present review, we summarized the structural and physiological characteristics of nodule senescence, as well as the corresponding detection methods, providing technical support for the identification of nodule senescence phenotype. We then outlined the key genes currently known to be involved in the regulation of nodule senescence, offering the molecular genetic information for breeding varieties with delayed nodule senescence. In addition, we reviewed various abiotic factors affecting nodule senescence, providing a theoretical basis for the interaction between molecular genetics and abiotic factors in the regulation of nodule senescence. Finally, we briefly prospected research foci of nodule senescence in the future.

Genome-wide survey of soybean papain-like cysteine proteases and their expression analysis in root nodule symbiosis

BACKGROUND

Plant papain-like cysteine proteases (PLCPs) are a large class of proteolytic enzymes and play important roles in root nodule symbiosis (RNS), while the whole-genome studies of PLCP family genes in legume are quite limited, and the roles of Glycine max PLCPs (GmPLCPs) in nodulation, nodule development and senescence are not fully understood.

RESULTS

In the present study, we identified 97 GmPLCPs and performed a genome-wide survey to explore the expansion of soybean PLCP family genes and their relationships to RNS. Nineteen paralogous pairs of genomic segments, consisting of 77 GmPLCPs, formed by whole-genome duplication (WGD) events were identified, showing a high degree of complexity in duplication. Phylogenetic analysis among different species showed that the lineage differentiation of GmPLCPs occurred after family expansion, and large tandem repeat segment were specifically in soybean. The expression patterns of GmPLCPs in symbiosis-related tissues and nodules identified RNS-related GmPLCPs and provided insights into their putative symbiotic functions in soybean. The symbiotic function analyses showed that a RNS-related GmPLCP gene (Glyma.04G190700) really participate in nodulation and nodule development.

CONCLUSIONS

Our findings improved our understanding of the functional diversity of legume PLCP family genes, and provided insights into the putative roles of the legume PLCPs in nodulation, nodule development and senescence.

Responses of mature symbiotic nodules to the whole-plant systemic nitrogen signaling

In symbiotic root nodules of legumes, terminally differentiated rhizobia fix atmospheric N2 producing an NH4+ influx that is assimilated by the plant. The plant, in return, provides photosynthates that fuel the symbiotic nitrogen acquisition. Mechanisms responsible for the adjustment of the symbiotic capacity to the plant N demand remain poorly understood. We have investigated the role of systemic signaling of whole-plant N demand on the mature N2-fixing nodules of the model symbiotic association Medicago truncatula/Sinorhizobium using split-root systems. The whole-plant N-satiety signaling rapidly triggers reductions of both N2 fixation and allocation of sugars to the nodule. These responses are associated with the induction of nodule senescence and the activation of plant defenses against microbes, as well as variations in sugars transport and nodule metabolism. The whole-plant N-deficit responses mirror these changes: a rapid increase of sucrose allocation in response to N-deficit is associated with a stimulation of nodule functioning and development resulting in nodule expansion in the long term. Physiological, transcriptomic, and metabolomic data together provide evidence for strong integration of symbiotic nodules into whole-plant nitrogen demand by systemic signaling and suggest roles for sugar allocation and hormones in the signaling mechanisms.

The rhizobial autotransporter determines the symbiotic nitrogen fixation activity of Lotus japonicus in a host-specific manner

Leguminous plants establish endosymbiotic associations with rhizobia and form root nodules in which the rhizobia fix atmospheric nitrogen. The host plant and intracellular rhizobia strictly control this symbiotic nitrogen fixation. We recently reported a Lotus japonicus Fix^- mutant, apn1 (aspartic peptidase nodule-induced 1), that impairs symbiotic nitrogen fixation. APN1 encodes a nodule-specific aspartic peptidase involved in the Fix^- phenotype in a rhizobial strain-specific manner. This host-strain specificity implies that some molecular interactions between host plant APN1 and rhizobial factors are required, although the biological function of APN1 in nodules and the mechanisms governing the interactions are unknown. To clarify how rhizobial factors are involved in strain-specific nitrogen fixation, we explored transposon mutants of Mesorhizobium loti strain TONO, which normally form Fix^- nodules on apn1 roots, and identified TONO mutants that formed Fix⁺ nodules on apn1 The identified causal gene encodes an autotransporter, part of a protein secretion system of Gram-negative bacteria. Expression of the autotransporter gene in M. loti strain MAFF3030399, which normally forms Fix⁺ nodules on apn1 roots, resulted in Fix^- nodules. The autotransporter of TONO functions to secrete a part of its own protein (a passenger domain) into extracellular spaces, and the recombinant APN1 protein cleaved the passenger protein in vitro. The M. loti autotransporter showed the activity to induce the genes involved in nodule senescence in a dose-dependent manner. Therefore, we conclude that the nodule-specific aspartic peptidase, APN1, suppresses negative effects of the rhizobial autotransporter in order to maintain effective symbiotic nitrogen fixation in root nodules.

Phytocystatins and their Potential Application in the Development of Drought Tolerance Plants in Soybeans (Glycine max L.)

Plant cystatins, also called phytocystatins constitute a family of specific cysteine protease inhibitors found in several monocots and dicots. In soybean, phytocystatins regulate several endogenous processes contributing immensely to this crop's tolerance to abiotic stress factors. Soybeans offer numerous nutritional, pharmaceutical and industrial benefits; however, their growth and yields is hampered by drought, which causes more than 10% yield losses recorded every harvest period worldwide. This review analyses the role of papain-like cysteine proteases and their inhibitors in soybean plant growth and development under drought stress. It also describes their localisation, regulation, target organs and tissues, and the overall impact of cystatins on generating drought tolerance soybean plants. These proteins have many functions that remain poorly characterized, particularly under abiotic stress. Although much information is available on the utilisation of proteases for industrial applications, very few reports have focused on the impact of proteases on plant stress responses. The exploitation of cystatins in plant engineering, as competitive proteases inhibitors is one of the means that will guarantee the continued utilisation of soybeans as an important oilseed crop.

Genome-Wide Identification and Characterization of Vacuolar Processing Enzyme Gene Family and Diverse Expression Under Stress in Apple (Malus × Domestic)

Vacuolar processing enzymes (VPEs) play an important role in stress resistance and development of plants. Despite their diverse roles, little information is available in apple (Malus × domestic). This study firstly presents the genome-wide identification of VPE family genes in apple, resulting in 20 family members those are unevenly distributed across six out of the 17 chromosomes. Phylogenetic analysis assigned these genes into four groups. Analysis of exon-intron junctions and motifs of each candidate gene revealed high levels of conservation within and between phylogenetic groups. Cis-element including w box, ABRE, LTR, and TC-rich repeats were found in promoters of MdVPEs. NCBI-GEO database shown that the expression of MdVPEs exhibited diverse patterns in different tissues as well as the infection of Pythium ultimum and Apple Stem Grooving Virus. Furthermore, qRT-PCR showed that MdVPE genes were responsive to salt, cadmium, low-temperature, and drought. Overexpression of MDP0000172014, which was strongly induced by salt and drought stress, significantly decreased Arabidopsis tolerance to salt stress. The genome-wide identification and characterization of MdVPEs in apple provided basic information for the potential utilization of MdVPEs in stress resistance.

Winged bean (Psophocarpus tetragonolobus (L.) DC.) for food and nutritional security: synthesis of past research and future direction

Winged bean is popularly known as "One Species Supermarket" for its nutrient-dense green pods, immature seeds, tubers, leaves, and mature seeds. This underutilised crop has potential beneficial traits related to its biological nitrogen-fixation to support low-input farming. Drawing from past knowledge, and based on current technologies, we propose a roadmap for research and development of winged bean for sustainable food systems. Reliance on a handful of "major" crops has led to decreased diversity in crop species, agricultural systems and human diets. To reverse this trend, we need to encourage the greater use of minor, "orphan", underutilised species. These could contribute to an increase in crop diversity within agricultural systems, to improve human diets, and to support more sustainable and resilient food production systems. Among these underutilised species, winged bean (Psophocarpus tetragonolobus) has long been proposed as a crop for expanded use particularly in the humid tropics. It is an herbaceous perennial legume of equatorial environments and has been identified as a rich source of protein, with most parts of the plant being edible when appropriately prepared. However, to date, limited progress in structured improvement programmes has restricted the expansion of winged bean beyond its traditional confines. In this paper, we discuss the reasons for this and recommend approaches for better use of its genetic resources and related Psophocarpus species in developing improved varieties. We review studies on the growth, phenology, nodulation and nitrogen-fixation activity, breeding programmes, and molecular analyses. We then discuss prospects for the crop based on the greater understanding that these studies have provided and considering modern plant-breeding technologies and approaches. We propose a more targeted and structured research approach to fulfil the potential of winged bean to contribute to food security.

Plant Vacuolar Processing Enzymes

Plant proteomes contain hundreds of proteases divided into different families based on evolutionary and functional relationship. In particular, plant cysteine proteases of the C1 (papain-like) and C13 (legumain-like) families play key roles in many physiological processes. The legumain-like proteases, also called vacuolar processing enzymes (VPEs), perform a multifunctional role in different plant organs and during different stages of plant development and death. VPEs are similar to animal caspases, and although caspase activity was identified in plants almost 40 years ago, there still remains much research to be done to gain a complete understanding of their various roles and functions in plants. Here we not only summarize the current existing knowledge of plant VPEs, including recent developments in the field, but also highlight the future prospective areas to be investigated to obtain a more detailed understanding of the role of VPEs in plants.

Uncovering Bax inhibitor-1 dual role in the legume-rhizobia symbiosis in common bean roots

Bax-inhibitor 1 (BI-1) is a cell death suppressor conserved in all eukaryotes that modulates cell death in response to abiotic stress and pathogen attack in plants. However, little is known about its role in the establishment of symbiotic interactions. Here, we demonstrate the functional relevance of an Arabidopsis thaliana BI-1 homolog (PvBI-1a) to symbiosis between the common bean (Phaseolus vulgaris) and Rhizobium tropici. We show that the changes in expression of PvBI-1a observed during early symbiosis resemble those of some defence response-related proteins. By using gain- and loss-of-function approaches, we demonstrate that the overexpression of PvBI-1a in the roots of common bean increases the number of rhizobial infection events (and therefore the final number of nodules per root), but induces the premature death of nodule cells, affecting their nitrogen fixation efficiency. Nodule morphological alterations are known to be associated with changes in the expression of genes tied to defence, autophagy, and vesicular trafficking. Results obtained in the present work suggest that BI-1 has a dual role in the regulation of programmed cell death during symbiosis, extending our understanding of its critical function in the modulation of host immunity while responding to beneficial microbes.

Drought-induced senescence of Medicago truncatula nodules involves serpin and ferritin to control proteolytic activity and iron levels

Drought is a major constraint for legume growth and yield. Senescence of nitrogen-fixing nodules is one of the early drought responses and may cause nutrient stress in addition to water stress in legumes. For nodule senescence to function as part of a drought-survival strategy, we propose that the intrinsically destructive senescence process must be tightly regulated. Medicago truncatula protease inhibitor and iron scavenger-encoding genes, possibly involved in controlling nodule senescence, were identified. RNA interference (RNAi) lines were constructed in which expression of a serpin or ferritins was knocked down. Both wild-type and RNAi lines were subjected to drought stress and nodule activity and plant physiological responses were measured. Drought caused M. truncatula to initiate nodule senescence before plant growth was affected and before an increase in papain-like proteolytic activity and free iron levels was apparent. Knock-down expression of serpin6 and ferritins caused increased protease activity, free iron levels, early nodule senescence and reduced plant growth. The results suggest that M. truncatula nodule-expressed serpin6 and ferritins mediate ordered drought-induced senescence by regulating papain-like cysteine protease activity and free iron levels. This strategy may allow the drought-stressed plants to benefit maximally from residual nitrogen fixation and nutrient recovery resulting from break down of macromolecules.

Identification of a type II cystatin in Fragaria chiloensis: A proteinase inhibitor differentially regulated during achene development and in response to biotic stress-related stimuli

The equilibrium between protein synthesis and degradation is key to maintaining efficiency in different physiological processes. The proteinase inhibitor cystatin regulates protease activities in different developmental and physiological contexts. Here we describe for the first time the identification and the biological function of the cysteine protease inhibitor cystatin of Fragaria chiloensis, FchCYS1. Based on primary sequence and 3D-structural homology modelling, FchCYS1 is a type II phytocystatin with high identity to other cystatins of the Fragaria genus. Both the papain-like and the legumain-like protease inhibitory domains are indeed functional, based on in vitro assays performed with Escherichia coli protein extracts containing recombinant FchCYS1. FchCYS1 is differentially-expressed in achenes of F. chiloensis fruits, with highest expression as the fruit reaches the ripened stage, suggesting a role in preventing degradation of storage proteins that will nourish the embryo during seed germination. Furthermore, FchCYS1 responds transcriptionally to the application of salicylic acid and to mechanical injury, strongly suggesting that FchCYS1 could be involved in the response against pathogen attack. Overall these results point to a role for FchCYS1 in diverse physiological processes in F. chiloensis.

Soil environment is a key driver of adaptation in Medicago truncatula: new insights from landscape genomics

Spatial differences in environmental selective pressures interact with the genomes of organisms, ultimately leading to local adaptation. Landscape genomics is an emergent research area that uncovers genome-environment associations, thus allowing researchers to identify candidate loci for adaptation to specific environmental variables. In the present study, we used latent factor mixed models (LFMMs) and Moran spectral outlier detection/randomization (MSOD-MSR) to identify candidate loci for adaptation to 10 environmental variables (climatic, soil and atmospheric) among 43 515 single nucleotide polymorphisms (SNPs) from 202 accessions of the model legume Medicago truncatula. Soil variables were associated with a large number of candidate loci identified through both LFMMs and MSOD-MSR. Genes tagged by candidate loci associated with drought and salinity are involved in the response to biotic and abiotic stresses, while those tagged by candidates associated with soil nitrogen and atmospheric nitrogen, participate in the legume-rhizobia symbiosis. Candidate SNPs identified through both LFMMs and MSOD-MSR explained up to 56% of variance in flowering traits. Our findings highlight the importance of soil in driving adaptation in the system and elucidate the basis of evolutionary potential of M. truncatula to respond to global climate change and anthropogenic disruption of the nitrogen cycle.

Papain-like cysteine proteases in Carica papaya: lineage-specific gene duplication and expansion

Background

Papain-like cysteine proteases (PLCPs), a large group of cysteine proteases structurally related to papain, play important roles in plant development, senescence, and defense responses. Papain, the first cysteine protease whose structure was determined by X-ray crystallography, plays a crucial role in protecting papaya from herbivorous insects. Except the four major PLCPs purified and characterized in papaya latex, the rest of the PLCPs in papaya genome are largely unknown.

Results

We identified 33 PLCP genes in papaya genome. Phylogenetic analysis clearly separated plant PLCP genes into nine subfamilies. PLCP genes are not equally distributed among the nine subfamilies and the number of PLCPs in each subfamily does not increase or decrease proportionally among the seven selected plant species. Papaya showed clear lineage-specific gene expansion in the subfamily III. Interestingly, all four major PLCPs purified from papaya latex, including papain, chymopapain, glycyl endopeptidase and caricain, were grouped into the lineage-specific expansion branch in the subfamily III. Mapping PLCP genes on chromosomes of five plant species revealed that lineage-specific expansions of PLCP genes were mostly derived from tandem duplications. We estimated divergence time of papaya PLCP genes of subfamily III. The major duplication events leading to lineage-specific expansion of papaya PLCP genes in subfamily III were estimated at 48 MYA, 34 MYA, and 16 MYA. The gene expression patterns of the papaya PLCP genes in different tissues were assessed by transcriptome sequencing and qRT-PCR. Most of the papaya PLCP genes of subfamily III expressed at high levels in leaf and green fruit tissues.

Conclusions

Tandem duplications played the dominant role in affecting copy number of PLCPs in plants. Significant variations in size of the PLCP subfamilies among species may reflect genetic adaptation of plant species to different environments. The lineage-specific expansion of papaya PLCPs of subfamily III might have been promoted by the continuous reciprocal selective effects of herbivore attack and plant defense.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4394-y) contains supplementary material, which is available to authorized users.

ARG1 Functions in the Physiological Adaptation of Undifferentiated Plant Cells to Spaceflight

Scientific access to spaceflight and especially the International Space Station has revealed that physiological adaptation to spaceflight is accompanied or enabled by changes in gene expression that significantly alter the transcriptome of cells in spaceflight. A wide range of experiments have shown that plant physiological adaptation to spaceflight involves gene expression changes that alter cell wall and other metabolisms. However, while transcriptome profiling aptly illuminates changes in gene expression that accompany spaceflight adaptation, mutation analysis is required to illuminate key elements required for that adaptation. Here we report how transcriptome profiling was used to gain insight into the spaceflight adaptation role of Altered response to gravity 1 (Arg1), a gene known to affect gravity responses in plants on Earth. The study compared expression profiles of cultured lines of Arabidopsis thaliana derived from wild-type (WT) cultivar Col-0 to profiles from a knock-out line deficient in the gene encoding ARG1 (ARG1 KO), both on the ground and in space. The cell lines were launched on SpaceX CRS-2 as part of the Cellular Expression Logic (CEL) experiment of the BRIC-17 spaceflight mission. The cultured cell lines were grown within 60 mm Petri plates in Petri Dish Fixation Units (PDFUs) that were housed within the Biological Research In Canisters (BRIC) hardware. Spaceflight samples were fixed on orbit. Differentially expressed genes were identified between the two environments (spaceflight and comparable ground controls) and the two genotypes (WT and ARG1 KO). Each genotype engaged unique genes during physiological adaptation to the spaceflight environment, with little overlap. Most of the genes altered in expression in spaceflight in WT cells were found to be Arg1-dependent, suggesting a major role for that gene in the physiological adaptation of undifferentiated cells to spaceflight. Key Words: ARG1-Spaceflight-Gene expression-Physiological adaptation-BRIC. Astrobiology 17, 1077-1111.

Cysteine proteases and wheat (Triticum aestivum L) under drought: A still greatly unexplored association

Bread wheat (Triticum aestivum L.) provides about 19% of global dietary energy. Environmental stress, such as drought, affects wheat growth causing premature plant senescence and ultimately plant death. A plant response to drought is an increase in protease-mediated proteolysis with rapid degradation of proteins required for metabolic processes. Among the plant proteases that are increased in their activity following stress, cysteine proteases are the best characterized. Very little is known about particular wheat cysteine protease sequences, their expression and also localization. The current knowledge on wheat cysteine proteases belonging to the five clans (CA, CD, CE, CF and CP) is outlined, in particular their expression and possible function under drought. The first successes in establishing an annotated wheat genome database are further highlighted which has allowed more detailed mining of cysteine proteases. We also share our thoughts on future research directions considering the growing availability of genomic resources of this very important food crop. Finally, we also outline future application of developed knowledge in transgenic wheat plants for environmental stress protection and also as senescence markers to monitor wheat growth under environmental stress conditions.

Use of CRISPR/Cas9 for Symbiotic Nitrogen Fixation Research in Legumes

Nitrogen-fixing rhizobia have established a symbiotic relationship with the legume family through more than 60 million years of evolution. Hundreds of legume host genes are involved in the SNF (symbiotic nitrogen fixation) process, such as recognition of the bacterial partners, nodulation signaling and nodule development, maintenance of highly efficient nitrogen fixation within nodules, regulation of nodule numbers, and nodule senescence. However, investigations of SNF-related gene functions and dissecting molecular mechanisms of the complicated signaling crosstalk on a genomic scale were significantly restricted by insufficient mutant resources of several representative model legumes. Targeted genome-editing technologies, including ZFNs, TALENs, and CRISPR-Cas systems, have been developed in recent years and rapidly revolutionized biological research in many fields. These technologies were also applied to legume plants, and significant progress has been made in the last several years. Here, we summarize the applications of these genome-editing technologies, especially CRISPR-Cas9, toward the study of SNF in legumes, which should greatly advance our understanding of the basic mechanisms underpinning the legume-rhizobia interactions and guide the engineering of the SNF pathway into nonlegume crops to reduce the dependence on the use of nitrogen fertilizers for sustainable development of modern agriculture.

RNA-Seq analysis of nodule development at five different developmental stages of soybean (Glycine max) inoculated with Bradyrhizobium japonicum strain 113-2

Nodule development directly affects nitrogen fixation efficiency during soybean growth. Although abundant genome-based information related to nodule development has been released and some studies have reported the molecular mechanisms that regulate nodule development, information on the way nodule genes operate in nodule development at different developmental stages of soybean is limited. In this report, notably different nodulation phenotypes in soybean roots inoculated with Bradyrhizobium japonicum strain 113-2 at five developmental stages (branching stage, flowering stage, fruiting stage, pod stage and harvest stage) were shown, and the expression of nodule genes at these five stages was assessed quantitatively using RNA-Seq. Ten comparisons were made between these developmental periods, and their differentially expressed genes were analysed. Some important genes were identified, primarily encoding symbiotic nitrogen fixation-related proteins, cysteine proteases, cystatins and cysteine-rich proteins, as well as proteins involving plant-pathogen interactions. There were no significant shifts in the distribution of most GO functional annotation terms and KEGG pathway enrichment terms between these five development stages. A cystatin Glyma18g12240 was firstly identified from our RNA-seq, and was likely to promote nodulation and delay nodule senescence. This study provides molecular material for further investigations into the mechanisms of nitrogen fixation at different soybean developmental stages.

Drought Stress Responses in Soybean Roots and Nodules

Drought is considered to be a major threat to soybean production worldwide and yet our current understanding of the effects of drought on soybean productively is largely based on studies on above-ground traits. Although the roots and root nodules are important sensors of drought, the responses of these crucial organs and their drought tolerance features remain poorly characterized. The symbiotic interaction between soybean and rhizobia facilitates atmospheric nitrogen fixation, a process that provides essential nitrogen to support plant growth and development. Symbiotic nitrogen fixation is important for sustainable agriculture, as it sustains plant growth on nitrogen-poor soils and limits fertilizer use for crop nitrogen nutrition. Recent developments have been made in our understanding of the drought impact on soybean root architecture and nodule traits, as well as underpinning transcriptome, proteome and also emerging metabolome information, with a view to improve the selection of more drought-tolerant soybean cultivars and rhizobia in the future. We conclude that the direct screening of root and nodule traits in the field as well as identification of genes, proteins and also metabolites involved in such traits will be essential in order to gain a better understanding of the regulation of root architecture, bacteroid development and lifespan in relation to drought tolerance in soybean.

Sequence comparison, molecular modeling, and network analysis predict structural diversity in cysteine proteases from the Cape sundew, Drosera capensis

Carnivorous plants represent a so far underexploited reservoir of novel proteases with potentially useful activities. Here we investigate 44 cysteine proteases from the Cape sundew, Drosera capensis, predicted from genomic DNA sequences. D. capensis has a large number of cysteine protease genes; analysis of their sequences reveals homologs of known plant proteases, some of which are predicted to have novel properties. Many functionally significant sequence and structural features are observed, including targeting signals and occluding loops. Several of the proteases contain a new type of granulin domain. Although active site residues are conserved, the sequence identity of these proteases to known proteins is moderate to low; therefore, comparative modeling with all-atom refinement and subsequent atomistic MD-simulation is used to predict their 3D structures. The structure prediction data, as well as analysis of protein structure networks, suggest multifarious variations on the papain-like cysteine protease structural theme. This in silico methodology provides a general framework for investigating a large pool of sequences that are potentially useful for biotechnology applications, enabling informed choices about which proteins to investigate in the laboratory.

•

44 new cysteine proteases from the carnivorous plant Drosera capensis are described.

•

Structure prediction and molecular dynamics simulation predict overall folds similar to papain.

•

Functionally significant sequence and structural features are observed, including targeting signals and occluding loops.

•

Several of the proteases contain a new type of granulin domain.

•

Protein structure networks reveal global differences in interactions among chemical groups.

Targeted expression of cystatin restores fertility in cysteine protease induced male sterile tobacco plants

Fertility restoration in male sterile plants is an essential requirement for their utilization in hybrid seed production. In an earlier investigation, we have demonstrated that the targeted expression of a cysteine protease in tapetal cell layer resulted in complete male sterility in tobacco transgenic plants. In the present investigation, we have used a cystatin gene, which encodes for a cysteine protease inhibitor, from a wild peanut, Arachis diogoi and developed a plant gene based restoration system for cysteine protease induced male sterile transgenic tobacco plants. We confirmed the interaction between the cysteine protease and a cystatin of the wild peanut, A. diogoi through in silico modeling and yeast two-hybrid assay. Pollen from primary transgenic tobacco plants expressing cystatin gene under the tapetum specific promoter- TA29 restored fertility on cysteine protease induced male sterile tobacco plants developed earlier. This has confirmed the in vivo interaction of cysteine protease and cystatin in the tapetal cells, and the inactivation of cysteine protease and modulation of its negative effects on pollen fertility. Both the cysteine protease and cystatin genes are of plant origin in contrast to the analogous barnase-barstar system that deploys genes of prokaryotic origin. Because of the deployment of genes of plant origin, this system might not face biosafety problems in developing hybrids in food crops.

Review: The future of cystatin engineering

Plant cystatins are naturally occurring protease inhibitors that prevent proteolysis by papain-like cysteine proteases. Their protective action against environmental stresses has been relatively well characterised. Still, there is a need to greatly improve both potency and specificity based on the current rather poor performance of cystatins in biotechnological applications. Research in creating more potent and specific cystatins, including amino acid substitutions in either conserved cystatin motifs and/or at variable amino acid sites, is reviewed. Existing gaps for better understanding of cystatin-protease interactions are further explored. Current knowledge on multi-cystatins or hybrid protease inhibitors involving cystatins as an additional option for cystatin engineering is further outlined along with the nuances of how cystatins with rather unusual amino acid sequences might actually help in cystatin engineering. Finally, future opportunities for application of cystatins are highlighted which include applications in genetically modified transgenic plants for environmental stress protection and also as nutraceuticals, as part of more nutritious food. Further opportunities might also include the possible management of diseases and disorders, often associated with lifestyle changes, and the most immediate and promising application which is inclusion into plant-based recombinant protein production platforms.

Search for Nodulation and Nodule Development-Related Cystatin Genes in the Genome of Soybean (Glycine max)

Nodulation, nodule development and senescence directly affects nitrogen fixation efficiency, and previous studies have shown that inhibition of some cysteine proteases delay nodule senescence, so their nature inhibitors, cystatin genes, are very important in nodulation, nodule development, and senescence. Although several cystatins are actively transcribed in soybean nodules, their exact roles and functional diversities in legume have not been well explored in genome-wide survey studies. In this report, we performed a genome-wide survey of cystatin family genes to explore their relationship to nodulation and nodule development in soybean and identified 20 cystatin genes that encode peptides with 97-245 amino acid residues, different isoelectric points (pI) and structure characteristics, and various putative plant regulatory elements in 3000 bp putative promoter fragments upstream of the 20 soybean cystatins in response to different abiotic/biotic stresses, hormone signals, and symbiosis signals. The expression profiles of these cystatin genes in soybean symbiosis with rhizobium strain Bradyrhizobium japonicum strain 113-2 revealed that 7 cystatin family genes play different roles in nodulation as well as nodule development and senescence. However, these genes were not root nodule symbiosis (RNS)-specific and did not encode special clade cystatin protein with structures related to nodulation and nodule development. Besides, only two of these soybean cystatins were not upregulated in symbiosis after ABA treatment. The functional analysis showed that a candidate gene Glyma.15G227500 (GmCYS16) was likely to play a positive role in soybean nodulation. Besides, evolutionary relationships analysis divided the cystatin genes from Arabidopsis thaliana, Nicotiana tabacum, rice, barley and four legume plants into three groups. Interestingly, Group A cystatins are special in legume plants, but only include one of the above-mentioned 7 cystatin genes related to nodulation and nodule development. Overall, our results provide useful information or clues for our understanding of the functional diversity of legume cystatin family proteins in soybean nodulation and nodule development and for finding nodule-specific cysteine proteases in soybean.

Cellular and Molecular Changes Associated with Onion Skin Formation Suggest Involvement of Programmed Cell Death

Skin formation of onion (Allium cepa L.) bulb involves scale desiccation accompanied by scale senescence, resulting in cell death and tissue browning. Understanding the mechanism of skin formation is essential to improving onion skin and bulb qualities. Although onion skin plays a crucial role in postharvest onion storage and shelf life, its formation is poorly understood. We investigated the mode of cell death in the outermost scales that are destined to form the onion skin. Surprisingly, fluorescein diacetate staining and scanning electron microscopy indicated that the outer scale desiccates from the inside out. This striking observation suggests that cell death in the outer scales, during skin formation, is an internal and organized process that does not derive only from air desiccation. DNA fragmentation, a known hallmark of programmed cell death (PCD), was revealed in the outer scales and gradually decreased toward the inner scales of the bulb. Transmission electron microscopy further revealed PCD-related structural alterations in the outer scales which were absent from the inner scales. De novo transcriptome assembly for three different scales: 1st (outer), 5th (intermediate) and 8th (inner) fleshy scales identified 2,542 differentially expressed genes among them. GO enrichment for cluster analysis revealed increasing metabolic processes in the outer senescent scale related to defense response, PCD processes, carbohydrate metabolism and flavonoid biosynthesis, whereas increased metabolism and developmental growth processes were identified in the inner scales. High expression levels of PCD-related genes were identified in the outer scale compared to the inner ones, highlighting the involvement of PCD in outer-skin development. These findings suggest that a program to form the dry protective skin exists and functions only in the outer scales of onion.

Potential use of phytocystatins in crop improvement, with a particular focus on legumes

Phytocystatins are a well-characterized class of naturally occurring protease inhibitors that function by preventing the catalysis of papain-like cysteine proteases. The action of cystatins in biotic stress resistance has been studied intensively, but relatively little is known about their functions in plant growth and defence responses to abiotic stresses, such as drought. Extreme weather events, such as drought and flooding, will have negative impacts on the yields of crop plants, particularly grain legumes. The concepts that changes in cellular protein content and composition are required for acclimation to different abiotic stresses, and that these adjustments are achieved through regulation of proteolysis, are widely accepted. However, the nature and regulation of the protein turnover machinery that underpins essential stress-induced cellular restructuring remain poorly characterized. Cysteine proteases are intrinsic to the genetic programmes that underpin plant development and senescence, but their functions in stress-induced senescence are not well defined. Transgenic plants including soybean that have been engineered to constitutively express phytocystatins show enhanced tolerance to a range of different abiotic stresses including drought, suggesting that manipulation of cysteine protease activities by altered phytocystatin expression in crop plants might be used to improve resilience and quality in the face of climate change.