Identification and expression analysis of the GDSL esterase/lipase family genes, and the characterization of SaGLIP8 in Sedum alfredii Hance under cadmium stress

Identification, evolution, and expression of GDSL-type Esterase/Lipase (GELP) gene family in three cotton species: a bioinformatic analysis

BACKGROUND

GDSL esterase/lipases (GELPs) play important roles in plant growth, development, and response to biotic and abiotic stresses. Presently, an extensive and in-depth analysis of GELP family genes in cotton is still not clear enough, which greatly limits the further understanding of cotton GELP function and regulatory mechanism.

RESULTS

A total of 389 GELP family genes were identified in three cotton species of Gossypium hirsutum (193), G. arboreum (97), and G. raimondii (99). These GELPs could be classified into three groups and eight subgroups, with the GELPs in same group to have similar gene structures and conserved motifs. Evolutionary event analysis showed that the GELP family genes tend to be diversified at the spatial dimension and certain conservative at the time dimension, with a trend of potential continuous expansion in the future. The orthologous or paralogous GELPs among different genomes/subgenomes indicated the inheritance from genome-wide duplication during polyploidization, and the paralogous GELPs were derived from chromosomal segment duplication or tandem replication. GELP genes in the A/D subgenome underwent at least three large-scale replication events in the evolutionary process during the period of 0.6-3.2 MYA, with two large-scale evolutionary events between 0.6-1.8 MYA that were associated with tetraploidization, and the large-scale duplication between 2.6-9.1 MYA that occurred during diploidization. The cotton GELPs indicated diverse expression patterns in tissue development, ovule and fiber growth, and in response to biotic and abiotic stresses, combining the existing cis-elements in the promoter regions, suggesting the GELPs involvements of functions to be diversification and of the mechanisms to be a hormone-mediated manner.

CONCLUSIONS

Our results provide a systematic and comprehensive understanding the function and regulatory mechanism of cotton GELP family, and offer an effective reference for in-depth genetic improvement utilization of cotton GELPs.

Phenotypic and Proteomic Insights into Differential Cadmium Accumulation in Maize Kernels

The contamination of agricultural soil with cadmium (Cd), a heavy metal, poses a significant environmental challenge, affecting crop growth, development, and human health. Previous studies have established the pivotal role of the ZmHMA3 gene, a P-type ATPase heavy metal transporter, in determining variable Cd accumulation in maize grains among 513 inbred lines. To decipher the molecular mechanism underlying mutation-induced phenotypic differences mediated by ZmHMA3, we conducted a quantitative tandem mass tag (TMT)-based proteomic analysis of immature maize kernels. This analysis aimed to identify differentially expressed proteins (DEPs) in wild-type B73 and ZmHMA3 null mutant under Cd stress. The findings demonstrated that ZmHMA3 accumulated higher levels of Cd compared to B73 when exposed to varying Cd concentrations in the soil. In comparison to soil with a low Cd concentration, B73 and ZmHMA3 exhibited 75 and 142 DEPs, respectively, with 24 common DEPs shared between them. ZmHMA3 showed a higher induction of upregulated genes related to Cd stress than B73. Amino sugar and nucleotide sugar metabolism was specifically enriched in B73, while phenylpropanoid biosynthesis, nitrogen metabolism, and glyoxylate and dicarboxylate metabolism appeared to play a more significant role in ZmHMA3. This study provides proteomics insights into unraveling the molecular mechanism underlying the differences in Cd accumulation in maize kernels.

Physiological and gene expression responses involved in teak (Tectona grandis L.) seedlings exposed to osmotic and salt stressors

BACKGROUND

Teak (Tectona grandis L.) is a forest tree having 2n = 2x = 36 diploid chromosomes. Plants are continually subjected to variety of abiotic stresses due to climate change, which alter their physiological processes and gene expression.

METHODS AND RESULTS

The current study sought to examine the physiological and differential gene expression of teak seedlings exposed to abiotic stresses (150 mM NaCl and 15% PEG-6000). Chlorophyll content, membrane stability index and relative water content were measured at 0, 2, 7 and 12 days after treatment. These parameters were initially numerically reduced, but they were significantly reduced during a longer period of treatment. Seedlings treated with 150 mM NaCl displayed more harmful effect on the plant than other treatments. The results showed that variety of stresses significantly affect the physiology of seedlings because they cause membrane damage, ROS generation, chlorophyll degradation, and reduction in water absorption. The gene expression of treated and control seedlings was also evaluated at 12 days after treatment. Ten stress-related genes were examined for their differential expression using RT-PCR under applied stress. The stress-treated seedlings' leaves showed an up-regulated expression of the genes MYB-3, HSP-1, BI-1 and CS-2.

