Genome-wide analysis reveals the evolution and structural features of WRINKLED1 in plants

Genome-wide association and RNA-seq analyses reveal a potential gene related to linolenic acid in soybean seeds

Linolenic acid (LA) has poor oxidative stability since it is a polyunsaturated fatty acid. Soybean oil has a high LA content and thus has poor oxidative stability. To identify candidate genes that affect the linolenic acid (LA) content in soybean seeds, a genome-wide association study (GWAS) was performed with 1,060 soybean cultivars collected in China between 2019-2021 and which LA content was measured using matrix-assisted laser desorption/ionization time-of-flight imaging mass spectrometry (MALDI-TOF IMS). A candidate gene, GmWRI14, encoding an APETALA2 (AP2)-type transcription factor, was detected by GWAS in cultivars from all three study years. Multiple sequence alignments showed that GmWRI14 belongs to the plant WRI1 family. The fatty acid contents of different soybean lines were evaluated in transgenic lines with a copy of GmWRI14, control lines without GmWRI14, and the gmwri14 mutant. MALDI-TOF IMS revealed that GmWRI14 transgenic soybeans had a lower LA content with a significant effect on seed size and shape, whereas gmwri14 mutants had a higher LA content. compared to control. The RNA-seq results showed that GmWRI14 suppresses GmFAD3s (GmFAD3B and GmFAD3C) and GmbZIP54 expression in soybean seeds, leading to decreased LA content. Based on the RNA-seq data, yeast one-hybrid (Y1H) and qRT-PCR were performed to confirm the transcriptional regulation of FAD3s by GmWRI14. Our results suggest that FAD3 is indirectly regulated by GmWRI14, representing a new molecular mechanism of fatty acid biosynthesis, in which GmWRI14 regulates LA content in soybean seeds.

GmWRI1c Increases Palmitic Acid Content to Regulate Seed Oil Content and Nodulation in Soybean (Glycine max)

Soybean (Glycine max) is an important oil crop, but the regulatory mechanisms underlying seed oil accumulation remain unclear. We identified a member of the GmWRI1s transcription factor family, GmWRI1c, that is involved in regulating soybean oil content and nodulation. Overexpression of GmWRI1c in soybean hairy roots increased the expression of genes involved in glycolysis and de novo lipogenesis, the proportion of palmitic acid (16:0), and the number of root nodules. The effect of GmWRI1c in increasing the number of root nodules via regulating the proportion of palmitic acid was confirmed in a recombinant inbred line (RIL) population. GmWRI1c shows abundant sequence diversity and has likely undergone artificial selection during domestication. An association analysis revealed a correlation between seed oil content and five linked natural variations (Hap1/Hap2) in the GmWRI1c promoter region. Natural variations in the GmWRI1c promoter were strongly associated with the GmWRI1c transcript level, with higher GmWRI1c transcript levels in lines carrying GmWRI1c^Hap1 than in those carrying GmWRI1c^Hap2. The effects of GmWRI1c alleles on seed oil content were confirmed in natural and RIL populations. We identified a favourable GmWRI1c allele that can be used to breed new varieties with increased seed oil content and nodulation.

Genome-wide characterization of C2H2 zinc-finger gene family provides insight into the mechanisms and evolution of the dehydration-rehydration responses in Physcomitrium and Arabidopsis

Dehydration tolerance is a vital factor for land plant evolution and world agricultural production. Numerous studies enlightened that the plant-specific C2H2-type zinc-finger proteins (C2H2-ZFPs) as master regulators played pivotal roles in the abiotic stress responses of plants. However, a comprehensive understanding of the evolution of C2H2-ZFPs in terrestrial plants and its regulatory mechanism in dehydration and rehydration response remains a mystery. In this study, the genome-wide identification of C2H2-ZFP genes revealed 549 homologs in the representatives of terrestrial plant lineages from liverwort to angiosperms. Based on the characteristics of the conserved C2H2-ZF domains, four major C2H2-ZF types (M-, Z-, Q-, and D-type) were identified in the C2H2-ZFPs, with the dominants of M-type in all selected species and followed by Z-type in non-seed plants and Q-type in seed plants, respectively. Phylogenetic analyses of the identified C2H2-ZFPs supported four major groups in the land plant representatives, among which the members from the desiccation-tolerant Physcomitrium patens and the dehydration-sensitive Arabidopsis thaliana displayed different topological relationships in the phylogenies reconstructed for a single species. C2H2-ZFPs clustered in the same subclades shared similar features in their conserved domains and gene structures. Approximately, 81% of the C2H2-ZFP promoters of all 549 identified C2H2-ZFPs harbored the conserved ABA-responsive elements (ABREs) and/or dehydration-responsive elements (DREs). Comparative transcriptomic analyses showed that 50 PpZFPs and 56 AtZFPs significantly changed their transcripts abundance. Interestingly, most of the dehydration- and rehydration-responsive PpZPFs and AtZFPs had been predicted to contain the ABRE and DRE elements in their promoter regions and with over half of which phylogenetically belonging to group III. The differences in the expression patterns of C2H2-ZFPs in responses to dehydration and rehydration between P. patens and A. thaliana reflected their different strategies to adapt to dehydration. The identified candidate PpZFPs were specifically induced by moderate dehydration and reached the peak transcript abundance in severe dehydration. Our study lays the foundations for further functional investigation of C2H2-ZFPs in dehydration responses from an evolutionary perspective in land plants. The findings will provide us with genetic resources and potential targets for drought tolerance breeding in crops and beyond.

