Genome-wide analysis of GDSL-type esterases/lipases in Arabidopsis

TBL38 atypical homogalacturonan-acetylesterase activity and cell wall microdomain localization in Arabidopsis seed mucilage secretory cells

Plant cell walls constitute complex polysaccharidic/proteinaceous networks whose biosynthesis and dynamics implicate several cell compartments. The synthesis and remodeling of homogalacturonan pectins involve Golgi-localized methylation/acetylation and subsequent cell wall-localized demethylation/deacetylation. So far, TRICHOME BIREFRINGENCE-LIKE (TBL) family members have been described as Golgi-localized acetyltransferases targeting diverse hemicelluloses or pectins. Using seed mucilage secretory cells (MSCs) from Arabidopsis thaliana, we demonstrate the atypical localization of TBL38 restricted to a cell wall microdomain. A tbl38 mutant displays an intriguing homogalacturonan immunological phenotype in this cell wall microdomain and in an MSC surface-enriched abrasion powder. Mass spectrometry oligosaccharide profiling of this fraction reveals an increased homogalacturonan acetylation phenotype. Finally, TBL38 displays pectin acetylesterase activity in vitro. These results indicate that TBL38 is an atypical cell wall-localized TBL that displays a homogalacturonan acetylesterase activity rather than a Golgi-localized acetyltransferase activity as observed in previously studied TBLs. TBL38 function during seed development is discussed.

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.

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis

Introduction

An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security.

Methods

In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season.

Results

The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation.

Discussion

The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Evaluation of the Structure-Function Relationship of SGNH Lipase from Streptomyces rimosus by Site-Directed Mutagenesis and Computational Approach

Streptomyces rimosus extracellular lipase (SrL) is a multifunctional hydrolase belonging to the SGNH family. Here site-directed mutagenesis (SDM) was used for the first time to investigate the functional significance of the conserved amino acid residues Ser10, Gly54, Asn82, Asn213, and His216 in the active site of SrL. The hydrolytic activity of SrL variants was determined using para-nitrophenyl (pNP) esters with C4, C8, and C16 fatty acid chains. Mutation of Ser10, Asn82, or His216, but not Gly54, to Ala abolished lipase activity for all substrates. In contrast, the Asn213Ala variant showed increased enzymatic activity for C8 and C16 pNP esters. Molecular dynamics (MD) simulations showed that the interactions between the long alkyl chain substrate (C16) and Ser10 and Asn82 were strongest in Asn213Ala SrL. In addition to Asn82, Gly54, and Ser10, several new constituents of the substrate binding site were recognized (Lys28, Ser53, Thr89, and Glu212), as well as strong electrostatic interactions between Lys28 and Glu212. In addition to the H bonds Ser10-His216 and His216-Ser214, Tyr11 interacted strongly with Ser10 and His216 in all complexes with an active enzyme form. A previously unknown strong H bond between the catalytically important Asn82 and Gly54 was uncovered, which stabilizes the substrate in an orientation suitable for the enzyme reaction.

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.

Identification and validation of seed dormancy loci and candidate genes and construction of regulatory networks by WGCNA in maize introgression lines

KEY MESSAGE

Seventeen PHS-QTLs and candidate genes were obtained, including eleven major loci, three under multiple environments and two with co-localization by the other mapping methods; The functions of three candidate genes were validated using mutants; nine target proteins and five networks were filtered by joint analysis of GWAS and WGCNA. Seed dormancy (SD) and pre-harvest sprouting (PHS) affect yield, as well as grain and hybrid quality in seed production. Therefore, identification of genetic and regulatory pathways underlying PHS and SD is key to gene function analysis, allelic variation mining and genetic improvement. In this study, 78,360 SNPs by SLAF-seq of 230 maize chromosome segment introgression lines (ILs), PHS under five environments were used to conduct GWAS (genome wide association study) (a threshold of 1/n), and seventeen unreported PHS QTLs were obtained, including eleven QTLs with PVE > 10% and three QTLs under multiple environments. Two QTL loci were co-located between the other two genetic mapping methods. Using differential gene expression analyses at two stages of grain development, gene functional analysis of Arabidopsis mutants, and gene functional analysis in the QTL region, seventeen PHS QTL-linked candidate genes were identified, and their five molecular regulatory networks constructed. Based on the Arabidopsis T-DNA mutations, three candidate genes were shown to regulate for SD and PHS. Meanwhile, using RNA-seq of grain development, the weighted correlation network analysis (WGCNA) was performed, deducing five regulatory pathways and target genes that regulate PHS and SD. Based on the conjoint analysis of GWAS and WGCNA, four pathways, nine target proteins and target genes were revealed, most of which regulate cell wall metabolism, cell proliferation and seed dehydration tolerance. This has important theoretical and practical significance for elucidating the genetic basis of maize PHS and SD, as well as mining of genetic resources and genetic improvement of traits.

Genome-wide analysis of the GDSL esterase/lipase family genes in Physcomitrium patens and the involvement of GELP31 in spore germination

