Genome-wide comparative analysis of papain-like cysteine protease family genes in castor bean and physic nut

Analysis of Carica papaya Informs Lineage-Specific Evolution of the Aquaporin (AQP) Family in Brassicales

Aquaporins (AQPs), a type of intrinsic membrane proteins that transport water and small solutes across biological membranes, play crucial roles in plant growth and development. This study presents a first genome-wide identification and comparative analysis of the AQP gene family in papaya (Carica papaya L.), an economically and nutritionally important fruit tree of tropical and subtropical regions. A total of 29 CpAQP genes were identified, which represent five subfamilies, i.e., nine plasma intrinsic membrane proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). Although the family is smaller than the 35 members reported in Arabidopsis, it is highly diverse, and the presence of CpXIP genes as well as orthologs in Moringa oleifera and Bretschneidera sinensis implies that the complete loss of the XIP subfamily in Arabidopsis is lineage-specific, sometime after its split with papaya but before Brassicaceae-Cleomaceae divergence. Reciprocal best hit-based sequence comparison of 530 AQPs and synteny analyses revealed that CpAQP genes belong to 29 out of 61 identified orthogroups, and lineage-specific evolution was frequently observed in Brassicales. Significantly, the well-characterized NIP3 group was completely lost; lineage-specific loss of the NIP8 group in Brassicaceae occurred sometime before the divergence with Cleomaceae, and lineage-specific loss of NIP2 and SIP3 groups in Brassicaceae occurred sometime after the split with Cleomaceae. In contrast to a predominant role of recent whole-genome duplications (WGDs) on the family expansion in B. sinensis, Tarenaya hassleriana, and Brassicaceae plants, no recent AQP repeats were identified in papaya, and ancient WGD repeats are mainly confined to the PIP subfamily. Subfamily even group-specific evolution was uncovered via comparing exon-intron structures, conserved motifs, the aromatic/arginine selectivity filter, and gene expression profiles. Moreover, down-regulation during fruit ripening and expression divergence of duplicated CpAQP genes were frequently observed in papaya. These findings will not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further studies of AQP genes in papaya and species beyond.

Identification and functional characterization of safflower cysteine protease 1 as negative regulator in response to low-temperature stress in transgenic Arabidopsis

We performed genome-wide and heterologous expression analysis of the safflower cysteine protease family and found that inhibition of CtCP1 expression enhanced plant cold resistance. Cysteine protease (CP) is mainly involved in plant senescence and stress responses. However, the molecular mechanism of endogenous cysteine protease inhibition in plant stress tolerance is yet unknown. Here, we report the discovery and functional characterization of a candidate CP1 gene from safflower. The conserved structural topology of CtCPs revealed important insights into their possible roles in plant growth and stress responses. The qRT-PCR results implied that most of CtCP genes were highly expressed at fading stage suggesting that they are most likely involved in senescence process. The CtCP1 expression was significantly induced at different time points under cold, NaCl, H₂O₂ and PEG stress, respectively. The in-vitro activity of heterologously expressed CtCP1 protein showed highest protease activity for casein and azocasein substrates. The expression and phenotypic data together with antioxidant activity and physiological indicators revealed that transgenic plants inhibited by CtCP1-anti showed higher tolerance to low temperature than WT and CtCP1-OE plants. Our findings demonstrated the discovery of a new Cysteine protease 1 gene that exerted a detrimental effect on transgenic Arabidopsis under low-temperature stress.

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.

"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.

Genome-Wide Identification and Comparison of Cysteine Proteases in the Pollen Coat and Other Tissues in Maize

Cysteine proteases, belonging to the C1-papain family, play a major role in plant growth and development, senescence, and immunity. There is evidence to suggest that pollen cysteine protease (CP) (ZmCP03) is involved in regulating the anther development and pollen formation in maize. However, there is no report on the genome-wide identification and comparison of CPs in the pollen coat and other tissues in maize. In this study, a total of 38 homologous genes of ZmCP03 in maize were identified. Subsequently, protein motifs, conserved domains, gene structures, and duplication patterns of 39 CPs are analyzed to explore their evolutionary relationship and potential functions. The cis-elements were identified in the upstream sequence of 39 CPs, especially those that are related to regulating growth and development and responding to environmental stresses and hormones. The expression patterns of these genes displayed remarked difference at a tissue or organ level in maize based on the available transcriptome data in the public database. Quantitative reverse transcription PCR (RT-qPCR) analysis showed that ZmCP03 was preferably expressed at a high level in maize pollen. Analyses by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and immunoblot, immunofluorescence and immunogold electron microscopy all validated the cellular localization of ZmCP03 in both the pollen coat and pollen cytoplasm. In addition, 142 CP genes from Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa) and cotton (Gossypium hirsutum), together with 39 maize CPs, were retrieved to analyze their evolution by comparing with orthologous genes. The results suggested that ZmCP03 was relatively conservative and stable during evolution. This study may provide a referential evidence on the function of ZmCP03 in pollen development and germination in maize.

