Genome-Wide Survey and Characterization of Fatty Acid Desaturase Gene Family in Brassica napus and Its Parental Species

Brassica napus cytochrome P450 superfamily: Origin from parental species and involvement in diseases resistance, abiotic stresses tolerance, and seed quality traits

Cytochromes P450 monooxygenases (CYP450s) constitute the largest enzymic protein family that is widely present in plants, animals, and microorganisms, participate in numerous metabolic pathways, and play diverse roles in development, metabolism, and defense. Rapeseed (Brassica napus) is an important oil crop worldwide and have many versions of reference genome. However, there is no systemically comparative genome-wide analysis of CYP450 family genes in rapeseed and its parental species B. rapa and B. oleracea. In this study, we identified 765, 293 and 437 CYP450 genes in B. napus, B. rapa and B. oleracea, respectively, which were unevenly located in A01-A10 and/or C01-C09 chromosomes in corresponding species. Phylogenetic relationship analysis indicated that 1745 CYP450 proteins from three Brassica species and Arabidopsis were divided into 4 groups. Whole genome duplication (WGD) or segmental duplication resulted in gene expansion of CYP450 family in three Brassica species. There were 33-83 SSR loci in CYP450 genes of three Brassica species, and numerous transcription factor binding sites were identified in their promoters. A total of 459-777 miRNAs were predicted to target 174-426 CYP450 genes in three Brassica species. Based on transcriptome data, BnCYP450s, BrCYP450s and BoCYP450s were differentially expressed in various tissues. There existed numerous BnCYP450 DEGs in response to pathogens and abiotic stresses. Besides, many BnCYP450 DEGs were involved in the regulation of important traits, such as seed germination, seed ALA content, and yellow-seed. The qRT-PCR experiment confirmed the transcriptome analysis results by validating two representative Sclerotinia-responsive BnCYP450 DEGs as an example. Three BnCYP450s genes (CYP707A1, CYP81F1, CYP81H1) might be regulated by seed-specific transcription factors BnTT1 and BnbZIP67 to participate in the development and metabolism of seed coat and embryo by undertaking related metabolic reactions.

Comparative Genome-Wide Identification of the Fatty Acid Desaturase Gene Family in Tea and Oil Tea

Camellia oil is valuable as an edible oil and serves as a base material for a range of high-value products. Camellia plants of significant economic importance, such as Camellia sinensis and Camellia oleifera, have been classified into sect. Thea and sect. Oleifera, respectively. Fatty acid desaturases play a crucial role in catalyzing the formation of double bonds at specific positions of fatty acid chains, leading to the production of unsaturated fatty acids and contributing to lipid synthesis. Comparative genomics results have revealed that expanded gene families in oil tea are enriched in functions related to lipid, fatty acid, and seed processes. To explore the function of the FAD gene family, a total of 82 FAD genes were identified in tea and oil tea. Transcriptome data showed the differential expression of the FAD gene family in mature seeds of tea tree and oil tea tree. Furthermore, the structural analysis and clustering of FAD proteins provided insights for the further exploration of the function of the FAD gene family and its role in lipid synthesis. Overall, these findings shed light on the role of the FAD gene family in Camellia plants and their involvement in lipid metabolism, as well as provide a reference for understanding their function in oil synthesis.

In Silico Identification and Characterization of Fatty Acid Desaturase (FAD) Genes in Argania spinosa L. Skeels: Implications for Oil Quality and Abiotic Stress

Fatty acid desaturase (FAD) is the key enzyme that leads to the formation of unsaturated fatty acids by introducing double bonds into hydrocarbon chains, and it plays a critical role in plant lipid metabolism. However, no data are available on enzyme-associated genes in argan trees. In addition, a candidate gene approach was adopted to identify and characterize the gene sequences of interest that are potentially involved in oil quality and abiotic stress. Based on phylogenetic analyses, 18 putative FAD genes of Argania spinosa L. (AsFAD) were identified and assigned to three subfamilies: stearoyl-ACP desaturase (SAD), Δ-12 desaturase (FAD2/FAD6), and Δ-15 desaturase (FAD3/FAD7). Furthermore, gene structure and motif analyses revealed a conserved exon-intron organization among FAD members belonging to the various oil crops studied, and they exhibited conserved motifs within each subfamily. In addition, the gene structure shows a wide variation in intron numbers, ranging from 0 to 8, with two highly conserved intron phases (0 and 1). The AsFAD and AsSAD subfamilies consist of three (H(X)2-4H, H(X)2-3HH, and H/Q (X)2-3HH) and two (EEN(K)RHG and DEKRHE) conserved histidine boxes, respectively. A set of primer pairs were designed for each FAD gene, and tested on DNA extracted from argan leaves, in which all amplicons of the expected size were produced. These findings of candidate genes in A spinosa L. will provide valuable knowledge that further enhances our understanding of the potential roles of FAD genes in the quality of oil and abiotic stress in the argan tree.

Integrated analysis of metabolome, transcriptome, and bioclimatic factors of Acer truncatum seeds reveals key candidate genes related to unsaturated fatty acid biosynthesis, and potentially optimal production area

BACKGROUND

Lipids found in plant seeds are essential for controlling seed dormancy, dispersal, and defenses against biotic and abiotic stress. Additionally, these lipids provide nutrition and energy and are therefore important to the human diet as edible oils. Acer truncatum, which belongs to the Aceaceae family, is widely cultivated around the world for its ornamental value. Further because its seed oil is rich in unsaturated fatty acids (UFAs)- i.e. α-linolenic acid (ALA) and nervonic acid (NA)- and because it has been validated as a new food resource in China, the importance of A. truncatum has greatly risen. However, it remains unknown how UFAs are biosynthesized during the growth season, to what extent environmental factors impact their content, and what areas are potentially optimal for their production.