In plants, the ability to produce hydrophobic substances that would provide protection from dehydration was required for the transition to land. This genome-wide investigation outlines the evolution of GDSL-type esterase/lipase (GELP) proteins in the moss Physcomitrium patens and suggests possible functions of some genes. GELP proteins play roles in the formation of hydrophobic polymers such as cutin and suberin that protect against dehydration and pathogen attack. GELP proteins are also implicated in processes such as pollen development and seed metabolism and germination. The P. patens GELP gene family comprises 48 genes and 14 pseudogenes. Phylogenetic analysis of all P. patens GELP sequences along with vascular plant GELP proteins with reported functions revealed that the P. patens genes clustered within previously identified A, B and C clades. A duplication model predicting the expansion of the GELP gene family within the P. patens lineage was constructed. Expression analysis combined with phylogenetic analysis suggested candidate genes for functions such as defence against pathogens, cutin metabolism, spore development and spore germination. The presence of relatively fewer GELP genes in P. patens may reduce the occurrence of functional redundancy that complicates the characterization of vascular plant GELP genes. Knockout lines of GELP31, which is highly expressed in sporophytes, were constructed. Gelp31 spores contained amorphous oil bodies and germinated late, suggesting (a) role(s) of GELP31 in lipid metabolism in spore development or germination. Future knockout studies of other candidate GELP genes will further elucidate the relationship between expansion of the family and the ability to withstand the harsh land environment.

Cosuppression of AtGELP22 and AtGELP23, two ubiquitinated target proteins of RING E3 ligase AtAIRP5, increases tolerance to drought stress in Arabidopsis

AtAIRP5 RING E3 ubiquitin ligase was recently identified as a positive regulator of the abscisic acid (ABA)-mediated drought stress response by stimulating the degradation of serine carboxypeptidase-like 1. Here, we identified GDSL-type esterase/lipase 22 (AtGELP22) and AtGELP23 as additional interacting partners of AtAIRP5. Yeast two-hybrid, pull-down, co-immunoprecipitation, and ubiquitination analyses verified that AtGELP22 and AtGELP23 are ubiquitinated target proteins of AtAIRP5. AtGELP22 and AtGELP23 were colocalized with AtAIRP5 to punctate-like structures in the cytosolic fraction, in which PYK10 and NAI2, two ER body marker proteins, are localized. T-DNA insertion atgelp22 and atgelp23 single knockout mutant plants showed phenotypes indistinguishable from those of wild-type plants under ABA treatment. In contrast, RNAi-mediated cosuppression of AtGELP22 and AtGELP23 resulted in hypersensitive ABA-mediated stomatal movements and higher tolerance to drought stress than that of the single mutant and wild-type plants. Taken together, our results suggest that the putative GDSL-type esterases/lipases AtGELP22 and AtGELP23 act as redundant negative regulators of the ABA-mediated drought stress response in Arabidopsis.

Transcriptome Analysis Reveals Brassinolide Signaling Pathway Control of Foxtail Millet Seedling Starch and Sucrose Metabolism under Freezing Stress, with Implications for Growth and Development

Low-temperature stress limits the growth and development of foxtail millet. Freezing stress caused by sudden temperature drops, such as late-spring coldness, often occurs in the seedling stage of foxtail millet. However, the ability and coping strategies of foxtail millet to cope with such stress are not clear. In the present study, we analyzed the self-regulatory mechanisms of freezing stress in foxtail millet. We conducted a physiological study on foxtail millet leaves at -4 °C for seven different durations (0, 2, 4, 6, 8, 10, and 12 h). Longer freezing time increased cell-membrane damage, relative conductance, and malondialdehyde content. This led to osmotic stress in the leaves, which triggered an increase in free proline, soluble sugar, and soluble protein contents. The increases in these substances helped to reduce the damage caused by stress. The activities of superoxide dismutase, peroxidase, and catalase increased reactive oxygen species (ROS) content. The optimal time point for the response to freezing stress was 8 h after exposure. The transcriptome analysis of samples held for 8 h at -4 °C revealed 6862 differentially expressed genes (DEGs), among which the majority are implicated in various pathways, including the starch and sucrose metabolic pathways, antioxidant enzyme pathways, brassinolide (BR) signaling pathway, and transcription factors, according to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. We investigated possible crosstalk between BR signals and other pathways and found that BR signaling molecules were induced in response to freezing stress. The beta-amylase (BAM) starch hydrolase signal was enhanced by the BR signal, resulting in the accelerated degradation of starch and the formation of sugars, which served as emerging ROS scavengers and osmoregulators to resist freezing stress. In conclusion, crosstalk between BR signal transduction, and both starch and sucrose metabolism under freezing stress provides a new perspective for improving freezing resistance in foxtail millet.

GDSL Esterase/Lipase GELP1 Involved in the Defense of Apple Leaves against Colletotrichum gloeosporioides Infection

GDSL esterases/lipases are a subclass of lipolytic enzymes that play critical roles in plant growth and development, stress response, and pathogen defense. However, the GDSL esterase/lipase genes involved in the pathogen response of apple remain to be identified and characterized. Thus, in this study, we aimed to analyze the phenotypic difference between the resistant variety, Fuji, and susceptible variety, Gala, during infection with C. gloeosporioides, screen for anti-disease-associated proteins in Fuji leaves, and elucidate the underlying mechanisms. The results showed that GDSL esterase/lipase protein GELP1 contributed to C. gloeosporioides infection defense in apple. During C. gloeosporioides infection, GELP1 expression was significantly upregulated in Fuji. Fuji leaves exhibited a highly resistant phenotype compared with Gala leaves. The formation of infection hyphae of C. gloeosporioides was inhibited in Fuji. Moreover, recombinant His:GELP1 protein suppressed hyphal formation during infection in vitro. Transient expression in Nicotiana benthamiana showed that GELP1-eGFP localized to the endoplasmic reticulum and chloroplasts. GELP1 overexpression in GL-3 plants increased resistance to C. gloeosporioides. MdWRKY15 expression was upregulated in the transgenic lines. Notably, GELP1 transcript levels were elevated in GL-3 after salicylic acid treatment. These results suggest that GELP1 increases apple resistance to C. gloeosporioides by indirectly regulating salicylic acid biosynthesis.