Overexpression of WRINKLED1 improves the weight and oil content in seeds of flax (Linum usitatissimum L.)

Seeds of flax (Linum usitatissimum L.) are highly rich in both oil and linolenic acid (LIN). It is crucial for flax agricultural production to identify positive regulators of fatty acid biosynthesis. In this study, we find that WRINKLED1 transcription factors play important positive roles during flax seed oil accumulation. Two WRINKLED1 genes, LuWRI1a and LuWRI1b, were cloned from flax, and LuWRI1a was found be expressed predominantly in developing seeds during maturation. Overexpression of LuWRI1a increased seed size, weight, and oil content in Arabidopsis and increased seed storage oil content in transgenic flax without affecting seed production or seed oil quality. The rise in oil content in transgenic flax seeds was primarily attributable to the increase in seed weight, according to a correlational analysis. Furthermore, overexpression or interference of LuWRI1a upregulated the expression of genes in the fatty acid biosynthesis pathway and LAFL genes, and the expression level of WRI1 was highly significantly positively associated between L1L, LEC1, and BCCP2. Our findings give a theoretical scientific foundation for the future application of genetic engineering to enhance the oil content of plant seeds.

Genome-Wide Association Study Identifies Candidate Genes Related to the Linoleic Acid Content in Soybean Seeds

Soybean (Glycine max (L.) Merrill) oil is a complex mixture of five fatty acids (palmitic, stearic, oleic, linoleic, and linolenic). The high content of linoleic acid (LA) contributes to the oil having poor oxidative stability. Therefore, soybean seed with a lower LA content is desirable. To investigate the genetic architecture of LA, we performed a genome-wide association study (GWAS) using 510 soybean cultivars collected from China. The phenotypic identification results showed that the content of LA varied from 36.22% to 72.18%. The GWAS analysis showed that there were 37 genes related to oleic acid content, with a contribution rate of 7%. The candidate gene Glyma.04G116500.1 (GmWRI14) on chromosome 4 was detected in three consecutive years. The GmWRI14 showed a negative correlation with the LA content and the correlation coefficient was -0.912. To test whether GmWRI14 can lead to a lower LA content in soybean, we introduced GmWRI14 into the soybean genome. Matrix-assisted laser desorption/ionization time-of-flight imaging mass spectrometry (MALDI-TOF IMS) showed that the overexpression of GmWRI14 leads to a lower LA content in soybean seeds. Meanwhile, RNA-seq verified that GmWRI14-overexpressed soybean lines showed a lower accumulation of GmFAD2-1A and GmFAD2-1B than control lines. Our results indicate that the down-regulation of the FAD2 gene triggered by the transcription factor GmWRI14 is the underlying mechanism reducing the LA level of seed. Our results provide novel insights into the genetic architecture of LA and pinpoint potential candidate genes for further in-depth studies.

Integrated Transcriptomic and Bioinformatics Analyses Reveal the Molecular Mechanisms for the Differences in Seed Oil and Starch Content Between Glycine max and Cicer arietinum

The seed oil and starch content of soybean are significantly different from that of chickpea. However, there are limited studies on its molecular mechanisms. To address this issue, we conducted integrated transcriptomic and bioinformatics analyses for species-specific genes and acyl-lipid-, starch-, and carbon metabolism-related genes. Among seven expressional patterns of soybean-specific genes, four were highly expressed at the middle- and late oil accumulation stages; these genes significantly enriched fatty acid synthesis and carbon metabolism, and along with common acetyl CoA carboxylase (ACCase) highly expressed at soybean middle seed development stage, common starch-degrading enzyme beta-amylase-5 (BAM5) was highly expressed at soybean early seed development stage and oil synthesis-related genes ACCase, KAS, KAR, ACP, and long-chain acyl-CoA synthetase (LACS) were co-expressed with WRI1, which may result in high seed oil content and low seed starch content in soybean. The common ADP-glucose pyrophosphorylase (AGPase) was highly expressed at chickpea middle seed development stage, along with more starch biosynthesis genes co-expressed with four-transcription-factor homologous genes in chickpea than in soybean, and the common WRI1 was not co-expressed with oil synthesis genes in chickpea, which may result in high seed starch content and low seed oil content in chickpea. The above results may be used to improve chickpea seed oil content in two ways. One is to edit CaWRI1 to co-express with oil synthesis-related genes, which may increase carbon metabolites flowing to oil synthesis, and another is to increase the expression levels of miRNA159 and miRNA319 to inhibit the expression of MYB33, which may downregulate starch synthesis-related genes, making more carbon metabolites flow into oil synthesis. Our study will provide a basis for future breeding efforts to increase the oil content of chickpea seeds.

