Multifunctionality and diversity of GDSL esterase/lipase gene family in rice (Oryza sativa L. japonica) genome: new insights from bioinformatics analysis

Genome-wide identification of the GDSL-type esterase/lipase protein (GELP) gene family in Ricinus communis and its transcriptional regulation during germination and seedling establishment under different abiotic stresses

GDSL-type esterase/lipase protein (GELP) genes are crucial in the specialized lipid metabolism, in the responses to abiotic stresses, and in the regulation of plant homeostasis. R. communis is an important oilseed crop species that can sustain growth and productivity when exposed to harsh environmental conditions. Herein, we raised the question of whether the GELP gene family could be involved in the acquisition of R. communis tolerance to abiotic stresses during seed germination and seedling establishment. Thus, we used bioinformatics and transcriptomics to characterize the R. communis GELP gene family. R. communis genome possesses 96 GELP genes that were characterized by extensive bioinformatics, including phylogenetic analysis, subcellular localization, exon-intron distribution, the analysis of regulatory cis-elements, tandem duplication, and physicochemical properties. Transcriptomics indicated that numerous RcGELP genes are readily responsive to high-temperature and salt stresses and might be potential candidates for genome editing techniques to develop abiotic stress-tolerant crops.

Understanding pyrethrin biosynthesis: toward and beyond natural pesticide overproduction

Pyrethrins are natural insecticides biosynthesised by Asteraceae plants, such as Tanacetum cinerariifolium and have a long history, dating back to ancient times. Pyrethrins are often used as low-persistence and safe insecticides to control household, horticultural, and agricultural insect pests. Despite its long history of use, pyrethrin biosynthesis remains a mystery, presenting a significant opportunity to improve yields and meet the growing demand for organic agriculture. To achieve this, both genetic modification and non-genetic methods, such as chemical activation and priming, are indispensable. Plants use pyrethrins as a defence against herbivores, but pyrethrin biosynthesis pathways are shared with plant hormones and signal molecules. Hence, the insight that pyrethrins may play broader roles than those traditionally expected is invaluable to advance the basic and applied sciences of pyrethrins.

A lipase from Lacticaseibacillus rhamnosus IDCC 3201 with thermostability and pH resistance for use as a detergent additive

Lipases are important biocatalysts and ubiquitous in plants, animals, and microorganisms. The high growth rates of microorganisms with low production costs have enabled the wide application of microbial lipases in detergent, food, and cosmetic industries. Herein, a novel lipase from Lacticaseibacillus rhamnosus IDCC 3201 (Lac-Rh) was isolated and its activity analyzed under a range of reaction conditions to evaluate its potential industrial application. The isolated Lac-Rh showed a molecular weight of 24 kDa and a maximum activity of 3438.5 ± 1.8 U/mg protein at 60 °C and pH 8. Additionally, Lac-Rh retained activity in alkaline conditions and in 10% v/v concentrations of organic solvents, including glycerol and acetone. Interestingly, after pre-incubation in the presence of multiple commercial detergents, Lac-Rh maintained over 80% of its activity and the stains from cotton were successfully removed under a simulated laundry  setting. Overall, the purified lipase from L. rhamnosus IDCC 3201 has potential for use as a detergent in industrial applications. KEY POINTS: • A novel lipase (Lac-Rh) was isolated from Lacticaseibacillus rhamnosus IDCC 3201 • Purified Lac-Rh exhibited its highest activity at a temperature of 60 °C and a pH of 8, respectively • Lac-Rh remains stable in commercial laundry detergent and enhances washing performance.

GDSL Lipase Gene HTA1 Negatively Regulates Heat Tolerance in Rice Seedlings by Regulating Reactive Oxygen Species Accumulation

High temperature is a significant environmental stress that limits plant growth and agricultural productivity. GDSL lipase is a hydrolytic enzyme with a conserved GDSL sequence at the N-terminus, which has various biological functions, such as participating in plant growth, development, lipid metabolism, and stress resistance. However, little is known about the function of the GDSL lipase gene in the heat tolerance of rice. Here, we characterized a lipase family protein coding gene HTA1, which was significantly induced by high temperature in rice. Rice seedlings in which the mutant hta1 was knocked out showed enhanced heat tolerance, whereas the overexpressing HTA1 showed more sensitivity to heat stress. Under heat stress, hta1 could reduce plant membrane damage and reactive oxygen species (ROS) levels and elevate the activity of antioxidant enzymes. Moreover, real-time quantitative PCR (RT-qPCR) analysis showed that mutant hta1 significantly activated gene expression in antioxidant enzymes, heat response, and defense. In conclusion, our results suggest that HTA1 negatively regulates heat stress tolerance by modulating the ROS accumulation and the expression of heat-responsive and defense-related genes in rice seedlings. This research will provide a valuable resource for utilizing HTA1 to improve crop heat tolerance.

OsGELP77, a QTL for broad-spectrum disease resistance and yield in rice, encodes a GDSL-type lipase

Lipids and lipid metabolites have essential roles in plant-pathogen interactions. GDSL-type lipases are involved in lipid metabolism modulating lipid homeostasis. Some plant GDSLs modulate lipid metabolism altering hormone signal transduction to regulate host-defence immunity. Here, we functionally characterized a rice lipase, OsGELP77, promoting both immunity and yield. OsGELP77 expression was induced by pathogen infection and jasmonic acid (JA) treatment. Overexpression of OsGELP77 enhanced rice resistance to both bacterial and fungal pathogens, while loss-of-function of osgelp77 showed susceptibility. OsGELP77 localizes to endoplasmic reticulum and is a functional lipase hydrolysing universal lipid substrates. Lipidomics analyses demonstrate that OsGELP77 is crucial for lipid metabolism and lipid-derived JA homeostasis. Genetic analyses confirm that OsGELP77-modulated resistance depends on JA signal transduction. Moreover, population genetic analyses indicate that OsGELP77 expression level is positively correlated with rice resistance against pathogens. Three haplotypes were classified based on nucleotide polymorphisms in the OsGELP77 promoter where OsGELP77Hap3 is an elite haplotype. Three OsGELP77 haplotypes are differentially distributed in wild and cultivated rice, while OsGELP77Hap3 has been broadly pyramided for hybrid rice development. Furthermore, quantitative trait locus (QTL) mapping and resistance evaluation of the constructed near-isogenic line validated OsGELP77, a QTL for broad-spectrum disease resistance. In addition, OsGELP77-modulated lipid metabolism promotes JA accumulation facilitating grain yield. Notably, the hub defence regulator OsWRKY45 acts upstream of OsGELP77 by initiating the JA-dependent signalling to trigger immunity. Together, OsGELP77, a QTL contributing to immunity and yield, is a candidate for breeding broad-spectrum resistant and high-yielding rice.

Molecular cloning, overexpression, characterization, and In silico modelling analysis of a novel GDSL autotransporter-dependent outer membrane lipase (OML) of Pseudomonas guariconensis

The outer membrane lipase (oml) gene, encoding a novel autotransporter-dependent lipase from Pseudomonas guariconensis, was cloned and sequenced. The oml gene has an open reading frame of 1866 bp. It encodes the 621 amino acid autotransporter-dependent GDSL lipase (OML), which has the highest sequence similarity (64.08 %) with the EstA of Pseudomonas aeruginosa (PDB:3kvn.1. A). OML was expressed and purified, which showed a purified band of approximately 70 kDa. The purified enzyme showed maximum activity at pH 9 and 40 °C. Substrate specificity studies and kinetic study by Lineweaver-Burk plot of purified OML showed Km of 1.27 mM and Vmax of 333.33 U/mL with p-nitrophenyl palmitate. The purified enzyme showed good stability in the presence of hexane, methanol, and ethanol, while the presence of the metal ion Mg2+ showed maximum lipase activity. Bioinformatics analysis supported the in vitro findings by predicting enzyme substrate specificity towards long-chain fatty acids and fatty acids with shorter chain lengths. The stability of the interaction of the protein-ligand complex (OML-ricinoleic acid) was confirmed using MDS and castor oil bioconversion using purified OML was confirmed using High-Performance Liquid Chromatography (HPLC).

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.