Novel mechanisms for the synthesis of important secondary metabolites in Ginkgo biloba seed revealed by multi-omics data

Although the detailed biosynthetic mechanism is still unclear, the unique secondary metabolites of Ginkgo biloba, including ginkgolic acids (GAs) and terpene trilactones, have attracted increasing attention for their potent medicinal, physiological and biochemical properties. In particular, GAs have shown great potential in the fields of antibacterial and insecticidal activities, making it urgent to elucidate their biosynthetic mechanism. In this study, we systematically revealed the landscape of metabolic-transcriptional regulation across continuous growth stages of G. biloba seeds (GBS) based on multi-omics mining and experimental verification, and successfully identified all major types of GAs and terpene trilactones along with more than a thousand kinds of other metabolites. The phenological changes and the essential gene families associated with these unique metabolites were analyzed in detail, and several potential regulatory factors were successfully identified based on co-expression association analysis. In addition, we unexpectedly found the close relationship between large introns and the biosynthesis of these secondary metabolites. These genes with large introns related to the synthesis of secondary metabolites showed higher gene expression and expression stability in different tissues or growth stages. Our results may provide a new perspective for the study of the regulatory mechanism of these unique secondary metabolites in GBS.

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.

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.

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.

Genetic Determinants of Fiber-Associated Traits in Flax Identified by Omics Data Integration

In this paper, we explore potential genetic factors in control of flax phenotypes associated with fiber by mining a collection of 306 flax accessions from the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. In total, 11 traits were assessed in the course of 3 successive years. A genome-wide association study was performed for each phenotype independently using six different single-locus models implemented in the GAPIT3 R package. Moreover, we applied a multivariate linear mixed model implemented in the GEMMA package to account for trait correlations and potential pleiotropic effects of polymorphisms. The analyses revealed a number of genomic variants associated with different fiber traits, implying the complex and polygenic control. All stable variants demonstrate a statistically significant allelic effect across all 3 years of the experiment. We tested the validity of the predicted variants using gene expression data available for the flax fiber studies. The results shed new light on the processes and pathways associated with the complex fiber traits, while the pinpointed candidate genes may be further used for marker-assisted selection.

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.

Gas Chromatography-Mass Spectrometry Metabolite Analysis Combined with Transcriptomic and Proteomic Provide New Insights into Revealing Cuticle Formation during Pepper Development

The cuticle plays an important role for the quality of pepper fruit. However, the molecular mechanism of cuticle formation in pepper fruit remains unclear. Our results showed that the wax was continuously accumulated during pepper development, while the cutin monomer first increased and then decreased. Hexadecanoic acid and 10,16-hydroxyhexadecanoic acid were the main components of wax and cutin, respectively. Combined with transcriptome and proteome, the formation patterns of wax and cutin polyester network for pepper cuticle was proposed. The 18 pairs of consistent expression genes and proteins involved in cuticle formation were revealed. Meanwhile, 12 key genes were screened from fatty acid biosynthesis, biosynthesis of unsaturated fatty acids, fatty acid elongation, cutin, suberine, and wax biosynthesis, glycerolipid metabolism, and transport pathway. This study would provide important candidate genes and theoretical basis for the molecular mechanism of cuticle formation, which is essential for the breeding of peppers.

bHLH010/089 Transcription Factors Control Pollen Wall Development via Specific Transcriptional and Metabolic Networks in Arabidopsis thaliana

The pollen wall is a specialized extracellular cell wall that protects male gametophytes from various environmental stresses and facilitates pollination. Here, we reported that bHLH010 and bHLH089 together are required for the development of the pollen wall by regulating their specific downstream transcriptional and metabolic networks. Both the exine and intine structures of bhlh010 bhlh089 pollen grains were severely defective. Further untargeted metabolomic and transcriptomic analyses revealed that the accumulation of pollen wall morphogenesis-related metabolites, including polysaccharides, glyceryl derivatives, and flavonols, were significantly changed, and the expression of such metabolic enzyme-encoding genes and transporter-encoding genes related to pollen wall morphogenesis was downregulated in bhlh010 bhlh089 mutants. Among these downstream target genes, CSLB03 is a novel target with no biological function being reported yet. We found that bHLH010 interacted with the two E-box sequences at the promoter of CSLB03 and directly activated the expression of CSLB03. The cslb03 mutant alleles showed bhlh010 bhlh089–like pollen developmental defects, with most of the pollen grains exhibiting defective pollen wall structures.

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.

Arabidopsis GELP7 functions as a plasma membrane-localized acetyl xylan esterase, and its overexpression improves saccharification efficiency

