Arabidopsis FAD2 gene encodes the enzyme that is essential for polyunsaturated lipid synthesis

Arabidopsis diacylglycerol acyltransferase1 mutants require fatty acid desaturation for normal seed development

Diacylglycerol acyltransferase1 (DGAT1) is the major enzyme that synthesizes triacylglycerols (TAG) during Arabidopsis seed development. Mutant dgat1 seeds possess low oil content in addition to a high polyunsaturated fatty acid (PUFA) composition. Two genes encoding endoplasmic reticulum localized desaturase enzymes, fatty acid desaturase2 (FAD2) and fatty acid desaturase3 (FAD3), were upregulated in both dgat1-1 and dgat1-2 developing seeds. Crosses between both dgat1 mutant alleles and fad2-1 failed to generate plants homozygous for both dgat1 and fad2. Reciprocal crosses with wild-type plants demonstrated that both male and female dgat1 fad2 gametophytes were viable. Siliques from DGAT1/dgat1-1 fad2-1/fad2-1 and dgat1-1/dgat1-1 FAD2/fad2-1 possessed abnormal looking seeds that were arrested in the torpedo growth stage. Approximately 25% of the seeds exhibited this arrested phenotype, genetically consistent with them possessing the double homozygous dgat1 fad2 genotype. In contrast, double homozygous dgat1-1 fad3-2 mutant plants were viable. Seeds from these plants possessed higher levels of 18:2 while their fatty acid content was lower than dgat1 mutant controls. The results are consistent with a model where in the absence of DGAT1 activity, desaturation of fatty acids by FAD2 becomes essential to provide PUFA substrates for phospholipid:diacylglycerol acyltransferase (PDAT) to synthesize TAG. In a dgat1 fad2 mutant, seed development is aborted because TAG is unable to be synthesized by either DGAT1 or PDAT.

Regulation of seed traits in soybean

Soybean (Glycine max) is an essential economic crop that provides vegetative oil and protein for humans, worldwide. Increasing soybean yield as well as improving seed quality is of great importance. Seed weight/size, oil and protein content are the three major traits determining seed quality, and seed weight also influences soybean yield. In recent years, the availability of soybean omics data and the development of related techniques have paved the way for better research on soybean functional genomics, providing a comprehensive understanding of gene functions. This review summarizes the regulatory genes that influence seed size/weight, oil content and protein content in soybean. We also provided a general overview of the pleiotropic effect for the genes in controlling seed traits and environmental stresses. Ultimately, it is expected that this review will be beneficial in breeding improved traits in soybean.

Comparative transcriptomic and evolutionary analysis of FAD-like genes of Brassica species revealed their role in fatty acid biosynthesis and stress tolerance

BACKGROUND

Fatty acid desaturases (FADs) are involved in regulating plant fatty acid composition by adding double bonds to growing hydrocarbon chain. Apart from regulating fatty acid composition FADs are of great importance, and are involved in stress responsiveness, plant development, and defense mechanisms. FADs have been extensively studied in crop plants, and are broadly classed into soluble and non-soluble fatty acids. However, FADs have not yet been characterized in Brassica carinata and its progenitors.

RESULTS

Here we have performed comparative genome-wide identification of FADs and have identified 131 soluble and 28 non-soluble FADs in allotetraploid B. carinata and its diploid parents. Most soluble FAD proteins are predicted to be resided in endomembrane system, whereas FAB proteins were found to be localized in chloroplast. Phylogenetic analysis classed the soluble and non-soluble FAD proteins into seven and four clusters, respectively. Positive type of selection seemed to be dominant in both FADs suggesting the impact of evolution on these gene families. Upstream regions of both FADs were enriched in stress related cis-regulatory elements and among them ABRE type of elements were in abundance. Comparative transcriptomic data analysis output highlighted that FADs expression reduced gradually in mature seed and embryonic tissues. Moreover, under heat stress during seed and embryo development seven genes remained up-regulated regardless of external stress. Three FADs were only induced under elevated temperature whereas five genes were upregulated under Xanthomonas campestris stress suggesting their involvement in abiotic and biotic stress response.

CONCLUSIONS

The current study provides insights into the evolution of FADs and their role in B. carinata under stress conditions. Moreover, the functional characterization of stress-related genes would exploit their utilization in future breeding programs of B. carinata and its progenitors.

A mutant cotton fatty acid desaturase 2-1d allele causes protein mistargeting and altered seed oil composition

BACKGROUND

Cotton (Gossypium sp.) has been cultivated for centuries for its spinnable fibers, but its seed oil also possesses untapped economic potential if, improvements could be made to its oleic acid content.

RESULTS

Previous studies, including those from our laboratory, identified pima accessions containing approximately doubled levels of seed oil oleic acid, compared to standard upland cottonseed oil. Here, the molecular properties of a fatty acid desaturase encoded by a mutant allele identified by genome sequencing in an earlier analysis were analyzed. The mutant sequence is predicted to encode a C-terminally truncated protein lacking nine residues, including a predicted endoplasmic reticulum membrane retrieval motif. We determined that the mutation was caused by a relatively recent movement of a Ty1/copia type retrotransposon that is not found associated with this desaturase gene in other sequenced cotton genomes. The mutant desaturase, along with its repaired isozyme and the wild-type A-subgenome homoeologous protein were expressed in transgenic yeast and stably transformed Arabidopsis plants. All full-length enzymes efficiently converted oleic acid to linoleic acid. The mutant desaturase protein produced only trace amounts of linoleic acid, and only when strongly overexpressed in yeast cells, indicating that the missing C-terminal amino acid residues are not strictly required for enzyme activity, yet are necessary for proper subcellular targeting to the endoplasmic reticulum membrane.

CONCLUSION

These results provide the biochemical underpinning that links a genetic lesion present in a limited group of South American pima cotton accessions and their rare seed oil oleic acid traits. Markers developed to the mutant desaturase allele are currently being used in breeding programs designed to introduce this trait into agronomic upland cotton varieties.

The genome of Orychophragmus violaceus provides genomic insights into the evolution of Brassicaceae polyploidization and its distinct traits

Orychophragmus violaceus, referred to as "eryuelan" (February orchid) in China, is an early-flowering ornamental plant. The high oil content and abundance of unsaturated fatty acids in O. violaceus seeds make it a potential high-quality oilseed crop. Here, we generated a whole-genome assembly for O. violaceus using Nanopore and Hi-C sequencing technologies. The assembled genome of O. violaceus was ∼1.3 Gb in size, with 12 pairs of chromosomes. Through investigation of ancestral genome evolution, we determined that the genome of O. violaceus experienced a tetraploidization event from a diploid progenitor with the translocated proto-Calepineae karyotype. Comparisons between the reconstructed subgenomes of O. violaceus identified indicators of subgenome dominance, indicating that subgenomes likely originated via allotetraploidy. O. violaceus was phylogenetically close to the Brassica genus, and tetraploidy in O. violaceus occurred approximately 8.57 million years ago, close in time to the whole-genome triplication of Brassica that likely arose via an intermediate tetraploid lineage. However, the tetraploidization in Orychophragmus was independent of the hexaploidization in Brassica, as evidenced by the results from detailed phylogenetic analyses and comparisons of the break and fusion points of ancestral genomic blocks. Moreover, identification of multi-copy genes regulating the production of high-quality oil highlighted the contributions of both tetraploidization and tandem duplication to functional innovation in O. violaceus. These findings provide novel insights into the polyploidization evolution of plant species and will promote both functional genomic studies and domestication/breeding efforts in O. violaceus.

The coexpression of two desaturases provides an optimized reduction of saturates in camelina oil

Reducing the saturate content of vegetable oils is key to increasing their utility and adoption as a feedstock for the production of biofuels. Expression of either the FAT5 16 : 0-CoA desaturase from Caenorhabditis elegans, or an engineered cyanobacterial 16 : 0/18 : 0-glycerolipid desaturase, DES9*, in seeds of Arabidopsis (Arabidopsis thaliana) substantially lowered oil saturates. However, because pathway fluxes and regulation of oil synthesis are known to differ across species, translating this transgene technology from the model plant to crop species requires additional investigation. In the work reported here, we found that high expression of FAT5 in seeds of camelina (Camelina sativa) provided only a moderate decrease in saturates, from 12.9% of total oil fatty acids in untransformed controls to 8.6%. Expression of DES9* reduced saturates to 4.6%, but compromised seed physiology and oil content. However, the coexpression of the two desaturases together cooperatively reduced saturates to only 4.0%, less than one-third of the level in the parental line, without compromising oil yield or seedling germination and establishment. Our successful lowering of oil saturates in camelina identifies strategies that can now be integrated with genetic engineering approaches that reduce polyunsaturates to provide optimized oil composition for biofuels in camelina and other oil seed crops.

Sequencing of Camelina neglecta, a diploid progenitor of the hexaploid oilseed Camelina sativa

