Subtilases - versatile tools for protein turnover, plant development, and interactions with the environment

Phytaspase Does Not Require Proteolytic Activity for Its Stress-Induced Internalization

Phytaspases differ from other members of the plant subtilisin-like protease family by having rare aspartate cleavage specificity and unusual localization dynamics. Phytaspases are secreted from healthy plant cells but are re-internalized upon perception of death-inducing stresses. Although proteolytic activity is required for the secretion of plant subtilases, its requirement for the retrograde transportation of phytaspases is currently unknown. To address this issue, we employed an approach to complement in trans the externalization of a prodomain-less form of Nicotiana tabacum phytaspase (NtPhyt) with the free prodomain in Nicotiana benthamiana leaf cells. Using this approach, the generation of the proteolytically active NtPhyt and its transport to the extracellular space at a level comparable to that of the native NtPhyt (synthesized as a canonical prodomain-containing precursor protein) were achieved. The application of this methodology to NtPhyt with a mutated catalytic Ser537 residue resulted in the secretion of the inactive, although processed (prodomain-free), protein as well. Notably, the externalized NtPhyt Ser537Ala mutant was still capable of retrograde transportation into plant cells upon the induction of oxidative stress. Our data thus indicate that the proteolytic activity of NtPhyt is dispensable for stress-induced retrograde transport of the enzyme.

Deciphering the regulatory networks involved in mild and severe salt stress responses in the roots of wild grapevine Vitis vinifera spp. sylvestris

Transcriptional regulatory networks are pivotal components of plant's response to salt stress. However, plant adaptation strategies varied as a function of stress intensity, which is mainly modulated by climate change. Here, we determined the gene regulatory networks based on transcription factor (TF) TF_gene co-expression, using two transcriptomic data sets generated from the salt-tolerant "Tebaba" roots either treated with 50 mM NaCl (mild stress) or 150 mM NaCl (severe stress). The analysis of these regulatory networks identified specific TFs as key regulatory hubs as evidenced by their multiple interactions with different target genes related to stress response. Indeed, under mild stress, NAC and bHLH TFs were identified as central hubs regulating nitrogen storage process. Moreover, HSF TFs were revealed as a regulatory hub regulating various aspects of cellular metabolism including flavonoid biosynthesis, protein processing, phenylpropanoid metabolism, galactose metabolism, and heat shock proteins. These processes are essentially linked to short-term acclimatization under mild salt stress. This was further consolidated by the protein-protein interaction (PPI) network analysis showing structural and plant growth adjustment. Conversely, under severe salt stress, dramatic metabolic changes were observed leading to novel TF members including MYB family as regulatory hubs controlling isoflavonoid biosynthesis, oxidative stress response, abscisic acid signaling pathway, and proteolysis. The PPI network analysis also revealed deeper stress defense changes aiming to restore plant metabolic homeostasis when facing severe salt stress. Overall, both the gene co-expression and PPI network provided valuable insights on key transcription factor hubs that can be employed as candidates for future genetic crop engineering programs.

Stomatal density suppressor PagSDD1 is a "generalist" gene that promotes plant growth and improves water use efficiency

The serine protease SDD1 regulates stomatal density, but its potential impact on plant vegetative growth is unclear. Our study reveals a substantial upregulation of SDD1 in triploid poplar apical buds and leaves, suggesting its possible role in their growth regulation. We cloned PagSDD1 from poplar 84 K (Populus alba × P. glandulosa) and found that overexpression in poplar, soybean, and lettuce led to decreased leaf stomatal density. Furthermore, PagSDD1 represses PagEPF1, PagEPF2, PagEPFL9, PagSPCH, PagMUTE, and PagFAMA expression. In contrast, PagSDD1 promotes the expression of its receptors, PagTMM and PagERECTA. PagSDD1-OE poplars showed stronger drought tolerance than wild-type poplars. Simultaneously, PagSDD1-OE poplar, soybean, and lettuce had vegetative growth advantages. RNA sequencing revealed a significant upregulation of genes PagLHCB2.1 and PagGRF5, correlating positively with photosynthetic rate, and PagCYCA3;4 and PagEXPA8 linked to cell division and differentiation in PagSDD1-OE poplars. This increase promoted leaf photosynthesis, boosted auxin and cytokinin accumulation, and enhanced vegetative growth. SDD1 overexpression can increase the biomass of poplar, soybean, and lettuce by approximately 70, 176, and 155 %, respectively, and increase the water use efficiency of poplar leaves by over 52 %, which is of great value for the molecular design and breeding of plants with growth and water-saving target traits.

Comprehensive transcriptome analyses of Fusarium-infected root xylem tissues to decipher genes involved in chickpea wilt resistance

Fusarium wilt is the most destructive soil-borne disease that poses a major threat to chickpea production. To comprehensively understand the interaction between chickpea and Fusarium oxysporum, the xylem-specific transcriptome analysis of wilt-resistant (WR315) and wilt-susceptible (JG62) genotypes at an early timepoint (4DPI) was investigated. Differential expression analysis showed that 1368 and 348 DEGs responded to pathogen infection in resistant and susceptible genotypes, respectively. Both genotypes showed transcriptional reprogramming in response to Foc2, but the responses in WR315 were more severe than in JG62. Results of the KEGG pathway analysis revealed that most of the DEGS in both genotypes with enrichment in metabolic pathways, secondary metabolite biosynthesis, plant hormone signal transduction, and carbon metabolism. Genes associated with defense-related metabolites synthesis such as thaumatin-like protein 1b, cysteine-rich receptor-like protein kinases, MLP-like proteins, polygalacturonase inhibitor 2-like, ethylene-responsive transcription factors, glycine-rich cell wall structural protein-like, beta-galactosidase-like, subtilisin-like protease, thioredoxin-like protein, chitin elicitor receptor kinase-like, proline transporter-like, non-specific lipid transfer protein and sugar transporter were mostly up-regulated in resistant as compared to susceptible genotypes. The results of this study provide disease resistance genes, which would be helpful in understanding the Foc resistance mechanism in chickpea.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03803-9.

A systematic analysis of the subtilase gene family and expression and subcellular localization investigation of anther-specific members in maize

Subtilases (SBTs), also known as Subtilisin-like serine proteases, are extracellular alkaline protease proteins. SBTs function in all stages of plant growth, development and stress responses. Maize (Zea mays L.) is a crop widely used worldwide as food, feed, and industrial materials. However, information about the members and their functions of the SBT proteins in maize is lacking. In this study, we identified 58 ZmSBT genes from the maize genome and conducted a comprehensive investigation of ZmSBTs by phylogenetic, gene duplication event, gene structure, and protein conserved motif analyses. The ZmSBT proteins were phylogenetically classified into seven groups, and collinearity analysis indicated that many ZmSBTs originate from tandem or segmental duplications. Structural and homolog protein comparison revealed ZmSBTs have conserved protein structures with reported subtilase proteins, suggesting the conserved functions. Further analysis showed that ZmSBTs are expressed in different tissues, and many are responses to specific abiotic stress. Analysis of the anther-specific ZmSBT genes showed their expression peaked at different developmental stages of maize anthers. Subcellular localization analysis of selected maize ZmSBTs showed they are located in different cellular compartments. The information provided in this study is valuable for further functional study of ZmSBTs.

Genome-Wide Identification, Phylogeny and Expression Analysis of Subtilisin (SBT) Gene Family under Wheat Biotic and Abiotic Stress

The subtilisin-like protease (SBT) family is widely known for its role in stress resistance to a number of stressors in different plant species, but is rarely studied in wheat. Subtilisin-like serine proteases (SBTs) are serine proteolytic enzymes that hydrolyze proteins into small peptides, which bind to receptors as signal molecules or ligands and participate in signal transduction. In this study, we identified 255 putative SBT genes from the wheat reference genome and then divided these into seven clades. Subsequently, we performed syntenic relation analysis, exon-intron organization, motif composition, and cis-element analysis. Further, expression analysis based on RNA-seq and tissue-specific expression patterns revealed that TaSBT gene family expression has multiple intrinsic functions during various abiotic and biotic stresses. Analysis of RNA-seq expression assays and further validation through qRT PCR suggested that some of the TaSBT genes have significant changes in expression levels during Pst interaction. TaSBT7, TaSBT26, TaSBT102, and TaSBT193 genes showed increasing expression levels during compatible and non-compatible interactions, while the expression levels of TaSBT111 and TaSBT213 showed a decreasing trend, indicating that these members of the wheat SBT gene family may have a role in wheat's defense against pathogens. In conclusion, these results expand our understanding of the SBT gene family, and provide a valuable reference for future research on the stress resistance function and comprehensive data of wheat SBT members.

