Identification of FAD2 and FAD3 genes in Brassica napus genome and development of allele-specific markers for high oleic and low linolenic acid contents

Modification of oleic acid (C18:1) and linolenic acid (C18:3) contents in seeds is one of the major goals for quality breeding after removal of erucic acid in oilseed rape (Brassica napus). The fatty acid desaturase genes FAD2 and FAD3 have been shown as the major genes for the control of C18:1 and C18:3 contents. However, the genome structure and locus distributions of the two gene families in amphidiploid B. napus are still not completely understood to date. In the present study, all copies of FAD2 and FAD3 genes in the A- and C-genome of B. napus and its two diploid progenitor species, Brassica rapa and Brassica oleracea, were identified through bioinformatic analysis and extensive molecular cloning. Two FAD2 genes exist in B. rapa and B. oleracea, and four copies of FAD2 genes exist in B. napus. Three and six copies of FAD3 genes were identified in diploid species and amphidiploid species, respectively. The genetic control of high C18:1 and low C18:3 contents in a double haploid population was investigated through mapping of the quantitative trait loci (QTL) for the traits and the molecular cloning of the underlying genes. One major QTL of BnaA.FAD2.a located on A5 chromosome was responsible for the high C18:1 content. A deleted mutation in the BnaA.FAD2.a locus was uncovered, which represented a previously unidentified allele for the high oleic variation in B. napus species. Two major QTLs on A4 and C4 chromosomes were found to be responsible for the low C18:3 content in the DH population as well as in SW Hickory. Furthermore, several single base pair changes in BnaA.FAD3.b and BnaC.FAD3.b were identified to cause the phenotype of low C18:3 content. Based on the results of genetic mapping and identified sequences, allele-specific markers were developed for FAD2 and FAD3 genes. Particularly, single-nucleotide amplified polymorphisms markers for FAD3 alleles were demonstrated to be a reliable type of SNP markers for unambiguous identification of genotypes with different content of C18:3 in amphidiploid B. napus.     

Characterization of the fan1 locus in soybean line A5 and development of molecular assays for high-throughput genotyping of FAD3 genes

Soybean is one of the most important oil crops worldwide, and reducing the linolenic acid content of soybean oil will provide increased stability of the oil to consumers and limit the amount of trans fat in processed foods. The linolenic content in soybean seed is controlled by three fatty acid desaturase (FAD) three enzymes, FAD3A, B, and C. The soybean lines with 1 % linolenic acid content which are widely used in breeding for reduced linolenic acid in the USA have mutations in each of the three FAD genes derived from lines A5 (deletion of FAD3A), A26, and A23 (missense mutations in FAD3B and C, respectively). Although soybean line A5 has been released for 30 years, the extent and definition of the deletion of the FAD3A gene has not been characterized, which has prevented researchers from designing robust molecular markers for effective marker-assisted selection (MAS). Using a PCR-based genomic strategy, we have identified a 6.4-kbp deletion of the FAD3A gene in A5 and developed a TaqMan detection assay by targeting the deletion junction in A5, which could be used to distinguish the homozygotes and heterozygotes of the gene. In addition, based on mutant single nucleotide polymorphisms in FAD3B and FAD3C identified in A26 and A23, respectively, we have also developed TaqMan assays for high-throughput MAS. The TaqMan assays have proven to be a very effective platform for detecting the mutant FAD3 alleles and thus will greatly facilitate high-throughput MAS for development of soybean lines with reduced linolenic acid content.     

A novel <Emphasis Type="Italic">FAD2</Emphasis>-<Emphasis Type="Italic">1 A</Emphasis> allele in a soybean plant introduction offers an alternate means to produce soybean seed oil with 85% oleic acid content