Camelina neglecta is a diploid species from the genus Camelina, which includes the versatile oilseed Camelina sativa. These species are closely related to Arabidopsis thaliana and the economically important Brassica crop species, making this genus a useful platform to dissect traits of agronomic importance while providing a tool to study the evolution of polyploids. A highly contiguous chromosome-level genome sequence of C. neglecta with an N50 size of 29.1 Mb was generated utilizing Pacific Biosciences (PacBio, Menlo Park, CA) long-read sequencing followed by chromosome conformation phasing. Comparison of the genome with that of C. sativa shows remarkable coincidence with subgenome 1 of the hexaploid, with only one major chromosomal rearrangement separating the two. Synonymous substitution rate analysis of the predicted 34 061 genes suggested subgenome 1 of C. sativa directly descended from C. neglecta around 1.2 mya. Higher functional divergence of genes in the hexaploid as evidenced by the greater number of unique orthogroups, and differential composition of resistant gene analogs, might suggest an immediate adaptation strategy after genome merger. The absence of genome bias in gene fractionation among the subgenomes of C. sativa in comparison with C. neglecta, and the complete lack of fractionation of meiosis-specific genes attests to the neopolyploid status of C. sativa. The assembled genome will provide a tool to further study genome evolution processes in the Camelina genus and potentially allow for the identification and exploitation of novel variation for Camelina crop improvement.

Metabolomics and Transcriptomics Provide Insights into Lipid Biosynthesis in the Embryos of Walnut (Juglans regia L.)

Walnut (Juglans regia L.) is an important woody oilseed tree species due to its commercial value. However, the regulation mechanism of walnut oil accumulation is still poorly understood, which restricted the breeding and genetic improvement of high-quality oil-bearing walnuts. In order to explore the metabolic mechanism that regulates the synthesis of walnut oil, we used transcriptome sequencing technology and metabolome technology to comprehensively analyze the key genes and metabolites involved in oil synthesis of the walnut embryo at 60, 90, and 120 days after pollination (DAP). The results showed that the oil and protein contents increased gradually during fruit development, comprising 69.61% and 18.32% of the fruit, respectively, during ripening. Conversely, the contents of soluble sugar and starch decreased gradually during fruit development, comprising 2.14% and 0.84%, respectively, during ripening. Transcriptome sequencing generated 40,631 unigenes across 9 cDNA libraries. We identified 51 and 25 candidate unigenes related to the biosynthesis of fatty acid and the biosynthesis of triacylglycerol (TAG), respectively. The expression levels of the genes encoding Acetyl-CoA carboxylase (ACCase), long-chain acyl-CoA synthetases (LACS), 3-oxoacyl-ACP synthase II (KASII), and glycerol-3-phosphate acyl transfer (GPAT) were upregulated at 60 DAP relative to the levels at 90 and 120 DAP, while the stearoyl-ACP-desaturase (SAD) and fatty acid desaturase 2 (FAD2) genes were highly abundantly expressed during all walnut developmental periods. We found that ABSCISIC ACID INSENSEITIVE3 (ABI3), WRINKLEDl (WRI1), LEAFY COTYLEDON1 (LEC1), and FUSCA3 (FUS3) may be key transcription factors involved in lipid synthesis. Additionally, the metabolomics analysis detected 706 metabolites derived from 18 samples, among which, 4 are implicated in the TAG synthesis, 2 in the glycolysis pathway, and 5 in the tricarboxylic acid cycle (TCA cycle) pathway. The combined analysis of the related genes and metabolites in TAG synthesis showed that phospholipid:diacylglycerol acyltransferase (PDAT) genes were highly abundantly expressed across walnut fruit developmental periods, and their downstream metabolite TAG gradually accumulated with the progression of fruit development. The FAD2 gene showed consistently higher expression during fruit development, and its downstream metabolites 18:2-PC and 18:3-PC gradually accumulated. The ACCase, LACS, SAD, FAD2, and PDAT genes may be crucial genes required for walnut oil synthesis. Our data will enrich public databases and provide new insights into functional genes related to lipid metabolism in walnut.

Do betaine lipids replace phosphatidylcholine as fatty acid editing hubs in microalgae?

Acyl editing refers to a deacylation and reacylation cycle on a lipid, which allows for fatty acid desaturation and modification prior to being removed and incorporated into other pools. Acyl editing is an important determinant of glycerolipid synthesis and has been well-characterized in land plants, thus this review begins with an overview of acyl editing in plants. Much less is known about acyl editing in algae, including the extent to which acyl editing impacts lipid synthesis and on which lipid substrate(s) it occurs. This review compares what is known about acyl editing on its major hub phosphatidylcholine (PC) in land plants with the evidence for acyl editing of betaine lipids such as diacylglyceryltrimethylhomoserine (DGTS), the structural analog that replaces PC in several species of microalgae. In land plants, PC is also known to be a major source of fatty acids and diacylglycerol (DAG) for synthesis of the neutral lipid triacylglycerol (TAG). We review the evidence that DGTS contributes substantially to TAG accumulation in algae as a source of fatty acids, but not as a precursor to DAG. We conclude with evidence of acyl editing on other membrane lipid substrates in plants and algae apart from PC or DGTS, and discuss future analyses to elucidate the role of DGTS and other betaine lipids in acyl editing in microalgae.

Indigenous and commercial isolates of arbuscular mycorrhizal fungi display differential effects in Pyrus betulaefolia roots and elicit divergent transcriptomic and metabolomic responses

BACKGROUND

Arbuscular mycorrhizal fungi (AMF) are beneficial soil fungi which can effectively help plants with acquisition of mineral nutrients and water and promote their growth and development. The effects of indigenous and commercial isolates of arbuscular mycorrhizal fungi on pear (Pyrus betulaefolia) trees, however, remains unclear.

METHODS

Trifolium repens was used to propagate indigenous AMF to simulate spore propagation in natural soils in three ways: 1. the collected soil was mixed with fine roots (R), 2. fine roots were removed from the collected soil (S), and 3. the collected soil was sterilized with 50 kGy ⁶⁰Co γ-radiation (CK). To study the effects of indigenous AMF on root growth and metabolism of pear trees, CK (sterilized soil from CK in T. repens mixed with sterilized standard soil), indigenous AMF (R, soil from R in T. repens mixed with sterilized standard soil; S, soil from S in T. repens mixed with sterilized standard soil), and two commercial AMF isolates (Rhizophagus intraradices(Ri) and Funneliformis mosseae (Fm)) inoculated in the media with pear roots. Effects on plant growth, root morphology, mineral nutrient accumulation, metabolite composition and abundance, and gene expression were analyzed.

RESULTS

AMF treatment significantly increased growth performance, and altered root morphology and mineral nutrient accumulation in this study, with the S treatment displaying overall better performance. In addition, indigenous AMF and commercial AMF isolates displayed common and divergent responses on metabolite and gene expression in pear roots. Compared with CK, most types of flavones, isoflavones, and carbohydrates decreased in the AMF treatment, whereas most types of fatty acids, amino acids, glycerolipids, and glycerophospholipids increased in response to the AMF treatments. Further, the relative abundance of amino acids, flavonoids and carbohydrates displayed different trends between indigenous and commercial AMF isolates. The Fm and S treatments altered gene expression in relation to root metabolism resulting in enriched fructose and mannose metabolism (ko00051), fatty acid biosynthesis (ko00061) and flavonoid biosynthesis (ko00941).

CONCLUSIONS

This study demonstrates that indigenous AMF and commercial AMF isolates elicited different effects in pear plants through divergent responses from gene transcription to metabolite accumulation.

LncRNAs elevate plant adaptation under low temperature by maintaining local chromatin landscape

Epigenetic regulation is one of the most precise and subtle ways of gene regulation, including DNA modification, histone modification, RNA modification, histone variants, chromatin remodeling, and long non-coding RNAs (lncRNAs). Chromatin modification is the most basic type of epigenetic regulation, which plays a key role in a myriad of developmental and physiological processes that have been thoroughly studied. These modifications are usually completed by a series of conserved chromatin modification complexes in eukaryotes. In recent years, a series of lncRNAs in organisms also have been described as having irreplaceable functions in biological environment adaptation, especially in biotic and abiotic stresses. Moreover, these molecules form a sophisticated regulatory network through mutual cross-regulation to achieve quantitative expression of key environmental response genes to external signals. For instance, the function of lncRNAs will directly or indirectly depend on the function of the chromatin modification complex. In this review, we mainly focus on chromatin modification, lncRNA, and their coordination mechanism to achieve the high adaptability of plants in low-temperature environments. We highlight recent findings and insights into lncRNA-mediated local chromatin environment changes during plant growth under low temperature via chromatin modification complexes, including target gene specificity for different lncRNA.

Transcriptome Analysis of Genes Involved in Fatty Acid and Lipid Biosynthesis in Developing Walnut (Juglans regia L.) Seed Kernels from Qinghai Plateau

Walnut (Juglans regia) is an important woody oil-bearing plant with high nutritional value. For better understanding of the underlying molecular mechanisms of its oil accumulation in the Qinghai Plateau, in this study we monitored walnut fruit development, and 15 cDNA libraries were constructed from walnut seed kernels collected at 72, 79, 93, 118 and 135 days after flowering (DAF). The candidate genes were identified using sequencing and expression analysis. The results showed that the oil content in the kernels increased dramatically in late July and reached the maximum value of 69% in mature seed. More than 90% of the oils were unsaturated fatty acids (UFAs) and linoleic acid (18:2) was the predominant UFA accumulated in mature seed. Differentially expressed genes (DEGs) in 15 KEGG pathways of lipid metabolism were detected. We identified 119 DEGs related to FA de novo biosynthesis (38 DEGs), FA elongation and desaturation (39 DEGs), triacylglycerol (TAG) assembly (24 DEGs), oil bodies (12 DEGs), and transcription factors (TFs, 6 DEGs). The abundantly expressed oleosins, caleosins and steroleosins may be important for timely energy reserve in oil bodies. Weighted gene coexpression network analysis (WGCNA) showed that AP2/ERF and bHLH were the key TFs, and were co-expressed with ACC1, α-CT, BCCP, MAT, KASII, LACS, FATA, and PDCT. Our transcriptome data will enrich public databases and provide new insights into functional genes related to the seed kernel lipid metabolism and oil accumulation in J. regia.