Genome-Wide Investigation and Co-Expression Network Analysis of SBT Family Gene in Gossypium

Subtilases (SBTs), which belong to the serine peptidases, control plant development by regulating cell wall properties and the activity of extracellular signaling molecules, and affect all stages of the life cycle, such as seed development and germination, and responses to biotic and abiotic environments. In this study, 146 Gossypium hirsutum, 138 Gossypium barbadense, 89 Gossypium arboreum and 84 Gossypium raimondii SBTs were identified and divided into six subfamilies. Cotton SBTs are unevenly distributed on chromosomes. Synteny analysis showed that the members of SBT1 and SBT4 were expanded in cotton compared to Arabidopsis thaliana. Co-expression network analysis showed that six Gossypium arboreum SBT gene family members were in a network, among which five SBT1 genes and their Gossypium hirsutum and Arabidopsis thaliana direct homologues were down-regulated by salt treatment, indicating that the co-expression network might share conserved functions. Through co-expression network and annotation analysis, these SBTs may be involved in the biological processes of auxin transport, ABA signal transduction, cell wall repair and root tissue development. In summary, this study provides valuable information for the study of SBT genes in cotton and excavates SBT genes in response to salt stress, which provides ideas for cotton breeding for salinity resistance.

Translational profile of coding and non-coding RNAs revealed by genome wide profiling of ribosome footprints in grapevine

Translation is a crucial process during plant growth and morphogenesis. In grapevine (Vitis vinifera L.), many transcripts can be detected by RNA sequencing; however, their translational regulation is still largely unknown, and a great number of translation products have not yet been identified. Here, ribosome footprint sequencing was carried out to reveal the translational profile of RNAs in grapevine. A total of 8291 detected transcripts were divided into four parts, including the coding, untranslated regions (UTR), intron, and intergenic regions, and the 26 nt ribosome-protected fragments (RPFs) showed a 3 nt periodic distribution. Furthermore, the predicted proteins were identified and classified by GO analysis. More importantly, 7 heat shock-binding proteins were found to be involved in molecular chaperone DNA J families participating in abiotic stress responses. These 7 proteins have different expression patterns in grape tissues; one of them was significantly upregulated by heat stress according to bioinformatics research and was identified as DNA JA6. The subcellular localization results showed that VvDNA JA6 and VvHSP70 were both localized on the cell membrane. Therefore, we speculate that DNA JA6 may interact with HSP70. In addition, overexpression of VvDNA JA6 and VvHSP70, reduced the malondialdehyde (MDA) content, improved the antioxidant enzyme activity of superoxide dismutase (SOD), catalase (CAT) and peroxidase (POD), increased the content of proline, an osmolyte substance, and affected the expression of the high-temperature marker genes VvHsfB1, VvHsfB2A, VvHsfC and VvHSP100. In summary, our study proved that VvDNA JA6 and the heat shock protein VvHSP70 play a positive role in the response to heat stress. This study lays a foundation for further exploring the balance between gene expression and protein translation in grapevine under heat stress.

Ecotin: A versatile protease inhibitor of bacteria and eukaryotes

Serine protease inhibitors are a large family of proteins involved in important pathways and processes, such as inflammatory responses and blood clotting. Most are characterized by a precise mode of action, thereby targeting a narrow range of protease substrates. However, the serine-protease inhibitor ecotin is able to inhibit a broad range of serine proteases that display a wide range of specificities. This specificity is driven by special structural features which allow unique flexibility upon binding to targets. Although frequently observed in many human/animal-associated bacteria, ecotin homologs may also be found in plant-associated taxa and environmental species. The purpose of this review is to provide an update on the biological importance, role in host-microbe interactions, and evolutionary relationship between ecotin orthologs isolated from Eukaryotic and Prokaryotic species across the Tree of Life.

Genome-Wide Identification and Analysis of the Maize Serine Peptidase S8 Family Genes in Response to Drought at Seedling Stage

Subtilisin-like proteases (subtilases) are found in almost all plant species and are involved in regulating various biotic and abiotic stresses. Although the literature on subtilases in different plant species is vast, the gene function of the serine peptidase S8 family and its maize subfamily is still unknown. Here, a bioinformatics analysis of this gene family was conducted by describing gene structure, conserved motifs, phylogenetic relationships, chromosomal distributions, gene duplications, and promoter cis-elements. In total, we identified 18 ZmSPS8 genes in maize, distributed on 7 chromosomes, and half of them were hydrophilic. Most of these proteins were located at the cell wall and had similar secondary and tertiary structures. Prediction of cis-regulatory elements in promoters illustrated that they were mainly associated with hormones and abiotic stress. Maize inbred lines B73, Zheng58, and Qi319 were used to analyze the spatial-temporal expression patterns of ZmSPS8 genes under drought treatment. Seedling drought results showed that Qi319 had the highest percent survival after 14 d of withholding irrigation, while B73 was the lowest. Leaf relative water content (LRWC) declined more rapidly in B73 and to lower values, and the nitrotetrazolium blue chloride (NBT) contents of leaves were higher in Qi319 than in the other inbreds. The qPCR results indicated that 6 serine peptidase S8 family genes were positively or negatively correlated with plant tolerance to drought stress. Our study provides a detailed analysis of the ZmSPS8s in the maize genome and finds a link between drought tolerance and the family gene expression, which was established by using different maize inbred lines.

A Multi-Year, Multi-Cultivar Approach to Differential Expression Analysis of High- and Low-Protein Soybean (Glycine max)

Soybean (Glycine max (L.) Merr.) is among the most valuable crops based on its nutritious seed protein and oil. Protein quality, evaluated as the ratio of glycinin (11S) to β-conglycinin (7S), can play a role in food and feed quality. To help uncover the underlying differences between high and low protein soybean varieties, we performed differential expression analysis on high and low total protein soybean varieties and high and low 11S soybean varieties grown in four locations across Eastern and Western Canada over three years (2018-2020). Simultaneously, ten individual differential expression datasets for high vs. low total protein soybeans and ten individual differential expression datasets for high vs. low 11S soybeans were assessed, for a total of 20 datasets. The top 15 most upregulated and the 15 most downregulated genes were extracted from each differential expression dataset and cross-examination was conducted to create shortlists of the most consistently differentially expressed genes. Shortlisted genes were assessed for gene ontology to gain a global appreciation of the commonly differentially expressed genes. Genes with roles in the lipid metabolic pathway and carbohydrate metabolic pathway were differentially expressed in high total protein and high 11S soybeans in comparison to their low total protein and low 11S counterparts. Expression differences were consistent between East and West locations with the exception of one, Glyma.03G054100. These data are important for uncovering the genes and biological pathways responsible for the difference in seed protein between high and low total protein or 11S cultivars.

Broad-spectrum resistance mechanism of serine protease Sp1 in Bacillus licheniformis W10 via dual comparative transcriptome analysis

Antagonistic microorganisms are considered to be the most promising biological controls for plant disease. However, they are still not as popular as chemical pesticides due to complex environmental factors in the field. It is urgent to exploit their potential genetic characteristics and excellent properties to develop biopesticides with antimicrobial substances as the main components. Here, the serine protease Sp1 isolated from the Bacillus licheniformis W10 strain was confirmed to have a broad antifungal and antibacterial spectrum. Sp1 treatment significantly inhibited fungal vegetative growth and damaged the structure of hyphae, in accordance with that caused by W10 strain. Furthermore, Sp1 could activate the systemic resistance of peach twigs, fruits and tobacco. Dual comparative transcriptome analysis uncovered how Sp1 resisted the plant pathogenic fungus Phomopsis amygdali and the potential molecular resistance mechanisms of tobacco. In PSp1 vs. P. amygdali, RNA-seq identified 150 differentially expressed genes (DEGs) that were upregulated and 209 DEGs that were downregulated. Further analysis found that Sp1 might act on the energy supply and cell wall structure to inhibit the development of P. amygdali. In TSp1 vs. Xanthi tobacco, RNA-seq identified that 5937 DEGs were upregulated and 2929 DEGs were downregulated. DEGs were enriched in the metabolic biosynthesis pathways of secondary metabolites, plant hormone signal transduction, plant–pathogen interactions, and MAPK signaling pathway–plant and further found that the genes of salicylic acid (SA) and jasmonic acid (JA) signaling pathways were highly expressed and the contents of SA and JA increased significantly, suggesting that systemic resistance induced by Sp1 shares features of SAR and ISR. In addition, Sp1 might induce the plant defense responses of tobacco. This study provides insights into the broad-spectrum resistance molecular mechanism of Sp1, which could be used as a potential biocontrol product.

A secreted fungal subtilase interferes with rice immunity via degradation of SUPPRESSOR OF G2 ALLELE OF skp1