CONCLUSION

Up-regulation of the genes confirmed the protective function of these genes in plants under abiotic stress. However, gene expression was affected by treatments, the extent of stress and the species of plant. This study came to the conclusion that physiological parameters could be utilized as marker indices to assess a tree's capability to withstand stress at seedling stage. The up-regulated genes will be further investigated and utilized to validate stress tolerance and susceptible teak seedlings.

Genome-wide exploration of the GDSL-type esterase/lipase gene family in rapeseed reveals several BnGELP proteins active during early seedling development

The Gly-Asp-Ser-Leu (GDSL)-type esterase/lipase proteins (GELP) are one of the most important families of lipolytic enzymes and play prominent roles in seed germination and early seedling establishment through mobilizing the lipids stored in seeds. However, there are no comprehensive studies systematically investigating the GELP gene family in Brassica napus (BnGELP), and their biological significance to these physiological processes are far from understood. In the present study, a total of 240 BnGELP genes were identified in B. napus cultivar "Zhongshuang 11" (ZS11), which is nearly 2.3-fold more GELP genes than in Arabidopsis thaliana. The BnGELP genes clustered into 5 clades based on phylogenetic analysis. Ten BnGELPs were identified through zymogram analysis of esterase activity followed by mass spectrometry, among which five clustered into the clade 5. Gene and protein architecture, gene expression, and cis-element analyses of BnGELP genes in clade 5 suggested that they may play different roles in different tissues and in response to different abiotic stresses. BnGELP99 and BnGELP159 were slightly induced by cold, which may be attributed to two low-temperature responsive cis-acting regulatory elements present in their promoters. An increased activity of esterase isozymes by cold was also observed, which may reflect other cold inducible esterases/lipases in addition to the ten identified BnGELPs. This study provides a systemic view of the BnGELP gene family and offers a strategy for researchers to identify candidate esterase/lipase genes responsible for lipid mobilization during seed germination and early seedling establishment.

The binding pocket properties were fundamental to functional diversification of the GDSL-type esterases/lipases gene family in cotton

Cotton is one of the most important crops in the world. GDSL-type esterases/lipases (GELPs) are widely present in all kingdoms and play an essential role in regulating plant growth, development, and responses to abiotic and biotic stresses. However, the molecular mechanisms underlying this functional diversity remain unclear. Here, based on the identification of the GELP gene family, we applied genetic evolution and molecular simulation techniques to explore molecular mechanisms in cotton species. A total of 1502 GELP genes were identified in 10 cotton species. Segmental duplication and differences in evolutionary rates are the leading causes of the increase in the number and diversity of GELP genes during evolution for ecological adaptation. Structural analysis revealed that the GELP family has high structural diversity. Moreover, molecular simulation studies have demonstrated significant differences in the properties of the binding pockets among cotton GELPs. In the process of adapting to the environment, GELPs not only have segmental duplication but also have different evolutionary rates, resulting in gene diversity. This diversity leads to significant differences in the 3D structure and binding pocket properties and, finally, to functional diversity. These findings provide a reference for further functional analyses of plant GELPs.

Genome wide identification of GDSL gene family explores a novel GhirGDSL26 gene enhancing drought stress tolerance in cotton

BACKGROUND

Current climate change scenarios are posing greater threats to the growth and development of plants. Thus, significant efforts are required that can mitigate the negative effects of drought on the cotton plant. GDSL esterase/lipases can offer an imperative role in plant development and stress tolerance. However, thesystematic and functional roles of the GDSL gene family, particularly in cotton under water deficit conditions have not yet been explored.

RESULTS