Genome-Wide Identification and Expression Analysis of Fatty Acid Desaturase (FAD) Genes in Camelina sativa (L.) Crantz

Camelina sativa (L.) Crantz is an indispensable oilseed crop, and its seeds contain many unsaturated fatty acids. FAD (fatty acid desaturase) regulates the synthesis of unsaturated fatty acids. In this research, we performed CsFAD gene family analysis and identified 24 CsFAD genes in Camelina, which were unevenly distributed on 14 of the 19 total chromosomes. Phylogenetic analysis showed that CsFAD includes four subfamilies, supported by the conserved structures and motifs of CsFAD genes. In addition, we investigated the expression patterns of the FAD family in the different tissues of Camelina. We found that CsFAD family genes were all expressed in the stem, and CsFAD2-2 was highly expressed in the early stage of seed development. Moreover, during low temperature (4 °C) stress, we identified that the expression level of CsFAD2-2 significantly changed. By observing the transient expression of CsFAD2-2 in Arabidopsis protoplasts, we found that CsFAD2-2 was located on the nucleus. Through the detection and analysis of fatty acids, we prove that CsFAD2-2 is involved in the synthesis of linolenic acid (C18:3). In conclusion, we identified CsFAD2-2 through the phylogenetic analysis of the CsFAD gene family and further determined the fatty acid content to find that CsFAD2-2 is involved in fatty acid synthesis in Camelina.

CRISPR/Cas9-mediated editing of double loci of BnFAD2 increased the seed oleic acid content of rapeseed (Brassica napus L.)

Seed oleic acid is an important quality trait sought in rapeseed breeding programs. Many methods exist to increase seed oleic acid content, such as the CRISPR/Cas9-mediated genome editing system, yet there is no report on seed oleic acid content improvement via this system's precise editing of the double loci of BnFAD2. Here, a precise CRISPR/Cas9-mediated genome editing of the encoded double loci (A5 and C5) of BnFAD2 was established. The results demonstrated high efficiency of regeneration and transformation, with the rapeseed genotype screened in ratios of 20.18% and 85.46%, respectively. The total editing efficiency was 64.35%, whereas the single locus- and double locus-edited ratios were 21.58% and 78.42%, respectively. The relative proportion of oleic acid with other fatty acids in seed oil of mutants was significantly higher for those that underwent the editing on A5 copy than that on C5 copy, but it was still less than 80%. For double locus-edited mutants, their relative proportion of oleic acid was more than 85% in the T₁ and T₄ generations. A comparison of the sequences between the double locus-edited mutants and reference showed that no transgenic border sequences were detected from the transformed vector. Analysis of the BnFAD2 sequence on A5 and C5 at the mutated locus of double loci mutants uncovered evidence for base deletion and insertion, and combination. Further, no editing issue of FAD2 on the copy of A1 was detected on the three targeted editing regions. Seed yield, yield component, oil content, and relative proportion of oleic acid between one selected double loci-edited mutant and wild type were also compared. These results showed that although the number of siliques per plant of the wild type was significantly higher than those of the mutant, the differences in seed yield and oil content were not significant between them, albeit with the mutant having a markedly higher relative proportion of oleic acid. Altogether, our results confirmed that the established CRISPR/Cas9-mediated genome editing of double loci (A5 and C5) of the BnFAD2 can precisely edit the targeted genes, thereby enhancing the seed oleic acid content to a far greater extent than can a single locus-editing system.

Insights into the Novel FAD2 Gene Regulating Oleic Acid Accumulation in Peanut Seeds with Different Maturity

AhFAD2 is a key enzyme catalyzing the conversion of oleic acid into linoleic acid. The high oleic acid characteristic of peanut mainly comes from the homozygous recessive mutation of AhFAD2A and AhFAD2B genes (aabb). However, even in high-oleic-acid varieties with the aabb genotype, the oleic acid content of seeds with different maturity varies significantly. Therefore, in addition to AhFAD2A and AhFAD2B, other FAD2 members or regulators may be involved in this process. Which FAD2 genes are involved in the regulatory processes associated with seed maturity is still unclear. In this study, four stable lines with different genotypes (AABB, aaBB, AAbb, and aabb) were used to analyze the contents of oleic acid and linoleic acid at different stages of seed development in peanut. Three new AhFAD2 genes (AhFAD2-7, AhFAD2-8, and AhFAD2-9) were cloned based on the whole-genome sequencing results of cultivated peanuts. All peanut FAD2 genes showed tissue preference in expression; however, only the expression level of AhFAD2-7 was positively correlated with the linoleic acid concentration in peanut seeds. These findings provide new insights into the regulation of oleic acid accumulation by maturity, and AhFAD2-7 plays an important role in the maturity dependent accumulation of oleic acid and linoleic acid in peanut.

Applications and prospects of genome editing in plant fatty acid and triacylglycerol biosynthesis

Triacylglycerol (TAG), which is a neutral lipid, has a structure in which three molecules of fatty acid (FA) are ester-bonded to one molecule of glycerol. TAG is important energy source for seed germination and seedling development in plants. Depending on the FA composition of the TAG, it is used as an edible oil or industrial material for cosmetics, soap, and lubricant. As the demand for plant oil is rising worldwide, either the type of FA must be changed or the total oil content of various plants must be increased. In this review, we discuss the regulation of FA metabolism by Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, a recent genome-editing technology applicable to various plants. The development of plants with higher levels of oleic acid or lower levels of very long-chain fatty acids (VLCFAs) in seeds are discussed. In addition, the current status of research on acyltransferases, phospholipases, TAG lipases, and TAG synthesis in vegetative tissues is described. Finally, strategies for the application of CRISPR/Cas9 in lipid metabolism studies are mentioned.

Molecular Approaches Reduce Saturates and Eliminate trans Fats in Food Oils

Vegetable oils composed of triacylglycerols (TAG) are a major source of calories in human diets. However, the fatty acid compositions of these oils are not ideal for human nutrition and the needs of the food industry. Saturated fatty acids contribute to health problems, while polyunsaturated fatty acids (PUFA) can become rancid upon storage or processing. In this review, we first summarize the pathways of fatty acid metabolism and TAG synthesis and detail the problems with the oil compositions of major crops. Then we describe how transgenic expression of desaturases and downregulation of the plastid FatB thioesterase have provided the means to lower oil saturates. The traditional solution to PUFA rancidity uses industrial chemistry to reduce PUFA content by partial hydrogenation, but this results in the production of trans fats that are even more unhealthy than saturated fats. We detail the discoveries in the biochemistry and molecular genetics of oil synthesis that provided the knowledge and tools to lower oil PUFA content by blocking their synthesis during seed development. Finally, we describe the successes in breeding and biotechnology that are giving us new, high-oleic, low PUFA varieties of soybean, canola and other oilseed crops.

Overexpression mutants reveal a role for a chloroplast MPD protein in regulation of reactive oxygen species during chilling in Arabidopsis

Reactive oxygen species (ROS) contribute to cellular damage in several different contexts, but their role during chilling damage is poorly defined. Chilling sensitivity both limits the distribution of plant species and causes devastating crop losses worldwide. Our screen of chilling-tolerant Arabidopsis (Arabidopsis thaliana) for mutants that suffer chilling damage identified a gene (At4g03410) encoding a chloroplast Mpv17_PMP22 protein, MPD1, with no previous connection to chilling. The chilling-sensitive mpd1-1 mutant is an overexpression allele that we successfully phenocopied by creating transgenic lines with a similar level of MPD1 overexpression. In mammals and yeast, MPD1 homologs are associated with ROS management. In chilling conditions, Arabidopsis overexpressing MPD1 accumulated H2O2 to higher levels than wild-type controls and exhibited stronger induction of ROS response genes. Paraquat application exacerbated chilling damage, confirming that the phenotype occurs due to ROS dysregulation. We conclude that at low temperature increased MPD1 expression results in increased ROS production, causing chilling damage. Our discovery of the effect of MPD1 overexpression on ROS production under chilling stress implies that investigation of the nine other members of the Mpv17_PMP22 family in Arabidopsis may lead to new discoveries regarding ROS signaling and management in plants.

Comparative Analysis and Structural Modeling of Elaeis oleifera FAD2, a Fatty Acid Desaturase Involved in Unsaturated Fatty Acid Composition of American Oil Palm

Simple Summary

Palm oil has become the world’s most important vegetable oil in terms of production quantity, and its overall demand is exponentially growing with the global population. The fatty acid composition and particularly the oleic/linoleic acid ratio are major factors influencing palm oil quality. In this study, we focused on FAD2, a fatty acid desaturase enzyme involved in the desaturation and conversion of oleic acid to linoleic acid in Elaeis oleifera, identified through in silico annotation analysis. Our phylogenetic and comparative studies revealed two SNP markers, SNP278 and SNP851, significantly correlated with the oleic/linoleic acid contents. Our study provides fundamental insights into the mechanism of fatty acids synthesis in oil palm and could support the application of molecular biology techniques to enhance the enzymatic activity and substrate affinity of EoFAD2.

Abstract