The alteration of fatty acid profiles in soybean to improve soybean oil quality has been a long-time goal of soybean researchers. Soybean oil with elevated oleic acid is desirable because this monounsaturated fatty acid improves the nutrition and oxidative stability of soybean oil compared to other oils. In the lipid biosynthetic pathway, the enzyme fatty acid desaturase 2 (FAD2) is responsible for the conversion of oleic acid precursors to linoleic acid precursors in developing soybean seeds. Two genes encoding FAD2-1A and FAD2-1B were identified to be expressed specifically in seeds during embryogenesis and have been considered to hold an important role in controlling the seed oleic acid content. A total of 22 soybean plant introduction (PI) lines identified to have an elevated oleic acid content were characterized for sequence mutations in the FAD 2-1A and FAD2-1B genes. PI 603452 was found to contain a deletion of a nucleotide in the second exon of FAD2-1A. These important SNPs were used in developing molecular marker genotyping assays. The assays appear to be a reliable and accurate tool to identify the FAD 2-1A and FAD2-1B genotype of wild-type and mutant plants. PI 603452 was subsequently crossed with PI 283327, a soybean line that has a mutation in FAD2-1B. Interestingly, soybean lines carrying both homozygous insertion/deletion mutation (indel) FAD2-1A alleles and mutant FAD2-1B alleles have an average of 82–86% oleic acid content, compared to 20% in conventional soybean, and low levels of linoleic and linolenic acids. The newly identified indel mutation in the FAD2-1A gene offers a simple method for the development of high oleic acid commercial soybean varieties.     

Two high linolenic mutants of Arabidopsis thaliana contain megabase-scale genome duplications encompassing the FAD3 locus

Understanding the quantitative control of fatty acid desaturation during the biosynthesis of seed storage oil has become a priority area for research, as a consequence of its importance for both human health and the substitution of mineral oil for industrial applications. We have analysed the genome structure of two mutants in Arabidopsis thaliana that show substantially elevated content of the omega‐3 polyunsaturated fatty acid linolenic acid in their seed oil. In one, rfc4, sequences totalling approximately 2 Mb from chromosome 2 have been duplicated and inserted into chromosome 3. In the other mutant, ife, chromosome 2 sequences totalling approximately 1.4 Mb have been duplicated and inserted into a linked position. In both cases, the duplications encompass the FAD3 locus, which encodes the linoleate desaturase responsible for the biosynthesis of linolenic acid for accumulation in seed storage oil. The results show that mutagens such as fast neutrons (used for the induction of rfc4) and T‐DNA (used for the induction of ife, which is not linked to the T‐DNA present in the line) can result in the duplication of very large genome segments. They also show that increasing the dosage of the FAD3‐containing genomic region results in an increase in the linolenic acid content of seed oil. Consequently, screening methods for duplication of FAD3 orthologues in oil crops may be an appropriate approach for the identification of germplasm for breeding varieties with increased proportions of linolenic acid in the oil that they produce.     

A Simple Allele-Specific PCR Assay for Detecting <Emphasis Type="Italic">FAD2</Emphasis> Alleles in Both A and B Genomes of the Cultivated Peanut for High-Oleate Trait Selection

In cultivated tetraploid peanut (2n = 4x = 40, AABB), the conversion of oleic acid to linoleic acid is mainly catalyzed by the Δ¹² fatty acid desaturase (FAD). Two homoeologous genes (FAD2A and FAD2B) encoding for the desaturase are located on the A and B genomes, respectively. Abolishing or reducing the desaturase activity by gene mutation can significantly increase the oleic acid/linoleic acid ratio. F435-derived high-oleate peanut cultivars contain two key mutations within the Δ¹² fatty acid desaturase gene which include a 1-bp substitution of G:C→A:T in the A genome and a 1-bp insertion of A:T in the B genome. Both of these mutations contribute to abolishing or reducing the desaturase activity, leading to accumulation of oleate versus linoleate. Currently, detection of FAD2 alleles can be achieved by a cleaved amplified polymorphic sequence marker for the A genome and a real-time polymerase chain reaction (PCR) marker for the B genome; however, detection of these key mutations has to use different assay platforms. Therefore, a simple PCR assay for detection of FAD2 alleles on both genomes was developed by designing allele-specific primers and altering PCR annealing temperatures. This assay was successfully used for detecting FAD2 alleles in peanut. Gas chromatography (GC) was used to determine fatty acid composition of PCR-assayed genotypes. The results from the PCR assay and GC analysis were consistent. This PCR assay is quick, reliable, economical, and easy to use. Implementation of this PCR assay will greatly enhance the efficiency of germplasm characterization and marker-assisted selection of high oleate in peanut.     

Eliminating expression of erucic acid-encoding loci allows the identification of “hidden” QTL contributing to oil quality fractions and oil content in <Emphasis Type="Italic">Brassica juncea</Emphasis> (Indian mustard)