RESULTS

In this study, transcriptome and metabolome of A. truncatum seeds at three representative developmental stages was used to find the accumulation patterns of all major FAs. Cumulatively, 966 metabolites and 87,343 unigenes were detected; the differential expressed unigenes and metabolites were compared between stages as follows: stage 1 vs. 2, stage 1 vs. 3, and stage 2 vs. 3 seeds, respectively. Moreover, 13 fatty acid desaturases (FADs) and 20 β-ketoacyl-CoA synthases (KCSs) were identified, among which the expression level of FAD3 (Cluster-7222.41455) and KCS20 (Cluster-7222.40643) were consistent with the metabolic results of ALA and NA, respectively. Upon analysis of the geographical origin-affected diversity from 17 various locations, we found significant variation in phenotypes and UFA content. Notably, in this study we found that 7 bioclimatic variables showed considerable influence on FAs contents in A. truncatum seeds oil, suggesting their significance as critical environmental parameters. Ultimately, we developed a model for potentially ecological suitable regions in China.

CONCLUSION

This study provides a comprehensive understanding of the relationship between metabolome and transcriptome in A. truncatum at various developmental stages of seeds and a new strategy to enhance seed FA content, especially ALA and NA. This is particularly significant in meeting the increasing demands for high-quality edible oil for human consumption. The study offers a scientific basis for A. truncatum's novel utilization as a woody vegetable oil rather than an ornamental plant, potentially expanding its cultivation worldwide.

Genome-Wide Characterization of Tomato FAD Gene Family and Expression Analysis under Abiotic Stresses

The fatty acid desaturase (FAD) gene family plays a crucial regulatory role in the resistance process of plant biomembranes. To understand the role of FADs in tomato growth and development, this study identified and analyzed the tomato FAD gene family based on bioinformatics analysis methods. In this study, 26 SlFADs were unevenly distributed on 10 chromosomes. Phylogenetic analysis showed that the SlFAD gene family was divided into six branches, and the exon-intron composition and conserved motifs of SlFADs clustered in the same branch were quite conservative. Several hormone and stress response elements in the SlFAD promoter suggest that the expression of SlFAD members is subject to complex regulation; the construction of a tomato FAD protein interaction network found that SlFAD proteins have apparent synergistic effects with SPA and GPAT proteins. qRT-PCR verification results show that SlFAD participates in the expression of tomato root, stem, and leaf tissues; SlFAD8 is mainly highly expressed in leaves; SlFAD9 plays a vital role in response to salt stress; and SlFAB5 regulates all stages of fruit development under the action of exogenous hormones. In summary, this study provides a basis for a systematic understanding of the SlFAD gene family. It provides a theoretical basis for in-depth research on the functional characteristics of tomato SlFAD genes.

A novel type of Brassica napus with higher stearic acid in seeds developed through genome editing of BnaSAD2 family

KEY MESSAGE

Modifications of multiple copies of the BnaSAD2 gene family with genomic editing technology result in higher stearic acid content in the seed of polyploidy rapeseed. Solid fats from vegetable oils are widely used in food processing industry. Accumulating data showed that stearic acid is more favorite as the major composite among the saturate fatty acids in solid fats in considerations of its effects on human health. Rapeseed is the third largest oil crop worldwide, and has potential to be manipulated to produce higher saturated fatty acids as raw materials of solid fats. Toward that end, we identified four SAD2 gene family members in B. napus genome and established spatiotemporal expression pattern of the BnaSAD2 members. Genomic editing technology was applied to mutate all the copies of BnaSAD2 in this allopolyploid species and mutants at multiple alleles were generated and characterized to understand the effect of each BnaSAD2 member on blocking desaturation of stearic acid. Mutations occurred at BnaSAD2.A3 resulted in more dramatic changes of fatty acid profile than ones on BnaSAD2.C3, BnaSAD2.A5 and BnaSAD2.C4. The content of stearic acid in mutant seeds with single locus increased dramatically with a range of 3.1-8.2%. Furthermore, combination of different mutated alleles of BnaSAD2 resulted in more dramatic changes in fatty acid profiles and the double mutant at BnaSAD2.A3 and BnaSAD2.C3 showed the most dramatic phenotypic changes compared with its single mutants and other double mutants, leading to 11.1% of stearic acid in the seeds. Our results demonstrated that the members of BnaSAD2 have differentiated in their efficacy as a Δ9-Stearoyl-ACP-Desaturase and provided valuable rapeseed germplasm for breeding high stearic rapeseed oil.

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.

Genome-wide analysis of fatty acid desaturase genes in chia (Salvia hispanica) reveals their crucial roles in cold response and seed oil formation

Chia (Salvia hispanica) is a functional food crop with high α-linolenic acid (ALA), the omega-3 essential fatty acid, but its worldwide plantation is limited by cold-intolerance and strict short-photoperiod flowering feature. Fatty acid desaturases (FADs) are responsible for seed oil accumulation, and play important roles in cold stress tolerance of plants. To date, there is no report on systemically genome-wide analysis of FAD genes in chia (ShiFADs). In this study, 31 ShiFAD genes were identified, 3 of which contained 2 alternative splicing transcripts, and they were located in 6 chromosomes of chia. Phylogenetic analysis classified the ShiFAD proteins into 7 groups, with conserved gene structure and MEME motifs within each group. Tandem and segmental duplications coursed the expansion of ShiFAD genes. Numerous cis-regulatory elements, including hormone response elements, growth and development elements, biotic/abiotic stress response elements, and transcription factor binding sites, were predicted in ShiFAD promoters. 24 miRNAs targeting ShiFAD genes were identified at whole-genome level. In total, 15 SSR loci were predicted in ShiFAD genes/promoters. RNA-seq data showed that ShiFAD genes were expressed in various organs with different levels. qRT-PCR detection revealed the inducibility of ShiSAD2 and ShiSAD7 in response to cold stress, and validated the seed-specific expression of ShiSAD11a. Yeast expression of ShiSAD11a confirmed the catalytic activity of its encoded protein, and its heterologous expression in Arabidopsis thaliana significantly increased seed oleic acid content. This work lays a foundation for molecular dissection of chia high-ALA trait and functional study of ShiFAD genes in cold tolerance.

