Acessibilidade / Reportar erro

Exploration of allelic diversity reveals a novel FAD2 (Oleate desaturase) gene in Brassica juncea

Abstract

Indian mustard (Brassica juncea) is an important contributor towards edible oil supply in India. Traditional Indian mustard varieties contain high proportion of 18C polyunsaturated fatty acids (linoleic and linolenic acids) and large amounts of long-chain monounsaturated FAs, mainly erucic acid, in seeds. Oleate desaturase (FAD2) regulates the composition of 18C PUFAs in cellular membranes and TAG in seed oil. The present study was conducted to gain insight into the allelic diversity of the FAD2 gene in Indian mustard. Analyses of cloned FAD2 genes of three Indian mustard varieties revealed a novel FAD2 gene that has a longer ORF (1167 bp) owing to insertions and several SNPs across its length that distinguish it from the more prevalent native FAD2 gene. Overall, the Indian mustard varieties possess three FAD2 alleles, but there is limited nucleotide diversity among members of each FAD2 type across varieties, suggesting narrow genetic diversity among the varieties examined.

Keywords:
FAD2; Indian mustard: Allelic diversity

INTRODUCTION

The composition of seed fatty acids (FAs) varies greatly depending on the plant species. Oil derived from vegetable sources is predominantly used for food-based applications, mainly for cooking purposes; therefore, the oxidative stability of the oil is one of the important properties. With the increase in unsaturation or introduction of double bonds in fatty acyl chains, the corresponding FAs become unstable at higher temperatures; therefore, the oxidative stability of oil is a critical consideration. Oilseed species typically contain five most common FAs, 16:0 (palmitic acid), 18:0 (stearic acid), 18:1 (oleic acid), and 18:3 (linolenic acid), in their seed oils; interestingly, these five FAs are also primary constituents of the bipolar lipid membrane of all plant cells (Kumar et al. 2006Kumar R, Wallis JG, Skidmore C, Browse J2006 A mutation in Arabidopsis cytochrome b5 reductase identified by high‐throughput screening differentially affects hydroxylation and desaturation. The Plant Journal 48:920-932) and hence are regarded as usual FAs. In addition to being stored in seeds as TAG (triacylglycerol) molecules, FAs and their derivatives are involved as components of the cellular membrane structure and precursors of signaling molecules that are associated with stress response and plant growth (Harwood 1996Harwood JL1996 Recent advances in the biosynthesis of plant fatty acids. Biochimica et Biophysica Acta (BBA)-Lipids and Lipid Metabolism 1301:7-56, Weber 2002Weber H2002 Fatty acid-derived signals in plants. Trends in Plant Science 7:217-224). In contrast to the usual FA-forming constituents of seed TAG in most cultivated oilseed species, the seeds of 300 plant species are reported to store FAs with hydroxy, epoxy or other functional groups or conjugated double bonds (Badami and Patil 1980Badami RC, Patil KB1980 Structure and occurrence of unusual fatty acids in minor seed oils. Progress in Lipid Research 19:119-153). Interestingly, these FAs are largely excluded from the lipid bilayer (Gurr 1980Gurr MI1980 The biosynthesis of triacylglycerols. In Stumpf PK (ed) The biochemistry of plants, lipids: Structure and function. Academic Press Inc., New York, p. 205-248), as their altered chemical structure has the potential to destabilize the fluidic nature of cellular membranes; hence, they are designated novel or unusual FAs.

The majority of polyunsaturated fatty acids (PUFAs) are synthesized through oleate desaturase, also known as fatty acid desaturase-2 (FAD2), which is a membrane-bound enzyme located in microsomes or the endoplasmic reticulum (ER). Moreover, FAD2 is responsible for more than 90% of PUFA synthesis in nonphotosynthetic tissues such as roots and in the developing seeds of oilseed crops, including Brassica (Okuley et al. 1994Okuley J, Lightner J, Feldmann K, Yadav N, Lark E, Browse J1994 Arabidopsis FAD2 gene encodes the enzyme that is essential for polyunsaturated lipid synthesis. The Plant Cell 6:147-158). Oleic acid (18:1) has one double bond, and by the catalytic action of the FAD2 enzyme, it is converted to linoleic acid (18:2) (two double bonds); by another round of desaturation, catalyzed by the FAD3 enzyme, 18:2 converts into linolenic acid (18:3) (three double bonds) (Browse et al. 1998Browse J, Spychalla J, Okuley J, Lightner J1998 Altering the fatty acid composition of vegetable oils. In Harwood JL (ed) Plant lipid biosynthesis: Fundamentals and agricultural applications. Cambridge University Press, United Kingdom, p. 131-153, Park and Kim 2022Park ME, Kim HU2022 Applications and prospects of genome editing in plant fatty acid and triacylglycerol biosynthesis. Frontiers in Plant Science 13:177-182). Interestingly, these three 18C unsaturated FAs constitute the majority of the total FAs found in the seed oil of many cultivated oilseed species. Among the three 18C unsaturated fatty acids, 18:3 is the least stable oxidatively and thus is utilized for drying oil. On the other hand, 18:1 is suitable for cooking purposes because it has greater oxidative stability at higher temperatures and is thus less prone to metabolizing into trans-fatty acids that cause various coronary diseases (Browse et al. 1998Browse J, Spychalla J, Okuley J, Lightner J1998 Altering the fatty acid composition of vegetable oils. In Harwood JL (ed) Plant lipid biosynthesis: Fundamentals and agricultural applications. Cambridge University Press, United Kingdom, p. 131-153). In addition to the abovementioned usual FAs, the seed oil of Brassica species also contains significant amounts of long-chain monounsaturated FAs, mainly erucic acid (C22:1), which are synthesized by the FAE1 enzyme (Millar and Kunst 1997Millar AA, Kunst L1997 Very‐long‐chain fatty acid biosynthesis is controlled through the expression and specificity of the condensing enzyme. The Plant Journal 12:121-131, Shi et al. 2017Shi J, Lang C, Wang F, Wu X, Liu R, Zheng T, Zhang D, Chen J, Wu G2017 Depressed expression of FAE1 and FAD2 genes modifies fatty acid profiles and storage compounds accumulation in Brassica napus seeds. Plant Science 263:177-182). Erucic acid in Brassica seed oil is a nutritionally undesirable constituent and is reported to be a health hazard. Vegetable oil from seeds of Brassica juncea (Indian mustard) contributes significantly to human consumption in the form of cooking oil in India. Indian mustard (Brassica juncea) constitutes the second largest oilseed crop after soybean, and Indian mustard is positioned as an important source of edible oil in the Indian subcontinent (Jat et al. 2019Jat RS, Singh VV, Sharma P, Rai PK2019 Oilseed brassica in India: Demand, supply, policy perspective and future potential. OCL 26:8).