The secret password: Cell death-inducing proteins in filamentous phytopathogens - As versatile tools to develop disease-resistant crops

Cell death-inducing proteins (CDIPs) are some of the secreted effector proteins manifested by filamentous oomycetes and fungal pathogens to invade the plant tissue and facilitate infection. Along with their involvement in different developmental processes and virulence, CDIPs play a crucial role in plant-pathogen interactions. As the name implies, CDIPs cause necrosis and trigger localised cell death in the infected host tissues by the accumulation of higher concentrations of hydrogen peroxide (H2O2), oxidative burst, accumulation of nitric oxide (NO), and electrolyte leakage. They also stimulate the biosynthesis of defense-related phytohormones such as salicylic acid (SA), jasmonic acid (JA), abscisic acid (ABA), and ethylene (ET), as well as the expression of pathogenesis-related (PR) genes that are important in disease resistance. Altogether, the interactions result in the hypersensitive response (HR) in the host plant, which might confer systemic acquired resistance (SAR) in some cases against a vast array of related and unrelated pathogens. The CDIPs, due to their capability of inducing host resistance, are thus unique among the array of proteins secreted by filamentous plant pathogens. More interestingly, a few transgenic plant lines have also been developed expressing the CDIPs with added resistance. Thus, CDIPs have opened an interesting hot area of research. The present study critically reviews the current knowledge of major types of CDIPs identified across filamentous phytopathogens and their modes of action in the last couple of years. This review also highlights the recent breakthrough technologies in studying plant-pathogen interactions as well as crop improvement by enhancing disease resistance through CDIPs.

Towards carotenoid biofortification in wheat: identification of XAT-7A1, a multicopy tandem gene responsible for carotenoid esterification in durum wheat

Yellow pigment content, mainly due to the accumulation of carotenoids, is a quality trait in durum wheat grain as it confers the bright yellow color to pasta preferred by consumers. Also, carotenoids are essential nutrients exerting important biological functions in human health. Consequently, biofortification strategies have been developed in many crops to increase carotenoid content. In this context, carotenoid esterification is emerging as a new breeding target for wheat biofortification, as carotenoid esters have been found to promote both carotenoid accumulation and stability. Until recently, no carotenoid esters have been identified in significant proportions in durum wheat grains, and interspecific breeding programs have been started to transfer esterification ability from common wheat and Hordeum chilense.In this work, XAT-7A1 is identified as the gene responsible for carotenoid esterification in durum wheat. Sequencing, copy number variation and mapping results show that XAT-7A1 is organized as tandem or proximal GDSL esterase/lipase copies in chromosome 7A. Three XAT-7A1 haplotypes are described: Type 1 copies, associated with high levels of carotenoid esters (diesters and monoesters) production and high expression in grain development; Type 2 copies, present in landraces with low levels of carotenoid esters (monoesters) or no esters; and Type 3 copies, without the signal peptide, resulting in zero-ester phenotypes.The identification of XAT-7A1 is a necessary step to make the carotenoid esterification ability available for durum and bread wheat breeding, which should be focused on the Type 1 XAT-7A1 haplotype, which may be assessed as a single gene since XAT-7A1 copies are inherited together.

Pan-transcriptomic Profiling Demarcates Serendipita Indica-Phosphorus Mediated Tolerance Mechanisms in Rice Exposed to Arsenic Toxicity

Inadvertent accumulation of arsenic (As) in rice (Oryza sativa L.) is a concern for people depending on it for their subsistence, as it verily causes epigenetic alterations across the genome as well as in specific cells. To ensure food safety, certain attempts have been made to nullify this highest health hazard encompassing physiological, chemical and biological methods. Albeit, the use of mycorrhizal association along with nutrient reinforcement strategy has not been explored yet. Mechanisms of response and resistance of two rice genotypes to As with or without phosphorus (P) nutrition and Serendipita indica (S. indica; S.i) colonization were explored by root transcriptome profiling in the present study. Results revealed that the resistant genotype had higher auxin content and root plasticity, which helped in keeping the As accumulation and P starvation response to a minimum under alone As stress. However, sufficient P supply and symbiotic relationship switched the energy resources towards plant's developmental aspects rather than excessive root proliferation. Higher As accumulating genotype (GD-6) displayed upregulation of ethylene signaling/biosynthesis, root stunting and senescence related genes under As toxicity. Antioxidant defense system and cytokinin biosynthesis/signaling of both genotypes were strengthened under As + S.i + P, while the upregulation of potassium (K) and zinc (Zn) transporters depicted underlying cross-talk with iron (Fe) and P. Differential expression of phosphate transporters, peroxidases and GSTs, metal detoxification/transport proteins, as well as phytohormonal metabolism were responsible for As detoxification. Taken together, S. indica symbiosis fortified with adequate P-fertilizer can prove to be effective in minimizing As acquisition and accumulation in rice plants.

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.

Identification and co-expression network analysis of plumule-preferentially expressed genes in Oryza sativa

BACKGROUND

The seedling establishment is controlled by the programmed expression of sets of genes at the specific tissues of seed, abundance and environment. Plumule is an important part of the seed embryo and expresses the suits of genes to exert distinct functions during seed germination. Although rice genomic resources are available and developed rapidly, thousands of transcripts have not previously been located in the plumule of rice.

OBJECTIVE

This study was performed to identify plumule-preferentially expressed (OsPluP) genes in rice and determine the expression profiles and functions of OsPluP genes.

METHODS

We identified the OsPluP genes through Affymetrix microarray data. Meanwhile, qRT-PCR was performed to validate the expression pattern, also found that OsPluP genes were regulated by dark/light treatment. The cis-acting regulatory elements were analyzed in the promoters' regions of OsPluP genes. The T-DNA mutant of the OsPluP seed was used to reveal the function in seed germination.

RESULTS

In this study, a genomic survey of OsPluP genes was performed, and we identified 88 OsPluP genes based on Affymetrix microarray data. The expression profiles of 88 OsPluP members in 24 representative tissues covering rice whole life cycle can be roughly classified into three major groups, suggesting functional divergence of OsPluP genes in seed germination. The microarray data, qRT-PCR, and promoter analysis results demonstrated that transcripts of more than half OsPluPs (54 genes) could be enhanced in the darkness and respond to phytohormone. Gene Ontology (GO)and Kyoto encyclopedia of genes and genomes (KEGG) analysis demonstrated that OsPluP and their co-expressed genes were highly enriched in fatty acid metabolism. Moreover, OsPluP82 T-DNA mutant seeds displayed short plumule length and storage lipid accumulation.

CONCLUSION

This study would enable the functions of OsPluP genes during seed germination and contribute to the goal of molecular regulatory networks that lay the foundation for further studies of seedling growth.

Bioinformatics-Based Screening Approach for the Identification and Characterization of Lipolytic Enzymes from the Marine Diatom Phaeodactylum tricornutum

Oleaginous diatoms accumulate lipids of biotechnological interest when exposed to nutrient stress conditions such as nitrogen starvation. While accumulation mechanisms are well-known and have been engineered to improve lipid production, degradation mechanisms remain poorly investigated in diatoms. Identifying lipid-degrading enzymes is the initial step to understanding the catabolic processes. In this study, an in silico screening of the genome of Phaeodactylum tricornutum led to the identification of 57 putative triacylglycerol lipases (EC 3.1.1.3) grouped in 4 families. Further analysis revealed the presence of conserved domains and catalytic residues of lipases. Physico-chemical characteristics and subcellular localization predictions highlighted that a majority of these putative proteins are hydrophilic and cytosolic, suggesting they could be recruited to lipid droplets directly from the cytosol. Among the 57 identified putative proteins, three lipases were identified as possibly involved in lipophagy due to a potential vacuolar localization. The expression of the mRNA corresponding to the 57 proteins was then searched in 3 transcriptomic datasets obtained under nitrogen starvation. Nine genes were highly regulated and were considered as encoding enzymes with a probable important function in lipid catabolism. A tertiary structure prediction of these nine candidates yielded eight functional 3D models. Among those, two downregulated enzymes, Phatr3_J54974 and Phatr3_EG00720, were highlighted as good targets for future functional genomics and purification studies to investigate their role in lipid degradation.

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.