Functional and Predictive Structural Characterization of WRINKLED2, A Unique Oil Biosynthesis Regulator in Avocado

WRINKLED1 (WRI1), a member of the APETALA2 (AP2) class of transcription factors regulates fatty acid biosynthesis and triacylglycerol (TAG) accumulation in plants. Among the four known Arabidopsis WRI1 paralogs, only WRI2 was unable to complement and restore fatty acid content in wri1-1 mutant seeds. Avocado (Persea americana) mesocarp, which accumulates 60-70% dry weight oil content, showed high expression levels for orthologs of WRI2, along with WRI1 and WRI3, during fruit development. While the role of WRI1 as a master regulator of oil biosynthesis is well-established, the function of WRI1 paralogs is poorly understood. Comprehensive and comparative in silico analyses of WRI1 paralogs from avocado (a basal angiosperm) with higher angiosperms Arabidopsis (dicot), maize (monocot) revealed distinct features. Predictive structural analyses of the WRI orthologs from these three species revealed the presence of AP2 domains and other highly conserved features, such as intrinsically disordered regions associated with predicted PEST motifs and phosphorylation sites. Additionally, avocado WRI proteins also contained distinct features that were absent in the nonfunctional Arabidopsis ortholog AtWRI2. Through transient expression assays, we demonstrated that both avocado WRI1 and WRI2 are functional and drive TAG accumulation in Nicotiana benthamiana leaves. We predict that the unique features and activities of ancestral PaWRI2 were likely lost in orthologous genes such as AtWRI2 during evolution and speciation, leading to at least partial loss of function in some higher eudicots. This study provides us with new targets to enhance oil biosynthesis in plants.

AtMYB92 enhances fatty acid synthesis and suberin deposition in leaves of Nicotiana benthamiana

Acyl lipids are important constituents of the plant cell. Depending on the cell type, requirements in acyl lipids vary greatly, implying a tight regulation of fatty acid and lipid metabolism. The discovery of the WRINKLED1 (WRI1) transcription factors, members of the AP2-EREBP (APETALA2-ethylene-responsive element binding protein) family, has emphasized the importance of transcriptional regulation for adapting the rate of acyl chain production to cell requirements. Here, we describe the identification of another activator of the fatty acid biosynthetic pathway, the Arabidopsis MYB92 transcription factor. This MYB and all the members of the subgroups S10 and S24 of MYB transcription factors can directly activate the promoter of BCCP2 that encodes a component of the fatty acid biosynthetic pathway. Two adjacent MYB cis-regulatory elements are essential for the binding and activation of the BCCP2 promoter by MYB92. Overexpression of MYB92 or WRI1 in Nicotiana benthamiana induces the expression of fatty acid biosynthetic genes but results in the accumulation of different types of acyl lipids. In the presence of WRI1, triacylglycerol biosynthetic enzymes coded by constitutively expressed genes efficiently channel the excess fatty acids toward reserve lipid accumulation. By contrast, MYB92 activates both fatty acid and suberin biosynthetic genes; hence, the remarkable increase in suberin monomers measured in leaves expressing MYB92. These results provide additional insight into the molecular mechanisms that control the biosynthesis of an important cell wall-associated acylglycerol polymer playing critical roles in plants.

C2H2 Zinc Finger Proteins: Master Regulators of Abiotic Stress Responses in Plants

Abiotic stresses such as drought and salinity are major environmental factors that limit crop yields. Unraveling the molecular mechanisms underlying abiotic stress resistance is crucial for improving crop performance and increasing productivity under adverse environmental conditions. Zinc finger proteins, comprising one of the largest transcription factor families, are known for their finger-like structure and their ability to bind Zn2+. Zinc finger proteins are categorized into nine subfamilies based on their conserved Cys and His motifs, including the Cys2/His2-type (C2H2), C3H, C3HC4, C2HC5, C4HC3, C2HC, C4, C6, and C8 subfamilies. Over the past two decades, much progress has been made in understanding the roles of C2H2 zinc finger proteins in plant growth, development, and stress signal transduction. In this review, we focus on recent progress in elucidating the structures, functions, and classifications of plant C2H2 zinc finger proteins and their roles in abiotic stress responses.