Comparative phylogenomic insights of KCS and ELO gene families in Brassica species indicate their role in seed development and stress responsiveness

Very long-chain fatty acids (VLCFAs) possess more than twenty carbon atoms and are the major components of seed storage oil, wax, and lipids. FAE (Fatty Acid Elongation) like genes take part in the biosynthesis of VLCFAs, growth regulation, and stress responses, and are further comprised of KCS (Ketoacyl-CoA synthase) and ELO (Elongation Defective Elongase) sub-gene families. The comparative genome-wide analysis and mode of evolution of KCS and ELO gene families have not been investigated in tetraploid Brassica carinata and its diploid progenitors. In this study, 53 KCS genes were identified in B. carinata compared to 32 and 33 KCS genes in B. nigra and B. oleracea respectively, which suggests that polyploidization might has impacted the fatty acid elongation process during Brassica evolution. Polyploidization has also increased the number of ELO genes in B. carinata (17) over its progenitors B. nigra (7) and B. oleracea (6). Based on comparative phylogenetics, KCS, and ELO proteins can be classified into eight and four major groups, respectively. The approximate date of divergence for duplicated KCS and ELO genes varied from 0.03 to 3.20 million years ago (MYA). Gene structure analysis indicated that the maximum number of genes were intron-less and remained conserved during evolution. The neutral type of selection seemed to be predominant in both KCS and ELO genes evolution. String-based protein-protein interaction analysis suggested that bZIP53, a transcription factor might be involved in the activation of transcription of ELO/KCS genes. The presence of biotic and abiotic stress-related cis-regulatory elements in the promoter region suggests that both KCS and ELO genes might also play their role in stress tolerance. The expression analysis of both gene family members reflect their preferential seed-specific expression, especially during the mature embryo development stage. Furthermore, some KCS and ELO genes were found to be specifically expressed under heat stress, phosphorus starvation, and Xanthomonas campestris infection. The current study provides a basis to understand the evolution of both KCS and ELO genes in fatty acid elongation and their role in stress tolerance.

Identification of Trehalose-6-Phosphate Synthase (TPS) Genes Associated with Both Source-/Sink-Related Yield Traits and Drought Response in Rapeseed (Brassica napus L.)

Trehalose-6-phosphate synthase (TPS) is an important enzyme for the synthesis of Trehalose-6-phosphate (T6P). In addition to being a signaling regulator of carbon allocation that improves crop yields, T6P also plays essential roles in desiccation tolerance. However, comprehensive studies, such as evolutionary analysis, expression analysis, and functional classification of the TPS family in rapeseed (Brassica napus L.) are lacking. Here, we identified 35 BnTPSs, 14 BoTPSs, and 17 BrTPSs in cruciferous plants, which were classified into three subfamilies. Phylogenetic and syntenic analysis of TPS genes in four cruciferous species indicated that only gene elimination occurred during their evolution. Combined phylogenetic, protein property, and expression analysis of the 35 BnTPSs suggested that changes in gene structures might have led to changes in their expression profiles and further functional differentiation during their evolution. In addition, we analyzed one set of transcriptome data from Zhongshuang11 (ZS11) and two sets of data from extreme materials associated with source-/sink-related yield traits and the drought response. The expression levels of four BnTPSs (BnTPS6, BnTPS8, BnTPS9, and BnTPS11) increased sharply after drought stress, and three differentially expressed genes (BnTPS1, BnTPS5, and BnTPS9) exhibited variable expression patterns among source and sink tissues between yield-related materials. Our findings provide a reference for fundamental studies of TPSs in rapeseed and a framework for future functional research of the roles of BnTPSs in both yield and drought resistance.

Genome-Wide Identification of Membrane-Bound Fatty Acid Desaturase Genes in Three Peanut Species and Their Expression in Arachis hypogaea during Drought Stress

As a crop irrigated primarily by rain, the quality and yield of peanuts are significantly limited by drought. To date, many studies have indicated that fatty acid desaturase (FAD) genes enhance plant tolerance to drought stresses. In this study, 16, 15, and 31 FADs were identified in Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively. All the FADs were divided into four subfamilies, which had relatively conserved gene structures, motifs, and domains. The synteny relationships and chromosomal position analysis showed that the FADs in subgenome pairs, A. duranensis-A. hypogaea (AA) and A. ipaensis-A. hypogaea (BB), were homologous, and their physical locations were consistent. The Ka/Ks results indicated that nine FAD genes underwent a purifying selection, and Ah|FAD3.2 experienced positive selection during tetraploid peanut speciation. Various cis-acting elements related to hormone signaling and stress responsiveness in promoters and the predicted miRNA targeting Ah|FADs suggested that these genes play crucial roles in drought tolerance. The expression profiles of Ah|FADs in 22 tissues and drought-tolerant and -sensitive cultivars under drought stress suggested that 4 and 6 FADs were putative genes related to oil accumulation and drought, respectively. These findings will help provide insight into the potential functional roles of the FAD genes, which may aid in dealing with plant drought stress.

Genome-Wide Characterization and Expression Analysis of Fatty acid Desaturase Gene Family in Poplar