Oil content and oil quality fractions (viz., oleic, linoleic and linolenic acid) are strongly influenced by the erucic acid pathway in oilseed Brassicas. Low levels of erucic acid in seed oil increases oleic acid content to nutritionally desirable levels, but also increases the linoleic and linolenic acid fractions and reduces oil content in Indian mustard (Brassica juncea). Analysis of phenotypic variability for oil quality fractions among a high-erucic Indian variety (Varuna), a low-erucic east-European variety (Heera) and a zero-erucic Indian variety (ZE-Varuna) developed by backcross breeding in this study indicated that lower levels of linoleic and linolenic acid in Varuna are due to substrate limitation caused by an active erucic acid pathway and not due to weaker alleles or enzyme limitation. To identify compensatory loci that could be used to increase oil content and maintain desirable levels of oil quality fractions under zero-erucic conditions, we performed Quantitative Trait Loci (QTL) mapping for the above traits on two independent F1 doubled haploid (F1DH) mapping populations developed from a cross between Varuna and Heera. One of the populations comprised plants segregating for erucic acid content (SE) and was used earlier for construction of a linkage map and QTL mapping of several yield-influencing traits in B. juncea. The second population consisted of zero-erucic acid individuals (ZE) for which, an Amplified Fragment Length Polymorphism (AFLP)-based framework linkage map was constructed in the present study. By QTL mapping for oil quality fractions and oil content in the ZE population, we detected novel loci contributing to the above traits. These loci did not co-localize with mapped locations of the fatty acid desaturase 2 (FAD2), fatty acid desaturase 3 (FAD3) or fatty acid elongase (FAE) genes unlike those of the SE population wherein major QTL were found to coincide with mapped locations of the FAE genes. Some of the new loci identified in the ZE population could be detected as ‘weak’ contributors (with LOD < 2.5) in the SE population in which their contribution to the traits was “masked” due to pleiotropic effects of erucic acid genes. The novel loci identified in this study could now be used to improve oil quality parameters and oil content in B. juncea under zero-erucic conditions.     

Mapping of the loci controlling oleic and linolenic acid contents and development of fad2 and fad3 allele-specific markers in canola (Brassica napus L.)

The quality of canola oil is determined by its constituent fatty acids such as oleic acid (C18:1), linoleic acid (C18:2) and linolenic acid (C18:3). Most canola cultivars normally produce oil with about 55–65% oleic acid and 8–12% linolenic acid. High concentrations of linolenic acid lead to oil instability and off-type flavor, while high levels of oleic acid increase oxidative stability and nutritional value of oil. Therefore, development of canola cultivars with increased oleic acid and reduced linolenic acid is highly desirable for canola oil quality. In this study, we have mapped one locus that has a major effect and one locus that has a minor effect for high oleic acid and two loci that have major effects for low linolenic acid in a doubled haploid population. The major locus for high C18:1 was proven to be the fatty acid desaturase-2 (fad2) gene and it is located on the linkage group N5; the minor locus is located on N1. One major QTL for C18:3 is the fatty acid desaturase-3 gene of the genome C (fad3c) and it is located on N14. The second major QTL resides on N4 and is the fad3a gene of the A genome. We have sequenced genomic clones of the fad2 and fad3c genes amplified from an EMS-induced mutant and a wild-type canola cultivar. A comparison of the mutant and wild-type allele sequences of the fad2 and fad3c genes revealed single nucleotide mutations in each of the genes. Detailed sequence analyses suggested mechanisms by which both the mutations can cause altered fatty acid content. Based on the sequence differences between the mutant and wild-type alleles, two single nucleotide polymorphism (SNP) markers, corresponding to the fad2 and fad3c gene mutations, were developed. These markers will be highly useful for direct selection of desirable fad2 and fad3c alleles during marker-assisted trait introgression and breeding of canola with high oleic and low linolenic acid.     

Combinations of mutant FAD2 and FAD3 genes to produce high oleic acid and low linolenic acid soybean oil

High oleic acid soybeans were produced by combining mutant FAD2-1A and FAD2-1B genes. Despite having a high oleic acid content, the linolenic acid content of these soybeans was in the range of 4-6 %, which may be high enough to cause oxidative instability of the oil. Therefore, a study was conducted to incorporate one or two mutant FAD3 genes into the high oleic acid background to further reduce the linolenic acid content. As a result, soybean lines with high oleic acid and low linolenic acid (HOLL) content were produced using different sources of mutant FAD2-1A genes. While oleic acid content of these HOLL lines was stable across two testing environments, the reduction of linolenic acid content varied depending on the number of mutant FAD3 genes combined with mutant FAD2-1 genes, on the severity of mutation in the FAD2-1A gene, and on the testing environment. Combination of two mutant FAD2-1 genes and one mutant FAD3 gene resulted in less than 2 % linolenic acid content in Portageville, Missouri (MO) while four mutant genes were needed to achieve the same linolenic acid in Columbia, MO. This study generated non-transgenic soybeans with the highest oleic acid content and lowest linolenic acid content reported to date, offering a unique alternative to produce a fatty acid profile similar to olive oil.     