Serine protease subtilase, found widely in both eukaryotes and prokaryotes, participates in various biological processes. However, how fungal subtilase regulates plant immunity is a major concern. Here, we identified a secreted fungal subtilase, UvPr1a, from the rice false smut (RFS) fungus Ustilaginoidea virens. We characterized UvPr1a as a virulence effector localized to the plant cytoplasm that inhibits plant cell death induced by Bax. Heterologous expression of UvPr1a in rice (Oryza sativa) enhanced plant susceptibility to rice pathogens. UvPr1a interacted with the important rice protein SUPPRESSOR OF G2 ALLELE OF skp1 (OsSGT1), a positive regulator of innate immunity against multiple rice pathogens, degrading OsSGT1 in a protease activity-dependent manner. Furthermore, host-induced gene silencing of UvPr1a compromised disease resistance of rice plants. Our work reveals a previously uncharacterized fungal virulence strategy in which a fungal pathogen secretes a subtilase to interfere with rice immunity through degradation of OsSGT1, thereby promoting infection. These genetic resources provide tools for introducing RFS resistance and further our understanding of plant-pathogen interactions.

Structure-based discovery of small molecules improving stability of human broadly-neutralizing anti-HIV antibody 2F5 in plant suspension cells

The production of biopharmaceuticals in engineered plant-based systems is a promising technology that has proven its suitability for the production of various recombinant glyco-proteins that are currently undergoing clinical trials. However, compared to mammalian cell lines, the productivity of plant-based systems still requires further improvement. A major obstacle is the proteolytic degradation of recombinant target proteins by endogenous plant proteases mainly from the subtilisin family of serine proteases. In the present study, the authors screened for putative small molecule inhibitors for subtilases that are secreted from tobacco BY-2 suspension cells using an in silico approach. The effectiveness of the substances identified in this screen was subsequently tested in degradation assays using the human broadly-neutralizing anti-HIV monoclonal antibody 2F5 (mAb2F5) and spent BY-2 culture medium as a model system. Among 16 putative inhibitors identified by in silico studies, three naphthalene sulfonic acid derivatives showed inhibitory activity in in vitro degradation assays and are similar to or even more effective than phenylmethylsulfonyl fluoride (PMSF), a classical inhibitor of serine proteases, which served as positive control.

Genome-Wide Association Study of Soybean Germplasm Derived From Canadian × Chinese Crosses to Mine for Novel Alleles to Improve Seed Yield and Seed Quality Traits

Genome-wide association study (GWAS) has emerged in the past decade as a viable tool for identifying beneficial alleles from a genomic diversity panel. In an ongoing effort to improve soybean [Glycine max (L.) Merr.], which is the third largest field crop in Canada, a GWAS was conducted to identify novel alleles underlying seed yield and seed quality and agronomic traits. The genomic panel consisted of 200 genotypes including lines derived from several generations of bi-parental crosses between modern Canadian × Chinese cultivars (CD-CH). The genomic diversity panel was field evaluated at two field locations in Ontario in 2019 and 2020. Genotyping-by-sequencing (GBS) was conducted and yielded almost 32 K high-quality SNPs. GWAS was conducted using Fixed and random model Circulating Probability Unification (FarmCPU) model on the following traits: seed yield, seed protein concentration, seed oil concentration, plant height, 100 seed weight, days to maturity, and lodging score that allowed to identify five QTL regions controlling seed yield and seed oil and protein content. A candidate gene search identified a putative gene for each of the three traits. The results of this GWAS study provide insight into potentially valuable genetic resources residing in Chinese modern cultivars that breeders may use to further improve soybean seed yield and seed quality traits.

Genome-Wide Investigation of SBT Family Genes in Pineapple and Functional Analysis of AcoSBT1.12 in Floral Transition

SBT (Subtilisin-like serine protease), a clan of serine proteolytic enzymes, plays a versatile role in plant growth and defense. Although SBT family genes have been obtained from studies of dicots such as Arabidopsis, little is known about the potential functions of SBT in the monocots. In this study, 54 pineapple SBT genes (AcoSBTs) were divided into six subfamilies and then identified to be experienced strong purifying selective pressure and distributed on 25 chromosomes unevenly. Cis-acting element analysis indicated that almost all AcoSBTs promoters contain light-responsive elements. Further, the expression pattern via RNA-seq data showed that different AcoSBTs were preferentially expressed in different above-ground tissues. Transient expression in tobacco showed that AcoSBT1.12 was located in the plasma membrane. Moreover, Transgenic Arabidopsis ectopically overexpressing AcoSBT1.12 exhibited delayed flowering time. In addition, under the guidance of bioinformatic prediction, we found that AcoSBT1.12 could interact with AcoCWF19L, AcoPUF2, AcoCwfJL, Aco012905, and AcoSZF1 by yeast-two hybrid (Y2H). In summary, this study provided valuable information on pineapple SBT genes and illuminated the biological function of AcoSBT1.12 in floral transition.

Transcriptome analysis reveals key genes associated with root-lesion nematode Pratylenchus thornei resistance in chickpea

The root-lesion nematode, Pratylenchus thornei, is one of the major plant-parasitic nematode species causing significant yield losses in chickpea (Cicer arietinum). In order to identify the underlying mechanisms of resistance to P. thornei, the transcriptomes of control and inoculated roots of three chickpea genotypes viz. D05253 > F3TMWR2AB001 (resistant advanced breeding line), PBA HatTrick (moderately resistant cultivar), and Kyabra (susceptible cultivar) were studied at 20 and 50 days post inoculation using the RNA-seq approach. On analyzing the 633.3 million reads generated, 962 differentially expressed genes (DEGs) were identified. Comparative analysis revealed that the majority of DEGs upregulated in the resistant genotype were downregulated in the moderately resistant and susceptible genotypes. Transcription factor families WRKY and bZIP were uniquely expressed in the resistant genotype. The genes Cysteine-rich receptor-like protein kinase 10, Protein lifeguard-like, Protein detoxification, Bidirectional sugar transporter Sugars Will Eventually be Exported Transporters1 (SWEET1), and Subtilisin-like protease were found to play cross-functional roles in the resistant chickpea genotype against P. thornei. The identified candidate genes for resistance to P. thornei in chickpea can be explored further to develop markers and accelerate the introgression of P. thornei resistance into elite chickpea cultivars.

Rhizospheric Communication through Mobile Genetic Element Transfers for the Regulation of Microbe-Plant Interactions

Simple Summary

Rhizosphere, where microbes and plants coexist, is a hotspot of mobile genetic element (MGE) transfers. It was suggested that ancient MGE transfers drove the evolution of both microbes and plants. On the other hand, recurrent MGE transfers regulate microbe-plant interaction and the adaptation of microbes and plants to the environment. The studies of MGE transfers in the rhizosphere provide useful information for the research on pathogenic/ beneficial microbe-plant interaction. In addition, MGE transfers between microbes and the influence by plant root exudates on such transfers provide useful information for the research on bioremediation.

Abstract

The transfer of mobile genetic elements (MGEs) has been known as a strategy adopted by organisms for survival and adaptation to the environment. The rhizosphere, where microbes and plants coexist, is a hotspot of MGE transfers. In this review, we discuss the classic mechanisms as well as novel mechanisms of MGE transfers in the rhizosphere. Both intra-kingdom and cross-kingdom MGE transfers will be addressed. MGE transfers could be ancient events which drove evolution or recurrent events which regulate adaptations. Recent findings on MGE transfers between plant and its interacting microbes suggest gene regulations brought forth by such transfers for symbiosis or defense mechanisms. In the natural environment, factors such as temperature and soil composition constantly influence the interactions among different parties in the rhizosphere. In this review, we will also address the effects of various environmental factors on MGE transfers in the rhizosphere. Besides environmental factors, plant root exudates also play a role in the regulation of MGE transfer among microbes in the rhizosphere. The potential use of microbes and plants for bioremediation will be discussed.

In Silico Identification of the Full Complement of Subtilase-Encoding Genes and Characterization of the Role of TaSBT1.7 in Resistance Against Stripe Rust in Wheat