In this study, 103, 103, 99, 198, 203, 239, 249, and 215 GDSL proteins were identified in eight cotton genomes i.e., Gossypium herbaceum (A1), Gossypium arboretum (A2), Gossypium raimondii (D5), Gossypium hirsutum (AD1), Gossypium barbadense (AD2), Gossypium tomentosum (AD3), Gossypium mustelinum (AD4), Gossypium darwinii (AD5), respectively. A total of 198 GDSL genes of Gossypium hirsutum were divided into eleven clades using phylogenetic analysis, and the number of GhirGDSL varied among different clades. The cis-elements analysis showed that GhirGDSL gene expression was mainly related to light, plant hormones, and variable tense environments. Combining the results of transcriptome and RT-qPCR, GhirGDSL26 (Gh_A01G1774), a highly up-regulated gene, was selected for further elucidating its tole in drought stress tolerance via estimating physiological and biochemical parameters. Heterologous expression of the GhirGDSL26 gene in Arabidopsis thaliana resulted in a higher germination and survival rates, longer root lengths, lower ion leakage and induced stress-responsive genes expression under drought stress. This further highlighted that overexpressed plants had a better drought tolerance as compared to the wildtype plants. Moreover, 3, 3'-diaminobenzidine (DAB) and Trypan staining results indicated reduced oxidative damage, less cell membrane damage, and lower ion leakage in overexpressed plants as compared to wild type. Silencing of GhirGDSL26 in cotton via VIGS resulting in a susceptible phenotype, higher MDA and H₂O₂ contents, lower SOD activity, and proline content.

CONCLUSION

Our results demonstrated that GhirGDSL26 plays a critical role in cotton drought stress tolerance. Current findings enrich our knowledge of GDSL genes in cotton and provide theoretical guidance and excellent gene resources for improving drought tolerance in cotton.

Genome-Wide Classification and Phylogenetic Analyses of the GDSL-Type Esterase/Lipase (GELP) Family in Flowering Plants

GDSL-type esterase/lipase (GELP) enzymes have key functions in plants, such as developmental processes, anther and pollen development, and responses to biotic and abiotic stresses. Genes that encode GELP belong to a complex and large gene family, ranging from tens to more than hundreds of members per plant species. To facilitate functional transfer between them, we conducted a genome-wide classification of GELP in 46 plant species. First, we applied an iterative phylogenetic method using a selected set of representative angiosperm genomes (three monocots and five dicots) and identified 10 main clusters, subdivided into 44 orthogroups (OGs). An expert curation for gene structures, orthogroup composition, and functional annotation was made based on a literature review. Then, using the HMM profiles as seeds, we expanded the classification to 46 plant species. Our results revealed the variable evolutionary dynamics between OGs in which some expanded, mostly through tandem duplications, while others were maintained as single copies. Among these, dicot-specific clusters and specific amplifications in monocots and wheat were characterized. This approach, by combining manual curation and automatic identification, was effective in characterizing a large gene family, allowing the establishment of a classification framework for gene function transfer and a better understanding of the evolutionary history of GELP.

Genome-Wide Analysis of the GDSL Genes in Pecan (Carya illinoensis K. Koch): Phylogeny, Structure, Promoter Cis-Elements, Co-Expression Networks, and Response to Salt Stresses

The Gly-Asp-Ser-Leu (GDSL)-lipase family is a large subfamily of lipolytic enzymes that plays an important role in plant growth and defense against environmental stress. However, little is known about their function in pecans (Carya illinoensis K. Koch). In this study, 87 CilGDSLs were identified and divided into 2 groups and 12 subgroups using phylogenetic analysis; members of the same sub-branch had conserved gene structure and motif composition. The majority of the genes had four introns and were composed of an α-helix and a β-strand. Subcellular localization analysis revealed that these genes were localized in the extracellular matrix, chloroplasts, cytoplasm, nucleus, vacuole, and endoplasmic reticulum, and were validated by transient expression in tobacco mesophyll cells. Furthermore, the analysis of the promoter cis-elements for the CilGDSLs revealed the presence of plant anaerobic induction regulatory, abscisic acid response, light response elements, jasmonic acid (JA) response elements, etc. The qRT-PCR analysis results in “Pawnee” with salt treatment showed that the CilGDSL42.93 (leaf) and CilGDSL39.88 (root) were highly expressed in different tissues. After salt stress treatment, isobaric tags for relative and absolute quantitation (iTRAQ) analysis revealed the presence of a total of ten GDSL proteins. Moreover, the weighted gene co-expression network analysis (WGCNA) showed that one set of co-expressed genes (module), primarily CilGDSL41.11, CilGDSL39.49, CilGDSL34.85, and CilGDSL41.01, was significantly associated with salt stress in leaf. In short, some of them were shown to be involved in plant defense against salt stress in this study.

Metagenomic Screening for Lipolytic Genes Reveals an Ecology-Clustered Distribution Pattern