Specific molecular marker of the genes controlling linolenic acid content in rapeseed

In rapeseed, which is an agronomically important oilseed, variation in the linolenic acid content of the oil has been obtained through chemical mutagenesis treatment. Conventional breeding of this quantitative trait, however requires specific molecular markers. By means of biochemical experiments, we have established that the induced variation in linolenic acid content is associated with the fad3 gene encoding the microsomal ¹⁵ desaturase. Using a pair of primers specific to this gene and a doubled haploid progeny derived from a low linolenic x high linolenic acid F₁hybrid, we have identified a polymorphism of the fad3 alleles between the low- and the high-linolenic acid genotypes. The structure exon/intron of the fad3 DNA sequence seems to be very similar to that of the Arabidopsis fad3 gene. The choice of the primer pair allows specific amplification of one of the two rapeseed fad3 genes. The value and contribution of specific markers to conventional plant breeding is discussed.     

DNA allelic variations at the loci putatively implicated in seed oil formation among Brassica oilseed cultivars

Knowledge about the genes implicated in lipid biosynthesis acquired from the model plant Arabidopsis is useful in understanding the formation of seed oil in Brassica oilseeds. In this paper, we report the screening of polymorphic markers at the loci putative for the seed oil formation between two geographically different genotypes: the Chinese cultivar Ningyou-7 and the European cultivar Tapidor. These primer pairs (150) were designed based on 75 Brassica genes that were Arabidopsis orthologues implicated in the oil formation. A total of 52 out of the 150 primer pairs associated with 47 of the 75 genes showed polymorphisms between the two genotypes. The type of polymorphisms that could be detected on capillary electrophoresis images and their respective visual futures are described. Further, we selected 34 polymorphic markers to scan allelic variations and found rich DNA polymorphisms among the 54 Brassica oilseed cultivars. On the average, each primer pair resulted in 5.6 alleles at the region that was covered. The correlation between the alleles and seed quality traits revealed that the alleles of BnFAD7 were related to the variation of linolenic acid (C_18:3) contents among the cultivars. The allele FAD7-ics11170 (3/4)-b that was significantly correlated with high linolenic acid content can be used as an efficient marker for the selection of breeding materials with high linolenic acid content.     

Genetic variation of six desaturase genes in flax and their impact on fatty acid composition

Flax (Linum usitatissimum L.) is one of the richest plant sources of omega-3 fatty acids praised for their health benefits. In this study, the extent of the genetic variability of genes encoding stearoyl-ACP desaturase (SAD), and fatty acid desaturase 2 (FAD2) and 3 (FAD3) was determined by sequencing the six paralogous genes from 120 flax accessions representing a broad range of germplasm including some EMS mutant lines. A total of 6 alleles for sad1 and sad2, 21 for fad2a, 5 for fad2b, 15 for fad3a and 18 for fad3b were identified. Deduced amino acid sequences of the alleles predicted 4, 2, 3, 4, 6 and 7 isoforms, respectively. Allele frequencies varied greatly across genes. Fad3a, with 110 SNPs and 19 indels, and fad3b, with 50 SNPs and 5 indels, showed the highest levels of genetic variations. While most of the SNPs and all the indels were silent mutations, both genes carried nonsense SNP mutations resulting in premature stop codons, a feature not observed in sad and fad2 genes. Some alleles and isoforms discovered in induced mutant lines were absent in the natural germplasm. Correlation of these genotypic data with fatty acid composition data of 120 flax accessions phenotyped in six field experiments revealed statistically significant effects of some of the SAD and FAD isoforms on fatty acid composition, oil content and iodine value. The novel allelic variants and isoforms identified for the six desaturases will be a resource for the development of oilseed flax with unique and useful fatty acid profiles.     

Functional complementation of a periila ω3 fatty acid desaturase under the seed-specific SeFAD2 promoter