Acetyl substitution on the xylan chain is critical for stable interaction with cellulose and other cell wall polymers in the secondary cell wall. Xylan acetylation pattern is governed by Golgi and extracellular localized acetyl xylan esterase (AXE). We investigated the role of Arabidopsis clade Id from the GDSL esterase/lipase or GELP family in polysaccharide deacetylation. The investigation of the AtGELP7 T-DNA mutant line showed a decrease in stem esterase activity and an increase in stem acetyl content. We further generated overexpressor AtGELP7 transgenic lines, and these lines showed an increase in AXE activity and a decrease in xylan acetylation compared to wild-type plants. Therefore, we have named this enzyme as AtAXE1. The subcellular localization and immunoblot studies showed that the AtAXE1 enzyme is secreted out, associated with the plasma membrane and involved in xylan de-esterification post-synthesis. The cellulose digestibility was improved in AtAXE1 overexpressor lines without pre-treatment, after alkali and xylanases pre-treatment. Furthermore, we have also established that the AtGELP7 gene is upregulated in the overexpressor line of AtMYB46, a secondary cell wall specific transcription factor. This transcriptional regulation can drive AtGELP7 or AtAXE1 to perform de-esterification of xylan in a tissue-specific manner. Overall, these data suggest that AtGELP7 overexpression in Arabidopsis reduces xylan acetylation and improves digestibility properties of polysaccharides of stem lignocellulosic biomass.

Reprogramming of sorghum proteome in response to sugarcane aphid infestation

Sugarcane aphid (SCA; Melanaphis sacchari Zehntner) is a key piercing-sucking pest of sorghum (Sorghum bicolor) that cause significant yield losses. While feeding on host plants, complex signaling networks are invoked from recognition of insect attack to induction of plant defenses. Consequently, these signaling networks lead to the production of insecticidal compounds or limited access of nutrients to insects. Previously, several studies were published on the transcriptomics analysis of sorghum in response to SCA infestation, but no information is available on the physiological changes of sorghum at the proteome level. We used the SCA resistant sorghum genotype SC265 for the global proteomics analysis after 1 and 7 days of SCA infestation using the TMT-plex technique. Peptides matching a total of 4211 proteins were identified and 158 proteins were differentially expressed at day 1 and 7. Overall, proteome profiling of SC265 after SCA infestation at days 1 and 7 revealed the suppression of plant defense-related proteins and upregulation of plant defense and signaling-related proteins, respectively. The plant defense responses based on proteome data were validated using electrical penetration graph (EPG) technique to observe changes in aphid feeding. Feeding behavior analyses revealed that SCA spent significantly longer time in phloem phase on SCA infested plants for day 1 and lesser time in day 7 SCA infested sorghum plants, compared to their respective control plants. Overall, our study provides insights into underlying mechanisms that contribute to sorghum resistance to SCA.

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.

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.

Assembled and annotated 26.5 Gbp coast redwood genome: a resource for estimating evolutionary adaptive potential and investigating hexaploid origin

Sequencing, assembly, and annotation of the 26.5 Gbp hexaploid genome of coast redwood (Sequoia sempervirens) was completed leading toward discovery of genes related to climate adaptation and investigation of the origin of the hexaploid genome. Deep-coverage short-read Illumina sequencing data from haploid tissue from a single seed were combined with long-read Oxford Nanopore Technologies sequencing data from diploid needle tissue to create an initial assembly, which was then scaffolded using proximity ligation data to produce a highly contiguous final assembly, SESE 2.1, with a scaffold N50 size of 44.9 Mbp. The assembly included several scaffolds that span entire chromosome arms, confirmed by the presence of telomere and centromere sequences on the ends of the scaffolds. The structural annotation produced 118,906 genes with 113 containing introns that exceed 500 Kbp in length and one reaching 2 Mb. Nearly 19 Gbp of the genome represented repetitive content with the vast majority characterized as long terminal repeats, with a 2.9:1 ratio of Copia to Gypsy elements that may aid in gene expression control. Comparison of coast redwood to other conifers revealed species-specific expansions for a plethora of abiotic and biotic stress response genes, including those involved in fungal disease resistance, detoxification, and physical injury/structural remodeling and others supporting flavonoid biosynthesis. Analysis of multiple genes that exist in triplicate in coast redwood but only once in its diploid relative, giant sequoia, supports a previous hypothesis that the hexaploidy is the result of autopolyploidy rather than any hybridizations with separate but closely related conifer species.

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.

Identification of heat responsive genes in pea stipules and anthers through transcriptional profiling

Field pea (Pisum sativum L.), a cool-season legume crop, is known for poor heat tolerance. Our previous work identified PR11-2 and PR11-90 as heat tolerant and susceptible lines in a recombinant inbred population. CDC Amarillo, a Canadian elite pea variety, was considered as another heat tolerant variety based on its similar field performance as PR11-2. This study aimed to characterize the differential transcription. Plants of these three varieties were stressed for 3 h at 38°C prior to self-pollination, and RNAs from heat stressed anthers and stipules on the same flowering node were extracted and sequenced via the Illumina NovaSeq platform for the characterization of heat responsive genes. In silico results were further validated by qPCR assay. Differentially expressed genes (DEGs) were identified at log2 |fold change (FC)| ≥ 2 between high temperature and control temperature, the three varieties shared 588 DEGs which were up-regulated and 220 genes which were down-regulated in anthers when subjected to heat treatment. In stipules, 879 DEGs (463/416 upregulation/downregulation) were consistent among varieties. The above heat-induced genes of the two plant organs were related to several biological processes i.e., response to heat, protein folding and DNA templated transcription. Ten gene ontology (GO) terms were over-represented in the consistently down-regulated DEGs of the two organs, and these terms were mainly related to cell wall macromolecule metabolism, lipid transport, lipid localization, and lipid metabolic processes. GO enrichment analysis on distinct DEGs of individual pea varieties suggested that heat affected biological processes were dynamic, and variety distinct responses provide insight into molecular mechanisms of heat-tolerance response. Several biological processes, e.g., cellular response to DNA damage stimulus in stipule, electron transport chain in anther that were only observed in heat induced PR11-2 and CDC Amarillo, and their relevance to field pea heat tolerance is worth further validation.