American oil palm (Elaeis oleifera) is an important source of dietary oil that could fulfill the increasing worldwide demand for cooking oil. Therefore, improving its production is crucial and could be realized through breeding and genetic engineering approaches aiming to obtain high-yielding varieties with improved oil content and quality. The fatty acid composition and particularly the oleic/linoleic acid ratio are major factors influencing oil quality. Our work focused on a fatty acid desaturase (FAD) enzyme involved in the desaturation and conversion of oleic acid to linoleic acid. Following the in silico identification and annotation of Elaeis oleifera FAD2, its molecular and structural features characterization was performed to better understand the mechanistic bases of its enzymatic activity. EoFAD2 is 1173 nucleotides long and encodes a protein of 390 amino acids that shares similarities with other FADs. Interestingly, the phylogenetic study showed three distinguished groups where EoFAD2 clustered among monocotyledonous taxa. EoFAD2 is a membrane-bound protein with five transmembrane domains presumably located in the endoplasmic reticulum. The homodimer organization model of EoFAD2 enzyme and substrates and respective substrate-binding residues were predicted and described. Moreover, the comparison between 24 FAD2 sequences from different species generated two interesting single-nucleotide polymorphisms (SNPs) associated with the oleic/linoleic acid contents.

Flavonoid deficiency disrupts redox homeostasis and terpenoid biosynthesis in glandular trichomes of tomato

Glandular trichomes (GTs) are epidermal structures that provide the first line of chemical defense against arthropod herbivores and other biotic threats. The most conspicuous structure on leaves of cultivated tomato (Solanum lycopersicum) is the type-VI GT (tVI-GT), which accumulates both flavonoids and volatile terpenoids. Although these classes of specialized metabolites are derived from distinct metabolic pathways, previous studies with a chalcone isomerase 1 (CHI1)-deficient mutant called anthocyanin free (af) showed that flavonoids are required for terpenoid accumulation in tVI-GTs. Here, we combined global transcriptomic and proteomic analyses of isolated trichomes as a starting point to show that the lack of CHI1 is associated with reduced levels of terpenoid biosynthetic transcripts and enzymes. The flavonoid deficiency in af trichomes also resulted in the upregulation of abiotic stress-responsive genes associated with DNA damage and repair. Several lines of biochemical and genetic evidence indicate that the terpenoid defect in af mutants is specific for the tVI-GT and is associated with the absence of bulk flavonoids rather than loss of CHI1 per se. A newly developed genome-scale model of metabolism in tomato tVI-GTs helped identify metabolic imbalances caused by the loss of flavonoid production. We provide evidence that flavonoid deficiency in this cell type leads to increased production of reactive oxygen species (ROS), which may impair terpenoid biosynthesis. Collectively, our findings support a role for flavonoids as ROS-scavenging antioxidants in GTs.

The genetic control of polyacetylenes involved in bitterness of carrots (Daucus carota L.): Identification of QTLs and candidate genes from the plant fatty acid metabolism

BACKGROUND

Falcarinol-type polyacetylenes (PAs) such as falcarinol (FaOH) and falcarindiol (FaDOH) are produced by several Apiaceae vegetables such as carrot, parsnip, celeriac and parsley. They are known for numerous biological functions and contribute to the undesirable bitter off-taste of carrots and their products. Despite their interesting biological functions, the genetic basis of their structural diversity and function is widely unknown. A better understanding of the genetics of the PA levels present in carrot roots might support breeding of carrot cultivars with tailored PA levels for food production or nutraceuticals.

RESULTS

A large carrot F₂ progeny derived from a cross of a cultivated inbred line with an inbred line derived from a Daucus carota ssp. commutatus accession rich in PAs was used for linkage mapping and quantitative trait locus (QTL) analysis. Ten QTLs for FaOH and FaDOH levels in roots were identified in the carrot genome. Major QTLs for FaOH and FaDOH with high LOD values of up to 40 were identified on chromosomes 4 and 9. To discover putative candidate genes from the plant fatty acid metabolism, we examined an extended version of the inventory of the carrot FATTY ACID DESATURASE2 (FAD2) gene family. Additionally, we used the carrot genome sequence for a first inventory of ECERIFERUM1 (CER1) genes possibly involved in PA biosynthesis. We identified genomic regions on different carrot chromosomes around the found QTLs that contain several FAD2 and CER1 genes within their 2-LOD confidence intervals. With regard to the major QTLs on chromosome 9 three putative CER1 decarbonylase gene models are proposed as candidate genes.

CONCLUSION

The present study increases the current knowledge on the genetics of PA accumulation in carrot roots. Our finding that carrot candidate genes from the fatty acid metabolism are significantly associated with major QTLs for both major PAs, will facilitate future functional gene studies and a further dissection of the genetic factors controlling PA accumulation. Characterization of such candidate genes will have a positive impact on carrot breeding programs aimed at both lowering or increasing PA concentrations in carrot roots.

Crop Quality Improvement Through Genome Editing Strategy

Good quality of crops has always been the most concerning aspect for breeders and consumers. However, crop quality is a complex trait affected by both the genetic systems and environmental factors, thus, it is difficult to improve through traditional breeding strategies. Recently, the CRISPR/Cas9 genome editing system, enabling efficiently targeted modification, has revolutionized the field of quality improvement in most crops. In this review, we briefly review the various genome editing ability of the CRISPR/Cas9 system, such as gene knockout, knock-in or replacement, base editing, prime editing, and gene expression regulation. In addition, we highlight the advances in crop quality improvement applying the CRISPR/Cas9 system in four main aspects: macronutrients, micronutrients, anti-nutritional factors and others. Finally, the potential challenges and future perspectives of genome editing in crop quality improvement is also discussed.

Plastidial wax ester biosynthesis as a tool to synthesize shorter and more saturated wax esters

BACKGROUND

Wax esters (WE) are neutral lipids that consist of a fatty alcohol esterified to a fatty acid. WE are valuable feedstocks in industry for producing lubricants, coatings, and cosmetics. They can be produced chemically from fossil fuel or plant-derived triacylglycerol. As fossil fuel resources are finite, the synthesis of WE in transgenic plants may serve as an alternative source. As chain length and desaturation of the alcohol and acyl moieties determine the physicochemical properties of WE and their field of application, tightly controlled and tailor-made WE synthesis in plants would be a sustainable, beneficial, and valuable commodity. Here, we report the expression of ten combinations of WE producing transgenes in Arabidopsis thaliana. In order to study their suitability for WE production in planta, we analyzed WE amount and composition in the transgenic plants.

RESULTS

The transgenes consisted of different combinations of a FATTY ACYL-COA/ACP REDUCTASE (FAR) and two WAX SYNTHASES/ACYL-COA:DIACYLGLYCEROL O-ACYLTRANSFERASES (WSD), namely WSD2 and WSD5 from the bacterium Marinobacter aquaeoleoi. We generated constructs with and without plastidial transit peptides to access distinct alcohol and acyl substrate pools within A. thaliana cells. We observed WE formation with plastid and cytosol-localized FAR and WSD in seeds. A comparative WE analysis revealed the production of shorter and more saturated WE by plastid-localized WE biosynthesis compared to cytosolic WE synthesis.

CONCLUSIONS

A shift of WE formation into seed plastids is a suitable approach for tailor-made WE production and can be used to synthesize WE that are mainly derived from mid- and long-chain saturated and monounsaturated substrates.

Genome-Wide Analysis of the Fatty Acid Desaturase Gene Family Reveals the Key Role of PfFAD3 in α-Linolenic Acid Biosynthesis in Perilla Seeds

Perilla (Perilla frutescens), a traditional medicinal and oilseed crop in Asia, contains extremely high levels of polyunsaturated α-linolenic acid (ALA) (up to 60.9%) in its seeds. ALA biosynthesis is a multistep process catalyzed by fatty acid desaturases (FADs), but the FAD gene family in perilla has not been systematically characterized. Here, we identified 42 PfFADs in the perilla genome and classified them into five subfamilies. Subfamily members of PfFADs had similar exon/intron structures, conserved domain sequences, subcellular localizations, and cis-regulatory elements in their promoter regions. PfFADs also possessed various expression patterns. PfFAD3.1 was highly expressed in the middle stage of seed development, whereas PfFAD7/8.3 and PfFAD7/8.5 were highly expressed in leaf and later stages of seed development, respectively. Phylogenetic analysis revealed that the evolutionary features coincided with the functionalization of different subfamilies of PUFA desaturase. Heterologous overexpression of PfFAD3.1 in Arabidopsis thaliana seeds increased ALA content by 17.68%–37.03%. These findings provided insights into the characteristics and functions of PfFAD genes in perilla.

A Tree Peony Trihelix Transcription Factor PrASIL1 Represses Seed Oil Accumulation

In many higher plants, seed oil accumulation is governed by complex multilevel regulatory networks including transcriptional regulation, which primarily affects fatty acid biosynthesis. Tree peony (Paeonia rockii), a perennial deciduous shrub endemic to China is notable for its seed oil that is abundant in unsaturated fatty acids. We discovered that a tree peony trihelix transcription factor, PrASIL1, localized in the nucleus, is expressed predominantly in developing seeds during maturation. Ectopic overexpression of PrASIL1 in Nicotiana benthamiana leaf tissue and Arabidopsis thaliana seeds significantly reduced total fatty acids and altered the fatty acid composition. These changes were in turn associated with the decreased expression of multitudinous genes involved in plastidial fatty acid synthesis and oil accumulation. Thus, we inferred that PrASIL1 is a critical transcription factor that represses oil accumulation by down-regulating numerous key genes during seed oil biosynthesis. In contrary, up-regulation of oil biosynthesis genes and a significant increase in total lipids and several major fatty acids were observed in PrASIL1-silenced tree peony leaves. Together, these results provide insights into the role of trihelix transcription factor PrASIL1 in controlling seed oil accumulation. PrASIL1 can be targeted potentially for oil enhancement in tree peony and other crops through gene manipulation.