Fatty acid desaturases (FADs) modulate carbon–carbon single bonds to form carbon–carbon double bonds in acyl chains, leading to unsaturated fatty acids (UFAs) that have vital roles in plant growth and development and their response to environmental stresses. In this study, a total of 23 Populus trichocarpaFAD (PtFAD) candidates were identified from the poplar genome and clustered into seven clades, including FAB2, FAD2, FAD3/7/8, FAD5, FAD6, DSD, and SLD. The exon–intron compositions and conserved motifs of the PtFADs, clustered into the same clade, were considerably conserved. It was found that segmental duplication events are predominantly attributable to the PtFAD gene family expansion. Several hormone- and stress-responsive elements in the PtFAD promoters implied that the expression of the PtFAD members was complicatedly regulated. A gene expression pattern analysis revealed that some PtFAD mRNA levels were significantly induced by abiotic stress. An interaction proteins and gene ontology (GO) analysis indicated that the PtFADs are closely associated with the UFAs biosynthesis. In addition, the UFA contents in poplars were significantly changed under drought and salt stresses, especially the ratio of linoleic and linolenic acids. The integration of the PtFAD expression patterns and UFA contents showed that the abiotic stress-induced PtFAD3/7/8 members mediating the conversion of linoleic and linolenic acids play vital roles in response to osmotic stress. This study highlights the profiles and functions of the PtFADs and identifies some valuable genes for forest improvements.

Genome-Wide Identification and Expression Analysis of nsLTP Gene Family in Rapeseed (Brassica napus) Reveals Their Critical Roles in Biotic and Abiotic Stress Responses

Non-specific lipid transfer proteins (nsLTPs) are small cysteine-rich basic proteins which play essential roles in plant growth, development and abiotic/biotic stress response. However, there is limited information about the nsLTP gene (BnLTP) family in rapeseed (Brassica napus). In this study, 283 BnLTP genes were identified in rapeseed, which were distributed randomly in 19 chromosomes of rapeseed. Phylogenetic analysis showed that BnLTP proteins were divided into seven groups. Exon/intron structure and MEME motifs both remained highly conserved in each BnLTP group. Segmental duplication and hybridization of rapeseed’s two sub-genomes mainly contributed to the expansion of the BnLTP gene family. Various potential cis-elements that respond to plant growth, development, biotic/abiotic stresses, and phytohormone signals existed in BnLTP gene promoters. Transcriptome analysis showed that BnLTP genes were expressed in various tissues/organs with different levels and were also involved in the response to heat, drought, NaCl, cold, IAA and ABA stresses, as well as the treatment of fungal pathogens (Sclerotinia sclerotiorum and Leptosphaeria maculans). The qRT-PCR assay validated the results of RNA-seq expression analysis of two top Sclerotinia-responsive BnLTP genes, BnLTP129 and BnLTP161. Moreover, batches of BnLTPs might be regulated by BnTT1 and BnbZIP67 to play roles in the development, metabolism or adaptability of the seed coat and embryo in rapeseed. This work provides an important basis for further functional study of the BnLTP genes in rapeseed quality improvement and stress resistance.

Three strategies of transgenic manipulation for crop improvement

Heterologous expression of exogenous genes, overexpression of endogenous genes, and suppressed expression of undesirable genes are the three strategies of transgenic manipulation for crop improvement. Up to 2020, most (227) of the singular transgenic events (265) of crops approved for commercial release worldwide have been developed by the first strategy. Thirty-eight of them have been transformed by synthetic sequences transcribing antisense or double-stranded RNAs and three by mutated copies for suppressed expression of undesirable genes (the third strategy). By the first and the third strategies, hundreds of transgenic events and thousands of varieties with significant improvement of resistance to herbicides and pesticides, as well as nutritional quality, have been developed and approved for commercial release. Their application has significantly decreased the use of synthetic pesticides and the cost of crop production and increased the yield of crops and the benefits to farmers. However, almost all the events overexpressing endogenous genes remain at the testing stage, except one for fertility restoration and another for pyramiding herbicide tolerance. The novel functions conferred by the heterologously expressing exogenous genes under the control of constitutive promoters are usually absent in the recipient crops themselves or perform in different pathways. However, the endogenous proteins encoded by the overexpressing endogenous genes are regulated in complex networks with functionally redundant and replaceable pathways and are difficult to confer the desirable phenotypes significantly. It is concluded that heterologous expression of exogenous genes and suppressed expression by RNA interference and clustered regularly interspaced short palindromic repeats-cas (CRISPR/Cas) of undesirable genes are superior to the overexpression of endogenous genes for transgenic improvement of crops.

Genome-Wide Identification and Functional Differentiation of Fatty Acid Desaturase Genes in Olea europaea L

Olive (Olea europaea L.) is a world-famous woody oil tree and popular for redundant unsaturated fatty acids. Fatty acid desaturase (FAD) genes are responsible for fatty acid desaturation and stress regulation but have not yet been identified in olive at the whole genome level. This study identified 40 and 27 FAD genes in the cultivated olive O. europaea cv. Farga and the wild olive O. europaea var. Sylvestris, respectively. Phylogenetic analysis showed that all the FAD genes could be classified into the soluble FAB2/SAD clade and membrane-bound clade, including ADS/FAD5, DES, FAD4, SLD, ω-6 and ω-3, with the high consistency of subcellular localization, motif composition and exon-intron organization in each group. FAD genes in olive showed the diverse functional differentiation in morphology of different tissues, fruit development and stress responses. Among them, OeFAB2.8 and OeFAD2.3 were up-regulated and OeADS.1, OeFAD4.1 and OeFAD8.2 were down-regulated under the wound, Verticillium dahliae and cold stresses. This study presents a comprehensive analysis of the FAD genes at the whole-genome level in olives and will provide guidance for the improvement of oil quality or stress tolerance of olive trees.