Genetic control of arbuscular mycorrhizal colonization by Rhizophagus intraradices in Helianthus annuus (L.)

Plant symbiosis with arbuscular mycorrhizal (AM) fungi provides many benefits, including increased nutrient uptake, drought tolerance, and belowground pathogen resistance. To develop a better understanding of the genetic architecture of mycorrhizal symbiosis, we conducted a genome-wide association study (GWAS) of this plant-fungal interaction in cultivated sunflower. A diversity panel of cultivated sunflower (Helianthus annuus L.) was phenotyped for root colonization under inoculation with the AM fungus Rhizophagus intraradices. Using a mixed linear model approach with a high-density genetic map, we identified genomic regions that are likely associated with R. intraradices colonization in sunflower. Additionally, we used a set of twelve diverse lines to assess the effect that inoculation with R. intraradices has on dried shoot biomass and macronutrient uptake. Colonization among lines in the mapping panel ranged from 0-70% and was not correlated with mycorrhizal growth response, shoot phosphorus response, or shoot potassium response among the Core 12 lines. Association mapping yielded three single-nucleotide polymorphisms (SNPs) that were significantly associated with R. intraradices colonization. This is the first study to use GWAS to identify genomic regions associated with AM colonization in an Asterid eudicot species. Three genes of interest identified from the regions containing these SNPs are likely related to plant defense.

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.

Genome-Wide Transcriptomic and Proteomic Exploration of Molecular Regulations in Quinoa Responses to Ethylene and Salt Stress

Quinoa (Chenopodiumquinoa Willd.), originated from the Andean region of South America, shows more significant salt tolerance than other crops. To reveal how the plant hormone ethylene is involved in the quinoa responses to salt stress, 4-week-old quinoa seedlings of 'NL-6' treated with water, sodium chloride (NaCl), and NaCl with ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) were collected and analyzed by transcriptional sequencing and tandem mass tag-based (TMT) quantitative proteomics. A total of 9672 proteins and 60,602 genes was identified. Among them, the genes encoding glutathione S-transferase (GST), peroxidase (POD), phosphate transporter (PT), glucan endonuclease (GLU), beta-galactosidase (BGAL), cellulose synthase (CES), trichome birefringence-like protein (TBL), glycine-rich cell wall structural protein (GRP), glucosyltransferase (GT), GDSL esterase/lipase (GELP), cytochrome P450 (CYP), and jasmonate-induced protein (JIP) were significantly differentially expressed. Further analysis suggested that the genes may mediate through osmotic adjustment, cell wall organization, reactive oxygen species (ROS) scavenging, and plant hormone signaling to take a part in the regulation of quinoa responses to ethylene and salt stress. Our results provide a strong foundation for exploration of the molecular mechanisms of quinoa responses to ethylene and salt stress.

RabA2b Overexpression Alters the Plasma-Membrane Proteome and Improves Drought Tolerance in Arabidopsis

Rab proteins are small GTPases that are important in the regulation of vesicle trafficking. Through data mining, we identified RabA2b to be stress responsive, though little is known about the involvement of RabA in plant responses to abiotic stresses. Analysis of the RabA2b native promoter showed strong activity during osmotic stress, which required the stress hormone Abscisic acid (ABA) and was restricted to the vasculature. Sequence analysis of the promoter region identified predicted binding motifs for several ABA-responsive transcription factors. We cloned RabA2b and overexpressed it in Arabidopsis. The resulting transgenic plants were strikingly drought resistant. The reduced water loss observed in detached leaves of the transgenic plants could not be explained by stomatal aperture or density, which was similar in all the genotypes. Subcellular localization studies detected strong colocalization between RabA2b and the plasma membrane (PM) marker PIP2. Further studies of the PM showed, for the first time, a distinguished alteration in the PM proteome as a result of RabA2b overexpression. Proteomic analysis of isolated PM fractions showed enrichment of stress-coping proteins as well as cell wall/cuticle modifiers in the transgenic lines. Finally, the cuticle permeability of transgenic leaves was significantly reduced compared to the wild type, suggesting that it plays a role in its drought resistant properties. Overall, these data provide new insights into the roles and modes of action of RabA2b during water stresses, and indicate that increased RabA2b mediated PM trafficking can affect the PM proteome and increase drought tolerance.

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.

Integrating a genome-wide association study with transcriptomic data to predict candidate genes and favourable haplotypes influencing Brassica napus seed phytate