Weighted Gene Co-Expression Network Analysis Reveals Hub Genes for Fuzz Development in Gossypium hirsutum

Fuzzless Gossypium hirsutum mutants are ideal materials for investigating cotton fiber initiation and development. In this study, we used the fuzzless G. hirsutum mutant Xinluzao 50 FLM as the research material and combined it with other fuzzless materials for verification by RNA sequencing to explore the gene expression patterns and differences between genes in upland cotton during the fuzz period. A gene ontology (GO) enrichment analysis showed that differentially expressed genes (DEGs) were mainly enriched in the metabolic process, microtubule binding, and other pathways. A weighted gene co-expression network analysis (WGCNA) showed that two modules of Xinluzao 50 and Xinluzao 50 FLM and four modules of CSS386 and Sicala V-2 were highly correlated with fuzz. We selected the hub gene with the highest KME value among the six modules and constructed an interaction network. In addition, we selected some genes with high KME values from the six modules that were highly associated with fuzz in the four materials and found 19 common differential genes produced by the four materials. These 19 genes are likely involved in the formation of fuzz in upland cotton. Several hub genes belong to the arabinogalactan protein and GDSL lipase, which play important roles in fiber development. According to the differences in expression level, 4 genes were selected from the 19 genes and tested for their expression level in some fuzzless materials. The modules, hub genes, and common genes identified in this study can provide new insights into the formation of fiber and fuzz, and provide a reference for molecular design breeding for the genetic improvement of cotton fiber.

Identification of multiple novel genetic mechanisms that regulate chilling tolerance in Arabidopsis

INTRODUCTION

Cold stress adversely affects the growth and development of plants and limits the geographical distribution of many plant species. Accumulation of spontaneous mutations shapes the adaptation of plant species to diverse climatic conditions.

METHODS

The genome-wide association study of the phenotypic variation gathered by a newly designed phenomic platform with the over six millions single nucleotide polymorphic (SNP) loci distributed across the genomes of 417 Arabidopsis natural variants collected from various geographical regions revealed 33 candidate cold responsive genes.

RESULTS

Investigation of at least two independent insertion mutants for 29 genes identified 16 chilling tolerance genes governing diverse genetic mechanisms. Five of these genes encode novel leucine-rich repeat domain-containing proteins including three nucleotide-binding site-leucine-rich repeat (NBS-LRR) proteins. Among the 16 identified chilling tolerance genes, ADS2 and ACD6 are the only two chilling tolerance genes identified earlier.

DISCUSSION

The 12.5% overlap between the genes identified in this genome-wide association study (GWAS) of natural variants with those discovered previously through forward and reverse genetic approaches suggests that chilling tolerance is a complex physiological process governed by a large number of genetic mechanisms.

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 Association Study (GWAS) and genome prediction of seedling salt tolerance in bread wheat (Triticum aestivum L.)

BACKGROUND

Salinity tolerance in wheat is imperative for improving crop genetic capacity in response to the expanding phenomenon of soil salinization. However, little is known about the genetic foundation underlying salinity tolerance at the seedling growth stage of wheat. Herein, a GWAS analysis was carried out by the random-SNP-effect mixed linear model (mrMLM) multi-locus model to uncover candidate genes responsible for salt tolerance at the seedling stage in 298 Iranian bread wheat accessions, including 208 landraces and 90 cultivars.

RESULTS

A total of 29 functional marker-trait associations (MTAs) were detected under salinity, 100 mM NaCl (sodium chloride). Of these, seven single nucleotide polymorphisms (SNPs) including rs54146, rs257, rs37983, rs18682, rs55629, rs15183, and rs63185 with R² ≥ 10% were found to be linked with relative water content, root fresh weight, root dry weight, root volume, shoot high, proline, and shoot potassium (K⁺), respectively. Further, a total of 27 candidate genes were functionally annotated to be involved in response to the saline environment. Most of these genes have key roles in photosynthesis, response to abscisic acid, cell redox homeostasis, sucrose and carbohydrate metabolism, ubiquitination, transmembrane transport, chromatin silencing, and some genes harbored unknown functions that all together may respond to salinity as a complex network. For genomic prediction (GP), the genomic best linear unbiased prediction (GBLUP) model reflected genetic effects better than both bayesian ridge regression (BRR) and ridge regression-best linear unbiased prediction (RRBLUP), suggesting GBLUP as a favorable tool for wheat genomic selection.

CONCLUSION

The SNPs and candidate genes identified in the current work can be used potentially for developing salt-tolerant varieties at the seedling growth stage by marker-assisted selection.

Transcriptome-wide modulation by Sargassum vulgare and Acanthophora spicifera extracts results in a prime-triggered plant signalling cascade in tomato and sweet pepper

Seaweed extracts (SWEs) are becoming integrated into crop production systems due to their multiple beneficial effects including growth promotion and induction of defence mechanisms. However, the comprehensive molecular mechanisms of these effects are yet to be elucidated. The current study investigated the transcriptomic changes induced by SWEs derived from Sargassum vulgare and Acanthophora spicifera on tomato and sweet pepper plants. Tomato and sweet pepper plants were subjected to foliar treatment with alkaline extracts prepared from the above seaweeds. Transcriptome changes in the plants were assessed 72 h after treatments using RNA sequencing. The treated plants were also analysed for defence enzyme activities, nutrient composition and phytohormonal profiles. The results showed the significant enrichment of genes associated with several growth and defence processes including photosynthesis, carbon and nitrogen metabolism, plant hormone signal transduction, plant-pathogen interaction, secondary metabolite metabolism, MAPK signalling and amino acid biosynthesis. Activities of defence enzymes were also significantly increased in SWE-treated plants. Plant nutrient profiling showed significant increases in calcium, potassium, nitrogen, sulphur, boron, copper, iron, manganese, zinc and phosphorous levels in SWE-treated plants. Furthermore, the levels of auxins, cytokinins and gibberellins were also significantly increased in the treated plants. The severity of bacterial leaf spot and early blight incidence in plants treated with SWE was significantly reduced, in addition to other effects like an increase in chlorophyll content, plant growth, and fruit yield. The results demonstrated the complex effect of S. vulgare and A. spicifera extracts on the plants' transcriptome and provided evidence of a strong role of these extracts in increasing plant growth responses while priming the plants against pathogenic attack simultaneously. The current study contributes to the understanding of the molecular mechanisms of SWEs in plants and helps their usage as a viable organic input for sustainable crop production.

A cold‐active esterase enhances mesophilic properties through Mn 2+ binding

A key aspect of adaptation to cold environments is the production of cold-active enzymes by psychrophilic organisms. These enzymes not only have high activity at low temperatures, but also exhibit remarkable structural flexibility and thermolability. In this context, the role of metal ions has been little explored, and the few available studies seem to suggest that metal binding counteracts structural flexibility. This article reports an investigation into the role of the binding of manganese ion (Mn²⁺ ) in the thermal adaptation of an esterase (M-Est) of the GDSx family, identified in the genome of the Antarctic bacterium Marinomonas sp. ef1. M-Est is specific for esters containing acetate groups and turned out to be a highly thermolabile cold-active enzyme, with a catalysis optimum temperature of 5 °C and a melting temperature of 31.7 °C. A combination of biochemical and computational analyses, including molecular dynamics simulations, revealed that M-Est binds Mn²⁺ ions via a single binding site located on the surface of the enzyme, close to the active site. Although the interaction between M-Est and Mn²⁺ induces only local conformational changes involving the active site, quite surprisingly they trigger an improvement in both thermal stability and catalytic efficiency under mild temperature conditions. These results, together with the conservation of the Mn²⁺ binding site among psychrophilic and psychrotolerant homologues, suggest that Mn²⁺ binding may be a useful, albeit atypical, strategy to mitigate the detrimental effects of temperature on true cold-active enzymes.

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.