Plant subtilases (SBTs) or subtilisin-like proteases comprise a very diverse family of serine peptidases that participates in a broad spectrum of biological functions. Despite increasing evidence for roles of SBTs in plant immunity in recent years, little is known about wheat (Triticum aestivum) SBTs (TaSBTs). Here, we identified 255 TaSBT genes from bread wheat using the latest version 2.0 of the reference genome sequence. The SBT family can be grouped into five clades, from TaSBT1 to TaSBT5, based on a phylogenetic tree constructed with deduced protein sequences. In silico protein-domain analysis revealed the existence of considerable sequence diversification of the TaSBT family which, together with the local clustered gene distribution, suggests that TaSBT genes have undergone extensive functional diversification. Among those TaSBT genes whose expression was altered by biotic factors, TaSBT1.7 was found to be induced in wheat leaves by chitin and flg22 elicitors, as well as six examined pathogens, implying a role for TaSBT1.7 in plant defense. Transient overexpression of TaSBT1.7 in Nicotiana benthamiana leaves resulted in necrotic cell death. Moreover, knocking down TaSBT1.7 in wheat using barley stripe mosaic virus-induced gene silencing compromised the hypersensitive response and resistance against Puccinia striiformis f. sp. tritici, the causal agent of wheat stripe rust. Taken together, this study defined the full complement of wheat SBT genes and provided evidence for a positive role of one particular member, TaSBT1.7, in the incompatible interaction between wheat and a stripe rust pathogen.

The semi-dwarfing gene Rht-dp from dwarf polish wheat (Triticum polonicum L.) is the "Green Revolution" gene Rht-B1b

Background

The wheat dwarfing gene increases lodging resistance, the grain number per spike and harvest index. Dwarf Polish wheat (Triticum polonicum L., 2n = 4x = 28, AABB, DPW), initially collected from Tulufan, Xinjiang, China, carries a semi-dwarfing gene Rht-dp on chromosome 4BS. However, Rht-dp and its dwarfing mechanism are unknown.

Results

Homologous cloning and mapping revealed that Rht-dp is the ‘Green Revolution’ gene Rht-B1b. A haplotype analysis in 59 tetraploid wheat accessions showed that Rht-B1b was only present in T. polonicum. Transcriptomic analysis of two pairs of near-isogenic lines (NILs) of DPW × Tall Polish wheat (Triticum polonicum L., 2n = 4x = 28, AABB, TPW) revealed 41 differentially expressed genes (DEGs) as potential dwarfism-related genes. Among them, 28 functionally annotated DEGs were classed into five sub-groups: hormone-related signalling transduction genes, transcription factor genes, cell wall structure-related genes, reactive oxygen-related genes, and nitrogen regulation-related genes.

Conclusions

These results indicated that Rht-dp is Rht-B1b, which regulates pathways related to hormones, reactive oxygen species, and nitrogen assimilation to modify the cell wall structure, and then limits cell wall loosening and inhibits cell elongation, thereby causing dwarfism in DPW.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07367-x.

System-wide characterization of subtilases reveals that subtilisin-like protease FgPrb1 of Fusarium graminearum regulates fungal development and virulence

Subtilases represent the second largest subfamily of serine proteases, and are important for various biological processes. However, the biological function of subtilases has not been systematically characterized in plant pathogens. In present study, 32 subtilases were identified in the genome of wheat scab fungus Fusarium graminearum, a devastating cereal plant pathogen. Deletion mutants of each subtilase were obtained and functionally characterized. Among them, the deletion of FgPrb1 resulted in greatly reduced virulence of F. graminearum. The regulatory mechanisms of FgPrb1 in virulence were investigated in details. Our results showed that the loss of FgPrb1 led to defects in deoxynivalenol (DON) production, responses to environmental stimuli, and lipid metabolism. Additionally, we found that FgPrb1 was involved in autophagy regulation. Taken together, the systematic functional characterization of subtilases showed that the FgPrb1 of F. graminearum is critical for plant infection by regulating multiple different cellular processes.

The fungal subtilase AsES elicits a PTI-like defence response in Arabidopsis thaliana plants independently of its enzymatic activity

Acremonium strictum elicitor subtilisin (AsES) is a 34-kDa serine-protease secreted by the strawberry fungal pathogen A. strictum. On AsES perception, a set of defence reactions is induced, both locally and systemically, in a wide variety of plant species and against pathogens of alternative lifestyles. However, it is not clear whether AsES proteolytic activity is required for triggering a defence response or if the protein itself acts as an elicitor. To investigate the necessity of the protease activity to activate the defence response, AsES coding sequences of the wild-type gene and a mutant on the active site (S226A) were cloned and expressed in Escherichia coli. Our data show that pretreatment of Arabidopsis plants with inactive proteins, i.e. inhibited with phenylmethylsulphonyl fluoride (PMSF) and mutant, resulted in an increased systemic resistance to Botrytis cinerea and expression of defence-related genes in a temporal manner that mimics the effect already reported for the native AsES protein. The data presented in this study indicate that the defence-eliciting property exhibited by AsES is not associated with its proteolytic activity. Moreover, the enhanced expression of some immune marker genes, seedling growth inhibition and the involvement of the co-receptor BAK1 observed in plants treated with AsES suggests that AsES is being recognized as a pathogen-associated molecular pattern by a leucine-rich repeat receptor. The understanding of the mechanism of action of AsES will contribute to the development of new breeding strategies to confer durable resistance in plants.

Cell Wall Proteome of Wheat Grain Endosperm and Outer Layers at Two Key Stages of Early Development

The cell wall is an important compartment in grain cells that fulfills both structural and functional roles. It has a dynamic structure that is constantly modified during development and in response to biotic and abiotic stresses. Non-structural cell wall proteins (CWPs) are key players in the remodeling of the cell wall during events that punctuate the plant life. Here, a subcellular and quantitative proteomic approach was carried out to identify CWPs possibly involved in changes in cell wall metabolism at two key stages of wheat grain development: the end of the cellularization step and the beginning of storage accumulation. Endosperm and outer layers of wheat grain were analyzed separately as they have different origins (maternal and seed) and functions in grains. Altogether, 734 proteins with predicted signal peptides were identified (CWPs). Functional annotation of CWPs pointed out a large number of proteins potentially involved in cell wall polysaccharide remodeling. In the grain outer layers, numerous proteins involved in cutin formation or lignin polymerization were found, while an unexpected abundance of proteins annotated as plant invertase/pectin methyl esterase inhibitors were identified in the endosperm. In addition, numerous CWPs were accumulating in the endosperm at the grain filling stage, thus revealing strong metabolic activities in the cell wall during endosperm cell differentiation, while protein accumulation was more intense at the earlier stage of development in outer layers. Altogether, our work gives important information on cell wall metabolism during early grain development in both parts of the grain, namely the endosperm and outer layers. The wheat cell wall proteome is the largest cell wall proteome of a monocot species found so far.

Involvement of Populus CLEL peptides in root development

As one of the major groups of small post-translationally modified peptides, the CLV3/EMBRYO SURROUNDING REGION-RELATED (CLE)-like (CLEL) peptide family has been reported to regulate root growth, lateral root development and plant gravitropic responses in Arabidopsis thaliana. In this study, we identified 12 CLEL genes in Populus trichocarpa and performed a comprehensive bioinformatics analysis on these genes. Among them, five P. trichocarpa CLELs (PtrCLELs) were revised with new gene models. All of these PtrCLEL proteins were structurally similar to the A. thaliana CLELs (AtCLELs), including an N-terminal signal peptide, a conserved C-terminal 13-amino-acid CLEL motif and a variable intermediate region. In silico and quantitative real-time PCR analyses showed that PtrCLELs were widely expressed in various tissues, including roots, leaves, buds and stems. Exogenous application of chemically synthesized PtrCLEL peptides resulted in wavy or curly roots and reduced lateral root formation in A. thaliana. Moreover, germinating Populus deltoides seedlings on a growth medium containing these peptides caused the roots to thicken and to form abnormal lateral roots, in many cases in clusters. Anatomical and histological changes in thickened roots were further investigated by treating Populus 717 cuttings with the PtrCLEL10 peptide. We observed that root thickening was mainly due to an increased number of cells in the epidermis, hypodermis and cortex. The results of our study suggested that PtrCLEL and AtCLEL genes encode proteins with similar protein structures, sequences of peptide motif and peptide activities on developing roots. The activities of PtrCLEL peptides in root development were species-dependent.

Genome wide association study reveals novel QTL for barley yellow dwarf virus resistance in wheat

Background

Barley yellow dwarf (BYD) is an important virus disease that causes significant reductions in wheat yield. For effective control of Barley yellow dwarf virus through breeding, the identification of genetic sources of resistance is key to success. In this study, 335 geographically diverse wheat accessions genotyped using an Illumina iSelect 90 K single nucleotide polymorphisms (SNPs) bead chip array were used to identify new sources of resistance to BYD in different environments.

Results

A genome-wide association study (GWAS) performed using all the generalised and mixed linkage models (GLM and MLM, respectively) identified a total of 36 significant marker-trait associations, four of which were consistently detected in the K model. These four novel quantitative trait loci (QTL) were identified on chromosomes 2A, 2B, 6A and 7A and associated with markers IWA3520, IWB24938, WB69770 and IWB57703, respectively. These four QTL showed an additive effect with the average visual symptom score of the lines containing resistance alleles of all four QTL being much lower than those with less favorable alleles. Several Chinese landraces, such as H-205 (Baimazha) and H-014 (Dahongmai) which have all four favorable alleles, showed consistently higher resistance in different field trials. None of them contained the previously described Bdv2, Bdv3 or Bdv4 genes for BYD resistance.