Phytate is the storage form of phosphorus in angiosperm seeds and plays vitally important roles during seed development. However, in crop plants phytate decreases bioavailability of seed-sourced mineral elements for humans, livestock and poultry, and contributes to phosphate-related water pollution. However, there is little knowledge about this trait in oilseed rape (Brassica napus). Here, a panel of 505 diverse B. napus accessions was screened in a genome-wide association study (GWAS) using 3.28 × 10⁶ single-nucleotide polymorphisms (SNPs). This identified 119 SNPs significantly associated with phytate concentration (PA_Conc) and phytate content (PA_Cont) and six candidate genes were identified. Of these, BnaA9.MRP5 represented the candidate gene for the significant SNP chrA09_5198034 (27 kb) for both PA_Cont and PA_Conc. Transcription of BnaA9.MRP5 in a low-phytate variety (LPA20) was significantly elevated compared with a high-phytate variety (HPA972). Association and haplotype analysis indicated that inbred lines carrying specific SNP haplotypes within BnaA9.MRP5 were associated with high- and low-phytate phenotypes. No significant differences in seed germination and seed yield were detected between low and high phytate cultivars examined. Candidate genes, favourable haplotypes and the low phytate varieties identified in this study will be useful for low-phytate breeding of B. napus.

Comparative Transcriptome and Weighted Gene Co-expression Network Analysis Identify Key Transcription Factors of Rosa chinensis 'Old Blush' After Exposure to a Gradual Drought Stress Followed by Recovery

Rose is one of the most fundamental ornamental crops, but its yield and quality are highly limited by drought. The key transcription factors (TFs) and co-expression networks during rose’s response to drought stress and recovery after drought stress are still limited. In this study, the transcriptomes of leaves of 2-year-old cutting seedlings of Rosa chinensis ‘Old Blush’ from three continuous droughted stages (30, 60, 90 days after full watering) and rewatering were analyzed using RNA sequencing. Weighted gene co-expression network analysis (WGCNA) was used to construct a co-expression network, which was associated with the physiological traits of drought response to discovering the hub TFs involved in drought response. More than 45 million high-quality clean reads were generated from the sample and used for comparison with the rose reference genome. A total of 46433 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that drought stress caused significant changes in signal transduction, plant hormones including ABA, auxin, brassinosteroid (BR), cytokinin, ethylene (ET), jasmonic acid (JA) and salicylic acid (SA), primary and secondary metabolism, and a certain degree of recovery after rewatering. Gene co-expression analysis identified 18 modules, in which four modules showed a high degree of correlation with physiological traits. In addition, 42 TFs including members of NACs, WRKYs, MYBs, AP2/ERFs, ARFs, and bHLHs with high connectivity in navajowhite1 and blue modules were screened. This study provides the transcriptome sequencing report of R. chinensis ‘Old Blush’ during drought stress and rewatering process. The study also identifies the response of candidate TFs to drought stress, providing guidelines for improving the drought tolerance of the rose through molecular breeding in the future.

Genetic dissection of natural variation in oilseed traits of camelina by whole-genome resequencing and QTL mapping

Camelina [Camelina sativa (L.) Crantz] is an oilseed crop in the Brassicaceae family that is currently being developed as a source of bioenergy and healthy fatty acids. To facilitate modern breeding efforts through marker-assisted selection and biotechnology, we evaluated genetic variation among a worldwide collection of 222 camelina accessions. We performed whole-genome resequencing to obtain single nucleotide polymorphism (SNP) markers and to analyze genomic diversity. We also conducted phenotypic field evaluations in two consecutive seasons for variations in key agronomic traits related to oilseed production such as seed size, oil content (OC), fatty acid composition, and flowering time. We determined the population structure of the camelina accessions using 161,301 SNPs. Further, we identified quantitative trait loci (QTL) and candidate genes controlling the above field-evaluated traits by genome-wide association studies (GWAS) complemented with linkage mapping using a recombinant inbred line (RIL) population. Characterization of the natural variation at the genome and phenotypic levels provides valuable resources to camelina genetic studies and crop improvement. The QTL and candidate genes should assist in breeding of advanced camelina varieties that can be integrated into the cropping systems for the production of high yield of oils of desired fatty acid composition.

Transcriptome analysis of Sonneratia caseolaris seedlings under chilling stress

Sonneratia caseolaris is a native mangrove species found in China. It is fast growing and highly adaptable for mangrove afforestation, but suffered great damage by chilling event once introduced to high latitude area. To understand the response mechanisms under chilling stress, physiological and transcriptomic analyses were conducted. The relative electrolyte conductivity, malondialdehyde (MDA) content, soluble sugar content and soluble protein content increased significantly under chilling stress. This indicated that S. caseolaris suffered great damage and increased the levels of osmoprotectants in response to the chilling stress. Gene expression comparison analysis of S. caseolaris leaves after 6 h of chilling stress was performed at the transcriptional scale using RNA-Seq. A total of 168,473 unigenes and 3,706 differentially expressed genes (DEGs) were identified. GO and KEGG enrichment analyses showed that the DEGs were mainly involved in carbohydrate metabolism, antioxidant enzyme, plant hormone signal transduction, and transcription factors (TFs). Sixteen genes associated with carbohydrate metabolism, antioxidant enzyme, phytohormones and TFs were selected for qRT-PCR verification, and they indicated that the transcriptome data were reliable. Our work provided a comprehensive review of the chilling response of S. caseolaris at both physiological and transcriptomic levels, which will prove useful for further studies on stress-responses in mangrove plants.