Identification of Two GDSL-Type Esterase/Lipase Genes Related to Tissue-Specific Lipolysis in Dendrobium catenatum by Multi-Omics Analysis

Dendrobium catenatum is an important herb and widely cultivated in China. GDSL-Type Esterase/Lipase proteins (GELPs) are widely distributed in plants and play crucial roles in stress responses, plant growth, and development. However, no identification or functional analysis of GELPs was reported in D. catenatum. This study identifies 52 GELPs in D. catenatum genome, which is classified into four groups by phylogenetic analysis. Four conservative blocks (Ser-Gly-Asn-His) are found in most GELP domains. Transcriptome analysis reveals the expression profiles of GELPs in different organs and flowering phases. Co-expression analysis of the transcriptome and lipidome identifies a GELP gene, Dca016600, that positively correlates with 23 lipids. The purified Dca016600 protein shows the optimum pH is active from 8.0 to 8.5, and the optimum temperature is active from 30 °C to 40 °C. The kinetic study provides V_max (233.43 μmol·min^-1·mg^-1) and K_m (1.49 mM) for substrate p-nitrophenyl palmitate (p-NPP). Integrated analysis of the transcriptome and proteome identifies a GELP gene, Dca005399, which is specially induced by freezing. Interestingly, Dca005399 shows high expression in symbiotic germination seeds and sepals. This study provides new insights into the function of D. catenatum GELPs in plant development and stress tolerance.

Mutation in BrGGL7 gene encoding a GDSL esterase / lipase causes male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

MutMap and KASP analyses revealed that the BrGGL7 gene is responsible for the male-sterile trait of ftms1 in Chinese cabbage, with functional verification in Arabidopsis. The application of a male-sterile line is an ideal approach of hybrid seed production in Chinese cabbage. In this study, we obtained a male-sterile mutant (ftms1) from the double haploid line 'FT' using ethyl methane sulfonate (EMS) mutagenesis. The mutant was completely sterile due to abnormal enlargement and vacuolization of the tapetum cells. A single recessive nuclear gene was found to control male sterility in the mutant, while MutMap and KASP analyses identified BraA05g022470.3C (BrGGL7), which encodes a GDSL esterase / lipase, as the candidate mutant gene. A single nucleotide substitution from C to T occurred within the domain of BrGGL7 in ftms1, resulting in premature translation termination in the fourth exon. Meanwhile, qRT-PCR analysis indicated that BrGGL7 was prominently expressed in the anthers, and expression was greater in the wild-type 'FT' than ftms1. Genetic complementation of the orthologous Arabidopsis ggl7 mutant further confirmed the role of BrGGL7 in pollen development. These findings suggest that BrGGL7 plays a fundamental role in pollen formation, providing important insight into the molecular mechanisms underlying male sterility in Chinese cabbage.

Transcriptome analysis of tea (Camellia sinensis) leaves in response to ammonium starvation and recovery

The tea plant is a kind of ammonium-preferring crop, but the mechanism whereby ammonium (NH₄ ⁺) regulate its growth is not well understood. The current study focused on the effects of NH₄ ⁺ on tea plants. Transcriptomic analysis was performed to investigate the early- and late-stage NH₄ ⁺ deprivation and resupply in tea plants shoots. Through short- and long-term NH₄ ⁺ deficiency, the dynamic response to NH₄ ⁺ stress was investigated. The most significant effects of NH₄ ⁺ deficiency were found to be on photosynthesis and gene ontology (GO) enrichment varied with the length of NH₄ ⁺ deprivation. Enriched KEGG pathways were also different when NH₄ ⁺ was resupplied at different concentrations which may indicate reasons for tolerance of high NH₄ ⁺ concentration. Using weighted gene co-expression network analysis (WGCNA), modules related to significant tea components, tea polyphenols and free amino acids, were identified. Hence, NH₄ ⁺ could be regarded as a signaling molecule with the response of catechins shown to be higher than that of amino acids. The current work represents a comprehensive transcriptomic analysis of plant responses to NH₄ ⁺ and reveals many potential genes regulated by NH₄ ⁺ in tea plants. Such findings may lead to improvements in nitrogen efficiency of tea plants.

Spatiotemporal Transcriptomic Atlas of Developing Embryos and Vegetative Tissues in Flax

Flax (Linum usitatissimum L.) is an important multipurpose crop widely grown for oil and fiber. Despite recent advances in genomics, detailed gene activities during the important reproductive phase of its development are not well defined. In this study, we employed high-throughput RNA-sequencing methods to generate in-depth transcriptome profiles of flax tissues with emphasis on the reproductive phases of five key stages of embryogenesis (globular embryo, heart embryo, torpedo embryo, cotyledon embryo, and mature embryo), mature seed, and vegetative tissues viz. ovary, anther, and root. These datasets were used to establish the co-expression networks covering 36 gene modules based on the expression patterns for each gene through weighted gene co-expression network analysis (WGCNA). Functional interrogation with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) of dominantly expressed genetic modules in tissues revealed pathways involved in the development of different tissues. Moreover, the essential genes in embryo development and synthesis of storage reserves were identified based on their dynamic expression patterns. Together, this comprehensive dataset for developing embryos, mature seeds and vegetative tissues provides new insights into molecular mechanisms of seed development with potential for flax crop improvement.

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.

Identification of Stable and Multiple Environment Interaction QTLs and Candidate Genes for Fiber Productive Traits Under Irrigated and Water Stress Conditions Using Intraspecific RILs of Gossypium hirsutum var. MCU5 X TCH1218

Cotton productivity under water-stressed conditions is controlled by multiple quantitative trait loci (QTL). Enhancement of these productivity traits under water deficit stress is crucial for the genetic improvement of upland cotton, Gossypium hirsutum. In the present study, we constructed a genetic map with 504 single nucleotide polymorphisms (SNPs) covering a total span length of 4,416 cM with an average inter-marker distance of 8.76 cM. A total of 181 intra-specific recombinant inbred lines (RILs) were derived from a cross between G. hirsutum var. MCU5 and TCH1218 were used. Although 2,457 polymorphic SNPs were detected between the parents using the CottonSNP50K assay, only 504 SNPs were found to be useful for the construction of the genetic map. In the SNP genotyping, a large number of SNPs showed either >20% missing data, duplication, or segregation distortion. However, the mapped SNPs of this study showed collinearity with the physical map of the reference genome (G. hirsutum var.TM-1), indicating that there was no chromosomal rearrangement within the studied mapping population. RILs were evaluated under multi-environments and seasons for which the phenotypic data were acquired. A total of 53 QTL controlling plant height (PH), number of sympodial branches, boll number (BN), and boll weight (BW) were dissected by QTL analysis under irrigated and water stress conditions. Additionally, it was found that nine QTL hot spots not only co-localized for more than one investigated trait but were also stable with major QTL, i.e., with > 10% of phenotypic variation. One QTL hotspot on chromosome 22 flanked by AX-182254626-AX-182264770 with a span length of 89.4 cM co-localized with seven major and stable QTL linked to a number of sympodial branches both under irrigated and water stress conditions. In addition, putative candidate genes associated with water stress in the QTL hotspots were identified. Besides, few QTL from the hotspots were previously reported across various genetic architects in cotton validating the potential applications of these identified QTL for cotton breeding and improvement. Thus, the major and stable QTL identified in the present study would improve the cotton productivity under water-limited environments through marker-assisted selection.

Transgenic manipulation of triacylglycerol biosynthetic enzymes in B. napus alters lipid-associated gene expression and lipid metabolism