Genome-Wide Identification and Expression Analysis of SNARE Genes in Brassica napus

SNAREs (soluble N-ethylmaleimide-sensitive factor attachment protein receptors) are central components that drive membrane fusion events during exocytosis and endocytosis and play important roles in different biological processes of plants. In this study, we identified 237 genes encoding SNARE family proteins in B. napus in silico at the whole-genome level. Phylogenetic analysis showed that BnaSNAREs could be classified into five groups (Q (a-, b-, c-, bc-) and R) like other plant SNAREs and clustered into twenty-five subclades. The gene structure and protein domain of each subclade were found to be highly conserved. In many subclades, BnaSNAREs are significantly expanded compared with the orthologous genes in Arabidopsis thaliana. BnaSNARE genes are expressed differentially in the leaves and roots of B. napus. RNA-seq data and RT-qPCR proved that some of the BnaSNAREs are involved in the plant response to S. sclerotiorum infection as well as treatments with toxin oxalic acid (OA) (a virulence factor often secreted by S. sclerotiorum) or abscisic acid (ABA), methyl jasmonate (MeJA), and salicylic acid (SA), which individually promote resistance to S. sclerotiorum. Moreover, the interacted proteins of BnaSNAREs contain some defense response-related proteins, which increases the evidence that BnaSNAREs are involved in plant immunity. We also found the co-expression of BnaSYP121/2s, BnaSNAPs, and BnaVAMP722/3s in B. napus due to S. sclerotiorum infection as well as the probable interaction among them.

Genome-Wide Analysis and Expression Profile of Nuclear Factor Y (NF-Y) Gene Family in Z. jujuba

Nuclear factor-Y (NF-Y) is an important transcription factor in the plant species, which potentially provides a higher level of functional diversity including for abiotic stress tolerance. The genome-wide study and expression analysis of NF-Y gene family in Ziziphus, an elite abiotic stress-tolerant species, assist bioprospecting of genes. Here, a total of 32 NF-Y (8 NF-YA, 15 NF-YB, and 9 NF-YC) genes were identified in genome-wide search of Z. jujuba genome. Physicochemical properties, cellular localization, gene structure, chromosomal location, and protein motifs were analyzed for structural and functional understanding. Identified 12 NF-Ys were responsible for the expansion of NF-Y gene family by tandem duplication in Z. jujuba. Phylogenetic and comparative physical mapping of Z. jujuba NF-Ys with its orthologs illustrated evolutionary and functional insights into NF-Y gene family. A total of 45 perfect microsatellites (20bp to 40bp) were extracted across the ZjNF-Y genes. The promoter and gene ontology study suggested that Z. jujuba NF-Y gene family is functionally diverse and could play a wide-ranging role in plant abiotic stress, development, and cellular processes. An expression study revealed that large numbers of the NF-Ys are differentially expressed in response to drought and salinity. The total 15 and 18 ZjNF-Y genes that are upregulated under drought and salinity stress, respectively, are the potential candidates for further functional analysis for development of climate-resilient crops. The present study established a base for understanding the role of NF-Ys in Z. jujuba under abiotic stress conditions and paved a way for further research.

Identification and Structure Analysis of KCS Family Genes Suggest Their Reponding to Regulate Fiber Development in Long-Staple Cotton Under Salt-Alkaline Stress

Plant 3-ketoacyl-CoA synthase (KCS) gene family catalyzed a β ketoacyl-CoA synthase, which was the rate-limiting enzyme for the synthesis of very long chain fatty acids (VLCFAs). Gossypium barbadense was well-known not only for high-quality fiber, which was perceived as a cultivated species of Gossypium. In this study, a total of 131 KCS genes were identified in four cotton species, there were 38, 44, 26, 23 KCS genes in the G. barbadense, the G. hirsutum, the G. arboreum and G. raimondii, respectively. The gene structure and expression pattern were analyzed. GBKCS genes were divided into six subgroups, the chromosome distribution of members of the family were mapped. The prediction of cis-acting elements of the GBKCS gene promoters suggested that the GBKCS genes may be involved in hormone signaling, defense and the stress response. Collinearity analysis on the KCS genes of the four cotton species were formulated. Tandem duplication played an indispensable role in the evolution of the KCS gene family. Specific expression analysis of 20 GBKCS genes indicated that GBKCS gene were widely expressed in the first 25 days of fiber development. Among them, GBKCS3, GBKCS8, GBKCS20, GBKCS34 were expressed at a high level in the initial long-term level of the G. barbadense fiber. This study established a foundation to further understanding of the evolution of KCS genes and analyze the function of GBKCS genes.

Genome-wide identification of FAD gene family and their contributions to the temperature stresses and mutualistic and parasitic fungi colonization responses in banana

Fatty acid desaturase (FAD) plays important roles in plant growth and development and plant defense processes. In this study, we identified 27 MaFAD genes from the banana genome. According to the amino acid sequence similarities, their encoded proteins could be classified into five subfamilies. This classification is consistently supported by their gene and protein structures, conserved motifs and subcellular localizations. Segmental duplication events were found to play predominant roles in the MaFAD gene family expansion. Thirty miRNAs targeting MaFADs were identified and many hormone- and stress-responsive cis-acting elements and transcription factor binding sites (TFBSs) were identified in their promoters, indicating that the MaFADs expression regulation was very complicated. Gene expression analysis showed that some MaFADs showed significant differential expression in response to high and low temperature. FocTR4 influenced greatly the expression of several MaFADs and greatly induced the fatty acid (FA) accumulations in roots. Although S. indica showed no significant influence on the expression of most MaFADs, it could greatly alleviate the influence of FocTR4 on several MaFADs and FA biosynthesis. Our study revealed that MaFADs contributed greatly to the responses of high and low temperature stresses and mutualistic and parasitic fungi colonization in banana.