Lipolytic enzymes are one of the most important enzyme types for application in various industrial processes. Despite the continuously increasing demand, only a small portion of the so far encountered lipolytic enzymes exhibit adequate stability and activities for biotechnological applications. To explore novel and/or extremophilic lipolytic enzymes, microbial consortia in two composts at thermophilic stage were analyzed using function-driven and sequence-based metagenomic approaches. Analysis of community composition by amplicon-based 16S rRNA genes and transcripts, and direct metagenome sequencing revealed that the communities of the compost samples were dominated by members of the phyla Actinobacteria, Proteobacteria, Firmicutes, Bacteroidetes, and Chloroflexi. Function-driven screening of the metagenomic libraries constructed from the two samples yielded 115 unique lipolytic enzymes. The family assignment of these enzymes was conducted by analyzing the phylogenetic relationship and generation of a protein sequence similarity network according to an integrated classification system. The sequence-based screening was performed by using a newly developed database, containing a set of profile Hidden Markov models, highly sensitive and specific for detection of lipolytic enzymes. By comparing the lipolytic enzymes identified through both approaches, we demonstrated that the activity-directed complements sequence-based detection, and vice versa. The sequence-based comparative analysis of lipolytic genes regarding diversity, function and taxonomic origin derived from 175 metagenomes indicated significant differences between habitats. Analysis of the prevalent and distinct microbial groups providing the lipolytic genes revealed characteristic patterns and groups driven by ecological factors. The here presented data suggests that the diversity and distribution of lipolytic genes in metagenomes of various habitats are largely constrained by ecological factors.

Comparative Transcriptome Analysis of Softening and Ripening-Related Genes in Kiwifruit Cultivars Treated with Ethylene

This work presents the transcriptome analysis of green ‘Hayward’ (Actinidia deliciosa) and gold ‘Haegeum’ (Actinidia chinensis) kiwifruit cultivars after treatment with ethylene for three days at 25 °C. Illumina high-throughput sequencing platform was used to sequence total mRNAs and the transcriptome gene set was constructed by de novo assembly. A total of 1287 and 1724 unigenes were differentially expressed during the comparison of ethylene treatment with control in green ‘Hayward’ and gold ‘Haegeum’, respectively. From the differentially expressed unigenes, 594 and 906 were upregulated, and 693 and 818 were downregulated in the green and gold kiwifruit cultivars, respectively, when treated with ethylene. We also identified a list of genes that were expressed commonly and exclusively in the green and gold kiwifruit cultivars treated with ethylene. Several genes were expressed differentially during the ripening of kiwifruits, and their cumulative effect brought about the softening- and ripening-related changes. This work also identified and categorized genes related to softening and other changes during ripening. Furthermore, the transcript levels of 12 selected representative genes from the differentially expressed genes (DEGs) identified in the transcriptome analysis were confirmed via quantitative real-time PCR (qRT-PCR) to validate the reliability of the expression profiles obtained from RNA-Seq. The data obtained from the present study will add to the information available on the molecular mechanisms of the effects of ethylene during the ripening of kiwifruits. This study will also provide resources for further studies of the genes related to ripening, helping kiwifruit breeders and postharvest technologists to improve ripening quality.

Comparative and Systematic Omics Revealed Low Cd Accumulation of Potato StMTP9 in Yeast: Suggesting a New Mechanism for Heavy Metal Detoxification

The metal tolerance protein (MTP) family is a very old family with evolutionary conservation and less specific amplification. It seems to retain the original functions of the ancestral genes and plays an important role in maintaining metal homeostasis in plant cells. We identified the potato MTP family members for the first time, the specific and conservative StMPTs were discovered by using systematic and comparative omics. To be surprised, members of the StMTP family seem to have mutated before the evolution of dicotyledon and monocotyledon, and even the loss of the entire subfamily (subfamily G6, G7). Interestingly, StMTP9 represents the conserved structure of the entire subfamily involved in toxic metal regulation. However, the gene structure and transmembrane domain of StMTP8 have undergone specific evolution, showing that the transmembrane domain (Motif13) located at the NH₂ terminal has been replaced by the signal peptide domain, so it was selected as the control gene of StMTP9. Through real-time fluorescence quantitative analysis of StMTPs under Cd and Zn stress, a co-expression network was constructed, and it was found that StMTP9 responded significantly to Cd stress, while StMTP8 did the opposite. What excites us is that by introducing StMTPs 8/9 into the ∆ycf1 yeast cadmium-sensitive mutant strain, the functional complementation experiment proved that StMTPs 8/9 can restore Cd tolerance. In particular, StMTP9 can greatly reduce the cadmium content in yeast cells, while StMTP8 cannot. These findings provide a reference for further research on the molecular mechanism of potato toxic metal accumulation.