Abscisic Acid Improves Linoleic Acid Accumulation Possibly by Promoting Expression of EgFAD2 and Other Fatty Acid Biosynthesis Genes in Oil Palm Mesocarp

Abscisic acid plays an important role in fruit development. However, the effect of ABA on fatty acid biosynthesis in oil palm is still unknown. In this study, ABA treatments (CK, A1-A4) were applied to oil palm fruit at 16 WAP (weeks after pollination), and fatty acids in the mesocarp at 24 WAP were analyzed by GC-MS. Results showed that linoleic acid content under treatment A2 (20 μM ABA) was significantly higher (slightly increased by 8.33%) than the control. Therefore, mesocarp samples of A2, and the control at 16, 20, and 24 WAP was sampled for RNA-Seq. KEGG pathway enrichment analysis showed that 43 genes were differentially expressed in the fatty acid biosynthesis pathway, of which expression of EgFAD2 (unigene 105050201) under 20 μM ABA treatment was 1.84-fold higher than in the control at 20 WAP. Further sequence analysis found that unigene 105050201 had more ABA-responsive elements (ABRE), complete conserved domains, and a C-terminal signaling motif among two FAD2 copies. Furthermore, WGCNA and correlation analysis showed co-expression of EgFAD2 (unigene 105050201) with transcription factors (TFs) (WRI1, AP2-EREBP, bZIP, bHLH, C2C2-Dof, MYB, NAC, and WRKY), ABA signaling genes (PYR, PP2C, SnRK, and ABI5), and other genes involved in fatty acid biosynthesis (FATA, FATB, LACS, SAD, Oleosins, and so on). These results indicated that ABA treatment promoted the expression of FAD2 and other genes involved in fatty acid biosynthesis, which possibly resulted in the accumulation of linoleic acid. This study will be helpful for understanding the possible mechanisms through which ABA affects fatty acid biosynthesis and their accumulation in the mesocarp of oil palm.

Quantitative proteomic and lipidomics analyses of high oil content GmDGAT1-2 transgenic soybean illustrate the regulatory mechanism of lipoxygenase and oleosin

KEY MESSAGE

Proteomic and lipidomics analyses of WT and GmDGAT1-2 transgenic soybeans showed that GmDGAT1-2 over-expression induced lipoxygenase down-regulatation and oleoin up-regulatation, which significantly changed the compositions and total fatty acid. The main goal of soybean breeding is to increase the oil content. Diacylglycerol acyltransferase (DGAT) is a key rate-limiting enzyme in fatty acid metabolism and may regulate oil content. Herein, 10 GmDGAT genes were isolated from soybean and transferred into wild-type (WT) Arabidopsis. The total fatty acid was 1.2 times higher in T3 GmDGAT1-2 transgenic Arabidopsis seeds than in WT. Therefore, GmDGAT1-2 was transferred into WT soybean (JACK), and four T3 transgenic soybean lines were obtained. The results of high-performance gas chromatography and Soxhlet extractor showed that, compared with those of JACK, oleic acid (18:1), and total fatty acid levels in transgenic soybean plants were much higher, but linoleic acid (18:2) was lower than WT. Palmitic acid (16:0), stearic acid (18:0), and linolenic acid (18:3) were not significantly different. For mechanistic studies, 436 differentially expressed proteins (DEPs) and 180 differentially expressed metabolites (DEMs) were identified between WT (JACK) and transgenic soybean pods using proteomic and lipidomics analyses. Four lipoxygenase proteins were down-regulated in linoleic acid metabolism while four oleosin proteins were up-regulated in the final oil formation. The results showed an increase in the total fatty acid and 18:1 composition, and a decrease in the 18:2 composition of fatty acid. Our study brings new insights into soybean genetic transformation and the deep study of molecular mechanism that changes the total fatty acid, 18:1, and 18:2 compositions in GmDGAT1-2 transgenic soybean.

Construction of Multiple Guide RNAs in CRISPR/Cas9 Vector Using Stepwise or Simultaneous Golden Gate Cloning: Case Study for Targeting the FAD2 and FATB Multigene in Soybean

CRISPR/Cas9 is a commonly used technique in reverse-genetics research to knock out a gene of interest. However, when targeting a multigene family or multiple genes, it is necessary to construct a vector with multiple single guide RNAs (sgRNAs) that can navigate the Cas9 protein to the target site. In this protocol, the Golden Gate cloning method was used to generate multiple sgRNAs in the Cas9 vector. The vectors used were pHEE401E_UBQ_Bar and pBAtC_tRNA, which employ a one-promoter/one-sgRNA and a polycistronic-tRNA-gRNA strategy, respectively. Golden Gate cloning was performed with type IIS restriction enzymes to generate gRNA polymers for vector inserts. Four sgRNAs containing the pHEE401E_UBQ_Bar vector and four to six sgRNAs containing the pBAtC_tRNA vector were constructed. In practice, we constructed multiple sgRNAs targeting multiple genes of FAD2 and FATB in soybean using this protocol. These three vectors were transformed into soybeans using the Agrobacterium-mediated method. Using deep sequencing, we confirmed that the T0 generation transgenic soybean was edited at various indel ratios in the predicted target regions of the FAD2 and FATB multigenes. This protocol is a specific guide that allows researchers to easily follow the cloning of multiple sgRNAs into commonly used CRISPR/Cas9 vectors for plants.

A genome variation map provides insights into the genetics of walnut adaptation and agronomic traits

BACKGROUND

Common walnut (Juglans regia L.) is one of the top four most consumed nuts in the world due to its health benefits and pleasant taste. Despite its economic importance, the evolutionary history and genetic control of its adaptation and agronomic traits remain largely unexplored.

RESULTS

We report a comprehensive walnut genomic variation map based on whole-genome resequencing of 815 walnut accessions. Evolutionary analyses suggest that Chinese J. regia diverged from J. sigillata with extensive hybridizations after the split of the two species. In contrast to annual crops, the genetic diversity and heterozygous deleterious mutations of Chinese common walnut trees have continued to increase during the improvement process. Selective sweep analyses identify 902 genes uniquely selected in the improved common walnut compared to its progenitor population. Five major-effect loci are identified to be involved in walnut adaptations to temperature, precipitation, and altitude. Genome-wide association studies reveal 27 genomic loci responsible for 18 important agronomic traits, among which JrFAD2 and JrANR are the potentially major-effect causative genes controlling linoleic acid content and color of the endopleura of the nut, respectively.

CONCLUSIONS

The largest genomic resource for walnuts to date has been generated and explored in this study, unveiling their evolutionary history and cracking the genetic code for agronomic traits and environmental adaptation of this economically crucial crop tree.

C-to-G Base Editing Enhances Oleic Acid Production by Generating Novel Alleles of FATTY ACID DESATURASE 2 in Plants

The demand for vegetable oil, which is mainly used for dietary purposes and cooking, is steadily increasing worldwide. It is often desirable to reduce unsaturation levels of fatty acids in order to increase storage stability and reduce trans-fat generation during cooking. Functional disruption of FATTY ACID DESATURASE 2 (FAD2) prevents the conversion of monounsaturated oleic acid to polyunsaturated linoleic acid, thereby enhancing the production of the desirable oleic acid. However, FAD2 null alleles, due to growth defects under stress conditions, are impractical for agronomical purposes. Here, we aimed to attenuate FAD2 activity in planta while avoiding adverse growth effects by introducing amino-acid substitutions using CRISPR base editors. In Arabidopsis, we applied the adenine base editor (ABE) and cytosine base editor (CBE) to induce semi-random base substitutions within several selected FAD2 coding regions. Isolation of base-edited fad2 alleles with higher oleic acid revealed that the CBE application induced C-to-T and/or C-to-G base substitutions within the targeted sequences, resulting in an alteration of the FAD2 enzyme activities; for example, fad2-144 with multiple C-to-G base substitutions showed less growth defects but with a significant increase in oleic acids by 3-fold higher than wild type. Our "proof-of-concept" approach suggests that equivalent alleles may be generated in vegetable oil crops via precision genome editing for practical cultivation. Our targeted semi-random strategy may serve as a new complementary platform for planta engineering of useful agronomic traits.

In Silico Analysis of Fatty Acid Desaturases Structures in Camelina sativa, and Functional Evaluation of Csafad7 and Csafad8 on Seed Oil Formation and Seed Morphology

Fatty acid desaturases add a second bond into a single bond of carbon atoms in fatty acid chains, resulting in an unsaturated bond between the two carbons. They are classified into soluble and membrane-bound desaturases, according to their structure, subcellular location, and function. The orthologous genes in Camelina sativa were identified and analyzed, and a total of 62 desaturase genes were identified. It was revealed that they had the common fatty acid desaturase domain, which has evolved separately, and the proteins of the same family also originated from the same ancestry. A mix of conserved, gained, or lost intron structure was obvious. Besides, conserved histidine motifs were found in each family, and transmembrane domains were exclusively revealed in the membrane-bound desaturases. The expression profile analysis of C. sativa desaturases revealed an increase in young leaves, seeds, and flowers. C. sativa ω3-fatty acid desaturases CsaFAD7 and CsaDAF8 were cloned and the subcellular localization analysis showed their location in the chloroplast. They were transferred into Arabidopsis thaliana to obtain transgenic lines. It was revealed that the ω3-fatty acid desaturase could increase the C18:3 level at the expense of C18:2, but decreases in oil content and seed weight, and wrinkled phenotypes were observed in transgenic CsaFAD7 lines, while no significant change was observed in transgenic CsaFAD8 lines in comparison to the wild-type. These findings gave insights into the characteristics of desaturase genes, which could provide an excellent basis for further investigation for C. sativa improvement, and overexpression of ω3-fatty acid desaturases in seeds could be useful in genetic engineering strategies, which are aimed at modifying the fatty acid composition of seed oil.