Interference with ACSL1 gene in bovine adipocytes: Transcriptome profiling of circRNA related to unsaturated fatty acid production

Long-chain acyl-CoA synthetase 1 (ACSL1) is a member of the acyl-CoA synthetase family that plays a vital role in lipid metabolism. We have previously shown that the ACSL1 gene regulates the composition of unsaturated fatty acids (UFAs) in bovine skeletal muscle, which in turn regulates the fatty acid synthesis and the generation of lipid droplets. Here, we used RNA-Seq to screen circRNAs that regulated the expression of ACSL1 gene and other UFA synthesis-related genes by RNA interference and noninterference in bovine adipocytes. The results of KEGG pathway analysis showed that the parental genes of differentially expressed (DE)-circRNAs were primarily enriched in the adipocytokine signaling pathway. The prediction results showed that novel_circ_0004855, novel_circ_0001507, novel_circ_0001731, novel_circ_0005276, novel_circ_0002060, novel_circ_0005405 and novel_circ_0004254 regulated UFA synthesis-related genes by interacting with the related miRNAs. These results could help expand our knowledge of the molecular mechanisms of circRNAs in the regulation of UFA synthesis in bovine adipocytes.

Genome-Wide Identification, Evolution, and Comparative Analysis of B-Box Genes in Brassica rapa, B. oleracea, and B. napus and Their Expression Profiling in B. rapa in Response to Multiple Hormones and Abiotic Stresses

The B-box zinc-finger transcription factors are important for plant growth, development, and various physiological processes such as photomorphogenesis, light signaling, and flowering, as well as for several biotic and abiotic stress responses. However, there is relatively little information available regarding Brassica B-box genes and their expression. In this study, we identified 51, 52, and 101 non-redundant genes encoding B-box proteins in Brassica rapa (BrBBX genes), B. oleracea (BoBBX genes), and B. napus (BnBBX genes), respectively. A whole-genome identification, characterization, and evolutionary analysis (synteny and orthology) of the B-box gene families in the diploid species B. rapa (A genome) and B. oleracea (C genome) and in the allotetraploid species B. napus (AC genome) revealed segmental duplications were the major contributors to the expansion of the BrassicaBBX gene families. The BrassicaBBX genes were classified into five subgroups according to phylogenetic relationships, gene structures, and conserved domains. Light-responsive cis-regulatory elements were detected in many of the BBX gene promoters. Additionally, BrBBX expression profiles in different tissues and in response to various abiotic stresses (heat, cold, salt, and drought) or hormones (abscisic acid, methyl jasmonate, and gibberellic acid) were analyzed by qRT-PCR. The data indicated that many B-box genes (e.g., BrBBX13, BrBBX15, and BrBBX17) may contribute to plant development and growth as well as abiotic stress tolerance. Overall, the identified BBX genes may be useful as functional genetic markers for multiple stress responses and plant developmental processes.

Based on the whole genome clarified the evolution and expression process of fatty acid desaturase genes in three soybeans

Soybean is an important oil crop cultivated worldwide. With the increasing global population crossed with growing challenging cultivation conditions, improving soybean breeding by selecting important traits is urgent needed. Genes coding for plant fatty acid desaturases (FADs) genes are major candidates for that, because they are involving in controlling fatty acid composition and holding membrane fluidity under abiotic stress. Here, 75 FADs were found in three soybean genomes, which were further classified into four sub-groups. Phylogenetic tree, gene structure, motif and promoter analysis showed that the FAD gene family was conserved in the three soybeans. In addition, the numbers of omega desaturase from Chinese cultivated varieties were significantly higher than those in Chinese wild soybean and ancient polyploid soybean, respectively. However, it was the opposite for the sphingolipid subfamily. These results indicated that each subfamily was subjected to different selection pressures during cultivation and domestication. As the extra genes of the subfamily were very close to other family members' positions on chromosomes, they should be produced by duplication. The cis-element analysis of FAD promoter sequences revealed that upstream sequences of FAD contained abundant light, hormone and abiotic stress responsive cis-elements, suggesting that the quality of soybean could be improved by regulating these stresses. Expression analysis of Chinese wild soybean under salt stress showed that GsDES1.1, GsDES1.2, GsFAD2.1 and GsSLD1 in leaves and GsSLD2, GsSLD5 and GsSLD6 in roots were not closely related to salt stress response. Therefore, we explored the significant role of conserved, duplicated and neofunctionalized FAD in the domestication of soybean, which contributes to the importance of soybean as a global oil crop.

Genome-wide analysis of Hsp70 and Hsp100 gene families in Ziziphus jujuba

The Ziziphus species are naturally tolerant to a range of abiotic stresses. Therefore, it is expected that they are an enriched source of genes conferring stress tolerance. Heat shock proteins (Hsps) play a significant role in plants in imparting tolerance against abiotic stress conditions. To get an insight into potential Hsp function in Ziziphus, we performed a genome-wide analysis and expression study of Hsp70 and Hsp100 gene families in Ziziphus jujuba. We identified 21 and 6 genes of the ZjHsp70 and ZjHsp100 families, respectively. Physiochemical properties, chromosomal location, gene structure, motifs, and protein domain organization were analysed for structural and functional characterization. We identified the contribution of tandem and segmental gene duplications in expansions of ZjHsp70s and ZjHsp100s in Z. jujuba. Promoter analysis suggested that ZjHsp70s and ZjHsp100s perform diverse functions related to abiotic stress. Furthermore, expression analyses revealed that most of the Z. jujuba Hsp genes are differentially expressed in response to heat, drought, and salinity stress. Our analyses suggested ZjHsp70-3, ZjHsp70-5, ZjHsp70-6, ZjHsp70-16, ZjHsp70-17, ZjHsp70-20, ZjHsp100-1, ZjHsp100-2, and ZjHsp100-3 are potential candidates for further functional analysis and with regard to breeding new more resilient strains. The present analysis laid the foundation for understanding the molecular mechanism of Hsps70 and Hsp100 gene families regulating abiotic stress tolerance in Z. jujuba.