The seed FA composition of Indian mustard primarily includes erucic acid (50-60%), linolenic acid (14-16%), linoleic acid (10-15%), oleic acid (10-15%) and palmitic acid (8-10%) (Suresha et al. 2012Suresha GS, Rai RD, Santha IM2012 Molecular cloning, expression analysis and growth temperature dependent regulation of a novel oleate desaturase gene (fad2) homologue from Brassica juncea. Australian Journal of Crop Science 6:296-308). In Indian mustard (Brassica juncea) varieties, the seed oleic acid content ranges from 14-15% in Pusa Agrani (Anand et al. 2010Anand A, Nagarajan S, Kishore N, Verma AP2010 Impact of high temperature at pod development stage on yield and quality of Brassica juncea cultivars under controlled conditions. Indian Journal of Agricultural Sciences 80:12) and 12-13% in Pusa Jagannath and 18-19% in Basanti (Chauhan and Kumar 2011Chauhan JS, Singh KH, Singh VV, Kumar S2011 Hundred years of rapeseed-mustard breeding in India: accomplishments and future strategies. Indian Journal of Agricultural Sciences 81:1093-109). However, Brassica juncea seed oil, which is significantly rich in oleic acid (18:1), preferably in the range of 70-80%, along with an optimal balance of linoleic (18:2) and linolenic (18:3) acids and a near-zero level of erucic acid (22:1), will be the most desirable cooking oil commodity due to its balanced nutritional composition (Dar et al. 2017Dar AA, Choudhury AR, Kancharla PK, Arumugam N2017 The FAD2 gene in plants: occurrence, regulation, and role. Frontiers in Plant Science 8:1789). Moreover, oils derived from vegetable sources rich in oleic acid (18:1) are helpful for reducing cholesterol, suppressing tumor formation, and protecting against inflammatory diseases (Yamaki et al. 2005Yamaki T, Nagamine I, Fukumoto K, Yano T, Miyahara M, Sakurai H2005 High oleic peanut oil modulates promotion stage in lung tumorigenesis of mice treated with methyl nitrosourea. Food Science and Technology Research 11:231-235). The human body is unable to synthesize the essential PUFAs linoleic (18:2) and linolenic acids (18:3) (Guan et al. 2012Guan LL, Wang YB, Shen H, Hou K, Xu YW, Wu W2012 Molecular cloning and expression analysis of genes encoding two microsomal oleate desaturases (FAD2) from safflower (Carthamus tinctorius L.). Plant Molecular Biology Reporter 30:139-148), and these two 18C PUFAs serve as precursors for long-chain PUFAs, such as EPA, DHA and ARA, which are important for maintaining a healthy body (Damude et al. 2006Damude HG, Zhang H, Farrall L, Ripp KG, Tomb JF, Hollerbach D, Yadav NS2006 Identification of bifunctional Δ12/ω3 fatty acid desaturases for improving the ratio of ω3 to ω6 fatty acids in microbes and plants. Proceedings of the National Academy of Sciences 103:9446-9451). There are several reports in which a reduction in PUFA levels in Brassica juncea was accomplished through genetic engineering techniques. For example, in the Australian Brassica juncea variety, cosuppression of fad2 resulted in an increase in oleic acid (18:1) content and simultaneous decreases in linoleic (18:2) and linolenic (18:3) acid levels (Stoutjesdijk et al. 2000Stoutjesdijk PA, Hurlestone C, Singh SP, Green AG2000 High-oleic acid Australian Brassica napus and B. juncea varieties produced by co-suppression of endogenous Δ12-desaturases. Biochemical Society Transactions 28:938-940). Efforts were also made to achieve desirable levels of 18C polyunsaturated fatty acids (FAs), namely, linoleic (18:2) and linolenic (18:3) acids, by antisensing the fad2 gene driven by a seed-specific promoter (napin) in the zero-erucic Indian mustard variety (Sivaraman et al. 2004Sivaraman I, Arumugam N, Sodhi YS, Gupta V, Mukhopadhyay A, Pradhan AK, Burma PK, Pental D2004 Development of high oleic and low linoleic acid transgenics in a zero erucic acid Brassica juncea L. (Indian mustard) line by antisense suppression of the fad2 gene. Molecular Breeding 13:365-375). Similarly, breeding successes have been achieved in developing nearly zero (less than 2%) erucic acid (EA) lines through introgression of the fae1 allele from the ZEM1 mustard variety (Australian near zero erucic acid B. juncea variety) in cultivated Indian mustard varieties (Yadava et al. 2014Yadava DK, Sujata V, Singh N, Prabhu KV, Yadav SK, Yadav MS, Dass R, Giri SC, Singh R2014 Variety of Pusa mustard 30. Indian Journal of Genetics and Plant Breeding 74:126-127).