Papain-Like Cysteine Protease Gene Family in Fig (Ficus carica L.): Genome-Wide Analysis and Expression Patterns

The papain-like cysteine proteases (PLCPs) are the most abundant family of cysteine proteases in plants, with essential roles in biotic/abiotic stress responses, growth and senescence. Papain, bromelain and ficin are widely used in food, medicine and other industries. In this study, 31 PLCP genes (FcPCLPs) were identified in the fig (Ficus carica L.) genome by HMM search and manual screening, and assigned to one of nine subfamilies based on gene structure and conserved motifs. SAG12 and RD21 were the largest subfamilies with 10 and 7 members, respectively. The FcPCLPs ranged from 1,128 to 5,075 bp in length, containing 1-10 introns, and the coding sequence ranged from 624 to 1,518 bp, encoding 207-505 amino acids. Subcellular localization analysis indicated that 24, 2, and 5 PLCP proteins were targeted to the lysosome/vacuole, cytoplasm and extracellular matrix, respectively. Promoter (2,000 bp upstream) analysis of FcPLCPs revealed a high number of plant hormone and low temperature response elements. RNA-seq revealed differential expression of 17 FcPLCPs in the inflorescence and receptacle, and RD21 subfamily members were the major PLCPs expressed in the fruit; 16 and 5 FcPLCPs responded significantly to ethylene and light, respectively. Proteome analyses revealed 18 and 5 PLCPs in the fruit cell soluble proteome and fruit latex, respectively. Ficins were the major PLCP in fig fruit, with decreased abundance in inflorescences, but increased abundance in receptacles of commercial-ripe fruit. FcRD21B/C and FcALP1 were aligned as the genes encoding the main ficin isoforms. Our study provides valuable multi-omics information on the FcPLCP family and lays the foundation for further functional studies.

Genes encoding light-harvesting chlorophyll a/b-binding proteins in papaya (Carica papaya L.) and insight into lineage-specific evolution in Brassicaceae

Light-harvesting chlorophyll a/b-binding (Lhc) proteins comprise a plant-specific superfamily involved in photosynthesis and stress responses. Despite their importance, little is known in papaya (Carica papaya), an economically important tree fruit crop as well as a species close to the model plant arabidopsis (Arabidopsis thaliana). This study reports a first genome-wide analysis of Lhc superfamily genes in papaya, and a total of 28 members that represent four defined families or 26 orthologous groups were identified from the papaya genome. The superfamily number is comparable to 28 or 27 reported in castor (Ricinus communis) and jatropha (Jatropha curcas), respectively, two Euphorbiaceous plants also without any recent whole-genome duplication (WGD), but relatively less than 35, 34, 32, 32, 37, 30 or 32 present in cassava (Manihot esculenta), arabidopsis, A. lyrata, A. halleri, Capsella rubella, C. grandiflora, and Eutrema salsugineum, respectively, representative species having experienced one or two recent WGDs. Local duplication was shown to play a predominant role in gene expansion in papaya, castor, and jatropha, which is only confined to the Lhcb1 group. By contrast, WGD plays a relatively more important role in cassava, arabidopsis, and other Brassicaceous plants. Further comparison of Brassicaceous plants revealed that loss of the SEP6 group in arabidopsis is lineage-specific, occurring sometime after papaya-arabidopsis divergence but before the radiation of Brassicaceous plants. Transcriptional profiling revealed a leaf-preferential expression pattern of most CpLhc superfamily genes and their transcript levels were markedly regulated by three abiotic stresses, i.e., mimicking drought, cold, and high salt. These findings not only facilitate further functional studies in papaya, but also improve our knowledge on lineage-specific evolution of this special gene superfamily in Brassicaceae.

Comparative modelling studies of fruit bromelain using molecular dynamics simulation