Conclusions

This study identified multiple novel QTL for BYD resistance and some resistant wheat genotypes. These will be useful for breeders to generate combinations with and/or without Bdv2 to achieve higher levels and more stable BYD resistance.

Molecular evidence for origin, diversification and ancient gene duplication of plant subtilases (SBTs)

Plant subtilases (SBTs) are a widely distributed family of serine proteases which participates in plant developmental processes and immune responses. Although SBTs are divided into seven subgroups in plants, their origin and evolution, particularly in green algae remain elusive. Here, we present a comprehensive large-scale evolutionary analysis of all subtilases. The plant subtilases SBT1-5 were found to be monophyletic, nested within a larger radiation of bacteria suggesting that they originated from bacteria by a single horizontal gene transfer (HGT) event. A group of bacterial subtilases comprising representatives from four phyla was identified as a sister group to SBT1-5. The phylogenetic analyses, based on evaluation of novel streptophyte algal genomes, suggested that the recipient of the HGT of bacterial subtilases was the common ancestor of Coleochaetophyceae, Zygnematophyceae and embryophytes. Following the HGT, the subtilase gene duplicated in the common ancestor and the two genes diversified into SBT2 and SBT1, 3–5 respectively. Comparative structural analysis of homology-modeled SBT2 proteins also showed their conservation from bacteria to embryophytes. Our study provides the first molecular evidence about the evolution of plant subtilases via HGT followed by a first gene duplication in the common ancestor of Coleochaetophyceae, Zygnematophyceae, and embryophytes, and subsequent expansion in embryophytes.

How Quercus ilex L. saplings face combined salt and ozone stress: a transcriptome analysis

BACKGROUND

Similar to other urban trees, holm oaks (Quercus ilex L.) provide a physiological, ecological and social service in the urban environment, since they remove atmospheric pollution. However, the urban environment has several abiotic factors that negatively influence plant life, which are further exacerbated due to climate change, especially in the Mediterranean area. Among these abiotic factors, increased uptake of Na + and Cl - usually occurs in trees in the urban ecosystem; moreover, an excess of the tropospheric ozone concentration in Mediterranean cities further affects plant growth and survival. Here, we produced and annotated a de novo leaf transcriptome of Q. ilex as well as transcripts over- or under-expressed after a single episode of O3 (80 nl l-1, 5 h), a salt treatment (150 mM for 15 days) or a combination of these treatments, mimicking a situation that plants commonly face, especially in urban environments.

RESULTS

Salinity dramatically changed the profile of expressed transcripts, while the short O3 pulse had less effect on the transcript profile. However, the short O3 pulse had a very strong effect in inducing over- or under-expression of some genes in plants coping with soil salinity. Many differentially regulated genes were related to stress sensing and signalling, cell wall remodelling, ROS sensing and scavenging, photosynthesis and to sugar and lipid metabolism. Most differentially expressed transcripts revealed here are in accordance with a previous report on Q. ilex at the physiological and biochemical levels, even though the expression profiles were overall more striking than those found at the biochemical and physiological levels.

CONCLUSIONS

We produced for the first time a reference transcriptome for Q. ilex, and performed gene expression analysis for this species when subjected to salt, ozone and a combination of the two. The comparison of gene expression between the combined salt + ozone treatment and salt or ozone alone showed that even though many differentially expressed genes overlap all treatments, combined stress triggered a unique response in terms of gene expression modification. The obtained results represent a useful tool for studies aiming to investigate the effects of environmental stresses in urban-adapted tree species.

Host Cell Proteome of Physcomitrella patens Harbors Proteases and Protease Inhibitors under Bioproduction Conditions

Host cell proteins are inevitable contaminants of biopharmaceuticals. Here, we performed detailed analyses of the host cell proteome of moss ( Physcomitrella patens) bioreactor supernatants using mass spectrometry and subsequent bioinformatics analysis. Distinguishing between the apparent secretome and intracellular contaminants, a complex extracellular proteolytic network including subtilisin-like proteases, metallo-proteases, and aspartic proteases was identified. Knockout of a subtilisin-like protease affected the overall extracellular proteolytic activity. Besides proteases, also secreted protease-inhibiting proteins such as serpins were identified. Further, we confirmed predicted cleavage sites of 40 endogenous signal peptides employing an N-terminomics approach. The present data provide novel aspects to optimize both product stability of recombinant biopharmaceuticals as well as their maturation along the secretory pathway. Data are available via ProteomeXchange with identifier PXD009517.

Mechanism of Salt-Induced Self-Compatibility Dissected by Comparative Proteomic Analysis in Brassica napus L

Self-incompatibility (SI) in plants genetically prevents self-fertilization to promote outcrossing and genetic diversity. Its hybrids in Brassica have been widely cultivated due to the propagation of SI lines by spraying a salt solution. We demonstrated that suppression of Brassica napus SI from edible salt solution treatment was ascribed to sodium chloride and independent of S haplotypes, but it did not obviously change the expression of SI-related genes. Using the isobaric tags for relative and absolute quantitation (iTRAQ) technique, we identified 885 differentially accumulated proteins (DAPs) in Brassica napus stigmas of un-pollinated (UP), pollinated with compatible pollen (PC), pollinated with incompatible pollen (PI), and pollinated with incompatible pollen after edible salt solution treatment (NA). Of the 307 DAPs in NA/UP, 134 were unique and 94 were shared only with PC/UP. In PC and NA, some salt stress protein species, such as glyoxalase I, were induced, and these protein species were likely to participate in the self-compatibility (SC) pathway. Most of the identified protein species were related to metabolic pathways, biosynthesis of secondary metabolites, ribosome, and so on. A systematic analysis implied that salt treatment-overcoming SI in B.napus was likely conferred by at least five different physiological mechanisms: (i) the use of Ca2+ as signal molecule; (ii) loosening of the cell wall to allow pollen tube penetration; (iii) synthesis of compatibility factor protein species for pollen tube growth; (iv) depolymerization of microtubule networks to facilitate pollen tube movement; and (v) inhibition of protein degradation pathways to restrain the SI response.

Proteases catalyzing vicilin cleavage in developing pea (Pisum sativum L.) seeds

Legume species differ in whether or not the 7S globulins stored in seeds undergo proteolytic processing during seed development, while preserving the bicupin structure and trimeric assembly necessary for accumulation and packing into protein storage vacuoles. Two such cleavage sites have been documented for the vicilins in pea cotyledons: one in the linker region between the two cupin domains, and another in an exposed loop in the C-terminal cupin. In this report, we explain the occurrence of vicilin cleavage in developing pea by showing that the storage vacuoles are already acidified before germination, in contrast to soybean and peanut where acidification occurs only after germination. We also show that the two cleavage reactions are catalyzed by two different proteases. The vicilin cleavage at the linker region was inhibited by AEBSF (4-(2-aminoethyl)benzenesulfonyl fluoride), indicative of a serine protease. The cleavage in the C-terminal cupin domain was sensitive to the sulfhydryl-reactive reagents p-chloromercuriphenylsulfonate and iodoacetate, but not to E-64 (N-[N-(L-3-transcarboxyirane-2-carbonyl)-l-leucyl]-agmatine), characteristic of the legumain class of cysteine proteases. During seed development, we found the predominant vicilin cleavage in this pea cultivar (Knight) to be at the site in the second cupin domain; but after germination, both sites were cleaved at about the same rate.

From structure to function - a family portrait of plant subtilases

Contents Summary 901 I. Introduction 901 II. Biochemistry and structure of plant SBTs 902 III. Phylogeny of plant SBTs and family organization 903 IV. Physiological roles of plant SBTs 905 V. Conclusions and outlook 911 Acknowledgements 912 References 912 SUMMARY: Subtilases (SBTs) are serine peptidases that are found in all three domains of life. As compared with homologs in other Eucarya, plant SBTs are more closely related to archaeal and bacterial SBTs, with which they share many biochemical and structural features. However, in the course of evolution, functional diversification led to the acquisition of novel, plant-specific functions, resulting in the present-day complexity of the plant SBT family. SBTs are much more numerous in plants than in any other organism, and include enzymes involved in general proteolysis as well as highly specific processing proteases. Most SBTs are targeted to the cell wall, where they contribute to the control of growth and development by regulating the properties of the cell wall and the activity of extracellular signaling molecules. Plant SBTs affect all stages of the life cycle as they contribute to embryogenesis, seed development and germination, cuticle formation and epidermal patterning, vascular development, programmed cell death, organ abscission, senescence, and plant responses to their biotic and abiotic environments. In this article we provide a comprehensive picture of SBT structure and function in plants.