Oilseed rape (Brassica napus) is an important crop that is cultivated for the oil (mainly triacylglycerol; TAG) it produces in its seeds. TAG synthesis is controlled mainly by key enzymes in the Kennedy pathway, such as glycerol 3-phosphate acyltransferase (GPAT), lysophosphatidate acyltransferase (LPAT) and diacylglycerol acyltransferase (DGAT) but can also be produced from phosphoglycerides such as phosphatidylcholine (PC) by the activity of the enzyme phospholipid: diacylglycerol acyltransferase (PDAT). To evaluate the potential for these enzymes to alter oil yields or composition, we analysed transgenic B. napus lines which overexpressed GPAT, LPAT or PDAT using heterologous transgenes from Arabidopsis and Nasturtium and examined lipid profiles and changes in gene expression in these lines compared to WT. Distinct changes in PC and TAG abundance and spatial distribution in embryonic tissues were observed in some of the transgenic lines, together with altered expression of genes involved generally in acyl-lipid metabolism. Overall our results show that up-regulation of these key enzymes differentially affects lipid composition and distribution as well as lipid-associated gene expression, providing important information which could be used to improve crop properties by metabolic engineering.

Plant GDSL Esterases/Lipases: Evolutionary, Physiological and Molecular Functions in Plant Development

GDSL esterases/lipases (GELPs), present throughout all living organisms, have been a very attractive research subject in plant science due mainly to constantly emerging properties and functions in plant growth and development under both normal and stressful conditions. This review summarizes the advances in research on plant GELPs in several model plants and crops, including Arabidopsis, rice, maize and tomato, while focusing on the roles of GELPs in regulating plant development and plant-environment interactions. In addition, the possible regulatory network and mechanisms of GELPs have been discussed.

Transcriptional regulation of cell growth and reprogramming of systemic response in wheat (Triticum turgidum subsp. durum) seedlings by Bacillus paralicheniformis TRQ65

Bacillus paralicheniformis TRQ65 reprograms the gene expression patterns associated with systemic response to potentially facilitate its colonization and stimulate cell growth and plant biomass. Plant growth-promoting rhizobacteria (PGPR) carry out numerous mechanisms that enhance growth in seedlings, such as nutrient solubilization, phytohormone production, biocontrol activity, and regulation of induced systemic resistance (ISR) and acquired systemic resistance (ASR). Bacillus paralicheniformis TRQ65 is a biological and plant growth-promoting bacterium isolated from wheat (Triticum turgidum subsp. durum) rhizosphere. In this study, we performed a transcriptomic analysis of wheat seedlings inoculated with the native rhizobacterium Bacillus paralicheniformis TRQ65 (1 × 10⁷ cells∙g ^-1 of soil) at early development stages (GS15). A morphometrical assay was carried out to confirm growth promotion and after the cultivation period, TRQ65 was re-isolated to define inoculum persistence. Inoculated seedlings showed a significant (P < 0.05) increase in shoot length (93.48%) and dry weight in both shoot (117.02%) and root (48.33%) tissues; also, the strain persisted in the soil at 1.4 × 10⁷ UFC∙g^-1 of soil. A total of 228 differentially expressed genes (DEGs) (FDR < 0.05 and |log2 fold change|≥ 1.3) were observed in response to TRQ65 inoculation, of which 185 were down-regulated and 43 were up-regulated. The transcriptional patterns were characterized by the regulation of multidimensional cell growth (ROS, Ca⁺² channel, and NADPH oxidases activity), suppression of defense mechanism (PR proteins, PDFs, ROS, transcription factors), induction of central stimuli receptors (RALF, WAK, MAPK), carbohydrate metabolism (invertase activity) and phytohormone-related transport (ABCG transporter and AAAP). These results suggest that B. paralicheniformis TRQ65 is a promising bioinoculant agent for increasing wheat growth and development by reprogramming ISR and ASR simultaneously, suppressing defense mechanisms and inducing central stimuli response.

Modulation of drought adversities in Vicia faba by the application of plant growth promoting rhizobacteria and biochar

Drought is the greatest threat to world food security, seen as the catalyst for the great famines of the past. Given that the world's water supply is limited, it is likely that future demand of food for increasing population will further exacerbate the drought effects. Therefore, the present study was aimed to investigate the effect of biochar and plant growth promoting rhizobacteria (PGPR) Sphingobacterium pakistanensis (NCCP246) and Cellulomonas pakistanensis (NCCP11) on agronomic and physiological attributes of Vicia faba two varieties Desi (V1) and Pulista (V2) under induced drought stress. The seeds were sown in earthen pots filled with 3 kg sand and soil (1:2), and biochar (0 and 5% w/w) in triplicate arranged in complete randomized design. Analysis of biochar possessed 0.49 g cm^-3 bulk density, 9.6 pH; 5.4 cmol kg^-1 cation exchange capacity, 3.64% organic carbon and EC 6.7 ds/m. Agronomic attributes including seed LAI, LAR, SVI, %PHSI and RWC were improved by 30.4-180.4%, 14.37-47.20%, 37.64-50.91%, 18.21-30.80, and 35.82-54.34% in both varieties by the co-application of biochar and PGPR. Stomatal physiology and epidermal vigor was successfully improved by the application of PGPR and biochar as analyzed by scanning electron microscopy (SEM). Photosynthetic pigments, flavonoids, phenols, proline and glycine betaine were amplified by 58.33-173.8%, 50.59-130.33%, 46.58-86.62%, 46.66-109.30%, 35.74-56.10%, and 21.96-77.22% in both varieties by the co-application of biochar and PGPR. So, the present work concluded that, combined application of biochar and PGPR could be an effective strategy to alleviate the adversities of drought in V. faba growing in drastic ecosystems.

Genome-Wide Identification of GDSL-Type Esterase/Lipase Gene Family in Dasypyrum villosum L. Reveals That DvGELP53 Is Related to BSMV Infection

GDSL-type esterase/lipase proteins (GELPs) characterized by a conserved GDSL motif at their N-terminus belong to the lipid hydrolysis enzyme superfamily. In plants, GELPs play an important role in plant growth, development and stress response. The studies of the identification and characterization of the GELP gene family in Triticeae have not been reported. In this study, 193 DvGELPs were identified in Dasypyrum villosum and classified into 11 groups (clade A–K) by means of phylogenetic analysis. Most DvGELPs contain only one GDSL domain, only four DvGELPs contain other domains besides the GDSL domain. Gene structure analysis indicated 35.2% DvGELP genes have four introns and five exons. In the promoter regions of the identified DvGELPs, we detected 4502 putative cis-elements, which were associated with plant hormones, plant growth, environmental stress and light responsiveness. Expression profiling revealed 36, 44 and 17 DvGELPs were highly expressed in the spike, the root and the grain, respectively. Further investigation of a root-specific expressing GELP, DvGELP53, indicated it was induced by a variety of biotic and abiotic stresses. The knockdown of DvGELP53 inhibited long-distance movement of BSMV in the tissue of D. villosum. This research provides a genome-wide glimpse of the D. villosum GELP genes and hints at the participation of DvGELP53 in the interaction between virus and plants.

Genome-wide identification of the tobacco GDSL family and apical meristem-specific expression conferred by the GDSL promoter

BACKGROUND

GDSL esterases/lipases are a large protein subfamily defined by the distinct GDSL motif, and play important roles in plant development and stress responses. However, few studies have reported on the role of GDSLs in the growth and development of axillary buds. This work aims to identify the GDSL family members in tobacco and explore whether the NtGDSL gene contributes to development of the axillary bud in tobacco.

RESULTS