Functional characterization of the fatty acid desaturase genes and seed-specific promoters is prerequisite for altering the unsaturated fatty acid content of oilseeds by genetic manipulation. The ω-6 fatty acid desaturase (FAD2) and ω-3 fatty acid desaturase (FAD3) catalyze extra-plastidial desaturation of oleic acid to linoleic acid and linoleic acid to linolenic acid, respectively. These are major constituents in seed storage oils. Here, we report the complementation of a perilla linoleic acid desaturase (PrFAD3) cDNA under the seed-specific sesame FAD2 (SeFAD2) promoter in the Arabidopsis fad3 mutant. PrFAD3 is functionally active and the SeFAD2 promoter is applicable for modifying fatty acid composition in developing seeds. Transient expression of the GUS gene under that promoter in the developing seeds and leaves of sesame, soybean, and corn via microprojectile bombardment indicated that the SeFAD2 promoter likely will be useful for altering the seed phenotypes of dicot and monocot crops.     

Mutant alleles of <Emphasis Type="Italic">FAD2-1A</Emphasis> and <Emphasis Type="Italic">FAD2-1B</Emphasis>combine to produce soybeans with the high oleic acid seed oil trait

Background  

The alteration of fatty acid profiles in soybean [Glycine max (L.) Merr.] to improve soybean oil quality is an important and evolving theme in soybean research to meet nutritional needs and industrial criteria in the modern market. Soybean oil with elevated oleic acid is desirable because this monounsaturated fatty acid improves the nutrition and oxidative stability of the oil. Commodity soybean oil typically contains 20% oleic acid and the target for high oleic acid soybean oil is approximately 80% of the oil; previous conventional plant breeding research to raise the oleic acid level to just 50-60% of the oil was hindered by the genetic complexity and environmental instability of the trait. The objective of this work was to create the high oleic acid trait in soybeans by identifying and combining mutations in two delta-twelve fatty acid desaturase genes, FAD2-1A and FAD2-1B.     

Genome-Wide Analysis of Fatty Acid Desaturases in Soybean (<Emphasis Type="Italic">Glycine max</Emphasis>)

Fatty acid desaturases can introduce double bonds into the hydrocarbon chains of fatty acids to produce unsaturated fatty acids. In the present study, 29 full-length desaturase genes were identified from soybean genome by a thorough annotation exercise. A comprehensive analysis was performed to characterize phylogeny, chromosomal locations, structures, conserved motifs, and expression patterns of those genes. The soybean genes were phylogenetically clustered into nine subfamilies with the Arabidopsis counterparts, FAB2, FAD2, FAD3, FAD5, FAD6, FAD7, FAD8, SLD1, and DES1. Twenty-nine desaturase genes were found to be distributed on at least 15 of the 20 soybean chromosomes. The gene structures and motif compositions were considerably conserved among the subfamilies. The majority of desaturase genes showed specific temporal and spatial expression patterns across different tissues and developmental stages based on microarray data analyses. The study may provide new insights into the origin and evolution of fatty acid biosynthesis pathways in higher plants. Additionally, the characterization of desaturases from soybean will lead to the identification of additional genes for genetic modification of plants to produce nutritionally important fatty acids.     

Microarray analysis of differentially expressed genes between <Emphasis Type="Italic">Brassica napus</Emphasis> strains with high- and low-oleic acid contents

An increase in oleic acid (C18:1) content is a desirable trait. Despite the critical roles of the two desaturases, FAD2 and FAD3, in the control of fatty acid desaturation, a dispute remains over whether inactivation of their genes alone is sufficient enough to generate the high-oleic trait. To address this question, we employed microarray technology to investigate the difference in gene expression profile between two different Brassica napus strains with high-C18:1 (71.71%) and low-C18:1 (55.6%) contents, respectively. Our study revealed 562 differentially expressed genes, of which 194 genes were up-regulated and 368 down-regulated. Based on the Gene Ontology classification, these genes were classified into 23 functional categories. Three of the up-regulated genes represent B. napus homologs of Arabidopsis genes encoding a cytosolic isoform of pyruvate kinase (AT3G55810), Δ9 acyl-lipid desaturase (AT1G06080, ADS1) and fatty acyl-ACP thioesterase B (AT1G08510), respectively. Conversely, the homologs of two Arabidopsis sequences encoding Δ9 acyl-lipid desaturase (AT2G31360, ADS2) and FAD3 desaturase (AT2G29980) were down-regulated in the high-oleic acid strain. Furthermore, 60 differentially expressed genes were classified as associated with relevant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Collectively, our results suggest that expressing the high-oleic acid trait may require a coordinated regulation of diverse regulatory and metabolic gene networks in addition to inactivation of the FAD2 and FAD3 genes in the oilseed. A set of the differentially expressed genes identified in this study will facilitate our efforts to tap the germplasms with the potential to express the high-oleic acid trait.     