Proteomic approach to understand the molecular physiology of symbiotic interaction between Piriformospora indica and Brassica napus

Many studies have been now focused on the promising approach of fungal endophytes to protect the plant from nutrient deficiency and environmental stresses along with better development and productivity. Quantitative and qualitative protein characteristics are regulated at genomic, transcriptomic, and posttranscriptional levels. Here, we used integrated in-depth proteome analyses to characterize the relationship between endophyte Piriformospora indica and Brassica napus plant highlighting its potential involvement in symbiosis and overall growth and development of the plant. An LC-MS/MS based label-free quantitative technique was used to evaluate the differential proteomics under P. indica treatment vs. control plants. In this study, 8,123 proteins were assessed, of which 46 showed significant abundance (34 downregulated and 12 upregulated) under high confidence conditions (p-value ≤ 0.05, fold change ≥2, confidence level 95%). Mapping of identified differentially expressed proteins with bioinformatics tools such as GO and KEGG pathway analysis showed significant enrichment of gene sets involves in metabolic processes, symbiotic signaling, stress/defense responses, energy production, nutrient acquisition, biosynthesis of essential metabolites. These proteins are responsible for root's architectural modification, cell remodeling, and cellular homeostasis during the symbiotic growth phase of plant's life. We tried to enhance our knowledge that how the biological pathways modulate during symbiosis?

Subtilisin-like proteases in plant defence: the past, the present and beyond

Subtilisin-like proteases (or subtilases) are a very diverse family of serine peptidases present in many organisms, but mostly in plants. With a broad spectrum of biological functions, ranging from protein turnover and plant development to interactions with the environment, subtilases have been gaining increasing attention with regard to their involvement in plant defence responses against the most diverse pathogens. Over the last 5 years, the number of published studies associating plant subtilases with pathogen resistance and plant immunity has increased tremendously. In addition, the observation of subtilases and serine protease inhibitors secreted by pathogens has also gained prominence. In this review, we focus on the active participation of subtilases in the interactions established by plants with the environment, highlighting their role in plant-pathogen communication.

Arabidopsis thaliana phytaspase: identification and peculiar properties

Phytaspases are plant cell death-related proteases of the subtilisin-like protease family that possess an unusual aspartate cleavage specificity. Although phytaspase activity is widespread in plants, phytaspase of Arabidopsis thaliana (L.) Heynh. has escaped detection and identification thus far. Here, we show that a single gene (At4 g10540) out of 56 A. thaliana subtilisin-like protease genes encodes a phytaspase. The recombinant phytaspase was overproduced in Nicotiana benthamiana Domin leaves, isolated, and its substrate specificity and properties were characterised. At pH 5.5, at physiological mildly acidic reaction conditions, the Arabidopsis phytaspase was shown to be strictly Asp-specific. The strongly preferred cleavage motifs of the enzyme out of a panel of synthetic peptide substrates were YVAD and IETD, while the VEID-based substrate preferred by the tobacco and rice phytaspases was almost completely resistant to hydrolysis. At neutral pH, however, the Arabidopsis phytaspase could hydrolyse peptide substrates after two additional amino acid residues, His and Phe, in addition to Asp. This observation may indicate that the repertoire of Arabidopsis phytaspase targets could possibly be regulated by the conditions of the cellular environment. Similar to tobacco and rice phytaspases, the Arabidopsis enzyme was shown to accumulate in the apoplast of epidermal leaf cells. However, in stomatal cells Arabidopsis phytaspase was observed inside the cells, possibly co-localising with vacuole. Our study thus demonstrates that the Arabidopsis phytaspase possesses both important similarities with and distinctions from the already known phytaspases, and is likely to be the most divergent member of the phytaspase family.

Overexpression of a SDD1-Like Gene From Wild Tomato Decreases Stomatal Density and Enhances Dehydration Avoidance in Arabidopsis and Cultivated Tomato

Stomata are microscopic valves formed by two guard cells flanking a pore, which are located on the epidermis of most aerial plant organs and are used for water and gas exchange between the plant and the atmosphere. The number, size and distribution of stomata are set during development in response to changing environmental conditions, allowing plants to minimize the impact of a stressful environment. In Arabidopsis, STOMATAL DENSITY AND DISTRIBUTION 1 (AtSDD1) negatively regulates stomatal density and optimizes transpiration and water use efficiency (WUE). Despite this, little is known about the function of AtSDD1 orthologs in crop species and their wild stress-tolerant relatives. In this study, SDD1-like from the stress-tolerant wild tomato Solanum chilense (SchSDD1-like) was identified through its close sequence relationship with SDD1-like from Solanum lycopersicum and AtSDD1. Both Solanum SDD1-like transcripts accumulated in high levels in young leaves, suggesting that they play a role in early leaf development. Arabidopsis sdd1-3 plants transformed with SchSDD1-like under a constitutive promoter showed a significant reduction in stomatal leaf density compared with untransformed sdd1-3 plants. Additionally, a leaf dehydration shock test demonstrated that the reduction in stomatal abundance of transgenic plants was sufficient to slow down dehydration. Overexpression of SchSDD1-like in cultivated tomato plants decreased the stomatal index and density of the cotyledons and leaves, and resulted in higher dehydration avoidance. Taken together, these results indicate that SchSDD1-like functions in a similar manner to AtSDD1 and suggest that Arabidopsis and tomatoes share this component of the stomatal development pathway that impinges on water status.

Post-translational maturation of IDA, a peptide signal controlling floral organ abscission in Arabidopsis

The abscission of sepals, petals and stamens in Arabidopsis flowers is controlled by a peptide signal called IDA (Inflorescence Deficient in Abscission). IDA belongs to the large group of small post-translationally modified signaling peptides that are synthesized as larger precursors and require proteolytic processing and specific side chain modifications for signal biogenesis. Using tissue-specific expression of proteinase inhibitors as a novel approach for loss-of-function analysis, we recently identified the peptidases responsible for IDA maturation within the large family of subtilisin-like proteinases (subtilases; SBTs). Further biochemical and physiological assays identified three SBTs (AtSBT5.2, AtSBT4.12, AtSBT4.13) that cleave the IDA precursor to generate the N-terminus of the mature peptide. The C-terminal processing enzyme(s) remain(s) to be identified. While proline hydroxylation was suggested as additional post-translational modification required for IDA maturation, hydroxylated and non-hydroxylated IDA peptides were found to be equally active in bioassays for abscission.

Proteolytic activities in cortex of apical parts of Vicia faba ssp. minor seedling roots during kinetin-induced programmed cell death

Programmed cell death (PCD) is a crucial process in plant development. In this paper, proteolytically related aspects of kinetin-induced PCD in cortex cells of Vicia faba ssp. minor seedlings were examined using morphological, fluorometric, spectrophotometric, and fluorescence microscopic analyses. Cell viability estimation after 46 μM kinetin treatment of seedling roots showed that the number of dying cortex cells increased with treatment duration, reaching maximum after 72 h. Weight of the apical root segments increased with time and was about 2.5-fold greater after 96 h, while the protein content remained unchanged, compared to the control. The total and cysteine-dependent proteolytic activities fluctuated during 1-96-h treatment, which was not accompanied by the changes in the protein amount, indicating that the absolute protein amounts decreased during kinetin-induced PCD. N-ethylmaleimide (NEM), phenylmethylsulfonyl fluoride (PMSF), and Z-Leu-Leu-Nva-H (MG115), the respective cysteine, serine, and proteasome inhibitors, suppressed kinetin-induced PCD. PMSF significantly decreased serine-dependent proteolytic activities without changing the amount of proteins, unlike NEM and MG115. More pronounced effect of PMSF over NEM indicated that in the root apical segments, the most important proteolytic activity during kinetin-induced PCD was that of serine proteases, while that of cysteine proteases may be important for protein degradation in the last phase of the process. Both NEM and PMSF inhibited apoptotic-like structure formation during kinetin-induced PCD. The level of caspase-3-like activity of β1 proteasome subunit increased after kinetin treatment. Addition of proteasome inhibitor MG-115 reduced the number of dying cells, suggesting that proteasomes might play an important role during kinetin-induced PCD.

Diverse Peptide Hormones Affecting Root Growth Identified in the Medicago truncatula Secreted Peptidome