One hundred fifty-nine GDSL esterase/lipase genes from cultivated tobacco (Nicotiana tabacum) were identified, and the dynamic changes in the expression levels of 93 of these genes in response to topping, as assessed using transcriptome data of topping-induced axillary shoots, were analysed. In total, 13 GDSL esterase/lipase genes responded with changes in expression level. To identify genes and promoters that drive the tissue-specific expression in tobacco apical and axillary buds, the expression patterns of these 13 genes were verified using qRT-PCR. GUS activity and a lethal gene expression pattern driven by the NtGDSL127 promoter in transgenic tobacco demonstrated that NtGDSL127 is specifically expressed in apical buds, axillary buds, and flowers. Three separate deletions in the NtGDSL127 promoter demonstrated that a minimum upstream segment of 235 bp from the translation start site can drive the tissue-specific expression in the apical meristem. Additionally, NtGDSL127 responded to phytohormones, providing strategies for improving tobacco breeding and growth.

CONCLUSION

We propose that in tobacco, the NtGDSL127 promoter directs expression specifically in the apical meristem and that expression is closely correlated with axillary bud development.

Expression Pattern and Functional Analyses of Arabidopsis Guard Cell-Enriched GDSL Lipases

There are more than 100 GDSL lipases in Arabidopsis, but only a few members have been functionally investigated. Moreover, no reports have ever given a comprehensive analysis of GDSLs in stomatal biology. Here, we systematically investigated the expression patterns of 19 putative Guard-cell-enriched GDSL Lipases (GGLs) at various developmental stages and in response to hormone and abiotic stress treatments. Gene expression analyses showed that these GGLs had diverse expression patterns. Fifteen GGLs were highly expressed in guard cells, with seven preferentially in guard cells. Most GGLs were localized in endoplasmic reticulum, and some were also localized in lipid droplets and nucleus. Some closely homologous GGLs exhibited similar expression patterns at various tissues and in response to hormone and abiotic stresses, or similar subcellular localization, suggesting the correlation of expression pattern and biological function, and the functional redundancy of GGLs in plant development and environmental adaptations. Further phenotypic identification of ggl mutants revealed that GGL7, GGL14, GGL22, and GGL26 played unique and redundant roles in stomatal dynamics, stomatal density and morphology, and plant water relation. The present study provides unique resources for functional insights into these GGLs to control stomatal dynamics and development, plant growth, and adaptation to the environment.

Novel Genomic Regions of Fusarium Wilt Resistance in Bottle Gourd [Lagenaria siceraria (Mol.) Standl.] Discovered in Genome-Wide Association Study

Fusarium wilt (FW) is a typical soil-borne disease that seriously affects the yield and fruit quality of bottle gourd. Thus, to improve resistance to FW in bottle gourd, the genetic mechanism underlying FW resistance needs to be explored. In this study, we performed a genome-wide association study (GWAS) based on 5,330 single-nucleotide polymorphisms (SNPs) and 89 bottle gourd accessions. The GWAS results revealed a total of 10 SNPs (P ≤ 0.01, -log₁₀ P ≥ 2.0) significantly associated with FW resistance that were detected in at least two environments (2019DI, 2020DI, and the average across the 2 years); these SNPs were located on chromosomes 1, 2, 3, 4, 8, and 9. Linkage disequilibrium (LD) block structure analysis predicted three potential candidate genes for FW resistance. Genes HG_GLEAN_10001030 and HG_GLEAN_10001042 were within the range of the mean LD block of the marker BGReSe_14202; gene HG_GLEAN_10011803 was 280 kb upstream of the marker BGReSe_00818. Real-time quantitative PCR (qRT-PCR) analysis showed that HG_GLEAN_10011803 was significantly up-regulated in FW-infected plants of YD-4, Yin-10, and Hanbi; HG_GLEAN_10001030 and HG_GLEAN_10001042 were specifically up-regulated in FW-infected plants of YD-4. Therefore, gene HG_GLEAN_10011803 is likely the major effect candidate gene for resistance against FW in bottle gourd. This work provides scientific evidence for the exploration of candidate gene and development of functional markers in FW-resistant bottle gourd breeding programs.

Proteometabolomic characterization of apical bud maturation in Pinus pinaster

Bud maturation is a physiological process that implies a set of morphophysiological changes that lead to the transition of growth patterns from young to mature. This transition defines tree growth and architecture, and in consequence traits such as biomass production and wood quality. In Pinus pinaster Aiton, a conifer of great timber value, bud maturation is closely related to polycyclism (multiple growth periods per year). This process causes a lack of apical dominance, and consequently increased branching that reduces its timber quality and value. However, despite its importance, little is known about bud maturation. In this work, proteomics and metabolomics were employed to study apical and basal sections of young and mature buds in P. pinaster. Proteins and metabolites in samples were described and quantified using (n)UPLC-LTQ-Orbitrap. The datasets were analyzed employing an integrative statistical approach, which allowed the determination of the interactions between proteins and metabolites and the different bud sections and ages. Specific dynamics of proteins and metabolites such as histones H3 and H4, ribosomal proteins L15 and L12, chaperonin TCP1, 14-3-3 protein gamma, gibberellins A1, A3 and A8, strigolactones and abscisic acid, involved in epigenetic regulation, proteome remodeling, hormonal signaling and abiotic stress pathways showed their potential role during bud maturation. Candidates and pathways were validated employing interaction databases and targeted transcriptomics. These results increase our understanding of the molecular processes behind bud maturation, a key step towards improving timber production and natural pine forests management in a future scenario of climate change. However, further studies are necessary using different P. pinaster populations that show contrasting wood quality and stress tolerance in order to generalize the results.

Cutin:cutin-acid endo-transacylase (CCT), a cuticle-remodelling enzyme activity in the plant epidermis

Cutin is a polyester matrix mainly composed of hydroxy-fatty acids that occurs in the cuticles of shoots and root-caps. The cuticle, of which cutin is a major component, protects the plant from biotic and abiotic stresses, and cutin has been postulated to constrain organ expansion. We propose that, to allow cutin restructuring, ester bonds in this net-like polymer can be transiently cleaved and then re-formed (transacylation). Here, using pea epicotyl epidermis as the main model, we first detected a cutin:cutin-fatty acid endo-transacylase (CCT) activity. In-situ assays used endogenous cutin as the donor substrate for endogenous enzymes; the exogenous acceptor substrate was a radiolabelled monomeric cutin-acid, 16-hydroxy-[3H]hexadecanoic acid (HHA). High-molecular-weight cutin became ester-bonded to intact [3H]HHA molecules, which thereby became unextractable except by ester-hydrolysing alkalis. In-situ CCT activity correlated with growth rate in Hylotelephium leaves and tomato fruits, suggesting a role in loosening the outer epidermal wall during organ growth. The only well-defined cutin transacylase in the apoplast, CUS1 (a tomato cutin synthase), when produced in transgenic tobacco, lacked CCT activity. This finding provides a reference for future CCT protein identification, which can adopt our sensitive enzyme assay to screen other CUS1-related enzymes.

GDSL lipase occluded stomatal pore 1 is required for wax biosynthesis and stomatal cuticular ledge formation

The plant leaf surface is coated with a waterproof cuticle layer. Cuticle facing the stomatal pore surface needs to be sculpted to form outer cuticular ledge (OCL) after stomatal maturation for efficient gas exchange. Here, we characterized the roles of Arabidopsis GDSL lipase, Occlusion of Stomatal Pore 1 (OSP1), in wax biosynthesis and stomatal OCL formation. OSP1 mutation results in significant reduction in leaf wax synthesis and occlusion of stomata, leading to increased epidermal permeability, decreased transpiration rate, and enhanced drought tolerance. We demonstrated that OSP1 activity is critical for its role in wax biosynthesis and stomatal function. In vitro enzymatic assays demonstrated that OSP1 possesses thioesterase activity, particularly on C22:0 and C26:0 acyl-CoAs. Genetic interaction analyses with CER1 (ECERIFERUM 1), CER3 (ECERIFERUM 3) and MAH1 (Mid-chain Alkane Hydroxylase 1) in wax biosynthesis and stomatal OCL formation showed that OSP1 may act upstream of CER3 in wax biosynthesis, and implicate that wax composition percentage changes and keeping ketones in a lower level play roles, at least partially, in forming stomatal ledges. Our findings provided insights into the molecular mechanism mediating wax biosynthesis and highlighted the link between wax biosynthesis and the process of stomatal OCL formation.