FAD2 is a regulatory gene that controls the synthesis and thus overall composition of downstream synthesized 18C PUFAs, viz. 18:2 and 18:3. The two 18C PUFAs, 18:2 (linoleic acid) and 18:3 (linolenic acid), constitute essential components of our diet, as humans lack the enzymes required for their biosynthesis (Guan et al. 2012Guan LL, Wang YB, Shen H, Hou K, Xu YW, Wu W2012 Molecular cloning and expression analysis of genes encoding two microsomal oleate desaturases (FAD2) from safflower (Carthamus tinctorius L.). Plant Molecular Biology Reporter 30:139-148). However, the increased number of double bonds in fatty acyl chains renders them susceptible to oxidation, particularly at high temperatures, and leads to the transformation of fatty acyl chains into trans fats (Browse et al. 1998Browse J, Spychalla J, Okuley J, Lightner J1998 Altering the fatty acid composition of vegetable oils. In Harwood JL (ed) Plant lipid biosynthesis: Fundamentals and agricultural applications. Cambridge University Press, United Kingdom, p. 131-153). Therefore, increased consumption of unsaturated fatty acids with trans double bonds enhances susceptibility to cardiovascular disease and metabolic syndrome in humans (Hirata et al. 2021Hirata Y, Takahashi M, Yamada Y, Matsui R, Inoue A, Ashida R, Noguchi T, Matsuzawa A2021 trans-Fatty acids promote p53-dependent apoptosis triggered by cisplatin-induced DNA interstrand crosslinks via the Nox-RIP1-ASK1-MAPK pathway. Scientific Reports 11:10350). Natural variability of low 18:2 and 18:3 levels has not been reported in Indian mustard (Brassica juncea) (Sivaraman et al. 2004Sivaraman I, Arumugam N, Sodhi YS, Gupta V, Mukhopadhyay A, Pradhan AK, Burma PK, Pental D2004 Development of high oleic and low linoleic acid transgenics in a zero erucic acid Brassica juncea L. (Indian mustard) line by antisense suppression of the fad2 gene. Molecular Breeding 13:365-375). In the present study, FAD2 genes were cloned from three Indian mustard varieties to aid in understanding the molecular nature of their alleles.

MATERIAL AND METHODS