Genome-wide characterization of the hyperaccumulator Sedum alfredii F-box family under cadmium stress

The F-box genes, which form one of the largest gene families in plants, are vital for plant growth, development and stress response. However, F-box gene family in Sedum alfredii remains unknown. Comprehensive studies addressing their function responding to cadmium stress is still limited. In the present study, 193 members of the F-box gene (SaFbox) family were identified, which were classified into nine subfamilies. Most of the SaFboxs had highly conserved domain and motif. Various functionally related cis-elements involved in plant growth regulation, stress and hormone responses were located in the upstream regions of SaFbox genes. RNA-sequencing and co-expression network analysis revealed that the identified SaFbox genes would be involved in Cd stress. Expression analysis of 16 hub genes confirmed their transcription level in different tissues. Four hub genes (SaFbox40, SaFbox51, SaFbox136 and SaFbox170) were heterologously expressed in a Cd-sensitive yeast cell to assess their effects on Cd tolerance. The transgenic yeast cells carrying SaFbox40, SaFbox51, SaFbox136, or SaFbox170 were more sensitive and accumulated more cadmium under Cd stress than empty vector transformed control cells. Our results performed a comprehensive analysis of Fboxs in S. alfredii and identified their potential roles in Cd stress response.

Genome-Wide Investigation and Expression Profiling Under Abiotic Stresses of a Soybean Unknown Function (DUF21) and Cystathionine-β-Synthase (CBS) Domain-Containing Protein Family

Cystathionine-β-synthase (CBS) domain-containing proteins (CDCPs) constitute a large family in plants, and members of this family have been implicated in a variety of biological processes. However, the precise functions and the underlying mechanisms of most members of this family in plants remain to be elucidated. CBSDUF proteins belong to the CDCP superfamily, which contains one domain of unknown function (DUF21) and an N terminus that is adjacent to two intracellular CBS domains. In this study, a comprehensive genome database analysis of soybean was performed to investigate the role(s) of these CBSDUFs and to explore their nomenclature, classification, chromosomal distribution, exon-intron organization, protein structure, and phylogenetic relationships; the analysis identified a total of 18 putative CBSDUF genes. Using specific protein domains and phylogenetic analysis, the CBSDUF gene family was subdivided into eight groups. The soybean CBSDUF genes showed an uneven distribution on 12 chromosomes of Glycine max. RNA-seq transcriptome data from different tissues in public databases revealed tissue-specific and differential expression profiles of the GmCBSDUFs, and qPCR analysis revealed that certain groups of soybean CBSDUFs are likely involved in specific stress responses. In addition, GmCBSDUF3 transgenic Arabidopsis was subjected to phenotypic analysis under NaCl, PEG, and ABA stress treatments. The overexpression of GmCBSDUF3 could enhance tolerance to drought and salt stress in Arabidopsis. This study presents a first comprehensive look at soybean CBSDUF proteins and provides valuable resources for functionally elucidating this protein subgroup within the CBS domain-containing protein family.

Comparative Transcriptome Analysis of the Molecular Mechanism of the Hairy Roots of Brassica campestris L. in Response to Cadmium Stress

Brassica campestris L., a hyperaccumulator of cadmium (Cd), is considered a candidate plant for efficient phytoremediation. The hairy roots of Brassica campestris L are chosen here as a model plant system to investigate the response mechanism of Brassica campestris L. to Cd stress. High-throughput sequencing technology is used to identify genes related to Cd tolerance. A total of 2394 differentially expressed genes (DEGs) are identified by RNA-Seq analysis, among which 1564 genes are up-regulated, and 830 genes are down-regulated. Data from the gene ontology (GO) analysis indicate that DEGs are mainly involved in metabolic processes. Glutathione metabolism, in which glutathione synthetase and glutathione S-transferase are closely related to Cd stress, is identified in the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. A Western blot shows that glutathione synthetase and glutathione S-transferase are involved in Cd tolerance. These results provide a preliminary understanding of the Cd tolerance mechanism of Brassica campestris L. and are, hence, of particular importance to the future development of an efficient phytoremediation process based on hairy root cultures, genetic modification, and the subsequent regeneration of the whole plant.