Metabolic Control of Gametophore Shoot Formation through Arginine in the Moss Physcomitrium patens

Shoot formation is accompanied by active cell proliferation and expansion, requiring that metabolic state adapts to developmental control. Despite the importance of such metabolic reprogramming, it remains unclear how development and metabolism are integrated. Here, we show that disruption of ANGUSTIFOLIA3 orthologs (PpAN3s) compromises gametophore shoot formation in the moss Physcomitrium patens due to defective cell proliferation and expansion. Trans-omics analysis reveals that the downstream activity of PpAN3 is linked to arginine metabolism. Elevating arginine level by chemical treatment leads to stunted gametophores and causes Ppan3 mutant-like transcriptional changes in the wild-type plant. Furthermore, ectopic expression of AtAN3 from Arabidopsis thaliana ameliorates the defective arginine metabolism and promotes gametophore formation in Ppan3 mutants. Together, these findings indicate that arginine metabolism is a key pathway associated with gametophore formation and provide evolutionary insights into the establishment of the shoot system in land plants through the integration of developmental and metabolic processes.

ABCG transporters export cutin precursors for the formation of the plant cuticle

The plant cuticle is deposited on the surface of primary plant organs, such as leaves, fruits, and floral organs, forming a diffusion barrier and protecting the plant against various abiotic and biotic stresses. Cutin, the structural polyester of the plant cuticle, is synthesized in the apoplast. Plasma-membrane-localized ATP-binding cassette (ABC) transporters of the G family have been hypothesized to export cutin precursors. Here, we characterize SlABCG42 of tomato representing an ortholog of AtABCG32 in Arabidopsis. SlABCG42 expression in Arabidopsis complements the cuticular deficiencies of the Arabidopsis pec1/abcg32 mutant. RNAi-dependent downregulation of both tomato genes encoding proteins highly homologous to AtABCG32 (SlABCG36 and SlABCG42) leads to reduced cutin deposition and formation of a thinner cuticle in tomato fruits. By using a tobacco (Nicotiana benthamiana) protoplast system, we show that AtABCG32 and SlABCG42 have an export activity for 10,16-dihydroxy hexadecanoyl-2-glycerol, a cutin precursor in vivo. Interestingly, also free ω-hydroxy hexadecanoic acid as well as hexadecanedioic acid were exported, furthering the research on the identification of cutin precursors in vivo and the respective mechanisms of their integration into the cutin polymer.

Plastidial and cytosolic thiol reductases participate in the control of stomatal functioning

Stomatal movements via the control of gas exchanges determine plant growth in relation to environmental stimuli through a complex signalling network involving reactive oxygen species that lead to post-translational modifications of Cys and Met residues, and alter protein activity and/or conformation. Thiol-reductases (TRs), which include thioredoxins, glutaredoxins (GRXs) and peroxiredoxins (PRXs), participate in signalling pathways through the control of Cys redox status in client proteins. Their involvement in stomatal functioning remains poorly characterized. By performing a mass spectrometry-based proteomic analysis, we show that numerous thiol reductases, like PRXs, are highly abundant in guard cells. When investigating various Arabidopsis mutants impaired in the expression of TR genes, no change in stomatal density and index was noticed. In optimal growth conditions, a line deficient in cytosolic NADPH-thioredoxin reductases displayed higher stomatal conductance and lower leaf temperature evaluated by thermal infrared imaging. In contrast, lines deficient in plastidial 2-CysPRXs or type-II GRXs exhibited compared to WT reduced conductance and warmer leaves in optimal conditions, and enhanced stomatal closure in epidermal peels treated with abscisic acid or hydrogen peroxide. Altogether, these data strongly support the contribution of thiol redox switches within the signalling network regulating guard cell movements and stomatal functioning.

Single-nucleus RNA and ATAC sequencing reveals the impact of chromatin accessibility on gene expression in Arabidopsis roots at the single-cell level

Similar to other complex organisms, plants consist of diverse and specialized cell types. The gain of unique biological functions of these different cell types is the consequence of the establishment of cell-type-specific transcriptional programs. As a necessary step in gaining a deeper understanding of the regulatory mechanisms controlling plant gene expression, we report the use of single-nucleus RNA sequencing (sNucRNA-seq) and single-nucleus assay for transposase accessible chromatin sequencing (sNucATAC-seq) technologies on Arabidopsis roots. The comparison of our single-nucleus transcriptomes to the published protoplast transcriptomes validated the use of nuclei as biological entities to establish plant cell-type-specific transcriptomes. Furthermore, our sNucRNA-seq results uncovered the transcriptomes of additional cell subtypes not identified by single-cell RNA-seq. Similar to our transcriptomic approach, the sNucATAC-seq approach led to the distribution of the Arabidopsis nuclei into distinct clusters, suggesting the differential accessibility of chromatin between groups of cells according to their identity. To reveal the impact of chromatin accessibility on gene expression, we integrated sNucRNA-seq and sNucATAC-seq data and demonstrated that cell-type-specific marker genes display cell-type-specific patterns of chromatin accessibility. Our data suggest that the differential chromatin accessibility is a critical mechanism to regulate gene activity at the cell-type level.

NtMYB12a acts downstream of sucrose to inhibit fatty acid accumulation by targeting lipoxygenase and SFAR genes in tobacco