Fruit bromelain is a cysteine protease accumulated in pineapple fruits. This proteolytic enzyme has received high demand for industrial and therapeutic applications. In this study, fruit bromelain sequences QIM61759, QIM61760 and QIM61761 were retrieved from the National Center for Biotechnology Information (NCBI) Genbank Database. The tertiary structure of fruit bromelain QIM61759, QIM61760 and QIM61761 was generated by using MODELLER. The result revealed that the local stereochemical quality of the generated models was improved by using multiple templates during modelling process. Moreover, by comparing with the available papain model, structural analysis provides an insight on how pro-peptide functions as a scaffold in fruit bromelain folding and contributing to inactivation of mature protein. The structural analysis also disclosed the similarities and differences between these models. Lastly, thermal stability of fruit bromelain was studied. Molecular dynamics simulation of fruit bromelain structures at several selected temperatures demonstrated how fruit bromelain responds to elevation of temperature.

Light-harvesting chlorophyll a/b-binding protein-coding genes in jatropha and the comparison with castor, cassava and arabidopsis

The Lhc (light-harvesting chlorophyll a/b-binding protein) superfamily represents a class of antennae proteins that play indispensable roles in capture of solar energy as well as photoprotection under stress conditions. Despite their importance, little information has been available beyond model plants. In this study, we presents a first genome-wide analysis of Lhc superfamily genes in jatropha (Jatropha curcas L., Euphorbiaceae), an oil-bearing plant for biodiesel purpose. A total of 27 members were identified from the jatropha genome, which were shown to distribute over nine out of the 11 chromosomes. The superfamily number is comparable to 28 present in castor (Ricinus communis, Euphorbiaceae), but relatively less than 35 in cassava (Manihot esculenta, Euphorbiaceae) and 34 in arabidopsis (Arabidopsis thaliana) that experienced one or two recent whole-genome duplications (WGDs), respectively. In contrast to a high number of paralogs present in cassava and arabidopsis, few duplicates were found in jatropha as observed in castor, corresponding to no recent WGD occurred in these two species. Nevertheless, 26 orthologous groups representing four defined families were found in jatropha, and nearly one-to-one orthologous relationship was observed between jatropha and castor. By contrast, a novel group named SEP6 was shown to have been lost in arabidopsis. Global transcriptome profiling revealed a predominant expression pattern of most JcLhc superfamily genes in green tissues, reflecting their key roles in photosynthesis. Moreover, their expression profiles upon hormones, drought, and salt stresses were also investigated. These findings not only improve our knowledge on species-specific evolution of the Lhc supergene family, but also provide valuable information for further studies in jatropha.

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.

Genomics analysis of genes encoding respiratory burst oxidase homologs (RBOHs) in jatropha and the comparison with castor bean

Respiratory burst oxidase homologs (RBOHs), which catalyze the production of superoxide from oxygen and NADPH, play key roles in plant growth and development, hormone signaling, and stress responses. Compared with extensive studies in model plants arabidopsis and rice, little is known about RBOHs in other species. This study presents a genome-wide analysis of Rboh family genes in jatropha (Jatropha curcas) as well as the comparison with castor bean (Ricinus communis), another economically important non-food oilseed crop of the Euphorbiaceae family. The family number of seven members identified from the jatropha genome is equal to that present in castor bean, and further phylogenetic analysis assigned these genes into seven groups named RBOHD, -C, -B, -E, -F, -N, and -H. In contrast to a high number of paralogs present in arabidopsis and rice that experienced several rounds of recent whole-genome duplications, no duplicate was identified in both jatropha and castor bean. Conserved synteny and one-to-one orthologous relationship were observed between jatropha and castor bean Rboh genes. Although exon-intron structures are usually highly conserved between orthologs, loss of certain introns was observed for JcRbohB, JcRbohD, and RcRbohN, supporting their divergence. Global gene expression profiling revealed diverse patterns of JcRbohs over various tissues. Moreover, expression patterns of JcRbohs during flower development as well as various stresses were also investigated. These findings will not only improve our knowledge on species-specific evolution of the Rboh gene family, but also provide valuable information for further functional analysis of Rboh genes in jatropha.

Genome-wide comparison reveals divergence of cassava and rubber aquaporin family genes after the recent whole-genome duplication

Background

Aquaporins (AQPs) are a class of integral membrane proteins that facilitate the passive transport of water and other small solutes across biological membranes. Despite their importance, little information is available in cassava (Manihot esculenta), a perennial shrub of the Euphorbiaceae family that serves the sixth major staple crop in the world.

Results