Molecular Basis of Plant Oil Biosynthesis: Insights Gained From Studying the WRINKLED1 Transcription Factor

Most plant species generate and store triacylglycerol (TAG) in their seeds, serving as a core supply of carbon and energy to support seedling development. Plant seed oils have a wide variety of applications, from being essential for human diets to serving as industrial renewable feedstock. WRINKLED1 (WRI1) transcription factor plays a central role in the transcriptional regulation of plant fatty acid biosynthesis. Since the discovery of Arabidopsis WRI1 gene (AtWRI1) in 2004, the function of WRI1 in plant oil biosynthesis has been studied intensively. In recent years, the identification of WRI1 co-regulators and deeper investigations of the structural features and molecular functions of WRI1 have advanced our understanding of the mechanism of the transcriptional regulation of plant oil biosynthesis. These advances also help pave the way for novel approaches that will better utilize WRI1 for bioengineering oil production in crops.

WRINKLED1 homologs highly and functionally express in oil-rich endosperms of oat and castor

Oat (Avena sativa) and castor (Ricinus communis) accumulate a large amount of lipids in their endosperms, however the molecular mechanism remains unknown. In this study, differences in oil regulators between oat and wheat (Triticum aestivum) as well as common features between oat and castor were tested by analyzing their transcriptomes with further q-PCR analysis. Results indicated that WRINKLED1 (WRI1) homologs and their target genes highly expressed in the endosperms of oat and castor, but not in the starchy endosperms of wheat. Expression pattern of WRI1s was in agreement with that of oil accumulation. Three AsWRI1s (AsWRI1a, AsWRI1b and AsWRI1c) and one RcWRI1 were identified in the endosperms of oat and castor, respectively. AsWRI1c lacks VYL motif, which is different from the other three WRI1s. Expressions of these four WRI1s all complemented the phenotypes of Arabidopsis wri1-1 mutant. Overexpression of these WRI1s in Arabidopsis and tobacco BY2 cells increased oil contents of seeds and total fatty acids of the cells, respectively. Moreover, this overexpression also resulted in up-regulations of WRI1 target genes, such as PKp-β1. Taken together, our results suggest that high and functional expression of WRI1 play a key role in the oil-rich endosperms and the VYL motif is dispensable for WRI1 function.

WRINKLED1, a "Master Regulator" in Transcriptional Control of Plant Oil Biosynthesis

A majority of plant species generate and accumulate triacylglycerol (TAG) in their seeds, which is the main resource of carbon and energy supporting the process of seedling development. Plant seed oils have broad ranges of uses, being not only important for human diets but also renewable feedstock of industrial applications. The WRINKLED1 (WRI1) transcription factor is vital for the transcriptional control of plant oil biosynthetic pathways. Since the identification of the Arabidopsis WRI1 gene (AtWRI1) fifteen years ago, tremendous progress has been made in understanding the functions of WRI1 at multiple levels, ranging from the identification of AtWRI1 target genes to location of the AtWRI1 binding motif, and from discovery of intrinsic structural disorder in WRI1 to fine-tuning of WRI1 modulation by post-translational modifications and protein-protein interactions. The expanding knowledge on the functional understanding of the WRI1 regulatory mechanism not only provides a clearer picture of transcriptional regulation of plant oil biosynthetic pathway, but also helps generate new strategies to better utilize WRI1 for developing novel oil crops.

Atypical Splicing Accompanied by Skipping Conserved Micro-exons Produces Unique WRINKLED1, An AP2 Domain Transcription Factor in Rice Plants

WRINKLED1 (WRI1), an AP2 domain transcription factor, is a master regulator of oil synthesis in plant seeds. Its closely related proteins (WRIs) are also involved in regulating the synthesis of fatty acids, which play a role in producing oils, membranes, and other important components in plants. We found two WRI1 genes, OsWRI1-1 and OsWRI1-2, and two additional WRI1 homologs, OsWRI3 and OsWRI4, in the rice genome. OsWRI1 was ubiquitously expressed in rice plants, including developing seeds. However, OsWRI3 was only significantly expressed in the leaf blade and OsWRI4 was not expressed at all. OsWRI1-1 contains amino acid sequence GCL instead of VYL, which is encoded by an independent 9-bp micro-exon that is conserved in many plant species. We found that the GCL sequence was produced by an atypical splicing accompanied by skipping of the micro-exon. Furthermore, OsWRI1-1 highly activates the transcription of the promoter for the biotin carboxyl transferase 2 gene in Arabidopsis, but its activity was reduced by amino acid replacement or deletion of the GCL sequence in a transient assay using Arabidopsis cells. Our results indicated that atypical splicing produced unique WRI1 in rice plants.