Development of B. carinata with super-high erucic acid content through interspecific hybridization

KEY MESSAGE

Disomic alien chromosome addition Brassica carinata lines with super-high erucic acid content were developed through interspecific hybridization with B. juncea and characterized using molecular, cytological and biochemical techniques. Brassica carinata [A.] Braun (BBCC, 2n = 34) is a climate-resilient oilseed. Its seed oil is high in erucic acid (> 40%), rendering it well suited for the production of biofuel and other bio-based applications. To enhance the competitiveness of B. carinata with high erucic B. napus (HEAR), lines with super-high erucic acid content were developed through interspecific hybridization. To this end, a fad2B null allele from Brassica juncea (AABB, 2n = 36) was introgressed into B. carinata, resulting in a B. carinata fad2B mutant with erucic acid levels of over 50%. Subsequently, the FAE allele from B. rapa spp. yellow sarson (AA, 2n = 20) was transferred to the fad2B B. carinata line, yielding lines with erucic acid contents of up to 57.9%. Molecular analysis using the Brassica 90 K Illumina Infinium™ SNP genotyping array identified these lines as disomic alien chromosome addition lines, with two extra A08 chromosomes containing the BrFAE gene. The alien chromosomes from B. rapa were clearly distinguished by molecular cytogenetics in one of the addition lines. Analysis of microspore-derived offspring and hybrids from crosses with a CMS B. carinata line showed that the transfer rate of the A08 chromosome into male gametes was over 98%, resulting in almost completely stable transmission of an A08 chromosome copy into the progeny. The increase in erucic acid levels was accompanied by changes in the proportions of other fatty acids depending on the genetic changes that were introduced in the interspecific hybrids, providing valuable insights into erucic acid metabolism in Brassica.

Specific Incorporation of Polyunsaturated Fatty Acids into the sn-2 Position of Phosphatidylglycerol Accelerates Photodamage to Photosystem II under Strong Light

Free fatty acids (FFAs) are generated by the reaction of lipases with membrane lipids. Generated polyunsaturated fatty acids (PUFAs) containing more than two double bonds have toxic effects in photosynthetic organisms. In the present study, we examined the effect of exogenous FFAs in the growth medium on the activity of photosystem II (PSII) under strong light in the cyanobacterium Synechocystis sp. PCC 6803 (Synechocystis). PUFAs but not monounsaturated fatty acids accelerated the rate of photodamage to PSII by inactivating electron transfer at the oxygen-evolving complex. Moreover, supplemented PUFAs were specifically incorporated into the sn-2 position of phosphatidylglycerol (PG), which usually contains C16 fatty acids at the sn-2 position in Synechocystis cells. The disruption of the gene for an acyl-ACP synthetase reduced the effect of PUFAs on the photoinhibition of PSII. Thus, the specific incorporation of PUFAs into PG molecules requires acyl-ACP synthetase and leads to an unstable PSII, thereby accelerating photodamage to PSII. Our results are a breakthrough into elucidating the molecular mechanism of the toxicity of PUFAs to photosynthetic organisms.

Characterization of fatty acid desaturases reveals stress-induced synthesis of C18 unsaturated fatty acids enriched in triacylglycerol in the oleaginous alga Chromochloris zofingiensis

BACKGROUND

The green microalga Chromochloris zofingiensis is capable of producing high levels of triacylglycerol rich in C18 unsaturated fatty acids (UFAs). FA desaturation degree is regulated by FA desaturases (FADs). Nevertheless, it remains largely unknown regarding what FADs are involved in FA desaturations and how these FADs collaborate to contribute to the high abundance of C18 UFAs in triacylglycerol in C. zofingiensis.

RESULTS

To address these issues, we firstly determined the transcription start sites of 11 putative membrane-bound FAD-coding genes (CzFADs) and updated their gene models. Functional validation of these CzFADs in yeast and cyanobacterial cells revealed that seven are bona fide FAD enzymes with distinct substrates. Combining the validated functions and predicted subcellular compartments of CzFADs and the FA profiles of C. zofingiensis, the FA desaturation pathways in this alga were reconstructed. Furthermore, a multifaceted lipidomic analysis by systematically integrating thin-layer chromatography, gas chromatography-mass spectrometry and liquid chromatography-mass spectrometry techniques was conducted, unraveling profiles of polar membrane lipids in C. zofingiensis and major desaturation steps occurring in these lipids. By correlating transcriptional patterns of CzFAD genes and changes of lipids upon abiotic stress conditions, our results highlighted collaboration of CzFADs for C18 UFA synthesis and supported that both de novo FA synthesis and membrane lipid remodeling contributed C18 UFAs to triacylglycerol for storage.

CONCLUSIONS

Taken together, our study for the first time elucidated the pathways of C18 FA desaturations and comprehensive profiles of polar membrane lipids in C. zofingiensis and shed light on collaboration of CzFADs for the synthesis and enrichment of C18 UFAs in triacylglycerol.

The chromosome-scale reference genome of safflower (Carthamus tinctorius) provides insights into linoleic acid and flavonoid biosynthesis

Safflower (Carthamus tinctorius L.), a member of the Asteraceae, is a popular crop due to its high linoleic acid (LA) and flavonoid (such as hydroxysafflor yellow A) contents. Here, we report the first high-quality genome assembly (contig N50 of 21.23 Mb) for the 12 pseudochromosomes of safflower using single-molecule real-time sequencing, Hi-C mapping technologies and a genetic linkage map. Phyloge nomic analysis showed that safflower diverged from artichoke (Cynara cardunculus) and sunflower (Helianthus annuus) approximately 30.7 and 60.5 million years ago, respectively. Comparative genomic analyses revealed that uniquely expanded gene families in safflower were enriched for those predicted to be involved in lipid metabolism and transport and abscisic acid signalling. Notably, the fatty acid desaturase 2 (FAD2) and chalcone synthase (CHS) families, which function in the LA and flavonoid biosynthesis pathways, respectively, were expanded via tandem duplications in safflower. CarFAD2-12 was specifically expressed in seeds and was vital for high-LA content in seeds, while tandemly duplicated CarFAD2 genes were up-regulated in ovaries compared to CarFAD2-12, which indicates regulatory divergence of FAD2 in seeds and ovaries. CarCHS1, CarCHS4 and tandem-duplicated CarCHS5˜CarCHS6, which were up-regulated compared to other CarCHS members at early stages, contribute to the accumulation of major flavonoids in flowers. In addition, our data reveal multiple alternative splicing events in gene families related to fatty acid and flavonoid biosynthesis. Together, these results provide a high-quality reference genome and evolutionary insights into the molecular basis of fatty acid and flavonoid biosynthesis in safflower.

Transcriptomic Analysis Reveals Key Genes Involved in Oil and Linoleic Acid Biosynthesis during Artemisia sphaerocephala Seed Development

Artemisia sphaerocephala seeds are rich in polysaccharides and linoleic acid (C18:2), which have been widely used as traditional medicine and to improve food quality. The accumulation patterns and molecular regulatory mechanisms of polysaccharides during A. sphaerocephala seed development have been studied. However, the related research on seed oil and C18:2 remain unclear. For this study, A. sphaerocephala seeds at seven different development stages at 10, 20, 30, 40, 50, 60, and 70 days after flowering (designated as S1~S7), respectively, were employed as experimental samples, the accumulation patterns of oil and fatty acids (FA) and the underlying molecular regulatory mechanisms were analyzed. The results revealed that oil content increased from 10.1% to 20.0% in the early stages of seed development (S1~S2), and up to 32.0% in mature seeds, of which C18:2 accounted for 80.6% of the total FA. FA and triacylglycerol biosynthesis-related genes jointly involved in the rapid accumulation of oil in S1~S2. Weighted gene co-expression network analysis showed that transcription factors FUS3 and bHLH played a critical role in the seed oil biosynthesis. The perfect harmonization of the high expression of FAD2 with the extremely low expression of FAD3 regulated the accumulation of C18:2. This study uncovered the gene involved in oil biosynthesis and molecular regulatory mechanisms of high C18:2 accumulation in A. sphaerocephala seeds; thus, advancing research into unsaturated fatty acid metabolism in plants while generating valuable genetic resources for optimal C18:2 breeding.

Increasing Monounsaturated Fatty Acid Contents in Hexaploid Camelina sativa Seed Oil by FAD2 Gene Knockout Using CRISPR-Cas9