Multigene families encoding diverse secreted peptide hormones play important roles in plant development. A need exists to efficiently elucidate the structures and post-translational-modifications of these difficult-to-isolate peptide hormones in planta so that their biological functions can be determined. A mass spectrometry and bioinformatics approach was developed to comprehensively analyze the secreted peptidome of Medicago hairy root cultures and xylem sap. We identified 759 spectra corresponding to the secreted products of twelve peptide hormones including four CEP (C-TERMINALLY ENCODED PEPTIDE), two CLE (CLV3/ENDOSPERM SURROUNDING REGION RELATED) and six XAP (XYLEM SAP ASSOCIATED PEPTIDE) peptides. The MtCEP1, MtCEP2, MtCEP5 and MtCEP8 peptides identified differed in post-translational-modifications. Most were hydroxylated at conserved proline residues but some MtCEP1 derivatives were tri-arabinosylated. In addition, many CEP peptides possessed unexpected N- and C-terminal extensions. The pattern of these extensions suggested roles for endo- and exoproteases in CEP peptide maturation. Longer than expected, hydroxylated and homogeneously modified mono- and tri-arabinosylated CEP peptides corresponding to their in vivo structures were chemically synthesized to probe the effect of these post-translational-modifications on function. The ability of CEP peptides to elevate root nodule number was increased by hydroxylation at key positions. MtCEP1 peptides with N-terminal extensions or with tri-arabinosylation modification, however, were unable to impart increased nodulation. The MtCLE5 and MtCLE17 peptides identified were of precise size, and inhibited main root growth and increased lateral root number. Six XAP peptides, each beginning with a conserved DY sulfation motif, were identified including MtXAP1a, MtXAP1b, MtXAP1c, MtXAP3, MtXAP5 and MtXAP7. MtXAP1a and MtXAP5 inhibited lateral root emergence. Transcriptional analyses demonstrated peptide hormone gene expression in the root vasculature and tip. Since hairy roots can be induced on many plants, their corresponding root cultures may represent ideal source materials to efficiently identify diverse peptide hormones in vivo in a broad range of species.

Identification of Proteases and Protease Inhibitors in Allergenic and Non-Allergenic Pollen

Pollen is one of the most common causes of allergy worldwide, making the study of their molecular composition crucial for the advancement of allergy research. Despite substantial efforts in this field, it is not yet clear why some plant pollens strongly provoke allergies while others do not. However, proteases and protease inhibitors from allergen sources are known to play an important role in the development of pollen allergies. In this study, we aim to uncover differences in the transcriptional pattern of proteases and protease inhibitors in Betula verrucosa and Pinus sylvestris pollen as models for high and low allergenic potential, respectively. We applied RNA sequencing to Betula verrucosa and Pinus sylvestris pollen. After de-novo assembly we derived general functional profiles of the protein coding transcripts. By utilization of domain based functional annotation we identified potential proteases and protease inhibitors and compared their expression in the two types of pollen. Functional profiles are highly similar between Betula verrucosa and Pinus sylvestris pollen. Both pollen contain proteases and inhibitors from 53 and 7 Pfam families, respectively. Some of the members comprised within those families are implicated in facilitating allergen entry, while others are known allergens themselves. Our work revealed several candidate proteins which, with further investigation, represent exciting new leads in elucidating the process behind allergic sensitization.

Evolutionary History of Subtilases in Land Plants and Their Involvement in Symbiotic Interactions

Subtilases, a family of proteases involved in a variety of developmental processes in land plants, are also involved in both mutualistic symbiosis and host-pathogen interactions in different angiosperm lineages. We examined the evolutionary history of subtilase genes across land plants through a phylogenetic analysis integrating amino acid sequence data from full genomes, transcriptomes, and characterized subtilases of 341 species of diverse green algae and land plants along with subtilases from 12 species of other eukaryotes, archaea, and bacteria. Our analysis reconstructs the subtilase gene phylogeny and identifies 11 new gene lineages, six of which have no previously characterized members. Two large, previously unnamed, subtilase gene lineages that diverged before the origin of angiosperms accounted for the majority of subtilases shown to be associated with symbiotic interactions. These lineages expanded through both whole-genome and tandem duplication, with differential neofunctionalization and subfunctionalization creating paralogs associated with different symbioses, including nodulation with nitrogen-fixing bacteria, arbuscular mycorrhizae, and pathogenesis in different plant clades. This study demonstrates for the first time that a key gene family involved in plant-microbe interactions proliferated in size and functional diversity before the explosive radiation of angiosperms.

Identification and expression analysis of 11 subtilase genes during natural and induced senescence of barley plants

Subtilases are one of the largest groups of the serine protease family and are involved in many aspects of plant development including senescence. In wheat, previous reports demonstrate an active participation of two senescence-induced subtilases, denominated P1 and P2, in nitrogen remobilization during whole plant senescence. The aim of the present study was to examine the participation of subtilases in senescence-associated proteolysis of barley leaves while comparing different senescence types. With this purpose, subtilase enzymatic activity, immunodetection with a heterologous antiserum and gene expression of 11 subtilase sequences identified in barley databases by homology to P1 were analyzed in barley leaves undergoing dark-induced or natural senescence at the vegetative or reproductive growth phase. Results showed that subtilase specific activity as well as two inmunoreactive bands representing putative subtilases increased in barley leaves submitted to natural and dark-induced senescence. Gene expression analysis showed that two of the eleven subtilase genes analyzed, HvSBT3 and HvSBT6, were up-regulated in all the senescence conditions tested while HvSBT2 was expressed and up-regulated only during dark-induced senescence. On the other hand, HvSBT1, HvSBT4 and HvSBT7 were down-regulated during senescence and two other subtilase genes (HvSBT10 and HvSBT11) showed no significant changes. The remaining subtilase genes were not detected. Results demonstrate an active participation of subtilases in protein degradation during dark-induced and natural leaf senescence of barley plants both at the vegetative and reproductive stage, and, based on their expression profile, postulate HvSBT3 and HvSBT6 as key components of senescence-associated proteolysis.

A novel subtilase inhibitor in plants shows structural and functional similarities to protease propeptides

The propeptides of subtilisin-like serine proteinases (subtilases, SBTs) serve dual functions as intramolecular chaperones that are required for enzyme folding and as inhibitors of the mature proteases. SBT propeptides are homologous to the I9 family of protease inhibitors that have only been described in fungi. Here we report the identification and characterization of subtilisin propeptide-like inhibitor 1 (SPI-1) from Arabidopsis thaliana Sequence similarity and the shared β-α-β-β-α-β core structure identified SPI-1 as a member of the I9 inhibitor family and as the first independent I9 inhibitor in higher eukaryotes. SPI-1 was characterized as a high-affinity, tight-binding inhibitor of Arabidopsis subtilase SBT4.13 with K_d and K_i values in the picomolar range. SPI-1 acted as a stable inhibitor of SBT4.13 over the physiologically relevant range of pH, and its inhibitory profile included many other SBTs from plants but not bovine chymotrypsin or bacterial subtilisin A. Upon binding to SBT4.13, the C-terminal extension of SPI-1 was proteolytically cleaved. The last four amino acids at the newly formed C terminus of SPI-1 matched both the cleavage specificity of SBT4.13 and the consensus sequence of Arabidopsis SBTs at the junction of the propeptide with the catalytic domain. The data suggest that the C terminus of SPI-1 acts as a competitive inhibitor of target proteases as it remains bound to the active site in a product-like manner. SPI-1 thus resembles SBT propeptides with respect to its mode of protease inhibition. However, in contrast to SBT propeptides, SPI-1 could not substitute as a folding assistant for SBT4.13.

Precursor processing for plant peptide hormone maturation by subtilisin-like serine proteinases

Peptide hormones that regulate plant growth and development are derived from larger precursor proteins by proteolytic processing. Our study addressed the role of subtilisin-like proteinases (SBTs) in this process. Using tissue-specific expression of proteinase inhibitors as a tool to overcome functional redundancy, we found that SBT activity was required for the maturation of IDA (INFLORESCENCE DEFICIENT IN ABSCISSION), a peptide signal for the abscission of floral organs in Arabidopsis We identified three SBTs that process the IDA precursor in vitro, and this processing was shown to be required for the formation of mIDA (the mature and bioactive form of IDA) as the endogenous signaling peptide in vivo. Hence, SBTs act as prohormone convertases in plants, and several functionally redundant SBTs contribute to signal biogenesis.

Revisiting Vitis vinifera Subtilase Gene Family: A Possible Role in Grapevine Resistance against Plasmopara viticola