The genomic DNA of three released Indian mustard varieties (Basanti, Pusa Jagannath, and Pusa Agrani) was isolated from developing seedlings according to the 2xCTAB protocol (Lukowitz et al. 2000Lukowitz W, Gillmor CS, Scheible WR2000 Positional cloning in Arabidopsis. Why it feels good to have a genome initiative working for you. Plant Physiology 123:795-806). The PCR (polymerase chain reaction) conditions for amplification were as follows: initial denaturation at 94 °C for 5 min; 30 cycles at 94 °C, 56 °C and 72 °C; and a final extension at 72 °C for 10 min. The primers (F-5’ ATGGGTGCAGGTGGAAGAATG 3’; R- 5’ TCATAACTTATTGTTGTACCAG 3’) were designed based on published sequence information of the FAD2 gene (NCBI accession numbers FJ696650 and EF639848) of oilseed Brassica juncea. The reaction mixture (20 µL) for PCR included 1 U of Ex Taq DNA polymerase (Takara, USA), 1x PCR buffer (Takara, USA), 0.5(M primers (IDT), 0.2 mM dNTP mixture (Takara, USA) and 1 µL of dissolved DNA as a template. The final volume was adjusted with nuclease-free water (Invitrogen, USA). Following PCR, the amplified products were electrophoresed on a 1% native agarose gel. The fragments of the expected size (~1.1 kb) were eluted from the gel by an extraction kit (Zymo Research, USA). The gel-purified DNA fragment was finally dissolved in the requisite volume of nuclease-free water and stored in a freezer until use.

The purified DNA fragment was subsequently cloned and inserted into the PCR2.1 vector (Invitrogen, USA) according to the protocol of the Topo T/A cloning kit (Invitrogen, USA). The positive clones, obtained after verification by PCR using both gene-specific and vector-borne primers, were subsequently subjected to plasmid isolation via a kit (Zymo Research, USA). The plasmids with inserts were subjected to bidirectional sequencing. Sequence analyses were carried out using vector NTI software (Invitrogen, USA) and the BLAST service of NCBI.

RESULTS AND DISCUSSION

Fatty acid desaturase (FAD) enzymes are resident enzymes of both plastids and the endoplasmic reticulum. FADs not only regulate seed FA composition but also helps in maintaining the bipolar nature of lipidic cellular membranes (Browse et al. 1998Browse J, Spychalla J, Okuley J, Lightner J1998 Altering the fatty acid composition of vegetable oils. In Harwood JL (ed) Plant lipid biosynthesis: Fundamentals and agricultural applications. Cambridge University Press, United Kingdom, p. 131-153). In addition, the degree of FA desaturation, mediated by FADs, is also associated with diverse abiotic stress responses and stress adaptation (Upchurch 2008Upchurch RG2008 Fatty acid unsaturation, mobilization, and regulation in the response of plants to stress. Biotechnology Letters 30:967-977, Singer et al. 2016Singer SD, Zou J, Weselake RJ2016 Abiotic factors influence plant storage lipid accumulation and composition. Plant Science 243:1-9, Xue et al. 2018Xue Y, Chen B, Wang R, Win AN, Li J, Chai Y2018 Genome-wide survey and characterization of fatty acid desaturase gene family in Brassica napus and its parental species. Applied Biochemistry and Biotechnology 184:582-598). FADs are also associated with maintaining phytohormone levels in plants. Modulation of PUFA content in hypomorphic fad2 mutants affects plant growth and development through the abscisic acid (ABA), salicylic acid (SA), and jasmonic acid (JA) pathways (Kachroo et al. 2003Kachroo A, Lapchyk L, Fukushige H, Hildebrand D, Klessig D, Kachroo P2003 Plastidial fatty acid signaling modulates salicylic acid-and jasmonic acid-mediated defense pathways in the Arabidopsis ssi2 mutant. The Plant Cell 15:2952-2965, Regente et al. 2008Regente M, Monzon GC, de la Canal L2008 Phospholipids are present in extracellular fluids of imbibing sunflower seeds and are modulated by hormonal treatments. Journal of Experimental Botany 59:553-562). Linolenic acid (18:3) is the precursor of jasmonate, which was corroborated by the fact that spraying either jasmonate or 18:3 on floral buds of triple mutants ofArabidopsis (fad3-2 fad7-2 fad8) resulted in fertility (Stintzi and Browse 2000Stintzi A, Browse J2000 The Arabidopsis male-sterile mutant, opr3, lacks the 12-oxophytodienoic acid reductase required for jasmonate synthesis. Proceedings of the National Academy of Sciences 97:10625-10630). Among the FADs, FAD2 synthesizes >90% of the total PUFAs in nonphotosynthetic tissues, such as seeds and roots (Miquel and Browse 1992Miquel M, Browse J1992 Arabidopsis mutants deficient in polyunsaturated fatty acid synthesis. Biochemical and genetic characterization of a plant oleoyl-phosphatidylcholine desaturase. Journal of Biological Chemistry 267:1502-1509). In seed storage, TAG and polyunsaturated fatty acids (PUFAs) are synthesized from phosphatidylcholine (PC), a membrane lipid (He et al. 2020He M, Qin CX, Wang X, Ding NZ2020 Plant unsaturated fatty acids: Biosynthesis and regulation. Frontiers in Plant Science 11:390). The desaturation of stearic acid (C18:0) to oleic acid (C18:1) is catalyzed by soluble stearoyl-ACP desaturase (SAD) housed inside the stroma of plastids. Once synthesized in plastids, oleic acid (18:1) is then transported out for further desaturation in the ER or microsomes. 18:1 undergoes the first round of desaturation by FAD2, after which it is converted to 18:2, and by another round of desaturation by FAD3, 18:2 is converted to 18:3 (Browse et al. 1998, Park and Kim 2022Park ME, Kim HU2022 Applications and prospects of genome editing in plant fatty acid and triacylglycerol biosynthesis. Frontiers in Plant Science 13:177-182) (Figure 1). The FAD2 gene harbors a single intron positioned at the 5’-untranslated region (UTR), which is evolutionarily conserved (Dar et al. 2017Dar AA, Choudhury AR, Kancharla PK, Arumugam N2017 The FAD2 gene in plants: occurrence, regulation, and role. Frontiers in Plant Science 8:1789, Cao et al. 2021Cao S, Zhang J, Cheng H, Aslam M, Lv H, Dong W, Hu A, Guo M, Liu Q, Qin Y2021 Identification and evolutionary analysis of FAD2 gene family in green plants. Tropical Plant Biology 14:239-250). The sequence variation of the cis-acting intron of this gene has been attributed to the primary regulation of desaturation by FAD2 in Arabidopsis (Menard et al. 2017Menard GN, Moreno JM, Bryant FM, Munoz-Azcarate O, Kelly AA, Hassani-Pak K, Kurup S, Eastmond PJ2017 Genome wide analysis of fatty acid desaturation and its response to temperature. Plant Physiology 173:1594-1605).

Figure 1
Simplified scheme of seed fatty acid biosynthesis in Brassica species.

In the present study, we explored the allelic nature of the FAD2 gene, which encodes an oleate desaturase/delta-12 desaturase or ω-6 desaturase enzyme, in three Indian mustard varieties (Basanti, Pusa Jagannath, and Pusa Agrani). Brassica juncea (AABB) is an amphidiploid derived via natural interspecific hybridization between Brassica rapa (AA, n=10) and Brassica nigra (BB, n=8) (Kang et al. 2021Kang L, Qian L, Zheng M, Chen L, Chen H, Yang L, You L, Yang B, Yan M, Gu Y, Wang T2021 Genomic insights into the origin, domestication, and diversification of Brassica juncea. Nature Genetics 53:1392-1402). Two earlier studies conducted independently (Singh et al. 1995Singh S, Van der Heide T, McKinney S, Green A1995 Nucleotide sequence of a cDNA (Accession No. X91139) from Brassica juncea encoding a microsomal omega-6 desaturase. Plant Physiology (Electronic Plant Gene Register PGR 95-107) 109:1498, Suresha et al. 2012Suresha GS, Rai RD, Santha IM2012 Molecular cloning, expression analysis and growth temperature dependent regulation of a novel oleate desaturase gene (fad2) homologue from Brassica juncea. Australian Journal of Crop Science 6:296-308) reported the occurrence of the 1155 bp long FAD2 gene in Brassica juncea. However, comparisons of the FAD2 cDNA sequence of B. jun-2 (EF639848) (Suresha et al. 2012) with that of the B. jun-1 (X91139) (Singh et al. 1995) gene of B. juncea revealed 77.69% and 95.09% similarity at both the nucleotide and amino acid levels, respectively. In the present investigation, two homologs of the FAD2 gene with 1155 bp long ORFs were identified in the Indian mustard variety (Pusa Jagannath). Based on the similarities of the 1155 bp long FAD2 gene sequence identified in the present study with those of the FAD2 genes of Brassica rapa (AA) and Brassica nigra (BB), the two homologs were designated FAD2-A and FAD2-B, respectively. Taken together, the above results point toward the natural occurrence of two homologs of the FAD2 gene, which is consistent with its origin from both the diploids B. rapa and B. nigra (Figure 2). In the Indian mustard variety Pusa Jagannath, two homologs of FAD2, i.e., FAD2-A and FAD2-B, were obtained; however, in the other two varieties (Basanti and Agrani), only one homolog, i.e., FAD2-A, was isolated (Figure 2).

Figure 2
Regions of nucleic acid alignment of FAD2 genes of Brassica juncea varieties (Basanti, Pusa Jagannath, and Pusa Agrani) showing major differences. The published FAD2 gene sequences of B. juncea [B. jun-1 (X91139) and B. jun-2 (EF639848)], B. nigra (HM138369) and B. rapa (JN859550) were used as references. [The NCBI accession number of Indian mustard in this study. Basanti (FAD2-N)- MN585114, Pusa Agrani (FAD2-N)- MN585115, Pusa Jagannath (FAD2-N)- MN585116, Pusa Jagannath (FAD2-A)- MN585117, Basanti (FAD2-B)- MN585119, Pusa Agrani (FAD2-B)- MN585118 and Pusa Jagannath (FAD2-B)- MN585120].

Figure 3
Alignment of FAD2 proteins (selected regions) of Brassica juncea varieties (Basanti, Pusa Jagannath and Pusa Agrani). The published FAD2 protein sequences of B. juncea [B. jun-1 (CAA62578) and B. jun-2 (ABR27357)], B. nigra (ADJ58018) and B. rapa (AFC41105) were used as references. The NCBI accession numbers of Indian mustard in this study include the FAD2-N proteins of Basanti (QGW48095), Pusa Agrani (QGW48096), and Pusa Jagannath (QGW48097); the FAD2-A protein of Pusa Jagannath (QGW48098); and the FAD2-B proteins of Basanti (QGW48100), Pusa Agrani (QGW48099) and Pusa Jagannath (QGW48101).

Brassica napus (AACC) is also an amphidiploid species, but unlike Brassica juncea, it originates from spontaneous hybridization of Brassica rapa (AA) and Brassica oleracea (CC) (Nagaharu and Nagaharu 1935Nagaharu UJ, Nagaharu NJ1935 Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. Japanese Journal of Botany 7:389-452). It was estimated that the Brassica napus genome contains four FAD2 genes (Scheffler et al. 1997Scheffler JA, Sharpe AG, Schmidt H, Sperling P, Parkin IA, Lühs W, Lydiate DJ, Heinz E1997 Desaturase multigene families of Brassica napus arose through genome duplication. Theoretical and Applied Genetics 94:583-591). Later, four homologs of the FAD2 gene were cloned (Yang et al. 2012Yang Q, Fan C, Guo Z, Qin J, Wu J, Li Q, Fu T, Zhou Y2012 Identification of FAD2 and FAD3 genes in Brassica napus genome and development of allele-specific markers for high oleic and low linolenic acid contents. Theoretical and Applied Genetics 125:715-729) and grouped into two clusters. The two homologs of each cluster were ~97% identical at the nucleotide level, but the nucleotide identity between homologs of the two clusters was less than 90% (Lee et al. 2013Lee KR, Sohn SI, Jung JH, Kim SH, Roh KH, Kim JB, Suh MC, Kim HU2013 Functional analysis and tissue-differential expression of four FAD2 genes in amphidiploid Brassica napus derived from Brassica rapa and Brassica oleracea. Gene 531:253-262). Like in Brassica napus, we also identified two molecularly different FAD2 genes (FAD2-A and -B) in the Pusa Jagannath variety of Brassica juncea (Indian mustard) that were 92% identical at the nucleotide level. However, we could not isolate the other homologs of each of the two divergent FAD2 gene types (A/B) in the three varieties. Sequence analysis of fewer clones could explain this lack of recovery, as we were unsuccessful in isolating the FAD2-A gene type in two varieties. Alternatively, the other homolog of each FAD2 type could be ~100% identical at the nucleotide level because of the low polymorphism level among Indian mustard varieties, which has been attributed to their narrow genetic base (Chauhan et al. 2011Chauhan JS, Singh KH, Singh VV, Kumar S2011 Hundred years of rapeseed-mustard breeding in India: accomplishments and future strategies. Indian Journal of Agricultural Sciences 81:1093-109).

However, interestingly, an altogether novel type of FAD2 gene with a 1167 bp long ORF was found in all three varieties. We designated the novel FAD2 gene FAD2-N. A comparison of the FAD2-N sequence among the three varieties revealed 8 SNPs in Basanti compared to the other two varieties (Pusa Jagannath and Pusa Agrani), but only three SNPs positioned at 368, 797, and 901 changed the amino acid composition to D123G, I266T and T301A, respectively. Among the two SNPs located at 96 and 242 in Pusa Jagannath, only one exhibited change in the amino acid H81R compared to the other two varieties. In Pusa Agrani, there was a lone SNP at 381, but this did not change the amino acid composition. When the FAD2-N sequences of three varieties were used as queries to explore the genome of Varuna, an oleiferous variety of Brassica juncea (Paritosh et al. 2021Paritosh K, Yadava SK, Singh P, Bhayana L, Mukhopadhyay A, Gupta V, Bisht NC, Zhang J, Kudrna DA, Copetti D, Wing RA, Lachagiri VBR, Pradhan AK, Pental D2021 A chromosome‐scale assembly of allotetraploid Brassica juncea (AABB) elucidates comparative architecture of the A and B genomes. Plant Biotechnology Journal 19:602-614) in the NCBI database, > 99% (on Chr. B1) nucleotide identity was detected. Interestingly, after searching against the genome of Tumida (http://brassicadb.cn/), a vegetable variety of Brassica juncea (Yang et al. 2016Yang J, Liu D, Wang X, Ji C, Cheng F, Liu B, Hu Z, Chen S, Pental D, Ju Y, Yao P2016 The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection. Nature Genetics 48:1225-1232), the nucleotide identity was also found to be ~ >99% (on Chr. 07). Sequence comparison of FAD2-N with FAD2 genes of oilseed Brassica species revealed the presence of three insertions in FAD2-N differing in nucleotide numbers (boxed) (Figure 2). In addition to having longer ORFs owing to insertions, FAD2-N also presented several SNPs across its length, which distinguishes it from the FAD2-A and FAD2-B of Indian mustard varieties and other Brassica species. Additionally, several pseudogenes of FAD2 were also identified in the varieties examined; these included internal stop codons that encode nonfunctional proteins (data not shown).

The translation of the novel FAD2-N ORF yielded 388 polypeptides. It is premature to draw conclusions about the functional nature of these proteins in plants, as further study is needed. Interestingly, the novel FAD2-N proteins of all three varieties were 99.8% identical at the nucleic acid (NA) level, while their translated protein showed 98.7% identity at the amino acid level, suggesting little structural or functional differences among them. Among the FAD2-B alleles of the three varieties, 23 SNPs were found in Pusa Agrani compared to Basanti and Pusa Jagannath. On the other hand, a single SNP (C596T) was found in Basanti as compared to Pusa Jagannath. Interestingly, except for one change in amino acid (S199L) in Basanti, the translated proteins of the three varieties were nearly identical. For the FAD2-B allele (nigra type), all three varieties shared 97.7% identity at the NA level and 99.7% identity at the amino acid level. Such a high level of amino acid identity suggested little or no difference in the functionality of the proteins. However, when the FAD2-B allele (nigra type) was compared to the FAD2-A allele (rapa type), both exhibited 90.8% identity at the NA level and 95.6% identity at the amino acid level, indicating allelic divergence. Overall, Indian mustard possess three different FAD2 alleles but limited nucleotide diversity among members of the FAD2-B type (B and N), indicating narrow genetic diversity among the varieties examined. Interestingly, among oilseed Brassica species, the FAD2-A/-B ORF is 1155 bp long in Brassica napus, Brassica juncea, Brassica carinata, and Brassica rapa. However, it is 1152 bp long in Brassica nigra and 1158 bp long in Sinapis alba. Arabidopsis, a close relative of Brassica, harbors a single FAD2 gene that is 1152 bp long, similar to that of Brassica nigra. Intriguingly, the amphidiploid Brassica juncea acquired a 1155 bp long ORF of the FAD2 gene even though one of the progenitors, Brassica nigra, possessed a 1152 bp long ORF of the FAD2 gene.

Alignment of translated proteins of FAD2 genes, including the novel gene (FAD2-N), which was cloned in the present study, revealed the presence of three conserved histidine boxes (red box) similar to those present in other Brassica species (Figure 3). Additionally, the FAD2 protein has the canonical amino acid sequence YNNKL at the 3’ end, which is essential for microsome/endoplasmic reticulum (ER) retention. The conserved histidine motifs are iron binding sites in the cytoplasmic domain of the FAD2 protein and are important for the interaction between Cytochrome b5 and ER desaturases (Okuley et al. 1994Okuley J, Lightner J, Feldmann K, Yadav N, Lark E, Browse J1994 Arabidopsis FAD2 gene encodes the enzyme that is essential for polyunsaturated lipid synthesis. The Plant Cell 6:147-158, Shanklin and Cahoon 1998Shanklin J, Cahoon EB1998 Desaturation and related modifications of fatty acids. Annual Review of Plant Biology 49:611-641). Unlike chloroplasts, where ferredoxin supplies electrons to desaturases, in microsomes, cytochrome b5 is the direct source of electrons for desaturation reactions (Kumar et al. 2006Kumar R, Wallis JG, Skidmore C, Browse J2006 A mutation in Arabidopsis cytochrome b5 reductase identified by high‐throughput screening differentially affects hydroxylation and desaturation. The Plant Journal 48:920-932, Kumar et al. 2012). Interestingly, stearoyl CoA desaturase, located in plastids, is a soluble FAD that contains two conserved histidine boxes, in contrast to membrane-bound FADs (FAD2 and FAD3), which reside in the ER and contain three conserved histidine boxes (Xue et al. 2018Xue Y, Chen B, Wang R, Win AN, Li J, Chai Y2018 Genome-wide survey and characterization of fatty acid desaturase gene family in Brassica napus and its parental species. Applied Biochemistry and Biotechnology 184:582-598). Furthermore, the translated FAD2 protein of B. nigra is having one amino acid less and Sinapis alba is having one additional amino acid in contrast to 384 amino acid long FAD2 proteins of four other Brassica species viz. B. napus, B. juncea, B. rapa and B. carinata respectively. There are several reports regarding genetic diversity among genotypes/accessions/varieties of Indian mustards. Many such reports have utilized SSR markers owing to their low cost, codominant nature, wide genome coverage and ease of experimentation (Sharma et al. 2020Sharma D, Nanjundan J, Singh L, Singh SP, Parmar N, Sujith Kumar MS, Singh KH, Mishra AK, Singh R, Verma KS, Thakur AK2020 Genetic diversity in leafy mustard (Brassica juncea var. rugosa) as revealed by agro-morphological traits and SSR markers. Physiology and Molecular Biology of Plants 26:2005-2018). However, it is interesting to note that the low polymorphism level among Indian mustard varieties has been attributed to their narrow genetic base (Chauhan et al. 2011Chauhan JS, Singh KH, Singh VV, Kumar S2011 Hundred years of rapeseed-mustard breeding in India: accomplishments and future strategies. Indian Journal of Agricultural Sciences 81:1093-109). Therefore, the near absence of nucleotide diversity among members of each FAD2 gene/allele type across varieties in the present study is not surprising. Additionally, since it encodes a vital enzyme (oleate desaturase) for plant growth and development, it is plausibly protected from random changes owing to selection pressure.

REFERENCES

  • Anand A, Nagarajan S, Kishore N, Verma AP2010 Impact of high temperature at pod development stage on yield and quality of Brassica juncea cultivars under controlled conditions. Indian Journal of Agricultural Sciences 80:12
  • Badami RC, Patil KB1980 Structure and occurrence of unusual fatty acids in minor seed oils. Progress in Lipid Research 19:119-153
  • Browse J, Spychalla J, Okuley J, Lightner J1998 Altering the fatty acid composition of vegetable oils. In Harwood JL (ed) Plant lipid biosynthesis: Fundamentals and agricultural applications. Cambridge University Press, United Kingdom, p. 131-153
  • Cao S, Zhang J, Cheng H, Aslam M, Lv H, Dong W, Hu A, Guo M, Liu Q, Qin Y2021 Identification and evolutionary analysis of FAD2 gene family in green plants. Tropical Plant Biology 14:239-250
  • Chauhan JS, Kumar S2011 Assessment of oil and seed meal quality parameters of rapeseed mustard group of crops. Indian Journal of Agricultural Sciences 81:2
  • Chauhan JS, Singh KH, Singh VV, Kumar S2011 Hundred years of rapeseed-mustard breeding in India: accomplishments and future strategies. Indian Journal of Agricultural Sciences 81:1093-109
  • Damude HG, Zhang H, Farrall L, Ripp KG, Tomb JF, Hollerbach D, Yadav NS2006 Identification of bifunctional Δ12/ω3 fatty acid desaturases for improving the ratio of ω3 to ω6 fatty acids in microbes and plants. Proceedings of the National Academy of Sciences 103:9446-9451
  • Dar AA, Choudhury AR, Kancharla PK, Arumugam N2017 The FAD2 gene in plants: occurrence, regulation, and role. Frontiers in Plant Science 8:1789
  • Guan LL, Wang YB, Shen H, Hou K, Xu YW, Wu W2012 Molecular cloning and expression analysis of genes encoding two microsomal oleate desaturases (FAD2) from safflower (Carthamus tinctorius L.). Plant Molecular Biology Reporter 30:139-148
  • Gurr MI1980 The biosynthesis of triacylglycerols. In Stumpf PK (ed) The biochemistry of plants, lipids: Structure and function. Academic Press Inc., New York, p. 205-248
  • Harwood JL1996 Recent advances in the biosynthesis of plant fatty acids. Biochimica et Biophysica Acta (BBA)-Lipids and Lipid Metabolism 1301:7-56
  • He M, Qin CX, Wang X, Ding NZ2020 Plant unsaturated fatty acids: Biosynthesis and regulation. Frontiers in Plant Science 11:390
  • Hirata Y, Takahashi M, Yamada Y, Matsui R, Inoue A, Ashida R, Noguchi T, Matsuzawa A2021 trans-Fatty acids promote p53-dependent apoptosis triggered by cisplatin-induced DNA interstrand crosslinks via the Nox-RIP1-ASK1-MAPK pathway. Scientific Reports 11:10350
  • Jat RS, Singh VV, Sharma P, Rai PK2019 Oilseed brassica in India: Demand, supply, policy perspective and future potential. OCL 26:8
  • Kachroo A, Lapchyk L, Fukushige H, Hildebrand D, Klessig D, Kachroo P2003 Plastidial fatty acid signaling modulates salicylic acid-and jasmonic acid-mediated defense pathways in the Arabidopsis ssi2 mutant. The Plant Cell 15:2952-2965
  • Kang L, Qian L, Zheng M, Chen L, Chen H, Yang L, You L, Yang B, Yan M, Gu Y, Wang T2021 Genomic insights into the origin, domestication, and diversification of Brassica juncea. Nature Genetics 53:1392-1402
  • Kumar R, Tran LS, Neelakandan AK, Nguyen HT2012 Higher plant cytochrome b5 polypeptides modulate fatty acid desaturation. PLoS One 7:e31370
  • Kumar R, Wallis JG, Skidmore C, Browse J2006 A mutation in Arabidopsis cytochrome b5 reductase identified by high‐throughput screening differentially affects hydroxylation and desaturation. The Plant Journal 48:920-932
  • Lee KR, Sohn SI, Jung JH, Kim SH, Roh KH, Kim JB, Suh MC, Kim HU2013 Functional analysis and tissue-differential expression of four FAD2 genes in amphidiploid Brassica napus derived from Brassica rapa and Brassica oleracea. Gene 531:253-262
  • Lukowitz W, Gillmor CS, Scheible WR2000 Positional cloning in Arabidopsis. Why it feels good to have a genome initiative working for you. Plant Physiology 123:795-806
  • Menard GN, Moreno JM, Bryant FM, Munoz-Azcarate O, Kelly AA, Hassani-Pak K, Kurup S, Eastmond PJ2017 Genome wide analysis of fatty acid desaturation and its response to temperature. Plant Physiology 173:1594-1605
  • Millar AA, Kunst L1997 Very‐long‐chain fatty acid biosynthesis is controlled through the expression and specificity of the condensing enzyme. The Plant Journal 12:121-131
  • Miquel M, Browse J1992 Arabidopsis mutants deficient in polyunsaturated fatty acid synthesis. Biochemical and genetic characterization of a plant oleoyl-phosphatidylcholine desaturase. Journal of Biological Chemistry 267:1502-1509
  • Nagaharu UJ, Nagaharu NJ1935 Genome analysis in Brassica with special reference to the experimental formation of B. napus and peculiar mode of fertilization. Japanese Journal of Botany 7:389-452
  • Okuley J, Lightner J, Feldmann K, Yadav N, Lark E, Browse J1994 Arabidopsis FAD2 gene encodes the enzyme that is essential for polyunsaturated lipid synthesis. The Plant Cell 6:147-158
  • Paritosh K, Yadava SK, Singh P, Bhayana L, Mukhopadhyay A, Gupta V, Bisht NC, Zhang J, Kudrna DA, Copetti D, Wing RA, Lachagiri VBR, Pradhan AK, Pental D2021 A chromosome‐scale assembly of allotetraploid Brassica juncea (AABB) elucidates comparative architecture of the A and B genomes. Plant Biotechnology Journal 19:602-614
  • Park ME, Kim HU2022 Applications and prospects of genome editing in plant fatty acid and triacylglycerol biosynthesis. Frontiers in Plant Science 13:177-182
  • Regente M, Monzon GC, de la Canal L2008 Phospholipids are present in extracellular fluids of imbibing sunflower seeds and are modulated by hormonal treatments. Journal of Experimental Botany 59:553-562
  • Scheffler JA, Sharpe AG, Schmidt H, Sperling P, Parkin IA, Lühs W, Lydiate DJ, Heinz E1997 Desaturase multigene families of Brassica napus arose through genome duplication. Theoretical and Applied Genetics 94:583-591
  • Shanklin J, Cahoon EB1998 Desaturation and related modifications of fatty acids. Annual Review of Plant Biology 49:611-641
  • Sharma D, Nanjundan J, Singh L, Singh SP, Parmar N, Sujith Kumar MS, Singh KH, Mishra AK, Singh R, Verma KS, Thakur AK2020 Genetic diversity in leafy mustard (Brassica juncea var. rugosa) as revealed by agro-morphological traits and SSR markers. Physiology and Molecular Biology of Plants 26:2005-2018
  • Shi J, Lang C, Wang F, Wu X, Liu R, Zheng T, Zhang D, Chen J, Wu G2017 Depressed expression of FAE1 and FAD2 genes modifies fatty acid profiles and storage compounds accumulation in Brassica napus seeds. Plant Science 263:177-182
  • Singer SD, Zou J, Weselake RJ2016 Abiotic factors influence plant storage lipid accumulation and composition. Plant Science 243:1-9
  • Singh S, Van der Heide T, McKinney S, Green A1995 Nucleotide sequence of a cDNA (Accession No. X91139) from Brassica juncea encoding a microsomal omega-6 desaturase. Plant Physiology (Electronic Plant Gene Register PGR 95-107) 109:1498
  • Sivaraman I, Arumugam N, Sodhi YS, Gupta V, Mukhopadhyay A, Pradhan AK, Burma PK, Pental D2004 Development of high oleic and low linoleic acid transgenics in a zero erucic acid Brassica juncea L. (Indian mustard) line by antisense suppression of the fad2 gene. Molecular Breeding 13:365-375
  • Stintzi A, Browse J2000 The Arabidopsis male-sterile mutant, opr3, lacks the 12-oxophytodienoic acid reductase required for jasmonate synthesis. Proceedings of the National Academy of Sciences 97:10625-10630
  • Stoutjesdijk PA, Hurlestone C, Singh SP, Green AG2000 High-oleic acid Australian Brassica napus and B. juncea varieties produced by co-suppression of endogenous Δ12-desaturases. Biochemical Society Transactions 28:938-940
  • Suresha GS, Rai RD, Santha IM2012 Molecular cloning, expression analysis and growth temperature dependent regulation of a novel oleate desaturase gene (fad2) homologue from Brassica juncea. Australian Journal of Crop Science 6:296-308
  • Upchurch RG2008 Fatty acid unsaturation, mobilization, and regulation in the response of plants to stress. Biotechnology Letters 30:967-977
  • Weber H2002 Fatty acid-derived signals in plants. Trends in Plant Science 7:217-224
  • Xue Y, Chen B, Wang R, Win AN, Li J, Chai Y2018 Genome-wide survey and characterization of fatty acid desaturase gene family in Brassica napus and its parental species. Applied Biochemistry and Biotechnology 184:582-598
  • Yadava DK, Sujata V, Singh N, Prabhu KV, Yadav SK, Yadav MS, Dass R, Giri SC, Singh R2014 Variety of Pusa mustard 30. Indian Journal of Genetics and Plant Breeding 74:126-127
  • Yamaki T, Nagamine I, Fukumoto K, Yano T, Miyahara M, Sakurai H2005 High oleic peanut oil modulates promotion stage in lung tumorigenesis of mice treated with methyl nitrosourea. Food Science and Technology Research 11:231-235
  • Yang J, Liu D, Wang X, Ji C, Cheng F, Liu B, Hu Z, Chen S, Pental D, Ju Y, Yao P2016 The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection. Nature Genetics 48:1225-1232
  • Yang Q, Fan C, Guo Z, Qin J, Wu J, Li Q, Fu T, Zhou Y2012 Identification of FAD2 and FAD3 genes in Brassica napus genome and development of allele-specific markers for high oleic and low linolenic acid contents. Theoretical and Applied Genetics 125:715-729

Publication Dates

  • Publication in this collection
    07 June 2024
  • Date of issue
    2024

History

  • Received
    28 Nov 2023
  • Accepted
    07 Feb 2024
  • Published
    20 Feb 2024
Crop Breeding and Applied Biotechnology Universidade Federal de Viçosa, Departamento de Fitotecnia, 36570-000 Viçosa - Minas Gerais/Brasil, Tel.: (55 31)3899-2611, Fax: (55 31)3899-2611 - Viçosa - MG - Brazil
E-mail: cbab@ufv.br