Seed oils are used as edible oils and increasingly also for industrial applications. Although high-oleic seed oil is preferred for industrial use, most seed oil is high in polyunsaturated fatty acids (PUFAs) and low in monounsaturated fatty acids (MUFAs) such as oleic acid. Oil from Camelina, an emerging oilseed crop with a high seed oil content and resistance to environmental stress, contains 60% PUFAs and 30% MUFAs. Hexaploid Camelina carries three homoeologs of FAD2, encoding fatty acid desaturase 2 (FAD2), which is responsible for the synthesis of linoleic acid from oleic acid. In this study, to increase the MUFA contents of Camelina seed oil, we generated CsFAD2 knockout plants via CRISPR-Cas9-mediated gene editing using the pRedU6fad2EcCas9 vector containing DsRed as a selection marker, the U6 promoter to drive a single guide RNA (sgRNA) covering the common region of the three CsFAD2 homoeologs, and an egg-cell-specific promoter to drive Cas9 expression. We analyzed CsFAD2 homoeolog-specific sequences by PCR using genomic DNA from transformed Camelina leaves. Knockout of all three pairs of FAD2 homoeologs led to a stunted bushy phenotype, but greatly enhanced MUFA levels (by 80%) in seeds. However, transformants with two pairs of CsFAD2 homoeologs knocked out but the other pair wild-type heterozygous showed normal growth and a seed MUFAs production increased up to 60%. These results provide a basis for the metabolic engineering of genes that affect growth in polyploid crops through genome editing.

Functional characterization and expression profile of microsomal FAD2 and FAD3 genes involved in linoleic and α-linolenic acid production in Leucas cephalotes

The genus Leucas belongs to Lamiaceae, and has attained more attention due to the presence of unusual allenic fatty acids called laballenic and phlomic acid in majority of its species. This genus has been known since traditional medicinal times and has numerous economical, nutritional, and industrial properties. So far genetic, molecular and biochemical analyses of lipid metabolism and fatty acid biosynthetic pathway in Leucas has not been reported. The objective of this study is to identify, isolate, analyze expression profiles, and functionally characterize the membrane-associated desaturases responsible for unsaturated fatty acid accumulation in Leucas cephalotes. Full-length LcFAD2 and LcFAD3 cDNAs were isolated and expressed in Saccharomyces cerevisiae BY4741 for functional characterization. Substrate feeding assay using S. cerevisiae confirmed that the LcFAD2 enzyme catalyzes desaturation of both palmitoleic (16:1∆⁹) and oleic (18:1∆⁹) acids to form palmitolinoleic (16:2∆^9,12) and linoleic (18:2∆^9,12) acids respectively. As a contrast, the heterologous activity of LcFAD2 enzyme in S. cerevisiae led to the synthesis of palmitolinoleic (16:2∆^9,12) acid, an unusual fatty acid that is not found naturally in Leucas cephalotes. While the LcFAD3 enzyme catalyzed linoleic acid (18:2∆^9,12) into α-linolenic acid (18:3∆^9,12,15). Furthermore, transcript abundance of LcFAD2 and LcFAD3 cDNAs were estimated from various plant parts such as roots, shoots, leaves, petals and developing seeds. Our results have shown that the differential transcriptional activity of LcFAD2 and LcFAD3 desaturase genes differs significantly in developing seeds, petals, leaves, stems, and roots of L. cephalotes. Furthermore, for the industrial production of these essential fatty acids, namely, linoleic and α-linolenic acid, FAD2 and FAD3 enzyme activity could be exploited from this upcoming significant oil plant, Leucas cephalotes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01016-z.

A robust and practical CRISPR/crRNA screening system for soybean cultivar editing using LbCpf1 ribonucleoproteins

KEY MESSAGE

Calli protoplasts isolated from three soybean cultivars are useful tools to evaluate guide RNAs for clustered regularly interspaced short palindromic repeats (CRISPR)-based precise gene editing. A type V CRISPR effector, LbCpf1(Cas12a) from Lachnospiraceae bacterium ND 2006, has been used for precision editing of the plant genome. We report that callus-derived protoplasts from three soybeans, including Glycine Max var. Williams 82 and two Korean cultivars (Kwangan and Daewon) represent efficient systems for the screening of active crRNA for CRISPR/LbCpf1. CRISPR/LbCpf1 ribonucleoproteins (RNPs) were delivered as complexes of purified endonucleases mixed with designed crRNA to simultaneously edit target genes of GlymaFAD2-1A and GlymaFAD2-1B transfected into three soybean protoplasts including genome-sequenced Williams 82 with cultivars, Kwangan and Daewon. Previously, we reported that nine crRNAs designed for LbCpf1 exhibited varying degrees of editing efficacy for two FAD2 genes. Among the nine crRNAs, the LbCpf1-crRNA3 complexes showed the highest efficiency in soybean cotyledon protoplasts. The new screening systems of callus protoplasts from three soybeans have been successfully used to transfect GFP-tagged markers and CRISPR/LbCpf1 RNPs. The callus protoplasts confirm that the LbCpf1-crRNA3 complex is an active crRNA for LbCpf1 to edit two FAD2 genes similar to cotyledon protoplasts. These results demonstrate that soybean callus protoplast-based CRISPR/crRNA selection is a new and practical tool to screen the efficacy of crRNAs and a prerequisite for progressive regeneration of the edited soybean.

Genetic Engineering of Lesquerella with Increased Ricinoleic Acid Content in Seed Oil

Seeds of castor (Ricinus communis) are enriched in oil with high levels of the industrially valuable fatty acid ricinoleic acid (18:1OH), but production of this plant is limited because of the cooccurrence of the ricin toxin in its seeds. Lesquerella (Physaria fendleri) is being developed as an alternative industrial oilseed because its seeds accumulate lesquerolic acid (20:1OH), an elongated form of 18:1OH in seed oil which lacks toxins. Synthesis of 20:1OH is through elongation of 18:1OH by a lesquerella elongase, PfKCS18. Oleic acid (18:1) is the substrate for 18:1OH synthesis, but it is also used by fatty acid desaturase 2 (FAD2) and FAD3 to sequentially produce linoleic and linolenic acids. To develop lesquerella that produces 18:1OH-rich seed oils such as castor, RNA interference sequences targeting KCS18, FAD2 and FAD3 were introduced to lesquerella to suppress the elongation and desaturation steps. Seeds from transgenic lines had increased 18:1OH to 1.1-26.6% compared with that of 0.4-0.6% in wild-type (WT) seeds. Multiple lines had reduced 18:1OH levels in the T₂ generation, including a top line with 18:1OH reduced from 26.7% to 19%. Transgenic lines also accumulated more 18:1 than that of WT, indicating that 18:1 is not efficiently used for 18:1OH synthesis and accumulation. Factors limiting 18:1OH accumulation and new targets for further increasing 18:1OH production are discussed. Our results provide insights into complex mechanisms of oil biosynthesis in lesquerella and show the biotechnological potential to tailor lesquerella seeds to produce castor-like industrial oil functionality.

Ecotopic over-expression of PoCHS from Paeonia ostii altered the fatty acids composition and content in Arabidopsis thaliana

Chalcone synthase (CHS) is the key enzyme in the flavonoid biosynthetic pathway and has been studied in many plants, but the function of the CHS gene has not been well characterized in Paeonia ostii. In this study, we obtained a CHS homolog gene from P. ostii, which possessed the putative conserved amino acids of chalcone synthase by multiple alignment analysis and demonstrated the highest expression in developing seeds. In vitro assays of the recombinant PoCHS protein confirmed enzymatic activity using malonyl-CoA and 4-coumaroyl-CoA as substrates, and the optimal pH and reaction temperature were 7.5 and 40 °C, respectively. Furthermore, ectopic over-expression of PoCHS in Arabidopsis up-regulated the expression levels of genes involved in seed development (ABI), glycolysis (PKp2, PDH-E1a, and SUS2/3), and especially fatty acid biosynthesis (BCCP2, CAC2, CDS2, FatA, and FAD3). This resulted in an increased unsaturated fatty acid content, especially α-linolenic acid, in transgenic Arabidopsis seeds. In this study, we examined the functions of CHS homolog of P. ostii and demonstrated its new function in seed fatty acid biosynthesis.

TILLING-by-Sequencing+ to Decipher Oil Biosynthesis Pathway in Soybeans: A New and Effective Platform for High-Throughput Gene Functional Analysis

Reverse genetic approaches have been widely applied to study gene function in crop species; however, these techniques, including gel-based TILLING, present low efficiency to characterize genes in soybeans due to genome complexity, gene duplication, and the presence of multiple gene family members that share high homology in their DNA sequence. Chemical mutagenesis emerges as a genetically modified-free strategy to produce large-scale soybean mutants for economically important traits improvement. The current study uses an optimized high-throughput TILLING by target capture sequencing technology, or TILLING-by-Sequencing+ (TbyS+), coupled with universal bioinformatic tools to identify population-wide mutations in soybeans. Four ethyl methanesulfonate mutagenized populations (4032 mutant families) have been screened for the presence of induced mutations in targeted genes. The mutation types and effects have been characterized for a total of 138 soybean genes involved in soybean seed composition, disease resistance, and many other quality traits. To test the efficiency of TbyS+ in complex genomes, we used soybeans as a model with a focus on three desaturase gene families, GmSACPD, GmFAD2, and GmFAD3, that are involved in the soybean fatty acid biosynthesis pathway. We successfully isolated mutants from all the six gene family members. Unsurprisingly, most of the characterized mutants showed significant changes either in their stearic, oleic, or linolenic acids. By using TbyS+, we discovered novel sources of soybean oil traits, including high saturated and monosaturated fatty acids in addition to low polyunsaturated fatty acid contents. This technology provides an unprecedented platform for highly effective screening of polyploid mutant populations and functional gene analysis. The obtained soybean mutants from this study can be used in subsequent soybean breeding programs for improved oil composition traits.

Dynamic transcriptome analysis identifies genes related to fatty acid biosynthesis in the seeds of Prunus pedunculata Pall

BACKGROUND