Genome-wide investigation and expression analysis of membrane-bound fatty acid desaturase genes under different biotic and abiotic stresses in sunflower (Helianthus annuus L.)

Membrane-bound fatty acid desaturase (FAD) gene family plays crucial roles in regulation of fatty acid (FA) compositions in plants. Sunflower (Helianthus annuus L.) is an important oilseed crop in the world; however, no comprehensive study on exploring the role of FAD family in relation to stress tolerance in sunflower has been performed yet. In this study, we identified 40 putative FAD genes in H. annuus (HaFAD), which were unevenly distributed across 13 of the total 17 chromosomes. Phylogenetic analysis indicated that HaFAD genes were divided into four subfamilies, as supported by highly conserved gene structures and motifs. Collinearity analysis showed that tandem duplication events played a crucial role in the expansion of HaFAD gene family. In addition, tissue-specific expression showed that 32 HaFAD genes were widely expressed in various tissues or organs of sunflower. Furthermore, qRT-PCR results revealed significant expression changes of HaFAD genes in response to abiotic (cadmium, drought) and biotic (Orobanche cumana) stresses, suggesting their important functions in response to different stresses. Therefore, our results provide insights into HaFAD gene family in response to different stresses, and some specific up-regulated genes such as HaFAD3.2, HaADS8, HaFAD2.1, and HaADS9 would be the potential candidate genes for the sunflower tolerance breeding.

Genome-Wide Analysis of Heat Shock Transcription Factors in Ziziphus jujuba Identifies Potential Candidates for Crop Improvement Under Abiotic Stress

Plant heat shock transcription factors (Hsfs) play a significant role in adoption under abiotic stress conditions by modulating the expression of several stress-responsive genes. Analysis of the Hsf gene family will serve to understand the molecular mechanism which is involved in response to abiotic stress. The Ziziphus species grows in warm and dry regions and is inherently tolerant to abiotic stress conditions; thus, Ziziphus is a highly enriched source of genes conferring abiotic stress tolerance. Therefore, the present study provides a comprehensive genome-wide analysis of the Hsf gene family in Z. jujuba. Identified 21 non-redundant Hsf genes were grouped into three major classes (classes A, B, and C) based on the phylogenetic analysis. Promoter and gene ontology analysis suggested that ZjHsfs perform diverse functions in response to abiotic stress conditions. Two paralogous pairs resulting from tandem gene duplication events were identified. Also, physio-chemical properties of chromosomal locations, gene structure, motifs, and protein domain organization of Hsfs were analyzed. Real-time PCR expression analyses revealed that most of the Z. jujuba Hsf genes are differentially expressed in response to heat stress. The analysis suggested ZjHsf-2, ZjHsf-3, ZjHsf-5, ZjHsf-7, ZjHsf-8, ZjHsf-10, ZjHsf-12, ZjHsf-17, and ZjHsf-18 were the outstanding candidate genes for imparting heat stress tolerance and for future functional analysis. The present analysis laid the foundation for understanding the molecular mechanism of the Hsf gene family regulating Z. jujuba development and tolerance to abiotic stress conditions.

Identification, evolution, expression, and docking studies of fatty acid desaturase genes in wheat (Triticum aestivum L.)

Backgrounds

Fatty acid desaturases (FADs) introduce a double bond into the fatty acids acyl chain resulting in unsaturated fatty acids that have essential roles in plant development and response to biotic and abiotic stresses. Wheat germ oil, one of the important by-products of wheat, can be a good alternative for edible oils with clinical advantages due to the high amount of unsaturated fatty acids. Therefore, we performed a genome-wide analysis of the wheat FAD gene family (TaFADs).

Results

68 FAD genes were identified from the wheat genome. Based on the phylogenetic analysis, wheat FADs clustered into five subfamilies, including FAB2, FAD2/FAD6, FAD4, DES/SLD, and FAD3/FAD7/FAD8. The TaFADs were distributed on chromosomes 2A-7B with 0 to 10 introns. The Ka/Ks ratio was less than one for most of the duplicated pair genes revealed that the function of the genes had been maintained during the evolution. Several cis-acting elements related to hormones and stresses in the TaFADs promoters indicated the role of these genes in plant development and responses to environmental stresses. Likewise, 72 SSRs and 91 miRNAs in 36 and 47 TaFADs have been identified. According to RNA-seq data analysis, the highest expression in all developmental stages and tissues was related to TaFAB2.5, TaFAB2.12, TaFAB2.15, TaFAB2.17, TaFAB2.20, TaFAD2.1, TaFAD2.6, and TaFAD2.8 genes while the highest expression in response to temperature stress was related to TaFAD2.6, TaFAD2.8, TaFAB2.15, TaFAB2.17, and TaFAB2.20. Furthermore, docking simulations revealed several residues in the active site of TaFAD2.6 and TaFAD2.8 in close contact with the docked oleic acid that could be useful in future site-directed mutagenesis studies to increase the catalytic efficiency of them and subsequently improve agronomic quality and tolerance of wheat against environmental stresses.

Conclusions

This study provides comprehensive information that can lead to the detection of candidate genes for wheat genetic modification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07199-1.