This study presents a genome-wide analysis of the AQP gene family in cassava. The family of 42 members in this species could be divided into five subfamilies based on phylogenetic analysis, i.e., 14 plasma membrane intrinsic proteins (PIPs), 13 tonoplast intrinsic proteins (TIPs), nine NOD26-like intrinsic proteins (NIPs), four X intrinsic proteins (XIPs), and two small basic intrinsic proteins (SIPs). Best-reciprocal-hit-based sequence comparison and synteny analysis revealed 34 orthologous groups (OGs) present in the Euphorbiaceae ancestor, and nearly one-to-one or two-to-one orthologous relationships were observed between cassava with rubber/physic nut, reflecting the occurrence of one so-called ρ recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber. In contrast to a predominant role of the ρ WGD on family expansion in rubber, cassava AQP duplicates were derived from the WGD as well as local duplication. Species-specific gene loss was also observed in cassava, which includes the entire NIP4 group and/or six OGs. Comparison of conserved motifs and gene expression profiles revealed divergence of paralogs in cassava as observed in rubber.

Conclusions

Our findings will not only improve our knowledge on family evolution in Euphorbiaceae, but also provide valuable information for further functional analysis of AQP genes in cassava and rubber.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5780-4) contains supplementary material, which is available to authorized users.

Genome-wide identification and comparative evolutionary analysis of the Dof transcription factor family in physic nut and castor bean

DNA-binding with one finger (Dof) proteins comprise a plant-specific transcription factor family involved in plant growth, development and stress responses. This study presents a genome-wide comparison of Dof family genes in physic nut (Jatropha curcas) and castor bean (Ricinus communis), two Euphorbiaceae plants that have not experienced any recent whole-genome duplication. A total of 25 or 24 Dof genes were identified from physic nut and castor genomes, respectively, where JcDof genes are distributed across nine out of 11 chromosomes. Phylogenetic analysis assigned these genes into nine groups representing four subfamilies, and 24 orthologous groups were also proposed based on comparison of physic nut, castor, Arabidopsis and rice Dofs. Conserved microsynteny was observed between physic nut and castor Dof-coding scaffolds, which allowed anchoring of 23 RcDof genes to nine physic nut chromosomes. In contrast to how no recent duplicate was present in castor, two tandem duplications and one gene loss were found in the Dof gene family of physic nut. Global transcriptome profiling revealed diverse patterns of Jc/RcDof genes over various tissues, and key Dof genes involved in flower development and stress response were also identified in physic nut. These findings provide valuable information for further studies of Dof genes in physic nut and castor.

Genome-Wide Identification of Papain-Like Cysteine Proteases in Gossypium hirsutum and Functional Characterization in Response to Verticillium dahliae

Cotton, a natural fiber producing crop of huge importance, is often prone to attack of Verticillium dahliae. Papain-like cysteine proteases (PLCPs) constitute a large family in plants and were proposed to involve in plant defense against pathogen attack in a number of studies. However, there is no detailed characterization of PLCP genes in cotton against infection of V. dahliae. In this study, we carried out a genome-wide analysis in cotton and identified seventy-eight PLCPs, which were divided into nine subfamilies based on their evolution phylogeny: RD21 (responsive to desiccation 21), CEP (cysteine endopeptidase), XCP (xylem cysteine peptidase), XBCP3 (xylem bark cysteine peptidase 3), THI, SAG12 (senescence-associated gene 12), RD19 (responsive to desiccation 19), ALP (aleurain-like protease) and CTB (cathepsin B-like). Genes in each subfamily exhibit a similar structure and motif composition. The expression patterns of these genes in different organs were examined, and subfamily RD21 was the most abundant in these families. Expression profiles under abiotic stress showed that thirty-five PLCP genes were induced by multiple stresses. Further transcriptome analysis showed that sixteen PLCP genes were up-regulated in response to V. dahliae in cotton. Among those, GhRD21-7 showed a higher transcription level than most other PLCP genes. Additionally, over-expression of GhRD21-7 led to enhanced resistance and RNAi lines were more susceptible to V. dahliae in cotton. Our results provide valuable information for future functional genomic studies of PLCP gene family in cotton.

Genome-wide identification and characterization of the Dof gene family in cassava (Manihot esculenta)

Dof proteins comprise a plant-specific transcription factor family involved in regulation of many plant processes. Despite their diverse roles, little information is available in cassava (Manihot esculenta, Euphorbiaceae), the sixth major staple crop in the world. This study presents the first genome-wide identification of Dof family genes in cassava, resulting in 45 family members that are unevenly distributed across 17 out of the 18 chromosomes. Phylogenetic analysis assigned these genes into eight groups representing four subfamilies. Synteny analysis and BRH (Best Reciprocal Hit) method were also adopted to study evolutionary relationships. Among 21 MeDof duplicates identified, the majority of them were derived from the recent whole-genome duplication (WGD), whereas only one from local duplication. Transcriptional profiling revealed diverse expression patterns among various tissues, and further comparison of conserved motifs and expression profiles revealed divergence of MeDof duplicates. These findings provide valuable information for functional analysis of Dof genes in cassava.