A GDSL lipase-like from Ipomoea batatas catalyzes efficient production of 3,5-diCQA when expressed in Pichia pastoris

The synthesis of 3,5-dicaffeoylquinic acid (3,5-DiCQA) has attracted the interest of many researchers for more than 30 years. Recently, enzymes belonging to the BAHD acyltransferase family were shown to mediate its synthesis, albeit with notably low efficiency. In this study, a new enzyme belonging to the GDSL lipase-like family was identified and proven to be able to transform chlorogenic acid (5-O-caffeoylquinic acid, 5-CQA, CGA) in 3,5-DiCQA with a conversion rate of more than 60%. The enzyme has been produced in different expression systems but has only been shown to be active when transiently synthesized in Nicotiana benthamiana or stably expressed in Pichia pastoris. The synthesis of the molecule could be performed in vitro but also by a bioconversion approach beginning from pure 5-CQA or from green coffee bean extract, thereby paving the road for producing it on an industrial scale.

Miguel et al. identify a new enzyme belonging to the GDSL lipase-like family that is involved in the final stage of transformation of 5-CQA into 3,5-diCQA. This enzyme is able to realize an efficient transformation by over 60%, making the transformation process a valuable technological tool that can be easily transferred on an industrial scale.

Influence of Maternal Habitat on Salt Tolerance During Germination and Growth in Zygophyllum coccineum

Zygophyllum coccineum is a facultative halophyte widely distributed in desert wadis and coastal areas in Egypt. Here, we investigated the influences of maternal habitat on tolerance to salt stress during germination and seedling growth under salinity (0, 100, 200, 400 mM NaCl) of three populations of Z. coccineum from a saline habitat (Manzala coast) and non-saline habitats (Wadi Houf and Wadi Asyuti). In all populations, seed germination started within two days in distilled water but germination indices were reduced significantly with salt level increase. Germination percentage was not significantly greater for seeds from non-saline habitats than for those from the saline habitat under moderate salinity (100, 200 mM NaCl), but only seeds from the saline habitat were able to germinate under high salt stress (400 mM NaCl). Germination recovery was greater for seeds from the saline habitat compared to non-saline populations. At the seedling level, the Manzala population showed the lowest inhibition of shoot length and leaf area under salinity (200 and 400 mM NaCl) compared to non-saline habitats. In the same context, the Manzala population had the maximum chlorophyll a content, superoxide dismutase and esterase activities under salinity compared to non-saline populations, but salinity had a non-significant effect on chlorophyll b between the three populations. Carotenoids were enhanced with the increase of salt levels in all populations. These results suggest the salt tolerance of Manzala population is derived from maternal salinity and adaptive plasticity of this species may play an important role in the wide distribution of Z. coccineum.

Elevating seed oil content in a polyploid crop by induced mutations in SEED FATTY ACID REDUCER genes

Plant-based oils are valuable agricultural products, and seed oil content (SOC) is the major yield component in oil crops. Increasing SOC has been successfully targeted through the selection and genetic modification of oil biosynthesis. The SOC in rapeseed declined during the seed maturation and eventually caused the final accumulated seed oil quantity. However, genes involved in oil degradation during seed maturity are not deeply studied so far. We performed a candidate gene association study using a worldwide collection of rapeseed germplasm. We identified SEED FATTY ACID REDUCER (SFAR) genes, which had a significant effect on SOC and fatty acid (FA) composition. SFAR genes belong to the GDSL lipases, and GDSL lipases have a broad range of functions in plants. After quantification of gene expression using RNA-seq and quantitative PCR, we used targeted (CRISPR-Cas mediated) and random (chemical) mutagenesis to modify turnover rates of seed oil in winter rapeseed. For the first time, we demonstrate significant increase of SOC in a crop after knocking out members of the BnSFAR4 and BnSFAR5 gene families without pleiotropic effects on seed germination, vigour and oil mobilization. Our results offer new perspectives for improving oil yield by targeted mutagenesis.

IRREGULAR POLLEN EXINE2 Encodes a GDSL Lipase Essential for Male Fertility in Maize

Anther cuticle and pollen exine are two physical barriers protecting plant reproductive cells against environmental stresses; defects in either often cause male sterility. Here, we report the characterization of a male-sterile mutant irregular pollen exine2 (ipe2) of maize (Zea mays), which displays shrunken anthers and no starch accumulation in mature pollen grains. We cloned the causal gene IPE2 and confirmed its role in male fertility in maize with a set of complementary experiments. IPE2 is specifically expressed in maize developing anthers during stages 8 to 9 and encodes an endoplasmic-reticulum-localized GDSL lipase. Dysfunction of IPE2 resulted in delayed degeneration of tapetum and middle layer, leading to defective formation of anther cuticle and pollen exine, and complete male sterility. Aliphatic metabolism was greatly altered, with the contents of lipid constituents, especially C16/C18 fatty acids and their derivatives, significantly reduced in ipe2 developing anthers. Our study elucidates GDSL function in anther and pollen development and provides a promising genetic resource for breeding hybrid maize.

Multi-tissue transcriptome analysis using hybrid-sequencing reveals potential genes and biological pathways associated with azadirachtin A biosynthesis in neem (azadirachta indica)

BACKGROUND

Azadirachtin A is a triterpenoid from neem tree exhibiting excellent activities against over 600 insect species in agriculture. The production of azadirachtin A depends on extraction from neem tissues, which is not an eco-friendly and sustainable process. The low yield and discontinuous supply of azadirachtin A impedes further applications. The biosynthetic pathway of azadirachtin A is still unknown and is the focus of our study.

RESULTS

We attempted to explore azadirachtin A biosynthetic pathway and identified the key genes involved by analyzing transcriptome data from five neem tissues through the hybrid-sequencing (Illumina HiSeq and Pacific Biosciences Single Molecule Real-Time (SMRT)) approach. Candidates were first screened by comparing the expression levels between the five tissues. After phylogenetic analysis, domain prediction, and molecular docking studies, 22 candidates encoding 2,3-oxidosqualene cyclase (OSC), alcohol dehydrogenase, cytochrome P450 (CYP450), acyltransferase, and esterase were proposed to be potential genes involved in azadirachtin A biosynthesis. Among them, two unigenes encoding homologs of MaOSC1 and MaCYP71CD2 were identified. A unigene encoding the complete homolog of MaCYP71BQ5 was reported. Accuracy of the assembly was verified by quantitative real-time PCR (qRT-PCR) and full-length PCR cloning.

CONCLUSIONS

By integrating and analyzing transcriptome data from hybrid-seq technology, 22 differentially expressed genes (DEGs) were finally selected as candidates involved in azadirachtin A pathway. The obtained reliable and accurate sequencing data provided important novel information for understanding neem genome. Our data shed new light on understanding the biosynthesis of other triterpenoids in neem trees and provides a reference for exploring other valuable natural product biosynthesis in plants.

GDSL esterase/lipases OsGELP34 and OsGELP110/OsGELP115 are essential for rice pollen development

Pollen exine contains complex biopolymers of aliphatic lipids and phenolics. Abnormal development of pollen exine often leads to plant sterility. Molecular mechanisms regulating exine formation have been studied extensively but remain ambiguous. Here we report the analyses of three GDSL esterase/lipase protein genes, OsGELP34, OsGELP110, and OsGELP115, for rice exine formation. OsGELP34 was identified by cloning of a male sterile mutant gene. OsGELP34 encodes an endoplasmic reticulum protein and was mainly expressed in anthers during pollen exine formation. osgelp34 mutant displayed abnormal exine and altered expression of a number of key genes required for pollen development. OsGELP110 was previously identified as a gene differentially expressed in meiotic anthers. OsGELP110 was most homologous to OsGELP115, and the two genes showed similar gene expression patterns. Both OsGELP110 and OsGELP115 proteins were localized in peroxisomes. Individual knockout of OsGELP110 and OsGELP115 did not affect the plant fertility, but double knockout of both genes altered the exine structure and rendered the plant male sterile. OsGELP34 is distant from OsGELP110 and OsGELP115 in sequence, and osgelp34 and osgelp110/osgelp115 mutants were different in anther morphology despite both were male sterile. These results suggested that OsGELP34 and OsGELP110/OsGELP115 catalyze different compounds for pollen exine development.