Significant enhancement of fatty acid composition in seeds of the allohexaploid,Camelina sativa,using CRISPR/Cas9 gene editing

The CRISPR/Cas9 nuclease system is a powerful and flexible tool for genome editing, and novel applications of this system are being developed rapidly. Here, we used CRISPR/Cas9 to target the FAD2 gene in Arabidopsis thaliana and in the closely related emerging oil seed plant, Camelina sativa, with the goal of improving seed oil composition. We successfully obtained Camelina seeds in which oleic acid content was increased from 16% to over 50% of the fatty acid composition. These increases were associated with significant decreases in the less desirable polyunsaturated fatty acids, linoleic acid (i.e. a decrease from ~16% to <4%) and linolenic acid (a decrease from ~35% to <10%). These changes result in oils that are superior on multiple levels: they are healthier, more oxidatively stable and better suited for production of certain commercial chemicals, including biofuels. As expected, A. thaliana T₂ and T₃ generation seeds exhibiting these types of altered fatty acid profiles were homozygous for disrupted FAD2 alleles. In the allohexaploid, Camelina, guide RNAs were designed that simultaneously targeted all three homoeologous FAD2 genes. This strategy that significantly enhanced oil composition in T₃ and T₄ generation Camelina seeds was associated with a combination of germ‐line mutations and somatic cell mutations in FAD2 genes in each of the three Camelina subgenomes.     

Selective gene dosage by CRISPR‐Cas9 genome editing in hexaploid Camelina sativa

In many plant species, gene dosage is an important cause of phenotype variation. Engineering gene dosage, particularly in polyploid genomes, would provide an efficient tool for plant breeding. The hexaploid oilseed crop Camelina sativa, which has three closely related expressed subgenomes, is an ideal species for investigation of the possibility of creating a large collection of combinatorial mutants. Selective, targeted mutagenesis of the three delta‐12‐desaturase (FAD2) genes was achieved by CRISPR‐Cas9 gene editing, leading to reduced levels of polyunsaturated fatty acids and increased accumulation of oleic acid in the oil. Analysis of mutations over four generations demonstrated the presence of a large variety of heritable mutations in the three isologous CsFAD2 genes. The different combinations of single, double and triple mutants in the T3 generation were isolated, and the complete loss‐of‐function mutants revealed the importance of delta‐12‐desaturation for Camelina development. Combinatorial association of different alleles for the three FAD2 loci provided a large diversity of Camelina lines with various lipid profiles, ranging from 10% to 62% oleic acid accumulation in the oil. The different allelic combinations allowed an unbiased analysis of gene dosage and function in this hexaploid species, but also provided a unique source of genetic variability for plant breeding.     

Cloning and expression of <Emphasis Type="Italic">Perilla frutescens FAD2</Emphasis> gene and polymorphism analysis among cultivars

?12 fatty acid desaturase (FAD2) is a key enzyme for linoleic acid and linolenic acid biosynthesis. Perilla frutescens is a special oil plant species with highest linolenic acid content. In this study, based on RACE, two alleles for one FAD2 gene were isolated from P. frutescens cultivar C2: the 3956 bp PfFAD2a and the 3959 bp PfFAD2b, both with a full-length cDNA of 1526 bp, and both encoding a 382aa basic protein. The alleles have identities of over 98%, and their encoded proteins differ only by substitution of a strongly similar residue. Saccharomyces cerevisiae heterologous expression suggested that PfFAD2a/b both encode a bio-functional FAD2 enzyme. Phylogenetic analyses indicated that PfFAD2 shows the highest homologies to FAD2 genes from dicots such as Boraginaceae and Burseraceae. PfFAD2a/b expressions are mainly restricted to developing seeds. PfFAD2a/b expression in the seedling leaf is upregulated by cold (4 °C) and repressed by heat (42 °C). Each of the eight cultivars contains two alleles for one PfFAD2 and 40 SNP sites are found. One allelic gene in cultivars C1 and P1 is pseudogene because of premature stop codon mutation in 5′ coding region. All other normal PfFAD2 genes/allelic genes encode identical or very similar proteins. PfFAD2a/b expression level in developing seeds also varies among the eight cultivars. This study provides systemic molecular and functional features of PfFAD2 and enables its application in the study of plant fatty acids traits.