MYB12 promotes flavonol biosynthesis in plants by targeting several early biosynthesis genes (EBGs) of this pathway. The transcriptions of these EBGs are also induced by sucrose signal. However, whether MYB12 is activated by sucrose signal and what the other roles MYB12 has in regulating plant metabolism are poorly understood. In this study, two NtMYB12 genes were cloned from Nicotiana tabacum. Both NtMYB12a and NtMYB12b are involved in regulating flavonoids biosynthesis in tobacco. NtMYB12a is further shown to inhibit the accumulation of fatty acid (FA) in tobacco leaves and seeds. Post-translational activation and chromatin immunoprecipitation assays demonstrate that NtMYB12a directly promotes the transcriptions of NtLOX6, NtLOX5, NtSFAR4 and NtGDSL2, which encode lipoxygenase (LOX) or SFAR enzymes catalyzing the degradation of FA. NtLOX6 and NtLOX5 are shown to prevent the accumulation of FA in the mature seeds and significantly reduced the percentage of polyunsaturated fatty acids (PUFAs) in tobacco. Sucrose stimulates the transcription of NtMYB12a, and loss function of NtMYB12a partially suppresses the decrease of FA content in tobacco seedlings caused by sucrose treatment. The regulation of sucrose on the expression of NtLOX6 and NtGDSL2 genes is mediated by NtMYB12a, whereas those of NtLOX5 and NtSFAR4 genes are independent of sucrose.

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.

Transcriptome analysis reveals mechanism of early ripening in Kyoho grape with hydrogen peroxide treatment

Background

In a previous study, the early ripening of Kyoho grape following H₂O₂ treatment was explored at the physiological level, but the mechanism by which H₂O₂ promotes ripening at the molecular level is unclear. To reveal the molecular mechanism, RNA-sequencing analysis was conducted on the different developmental stages of Kyoho berry treated with H₂O₂.

Results

In the comparison of treatment and control groups, 406 genes were up-regulated and 683 were down-regulated. Time course sequencing (TCseq) analysis showed that the expression patterns of most of the genes were similar between the treatment and control, except for some genes related to chlorophyll binding and photosynthesis. Differential expression analysis and the weighted gene co-expression network were used to screen significantly differentially expressed genes and hub genes associated with oxidative stress (heat shock protein, HSP), cell wall deacetylation (GDSL esterase/lipase, GDSL), cell wall degradation (xyloglucan endotransglucosylase/ hydrolase, XTH), and photosynthesis (chlorophyll a-b binding protein, CAB1). Gene expression was verified with RT-qPCR, and the results were largely consistent with those of RNA sequencing.

Conclusions

The RNA-sequencing analysis indicated that H₂O₂ treatment promoted the early ripening of Kyoho berry by affecting the expression levels of HSP, GDSL, XTH, and CAB1 and- photosynthesis- pathways.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07180-y.

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.

Evolutionary analysis of the Moringa oleifera genome reveals a recent burst of plastid to nucleus gene duplications

It is necessary to identify suitable alternative crops to ensure the nutritional demands of a growing global population. The genome of Moringa oleifera, a fast-growing drought-tolerant orphan crop with highly valuable agronomical, nutritional and pharmaceutical properties, has recently been reported. We model here gene family evolution in Moringa as compared with ten other flowering plant species. Despite the reduced number of genes in the compact Moringa genome, 101 gene families, grouping 957 genes, were found as significantly expanded. Expanded families were highly enriched for chloroplastidic and photosynthetic functions. Indeed, almost half of the genes belonging to Moringa expanded families grouped with their Arabidopsis thaliana plastid encoded orthologs. Microsynteny analysis together with modeling the distribution of synonymous substitutions rates, supported most plastid duplicated genes originated recently through a burst of simultaneous insertions of large regions of plastid DNA into the nuclear genome. These, together with abundant short insertions of plastid DNA, contributed to the occurrence of massive amounts of plastid DNA in the Moringa nuclear genome, representing 4.71%, the largest reported so far. Our study provides key genetic resources for future breeding programs and highlights the potential of plastid DNA to impact the structure and function of nuclear genes and genomes.

Molecular basis for neofunctionalization of duplicated E3 ubiquitin ligases underlying adaptation to drought tolerance in Arabidopsis thaliana

Multigene families in plants expanded from ancestral genes via gene duplication mechanisms constitute a significant fraction of the coding genome. Although most duplicated genes are lost over time, many are retained in the genome. Clusters of tandemly arrayed genes are commonly found in the plant genome where they can promote expansion of gene families. In the present study, promoter fusion to the GUS reporter gene was used to examine the promoter architecture of duplicated E3 ligase genes that are part of group C in the Arabidopsis thaliana ATL family. Acquisition of gene expression by AtATL78, possibly generated from defective AtATL81 expression, is described. AtATL78 expression was purportedly enhanced by insertion of a TATA box within the core promoter region after a short tandem duplication that occurred during evolution of Brassicaceae lineages. This gene is associated with an adaptation to drought tolerance of A. thaliana. These findings also suggest duplicated genes could serve as a reservoir of tacit genetic information, and expression of these duplicated genes is activated upon acquisition of core promoter sequences. Remarkably, drought transcriptome profiling in response to rehydration suggests that ATL78-dependent gene expression predominantly affects genes with root-specific activities.