Prunus pedunculata Pall, the deciduous shrub of Amygdalus subgenus in Rosaceae, is a new kind of desert oil-bearing tree. It has a long story of being planted in the West and North of China for sand fixation and desert control. In addition, the seeds of P. pedunculata are rich of oil, especially the monounsaturated fatty acid and polyunsaturated fatty acid. However, little is known about the molecular mechanisms of oil accumulation during the seed development of P. pedunculata.

RESULTS

The seeds of P. pedunculata from three independent plants at 10, 18, 24, 31, 39, 45, 59 and 73 days after flowering (DAF) were obtained and the oil compositions were evaluated. It showed that oleic acid was the dominant type of oil content in the mature seeds (from 32.724% at 10DAF to 72.06% at 73DAF). Next, transcriptome sequencing for the developing seeds produced 988.795 million high quality reads and TRINITY assembled 326,271 genes for the first transcriptome for P. pedunculata. After the assembled transcriptome was evaluated by BUSCO with 85.9% completeness, we identified 195,342, 109,850 and 121,897 P. pedunculata genes aligned to NR, GO and KEGG pathway databases, respectively. Then, we predicted 23,229 likely proteins from the assembled transcriptome and identified 1917 signal peptides and 5512 transmembrane related proteins. In the developing seeds we detected 91,362 genes (average FPKM > 5) and correlation analysis indicated three possible development stages - early (10 ~ 24DAF), middle (31 ~ 45DAF) and late (59 ~ 73DAF). We next analyzed the differentially expressed genes (DEGs) in the developing seeds. Interestingly, compared to 10DAF the number of DEGs was increased from 4406 in 18DAF to 27,623 in 73DAF. Based on the gene annotation, we identified 753, 33, 8 and 645 DEGs related to the fatty acid biosynthesis, lipid biosynthesis, oil body and transcription factors. Notably, GPAT, DGD1, LACS2, UBC and RINO were highly expressed at the early development stage, ω6-FAD, SAD, ACP, ACCA and AHG1 were highly expressed at the middle development stage, and LACS6, DGD1, ACAT1, AGPAT, WSD1, EGY2 and oleosin genes were highly expressed at the late development stage.

CONCLUSIONS

This is the first time to study the developing seed transcriptome of P. pedunculata and our findings will provide a valuable resource for future studies. More importantly, it will improve our understanding of molecular mechanisms of oil accumulation in P. pedunculata.

Analyzing Mass Spectrometry Imaging Data of 13C-Labeled Phospholipids in Camelina sativa and Thlaspi arvense (Pennycress) Embryos

The combination of ¹³C-isotopic labeling and mass spectrometry imaging (MSI) offers an approach to analyze metabolic flux in situ. However, combining isotopic labeling and MSI presents technical challenges ranging from sample preparation, label incorporation, data collection, and analysis. Isotopic labeling and MSI individually create large, complex data sets, and this is compounded when both methods are combined. Therefore, analyzing isotopically labeled MSI data requires streamlined procedures to support biologically meaningful interpretations. Using currently available software and techniques, here we describe a workflow to analyze ¹³C-labeled isotopologues of the membrane lipid and storage oil lipid intermediate―phosphatidylcholine (PC). Our results with embryos of the oilseed crops, Camelina sativa and Thlaspi arvense (pennycress), demonstrated greater ¹³C-isotopic labeling in the cotyledons of developing embryos compared with the embryonic axis. Greater isotopic enrichment in PC molecular species with more saturated and longer chain fatty acids suggest different flux patterns related to fatty acid desaturation and elongation pathways. The ability to evaluate MSI data of isotopically labeled plant embryos will facilitate the potential to investigate spatial aspects of metabolic flux in situ.

Identification of genes associated with the biosynthesis of unsaturated fatty acid and oil accumulation in herbaceous peony 'Hangshao' (Paeonia lactiflora 'Hangshao') seeds based on transcriptome analysis

Background

Paeonia lactiflora ‘Hangshao’ is widely cultivated in China as a traditional Chinese medicine ‘Radix Paeoniae Alba’. Due to the abundant unsaturated fatty acids in its seed, it can also be regarded as a new oilseed plant. However, the process of the biosynthesis of unsaturated fatty acids in it has remained unknown. Therefore, transcriptome analysis is helpful to better understand the underlying molecular mechanisms.

Results

Five main fatty acids were detected, including stearic acid, palmitic acid, oleic acid, linoleic acid and α-linolenic acid, and their absolute contents first increased and then decreased during seed development. A total of 150,156 unigenes were obtained by transcriptome sequencing. There were 15,005 unigenes annotated in the seven functional databases, including NR, NT, GO, KOG, KEGG, Swiss-Prot and InterPro. Based on the KEGG database, 1766 unigenes were annotated in the lipid metabolism. There were 4635, 12,304, and 18,291 DEGs in Group I (60 vs 30 DAF), Group II (90 vs 60 DAF) and Group III (90 vs 30 DAF), respectively. A total of 1480 DEGs were detected in the intersection of the three groups. In 14 KEGG pathways of lipid metabolism, 503 DEGs were found, belonging to 111 enzymes. We screened out 123 DEGs involved in fatty acid biosynthesis (39 DEGs), fatty acid elongation (33 DEGs), biosynthesis of unsaturated fatty acid (24 DEGs), TAG assembly (17 DEGs) and lipid storage (10 DEGs). Furthermore, qRT-PCR was used to analyze the expression patterns of 16 genes, including BBCP, BC, MCAT, KASIII, KASII, FATA, FATB, KCR, SAD, FAD2, FAD3, FAD7, GPAT, DGAT, OLE and CLO, most of which showed the highest expression at 45 DAF, except for DGAT, OLE and CLO, which showed the highest expression at 75 DAF.

Conclusions

We predicted that MCAT, KASIII, FATA, SAD, FAD2, FAD3, DGAT and OLE were the key genes in the unsaturated fatty acid biosynthesis and oil accumulation in herbaceous peony seed. This study provides the first comprehensive genomic resources characterizing herbaceous peony seed gene expression at the transcriptional level. These data lay the foundation for elucidating the molecular mechanisms of fatty acid biosynthesis and oil accumulation for herbaceous peony.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07339-7.

Transcriptome Analysis and Identification of Lipid Genes in Physaria lindheimeri, a Genetic Resource for Hydroxy Fatty Acids in Seed Oil

Hydroxy fatty acids (HFAs) have numerous industrial applications but are absent in most vegetable oils. Physaria lindheimeri accumulating 85% HFA in its seed oil makes it a valuable resource for engineering oilseed crops for HFA production. To discover lipid genes involved in HFA synthesis in P. lindheimeri, transcripts from developing seeds at various stages, as well as leaf and flower buds, were sequenced. Ninety-seven percent clean reads from 552,614,582 raw reads were assembled to 129,633 contigs (or transcripts) which represented 85,948 unique genes. Gene Ontology analysis indicated that 60% of the contigs matched proteins involved in biological process, cellular component or molecular function, while the remaining matched unknown proteins. We identified 42 P. lindheimeri genes involved in fatty acid and seed oil biosynthesis, and 39 of them shared 78-100% nucleotide identity with Arabidopsis orthologs. We manually annotated 16 key genes and 14 of them contained full-length protein sequences, indicating high coverage of clean reads to the assembled contigs. A detailed profiling of the 16 genes revealed various spatial and temporal expression patterns. The further comparison of their protein sequences uncovered amino acids conserved among HFA-producing species, but these varied among non-HFA-producing species. Our findings provide essential information for basic and applied research on HFA biosynthesis.

CRISPR/Cas9-Induced fad2 and rod1 Mutations Stacked With fae1 Confer High Oleic Acid Seed Oil in Pennycress (Thlaspi arvense L.)

Pennycress (Thlaspi arvense L.) is being domesticated as an oilseed cash cover crop to be grown in the off-season throughout temperate regions of the world. With its diploid genome and ease of directed mutagenesis using molecular approaches, pennycress seed oil composition can be rapidly tailored for a plethora of food, feed, oleochemical and fuel uses. Here, we utilized Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 technology to produce knockout mutations in the FATTY ACID DESATURASE2 (FAD2) and REDUCED OLEATE DESATURATION1 (ROD1) genes to increase oleic acid content. High oleic acid (18:1) oil is valued for its oxidative stability that is superior to the polyunsaturated fatty acids (PUFAs) linoleic (18:2) and linolenic (18:3), and better cold flow properties than the very long chain fatty acid (VLCFA) erucic (22:1). When combined with a FATTY ACID ELONGATION1 (fae1) knockout mutation, fad2 fae1 and rod1 fae1 double mutants produced ∼90% and ∼60% oleic acid in seed oil, respectively, with PUFAs in fad2 fae1 as well as fad2 single mutants reduced to less than 5%. MALDI-MS spatial imaging analyses of phosphatidylcholine (PC) and triacylglycerol (TAG) molecular species in wild-type pennycress embryo sections from mature seeds revealed that erucic acid is highly enriched in cotyledons which serve as storage organs, suggestive of a role in providing energy for the germinating seedling. In contrast, PUFA-containing TAGs are enriched in the embryonic axis, which may be utilized for cellular membrane expansion during seed germination and seedling emergence. Under standard growth chamber conditions, rod1 fae1 plants grew like wild type whereas fad2 single and fad2 fae1 double mutant plants exhibited delayed growth and overall reduced heights and seed yields, suggesting that reducing PUFAs below a threshold in pennycress had negative physiological effects. Taken together, our results suggest that combinatorial knockout of ROD1 and FAE1 may be a viable route to commercially increase oleic acid content in pennycress seed oil whereas mutations in FAD2 will likely require at least partial function to avoid fitness trade-offs.