Subtilisin-like proteases, also known as subtilases, are a very diverse family of serine peptidases present in many organisms. In grapevine, there are hints of the involvement of subtilases in defense mechanisms, but their role is not yet understood. The first characterization of the subtilase gene family was performed in 2014. However, simultaneously, the grapevine genome was re-annotated and several sequences were re-annotated or retrieved. We have performed a re-characterization of this family in grapevine and identified 82 genes coding for 97 putative proteins, as result of alternative splicing. All the subtilases identified present the characteristic S8 peptidase domain and the majority of them also have a pro-domain I9 inhibitor, a protease-associated (PA) domain, and a signal peptide for targeting to the secretory pathway. Phylogenetic studies revealed six subtilase groups denominated VvSBT1 to VvSBT6. As several evidences have highlighted the participation of plant subtilases in response to biotic stimulus, we have investigated subtilase participation in grapevine resistance to Plasmopara viticola, the causative agent of downy mildew. Fourteen grapevine subtilases presenting either high homology to P69C from tomato, SBT3.3 from Arabidopsis thaliana or located near the Resistance to P. viticola (RPV) locus were selected. Expression studies were conducted in the grapevine-P. viticola pathosystem with resistant and susceptible cultivars. Our results may indicate that some of grapevine subtilisins are potentially participating in the defense response against this biotrophic oomycete.

Plant Proteases Involved in Regulated Cell Death

Each plant genome encodes hundreds of proteolytic enzymes. These enzymes can be divided into five distinct classes: cysteine-, serine-, aspartic-, threonine-, and metalloproteinases. Despite the differences in their structural properties and activities, members of all of these classes in plants are involved in the processes of regulated cell death - a basic feature of eukaryotic organisms. Regulated cell death in plants is an indispensable mechanism supporting plant development, survival, stress responses, and defense against pathogens. This review summarizes recent advances in studies of plant proteolytic enzymes functioning in the initiation and execution of distinct types of regulated cell death.

A non canonical subtilase attenuates the transcriptional activation of defence responses in Arabidopsis thaliana

Proteases play crucial physiological functions in all organisms by controlling the lifetime of proteins. Here, we identified an atypical protease of the subtilase family [SBT5.2(b)] that attenuates the transcriptional activation of plant defence independently of its protease activity. The SBT5.2 gene produces two distinct transcripts encoding a canonical secreted subtilase [SBT5.2(a)] and an intracellular protein [SBT5.2(b)]. Concomitant to SBT5.2(a) downregulation, SBT5.2(b) expression is induced after bacterial inoculation. SBT5.2(b) localizes to endosomes where it interacts with and retains the defence-related transcription factor MYB30. Nuclear exclusion of MYB30 results in its reduced transcriptional activation and, thus, suppressed resistance. sbt5.2 mutants, with abolished SBT5.2(a) and SBT5.2(b) expression, display enhanced defence that is suppressed in a myb30 mutant background. Moreover, overexpression of SBT5.2(b), but not SBT5.2(a), in sbt5.2 plants reverts the phenotypes displayed by sbt5.2 mutants. Overall, we uncover a regulatory mode of the transcriptional activation of defence responses previously undescribed in eukaryotes.

A non canonical subtilase attenuates the transcriptional activation of defence responses in Arabidopsis thaliana

Proteases play crucial physiological functions in all organisms by controlling the lifetime of proteins. Here, we identified an atypical protease of the subtilase family [SBT5.2(b)] that attenuates the transcriptional activation of plant defence independently of its protease activity. The SBT5.2 gene produces two distinct transcripts encoding a canonical secreted subtilase [SBT5.2(a)] and an intracellular protein [SBT5.2(b)]. Concomitant to SBT5.2(a) downregulation, SBT5.2(b) expression is induced after bacterial inoculation. SBT5.2(b) localizes to endosomes where it interacts with and retains the defence-related transcription factor MYB30. Nuclear exclusion of MYB30 results in its reduced transcriptional activation and, thus, suppressed resistance. sbt5.2 mutants, with abolished SBT5.2(a) and SBT5.2(b) expression, display enhanced defence that is suppressed in a myb30 mutant background. Moreover, overexpression of SBT5.2(b), but not SBT5.2(a), in sbt5.2 plants reverts the phenotypes displayed by sbt5.2 mutants. Overall, we uncover a regulatory mode of the transcriptional activation of defence responses previously undescribed in eukaryotes.

DOI:http://dx.doi.org/10.7554/eLife.19755.001

Like animals, plants have evolved numerous ways to protect themselves from disease. When a plant detects an invading microbe, it massively changes which genes it expresses to establish a defensive response. This is possible thanks to the action of a type of protein, named transcription factors, which are able to bind to DNA in the cell nucleus and regulate gene expression. However, triggering such a response comes at a cost, and so plants must keep their defensive response in check such that they can allocate resources in a balanced way.

In the model plant Arabidopsis, a protein named MYB30 is one transcription factor that is able to promote disease resistance. Previous research identified some proteins that can reduce the activity of this transcription factor to avoid triggering a response when it is not needed, for example, when no infectious microbes are present. However, it was likely that other proteins were also involved in the process.

Now, Serrano et al. report that an enzyme called SBT5.2 is an additional negative regulator of MYB30 activity. SBT5.2 belongs to a family of protein-degrading enzymes called subtilases, which are typically localized outside cells. As such, it was unclear how SBT5.2 could interact and regulate a transcription factor that is found inside the nucleus of plant cells.

Nevertheless, Serrano et al. found that the gene that encodes SBT5.2 actually gives rise to two distinct proteins. The first is a classical subtilase that is indeed located outside of the cell, and so cannot interact with MYB30 and does not affect its activity. The second protein is an atypical subtilase that localises to bubble-like compartments called vesicles within the cell and is able to highjack MYB30 on its way to the nucleus. When the atypical subtilase interacts with MYB30 at vesicles, it stops MYB30 from entering the nucleus. As a result, MYB30 cannot bind to the DNA nor activate its target genes. This means that the defensive response that normally depends on MYB30 is weakened.

The work of Serrano et al. uncovers a new way to regulate the expression of defence-related genes. Further unravelling the molecular mechanisms involved in the fine-tuning of gene expression represents a challenging task for future research.

DOI:http://dx.doi.org/10.7554/eLife.19755.002

Proteases from Canavalia ensiformis: Active and Thermostable Enzymes with Potential of Application in Biotechnology

Extracts of leaves, seeds, roots, and stem from a tropical legume, C. ensiformis, were prepared employing buffers and detergent in aqueous solution. Leaf extracts had the highest protein content and the most pronounced peptidase activity with optimal pH in the neutral to alkaline range. All extracts exhibited peaks of activity at various pH values, suggesting the presence of distinctive classes of proteases. N-α-Tosyl-L-arginine methyl ester hydrolysis was maximal at 30°C to 60°C and peptidase activity from all extracts presented very good thermal stability after 24 h incubation at 70°C. C. ensiformis proteases exhibited molecular masses of about 200–57, 40–37, and 20–15 kDa by SDS-PAGE analysis. These enzymes cleaved hemoglobin, bovine serum albumin, casein, and gelatin at different levels. Serine and metalloproteases are the major proteases in C. ensiformis extracts, modulated by divalent cations, stable at 1% of surfactant Triton X-100 and at different concentrations of the reducing agent β-mercaptoethanol. Thus, C. ensiformis expresses a particular set of proteases in distinctive organs with high activity and stability, making this legume an important source of proteases with biotechnological potential.

The SBT6.1 subtilase processes the GOLVEN1 peptide controlling cell elongation

The GOLVEN (GLV) gene family encode small secreted peptides involved in important plant developmental programs. Little is known about the factors required for the production of the mature bioactive GLV peptides. Through a genetic suppressor screen in Arabidopsis thaliana, two related subtilase genes, AtSBT6.1 and AtSBT6.2, were identified that are necessary for GLV1 activity. Root and hypocotyl GLV1 overexpression phenotypes were suppressed by mutations in either of the subtilase genes. Synthetic GLV-derived peptides were cleaved in vitro by the affinity-purified SBT6.1 catalytic enzyme, confirming that the GLV1 precursor is a direct subtilase substrate, and the elimination of the in vitro subtilase recognition sites through alanine substitution suppressed the GLV1 gain-of-function phenotype in vivo Furthermore, the protease inhibitor Serpin1 bound to SBT6.1 and inhibited the cleavage of GLV1 precursors by the protease. GLV1 and its homolog GLV2 were expressed in the outer cell layers of the hypocotyl, preferentially in regions of rapid cell elongation. In agreement with the SBT6 role in GLV precursor processing, both null mutants for sbt6.1 and sbt6.2 and the Serpin1 overexpression plants had shorter hypocotyls. The biosynthesis of the GLV signaling peptides required subtilase activity and might be regulated by specific protease inhibitors. The data fit with a model in which the GLV1 signaling pathway participates in the regulation of hypocotyl cell elongation, is controlled by SBT6 subtilases, and is modulated locally by the Serpin1 protease inhibitor.