Transformation and Characterization of Δ12-Fatty Acid Acetylenase and Δ12-Oleate Desaturase Potentially Involved in the Polyacetylene Biosynthetic Pathway from Bidens pilosa

Bidens pilosa is commonly used as an herbal tea component or traditional medicine for treating several diseases, including diabetes. Polyacetylenes have two or more carbon–carbon triple bonds or alkynyl functional groups and are mainly derived from fatty acid and polyketide precursors. Here, we report the cloning of full-length cDNAs that encode Δ12-fatty acid acetylenase (designated BPFAA) and Δ12-oleate desaturase (designated BPOD) from B. pilosa, which we predicted to play a role in the polyacetylene biosynthetic pathway. Subsequently, expression vectors carrying BPFAA or BPOD were constructed and transformed into B. pilosa via the Agrobacterium-mediated method. Genomic PCR analysis confirmed the presence of transgenes and selection marker genes in the obtained transgenic lines. The copy numbers of transgenes in transgenic lines were determined by Southern blot analysis. Furthermore, 4–5 FAA genes and 2–3 OD genes were detected in wild-type (WT) plants. Quantitative real time-PCR revealed that some transgenic lines had higher expression levels than WT. Western blot analysis revealed OD protein expression in the selected transformants. High-performance liquid chromatography profiling was used to analyze the seven index polyacetylenic compounds, and fluctuation patterns were found.

The Use of Genetic and Gene Technologies in Shaping Modern Rapeseed Cultivars (Brassica napus L.)

Since their domestication, Brassica oilseed species have undergone progressive transformation allied with the development of breeding and molecular technologies. The canola (Brassica napus) crop has rapidly expanded globally in the last 30 years with intensive innovations in canola varieties, providing for a wider range of markets apart from the food industry. The breeding efforts of B. napus, the main source of canola oil and canola meal, have been mainly focused on improving seed yield, oil quality, and meal quality along with disease resistance, abiotic stress tolerance, and herbicide resistance. The revolution in genetics and gene technologies, including genetic mapping, molecular markers, genomic tools, and gene technology, especially gene editing tools, has allowed an understanding of the complex genetic makeup and gene functions in the major bioprocesses of the Brassicales, especially Brassica oil crops. Here, we provide an overview on the contributions of these technologies in improving the major traits of B. napus and discuss their potential use to accomplish new improvement targets.

Genome-wide mining and comparative analysis of fatty acid elongase gene family in Brassica napus and its progenitors

Very long chain fatty acids (VLCFAs) that are structural components of cell membrane lipid, cuticular waxes and seed oil, play crucial roles in plant growth, development and stress response. Fatty acid elongases (FAEs) comprising KCS and ELO, are key enzymes for VLCFA biosynthesis in plants. Although reference genomes of Brassica napus and its parental speices both have been sequenced, whole-genome analysis of FAE gene family in these Brassica speices is not reported. Here, 58, 33 and 30 KCS genes were identified in B. napus, B. rapa and B. oleracea genomes, respectively, whereas 14, 6 and 8 members were obtained for ELO genes. These KCS genes were unevenly located in 37 chromosomes and 3 scaffolds of 3 Brassica species, while these ELO genes were mapped to 19 chromosomes. The KCS and ELO proteins were divided into 8 and 4 subclasses, respectively. Gene structure and protein motifs remained highly conserved in each KCS or ELO subclass. Most promoters of KCS and ELO genes harbored various plant growth-, phytohormone-, and stress response-related cis-acting elements. 20 SSR loci existed in the KCS and ELO genes/promoters. The whole-genome duplication and segmental duplication mainly contributed to expansion of KCS and ELO genes in these genomes. Transcriptome analysis showed that KCS and ELO genes in 3 Brassica species were expressed in various tissues/organs with different levels, whereas 1 BnELO gene and 6 BnKCS genes might be pathogen-responsive genes. The qRT-PCR assay showed that BnKCS22 and BnELO04 responded to various phytohormone treatments and abiotic stresses. This work lays the foundation for further function identification of KCS and ELO genes in B. napus and its progenitors.

Influence of nitrogen availability on biomass, lipid production, fatty acid profile, and the expression of fatty acid desaturase genes in Messastrum gracile SE-MC4

In this study, the effects of limited and excess nitrate on biomass, lipid production, and fatty acid profile in Messastrum gracile SE-MC4 were determined. The expression of fatty acid desaturase genes, namely stearoyl-ACP desaturase (SAD), omega-6 fatty acid desaturase (ω-6 FAD), omega-3 fatty acid desaturase isoform 1 (ω-3 FADi1), and omega-3 fatty acid desaturase isoform 2 (ω-3 FADi2) was also assessed. It was found that nitrate limitation generally increased the total oil, α-linolenic acid (C18:3n3) and total polyunsaturated fatty acid (PUFA) contents in M. gracile. The reduction of nitrate concentration from 1.76 to 0.11 mM increased the total oil content significantly from 32.5 to 41.85% (dry weight). Palmitic (C16:0) and oleic (C18:1) acids as the predominant fatty acids in this microalgae constituted between 82 and 87% of the total oil content and were relatively consistent throughout all nitrate concentrations tested. The expression of SAD, ω-6 FAD, and ω-3 FADi2 genes increased under nitrate limitation, especially at 0.11 mM nitrate. The ω-3 FADi1 demonstrated a binary up-regulation pattern of expression under both nitrate-deficient (0.11 mM) and -excess (3.55 mM) conditions. Thus, findings from this study suggested that limited or excess nitrate could be used as part of a cultivation strategy to increase oil and PUFA content following media optimisation and more efficient culture methodology. Data obtained from the expression of desaturase genes would provide valuable insights into their roles under excess and limited nitrate conditions in M. gracile, potentially paving the way for future genetic modifications.