Functional Omics Identifies Serine Hydrolases That Mobilize Storage Lipids during Rice Seed Germination

Elucidating proteolipidome dynamics is crucial for understanding the roles of these molecules in plant physiology and disease. Sequence-based functional annotation of the protein is inadequate, since protein activities depend on posttranslational modification. In this study, we applied a gel-free activity-based protein profiling approach to unravel the active lipases, including other Serine hydrolases (SHs), expressed during seed germination in rice (Oryza sativa). We successfully mapped the active sites of 43 active SHs encompassing lipases/esterases, GDSL lipases, proteases, Ser carboxypeptidases, ABHD protein, pectin acetylesterase, and other SHs. The mRNA expression levels of those genes encoding the identified SHs were monitored using microarray analysis. The lipidome analysis revealed distinct patterns of molecular species distribution in individual lipid classes and displayed the metabolic connections between lipid mobilization and rice seedling growth. Changes in the mobilization of storage lipids and their molecular species remodeling were correlated with the expression of the identified lipases and their lipase activity in a time-dependent manner. The physiological significance of the identified SHs was explored during biotic stress with Fusarium verticillioides infection. The fungal infection significantly reduced lipase activity and lipid mobilization, thus impairing the rice seedling. Collectively, our data demonstrate application of the functional proteome strategy along with the shotgun lipidome approach for the identification of active SHs, and thus for deciphering the role of lipid homeostasis during rice seed germination.

Activity-based protein profiling of rice (Oryza sativa L.) bran serine hydrolases

Rice bran is an underutilized agricultural by-product with economic importance. The unique phytochemicals and fatty acid compositions of bran have been targeted for nutraceutical development. The endogenous lipases and hydrolases are responsible for the rapid deterioration of rice bran. Hence, we attempted to provide the first comprehensive profiling of active serine hydrolases (SHs) present in rice bran proteome by activity-based protein profiling (ABPP) strategy. The active site-directed fluorophosphonate probe (rhodamine and biotin-conjugated) was used for the detection and identification of active SHs. ABPP revealed 55 uncharacterized active-SHs and are representing five different known enzyme families. Based on motif and domain analyses, one of the uncharacterized and miss annotated SHs (Os12Ssp, storage protein) was selected for biochemical characterization by overexpressing in yeast. The purified recombinant protein authenticated the serine protease activity in time and protein-dependent studies. Os12Ssp exhibited the maximum activity at a pH between 7.0 and 8.0. The protease activity was inhibited by the covalent serine protease inhibitor, which suggests that the ABPP approach is indeed reliable than the sequence-based annotations. Collectively, the comprehensive knowledge generated from this study would be useful in expanding the current understanding of rice bran SHs and paves the way for better utilization/stabilization of rice bran.

The DROOPING LEAF (DR) gene encoding GDSL esterase is involved in silica deposition in rice (Oryza sativa L.)

Leaf morphology is one of the most important agronomic traits in rice breeding because of its contribution to crop yield. The drooping leaf (dr) mutant was developed from the Ilpum rice cultivar by ethyl methanesulfonate (EMS) mutagenesis. Compared with the wild type, dr plants exhibited drooping leaves accompanied by a small midrib, short panicle, and reduced plant height. The phenotype of the dr plant was caused by a mutation within a single recessive gene on chromosome 2, dr (LOC_Os02g15230), which encodes a GDSL esterase. Analysis of wild-type and dr sequences revealed that the dr allele carried a single nucleotide substitution, glycine to aspartic acid. RNAi targeted to LOC_Os02g15230 produced same phenotypes to the dr mutation, confirming LOC_Os02g15230 as the dr gene. Microscopic observations and plant nutrient analysis of SiO₂ revealed that silica was less abundant in dr leaves than in wild-type leaves. This study suggests that the dr gene is involved in the regulation of silica deposition and that disruption of silica processes lead to drooping leaf phenotypes.

Physiological and genome-wide RNA-sequencing analyses identify candidate genes in a nitrogen-use efficient potato cv. Kufri Gaurav

Nitrogen (N) is an important nutrient for plant growth. However, its excess application leads to environmental damage. Hence, improving nitrogen use efficiency (NUE) of plant is one of the plausible options to solve the problems. Aim of this study was to identify candidate genes involved in enhancing NUE in potato cv. Kufri Gaurav (N efficient). Plants were grown in aeroponic with two contrasting N regimes (low N: 0.75 mM, and high N: 7.5 mM). Higher NUE in Kufri Gaurav was observed in low N based on the parameters like NUE, NUpE (N uptake efficiency), NUtE (N utilization efficiency) and AgNUE (agronomic NUE). Further, global gene expression profiles in root, leaf and stolon tissues were analyzed by RNA-sequencing using Ion Proton™ System. Quality data (≥Q20) of 2.04-2.73 Gb per sample were mapped with the potato genome. Statistically significant (P ≤ 0.05) differentially expressed genes (DEGs) were identified such as 176 (up-regulated) and 30 (down-regulated) in leaves, 39 (up-regulated) and 105 (down-regulated) in roots, and 81 (up-regulated) and 694 (down-regulated) in stolons. The gene ontology (GO) terms like metabolic process, cellular process and catalytic activity were predominant. Our RT-qPCR analysis confirmed the gene expression profiles of RNA-seq. Overall, we identified candidate genes associated with improving NUE such as superoxide dismutase, GDSL esterase lipase, probable phosphatase 2C, high affinity nitrate transporters, sugar transporter, proline rich proteins, transcription factors (VQ motif, SPX domain, bHLH) etc. Our findings suggest that these candidate genes probably play crucial roles in enhancing NUE in potato.

Genome-Wide Association Study Dissects the Genetic Architecture of Maize Husk Tightness

The husk is a leafy outer tissue that encloses a maize ear. Previously, we identified the optimum husk structure by measuring the husk length, husk layer number, husk thickness and husk width. Husk tightness (HTI) is a combined trait based on the above four husk measurements. Unveiling the genetic basis of HTI will aid in guiding the genetic improvement of maize for mechanical harvesting and for protecting the ear from pest damage and pathogen infection. Here, we used a maize associate population of 508 inbred lines with tropical, subtropical and temperate backgrounds to analyze the genetic architecture of HTI. Evaluating the phenotypic diversity in three different environments showed that HTI exhibited broad natural variations and a moderate heritability level of 0.41. A diversity analysis indicated that the inbred lines having a temperate background were more loosely related than those having a tropical or subtropical background. HTI showed significant negative correlations with husk thickness and width, which indicates that thicker and wider husks wrapped the ear tighter than thinner and slimmer husks. Combining husk traits with ∼1.25 million single nucleotide polymorphisms in a genome-wide association study revealed 27 variants that were significantly associated with HTI above the threshold of P < 7.26 × 10^-6. We found 27 candidate genes for HTI that may participate in (1) husk senescence involving lipid peroxidation (GRMZM2G017616) and programmed cell death (GRMZM2G168898 and GRMZM2G035045); (2) husk morphogenesis involving cell division (GRMZM5G869246) and cell wall architecture (GRMZM2G319798); and (3) cell signal transduction involving protein phosphorylation (GRMZM2G149277 and GRMZM2G004207) and the ABSISIC ACID INSENSITIVE3/VIVIPAROUS1 transcription factor (GRMZM2G088427). These results provide useful information for understanding the genetic basis of husk development. Further studies of identified candidate genes will help elucidate the molecular pathways that regulate HTI in maize.