Original Research ARTICLE
QTL Landscape for Oil Content in Brassica juncea: Analysis in Multiple Bi-Parental Populations in High and “0” Erucic Background
- 1Department of Genetics, University of Delhi South Campus, New Delhi, India
- 2Centre for Genetic Manipulation of Crop Plants, University of Delhi South Campus, New Delhi, India
Increasing oil content in oilseed mustard (Brassica juncea) is a major breeding objective—more so, in the lines that have “0” erucic acid content (< 2% of the seed oil) as earlier studies have shown negative pleiotropic effect of erucic acid loci on the oil content, both in oilseed mustard and rapeseed. We report here QTL analysis of oil content in eight different mapping populations involving seven different parents—including a high oil content line J8 (~49%). The parental lines of the mapping populations contained wide variation in oil content and erucic acid content. The eight mapping populations were categorized into two sets—five populations with individuals segregating for erucic acid (SE populations) and the remaining three with zero erucic acid segregants (ZE populations). Meta-analysis of QTL mapped in individual SE populations identified nine significant C-QTL, with two of these merging most of the major oil QTL that colocalized with the erucic acid loci on the linkage groups A08 and B07. QTL analysis of oil content in ZE populations revealed a change in the landscape of the oil QTL compared to the SE populations, in terms of altered allelic effects and phenotypic variance explained by ZE QTL at the “common” QTL and observation of “novel” QTL in the ZE background. The important loci contributing to oil content variation, identified in the present study could be used in the breeding programmes for increasing the oil content in high erucic and “0” erucic backgrounds.
Oilseed mustard, Brassica juncea (AABB), is an important oilseed crop of the south Asian region. In India alone, the crop is grown on around six million hectares of land (Chauhan et al., 2011) under low moisture regimes during the winter growing season. Breeding for high oil content and improved oil quality are two key breeding objectives in oilseed mustard. The quality of oil in Brassica species, including mustard, is determined by the relative proportions of different fatty acid fractions—saturated fatty acids such as palmitic (C16:0), stearic (C18:0) and unsaturated fractions such as oleic (C18:1), eicosenoic (C20:1), erucic (C22:1), linoleic (C18:2), and linolenic acid (C18:3). Erucic acid (C22:1), a long-chain monounsaturated fatty acid is a major fraction (~50%) of the seed oil in Brassica species. Presence of high erucic acid content has been shown to be anti-nutritional (Ackman et al., 1977). Breeding for “0” erucic acid content (< 2% of the seed oil) has been a major breeding objective in mustard. Almost all the pureline varieties and hybrids of rapeseed (Brassica napus, AACC) grown in north America, Europe, China, and many other parts of the world are “0” erucic (ZE). However, development of productive ZE mustard is still an unachieved objective even though experimental lines with “0” erucic content have been available for quite some time (Kirk and Oram, 1981). One of the major reasons for the low acceptability of ZE mustard is its lower seed oil content as compared to high erucic (HE) lines.
We had earlier mapped the fatty acid elongase (FAE1) genes involved in erucic acid biosynthesis in B. juncea to LGs A08 (FAE1.1) and B07 (FAE1.2) by using a F1DH mapping population derived from a cross of an Indian gene pool line Varuna (HE) with an exotic line Heera (ZE) (Gupta et al., 2004). Transfer of these two ZE loci to the extensively grown mega-variety Varuna through marker-assisted backcross breeding led to the development of ZE-Varuna; however, the oil content in ZE-Varuna was reduced to ~36% as compared to ~43% in HE-Varuna (Jagannath et al., 2011). Many other studies have also reported such a kind of negative association between erucic acid content and overall seed oil content, suggesting a negative pleiotropic effect of ZE trait on the oil content both in rapeseed (Ecke et al., 1995; Burns et al., 2003) and mustard (Mahmood et al., 2006; Jagannath et al., 2011). Breeding for high oil content concomitant with ZE trait would require identification of compensatory loci that can contribute to high oil content in ZE background (Pradhan and Pental, 2011).
A number of studies have mapped QTL for seed oil content in mustard (Cheung et al., 1998; Mahmood et al., 2006; Ramchiary et al., 2007; Yadava et al., 2012) and rapeseed (Ecke et al., 1995; Delourme et al., 2006; Qiu et al., 2006; Cao et al., 2010; Sun et al., 2012; Zhao et al., 2012; Wang et al., 2013; Jiang et al., 2014; Fu et al., 2017). Most of the studies have used a single bi-parental population with small parental difference for the oil content. Only two studies in B. napus, Sun et al. (2012) and Wang et al. (2013), have used parents with more than 10% difference in oil content. Recently, Fu et al. (2017) undertook QTL analysis of oil content by association mapping using 142 rapeseed breeding lines with diverse oil content ranging from 37.39 to 55.28%. However, none of these studies have identified and highlighted QTL that would be useful in breeding for oil content in ZE background. We have earlier mapped QTL for oil content in two bi-parental populations derived from the F1 of a cross between Varuna and Heera, one subset segregating for erucic acid loci (SE population) and the other subset consisting only of ZE individuals (Jagannath et al., 2011). However, for broader MAS (marker-assisted selection) application, consensus QTL identified through consolidation of QTL data from multiple bi-parental mapping populations by a joint and integrative analysis would be more pertinent than the QTL identified in a single bi-parental population. Moreover, QTL data derived from multiple populations with wider genetic variability for the trait is more informative than using a single bi-parental population with large parental contrast as the former would allow identification of more number of alleles with their effects and assist breeders in marker assisted transfer of the trait more effectively.
Therefore, a major objective of the present study was to define the QTL involved in seed oil content variation in B. juncea using multiple bi-parental mapping populations. Another objective was to identify specific loci operative in a ZE acid background. The parents used in the study encapsulate maximum genetic variability for the seed oil content in B. juncea. Through the use of two sets of bi-parental mapping populations—one set segregating for erucic acid (named as SE populations) and the other set consisting of zero erucic acid segregants (named as ZE populations), the QTL analysis identified pleiotropic QTL for seed oil content in B. juncea. The study also revealed a change in the genetic landscape of QTL and identified “novel” QTL in ZE populations which otherwise remained masked due to the influence of erucic acid in the SE populations. Based on a detailed genetic analysis of seed oil content, a relevant breeding strategy for the enhancement of oil content in ZE mustard has been discussed.
Materials and Methods
Five bi-parental F1DH mapping populations derived from seven purelines of B. juncea were used in the study. The seven parental lines contained wide variation for both seed oil content and erucic acid, including a pureline J8 with an oil content of ~49%—highest reported so far in oilseed mustard. The details of these parental lines are given in Table 1. These mapping populations were developed by microspore culture following an established protocol (Mukhopadhyay et al., 2007). Eight different mapping populations were sourced from these five bi-parental source DH populations (Table 2). Of these eight populations, five populations consisted of individuals segregating for the erucic acid trait (called SE populations) and were named as EJ8A8B7, DEB7, EPJA8B7, TDA8, and VHA8B7, the remaining three mapping populations consisted of only “0” erucic acid segregants (called ZE populations) and were named as EJ8Z, DEZ, and VHZ (Table 2). As all the five F1DH source populations segregate for one or both the erucic acid genes, five SE populations were constructed (EJ8A8B7, DEB7, EPJA8B7, TDA8, and VHA8B7) from their respective F1DH source population by random sampling. On the other hand, for the construction of three ZE populations (EJ8Z, DEZ, and VHZ), the selfed seeds of respective F1DH source population were initially analyzed by GC to identify the zero erucic acid segregants and the mapping populations were constructed by random sampling from these zero erucic acid segregants. Detailed genesis of these populations is shown in Supplementary Figure 1. Genotype of parental lines for two erucic acid genes in LG A08 and B07 was ascertained by FAE1 gene markers following the protocol of Gupta et al. (2004).
Table 1. Phenotypic variation for oil content and erucic acid in seven parental lines used for developing the five bi-parental mapping populations.
Table 2. Phenotypic performance of the parents and eight bi-parental mapping populations of B. juncea for seed oil content at different locations.
Field Experiments and Trait Measurements
The oil content and erucic acid content data for EJ8A8B7/Z, DEB7/Z, and EPJA8B7 populations were estimated from the seeds by growing all the bi-parental mapping populations either in three different locations or a single location over three years during the mustard growing season (Table 2). A randomized block design with three replications for each entry was adopted for field experiments of all the mapping populations. Each one of the mapping population and their respective parental lines were planted in three rows with a row length of 3 m per replication. On maturity, open pollinated seeds of five competitive plants from the middle row of each replication were pooled, desiccated and used for the estimation of the seed oil content (expressed in percentage) by NIRS (Mika et al., 2003). Erucic acid content was estimated by GC using a Perkin Elmer AutosystemXL GC with an Elite-FFAP column (3 mm length × 0.25 mm ID; 0.25 μm thick film) from the selfed seeds following an established protocol (Thies, 1971). The trait data for VHA8B7/Z and TDA8 populations were retrieved from earlier reports (Gupta et al., 2004; Ramchiary et al., 2007; Jagannath et al., 2011; Yadava et al., 2012). All the estimations were done in triplicate and the mean value was used as the trait value.
Linkage Map Construction and QTL Analysis
For QTL mapping, eight bi-parental linkage maps were developed using the five SE and three ZE mapping populations (Table 3). Genotyping of the mapping populations was done by using various marker systems, such as IP (Intron Polymorphism) (Panjabi et al., 2008), genic SSRs (Dhaka et al., 2017), genomic SSRs (Padmaja et al., 2005; Kim et al., 2009; Xu et al., 2010), SNPs (Paritosh et al., 2014), RFLP and AFLP markers (Pradhan et al., 2003). The linkage maps were constructed using JoinMap 4.0 (Van Ooijen, 2006) following the parameters used by Pradhan et al. (2003). The linkage groups were established with a minimum LOD threshold of 7.0, and the Kosambi function was used to convert recombination frequencies to cM distances. Of the eight bi-parental maps used in the study, three maps namely, EJ8A8B7, EJ8Z, and DEZ were constructed in the present study (Table 3). EJ8A8B7 map consists of 388 markers and covers a total map length of 1125.6 cM. EJ8Z and DEZ maps consist of 470 and 650 markers and cover map lengths of 1244.0 and 1680.2 cM, respectively. The remaining five bi-parental maps that have been reported earlier (Table 3), were updated with some SSR markers to enhance the genome coverage. Detailed map characteristics of all the eight bi-parental linkage maps used for the QTL analysis are given in Table 3.
Composite interval mapping (CIM) in Windows QTL Cartographer 2.5 (Wang et al., 2012), with a window size of 10 cM and a walk speed of 1 cM with five control markers, was performed for QTL analysis. For declaring the presence of a QTL, genome-wide threshold values (P = 0.05) were estimated from 1,000 permutations of trait data across all genetic intervals (Churchill and Doerge, 1994; Doerge and Churchill, 1996). QTL detected with more than 10% of the phenotypic variance (R2) were considered as major QTL. In addition, to correlate structures from among the two quantitative traits—erucic acid and oil content in seeds, we used multiple-trait mapping (Mt-CIM) implemented in Windows QTL Cartographer 2.5 to evaluate the pairs of variables jointly (Jiang and Zeng, 1995). This analysis was based on composite interval mapping similar to the procedure adopted for QTL analysis of single traits.
QTL meta-analysis was performed using BioMercator 2.1 (Arcade et al., 2004). Initially, one integrated map (INT map) was developed incorporating the eight bi-parental maps through iterative projection of genetic maps of both the SE and ZE populations on to the VHA8B7 map (reference map) using the ‘Map projection' command. The markers mapped on the integrated map were assigned to the Arabidopsis genomic blocks (Schranz et al., 2006) based on their similarity with B. rapa genes available in the Brassica database (http://brassicadb.org/brad), as obtained through BLASTN search.
For QTL meta-analysis, the identified QTL in the respective populations were projected on to the INT map and consensus QTL were discerned using a two-step approach. In the first step, QTL obtained from the five SE populations were integrated to identify consensus QTL that included the merged QTL from different mapping populations into a single QTL and was named as a C-QTL. Those QTL which did not merge were designated as “population-specific” QTL. In the second step, QTL identified from both SE and ZE populations were pooled together to identify the “common” QTL between the SE and ZE populations. The name of “common” QTL started with the trait name followed by LG name and the serial number of the QTL on that LG (e.g., Oil-A2-1).
The allelic effects of parental lines at these QTL merging into a C-QTL and “common” QTL were estimated following the procedure of Chardon et al. (2014), followed by elucidation of the allelic status of the parental lines using Heera/EH-2 as the common parent (Rout et al., 2015).
Trait Variation in Parents and Mapping Populations
Trait values for oil content and erucic acid content among the seven parental lines used in the study for genetic mapping of oil content trait are given in Table 1. These parental lines covered maximum phenotypic variability reported for the two traits in B. juncea. The phenotypic variation for the oil content ranged from ~36.9% in EH-2 to ~49.3% in J8. For the erucic acid content, two parents Heera and EH-2, both ZE lines, harbor the recessive alleles for the two erucic acid genes on the LGs A08 and B07. The lines, Varuna, J8, Pusa Jaikisan and TM-4 with high erucic acid content (>47.0%) contained both the dominant alleles for high erucic acid content on the LGs A08 and B07. The parental line, Donskaja-IV with intermediate levels of erucic acid (~20.5%) is homozygous dominant for the LG B07 locus and homozygous recessive for the LG A08 locus of the erucic acid trait (Table 1).
Frequency distributions of oil content showed continuous variation in all the eight mapping populations in all the locations, signifying the quantitative nature of the trait (Supplementary Figure 2). Transgressive segregation in both the directions was apparent in four of the five SE populations (VHA8B7, EPJA8B7, DEB7, and TDA8) whereas in the EJ8A8B7 population, transgressive variation beyond the high value parent was not observed. It can be concluded that the parental line J8 with the highest trait value has accumulated maximum number of positive alleles for oil content. A very high positive correlation (significant at 1% level) between erucic acid content and seed oil content was observed in the five SE populations in all the environments (Supplementary Table 1). The correlation coefficients ranged from 0.77 in EJA8B7 population to 0.39 in DEB7 population.
All the three ZE populations showed reduction in their mean oil content as compared to the respective SE populations (Table 2). ZE populations did not show any transgressive segregation for oil content beyond the contents in the high value parents in all the environments (Supplementary Figure 2).
QTL Analysis of Seed Erucic Acid Content (Eru) in SE Populations
QTL analysis of erucic acid content in EJ8A8B7, EPJA8B7, and VHA8B7 populations, identified major QTL on the LGs A08 and B07 in all the three environments co-mapping to the two FAE1 loci (Supplementary Table 2). The average phenotypic variance (R2) explained by these QTL in these three populations varied from 87.5% (in EPJA8B7 population) to 94.0% (in VHA8B7 population). A few minor loci explaining 1.0 to 2.1% of phenotypic variance were also identified in the EJ8A8B7and EPJA8B7 populations.
QTL analysis of erucic acid in the other two SE populations (DEB7 and TDA8) wherein only one locus was segregating for erucic acid, detected major QTL for the trait in all the environments co-mapping with the respective FAE1 locus. In DEB7 population segregating for erucic acid on LG B07, the average phenotypic variance explained by the major QTL was 90.7% whereas in the TDA8 population segregating for erucic acid on LG A08, the average phenotypic variance explained by the major QTL was 46.6%. In addition, in both the populations some minor QTL were also detected (Supplementary Table 2).
QTL Analysis of Seed Oil Content (Oil) in SE Populations
Mapping of the oil content in the EJ8A8B7 population showing maximum parental contrast for the oil content, detected a total of 15 Oil QTL from three environments (Supplementary Table 3). These QTL were distributed over two A sub-genome LGs (A07 and A08) and six B sub-genome LGs (B01, B02, B03, B04, B06, and B07). The phenotypic variance (R2) explained by the individual QTL varied from 2.5 to 43.3% with additive effects ranging from 0.4 to 1.8%. Six QTL were detected as major effect QTL, three each on LGs A08 (mean R2 = 34.0%) and B07 (mean R2 = 12.9%), all co-mapping with the two FAE1 loci on the corresponding LGs—A08 and B07. Remaining nine QTL were detected as minor QTL with phenotypic variance explained ranging from 2.5 to 7.7%. All the QTL, except Oil-B2-1-EJ and Oil-B2-2-EJ, had the trait enhancing allele from the high oil parent J8.
The EPJA8B7 population detected the highest number of 23 Oil QTL from three environments (Supplementary Table 3). These QTL were distributed over nine LGs, four in the A sub-genome and five in the B sub-genome. The phenotypic variance explained by these 23 Oil QTL ranged from 2.6 to 35.2% wherein Pusa Jaikisan and EH-2 contributed the positive alleles at 11 and 12 QTL, respectively. Of these 23 QTL, three QTL that co-mapped to the erucic acid QTL region on LG A08 were detected as major effect QTL (mean R2 = 33.2%). All other QTL including the three QTL in erucic acid region on LG B07 were detected as minor QTL.
The Oil QTL analysis in VHA8B7 population detected a total of 12 QTL of which eight QTL were detected on LGs A07, A08, A09, and A10, and four QTL on the LGs B07 and B08 (Supplementary Table 3). Eight QTL had trait enhancing allele from the high-value parent Varuna and four QTL from the low-value parent Heera. Five of the seven major QTL were co-mapped to the erucic acid QTL regions on the LGs A08 and B07.
The DEB7 population, segregating only for the erucic acid locus on LG B07 (Table 1), identified a total of 13 Oil QTL distributed over eight linkage groups, four LGs each of the A and B sub-genomes (Supplementary Table 3). The phenotypic variance explained by these QTL ranged from 4.7 to 23.0%. In six of the 13 QTL, the positive allele was contributed by Donskaja-IV, whereas EH-2 contributed the remaining seven QTL. Three major QTL were detected on LG B07, all of which co-mapped to the erucic acid QTL region. On LG A08, a minor QTL (Oil-A8-1-DE) was detected from a single location that mapped at a distance of ~25 cM from the position of the FAE1.1 gene.
The TDA8 population showing segregation for erucic acid locus on LG A08 had dominant allele in both the parents for the erucic acid locus on LG B07 (Table 1). QTL analysis identified nine Oil QTL from three environments (Supplementary Table 3). These QTL were distributed over three LGs of the A sub-genome and four LGs of the B sub-genome. The phenotypic variance explained by the QTL ranged from 7.2 to 22.4%. Of the nine QTL detected in the population, seven had acquired the positive allele from TM-4, while two had positive allele from Donskaja-IV. Six QTL were detected as major QTL, of which two QTL overlapped with the erucic acid QTL containing the erucic acid gene on LG A08 (Supplementary Table 3).
Oil QTL analysis from five bi-parental SE populations identified a total of 72 QTL, of which 40 were derived from the parental lines with high trait value, and the remaining 32 had the trait enhancing allele from the low oil content parents. The identified QTL were found to be distributed over 16 of the 18 LGs of Brassica juncea, except LGs A06 and B05, with almost equal contribution from both the A and B sub-genomes showing a wide distribution of loci influencing oil content. A total of 25 QTL were detected as major QTL with R2 greater than 10%, of which 19 QTL overlapped with erucic acid QTL regions including six high effect QTL explaining ≥25% of the phenotypic variance of oil content (Supplementary Table 3).
Meta-QTL Analysis of Oil Content in SE Populations
An integrated (INT) map was developed by merging the markers from eight bi-parental maps and contained a total of 3510 markers including 974 Intron Polymorphism, 917 SNP, 535 genic SSR, 428 genomic SSR, 645 AFLP, and 11 RFLP markers covering a genetic length of 2002.6 cM (Supplementary Table 4). About 70% of the mapped markers could be assigned Arabidopsis thaliana gene IDs and the corresponding genomic blocks of Arabidopsis (Schranz et al., 2006).
The 72 Oil QTL identified from the five SE mapping populations were distributed on the 16 LGs of B. juncea. QTL meta-analysis was undertaken only for 65 QTL that were distributed over 12 LGs, while excluding the remaining seven QTL which were detected only in one of the five bi-parental mapping populations. Meta-analysis resulted in nine C-QTL (Supplementary Table 4) distributed over nine LGs by merging a total of 47 Oil QTL mapped in the individual populations (Table 4, Supplementary Figure 3). The number of component QTL that formed a C-QTL varied from two to 12. No C-QTL was formed by merging component QTL from all the five SE populations. Seven out of the nine C-QTL merged the component QTL from two populations, each. Wherein, the remaining two C-QTL were formed by merging component QTL from four populations, each (Table 4).
Table 4. Details of consensus QTL (Oil-C) detected by merging the QTL data of five SE populations (EJ8A8B7, EPJA8B7, VHA8B7, DEB7, and TDA8) through QTL meta-analysis.
The two C-QTL—Oil-C-A8-1 and Oil-C-B7-1 obtained from four populations, merged 11 and 12 component QTL, respectively. These two C-QTL collectively accommodated 19 of the 25 major component Oil QTL, and co-mapped with the two erucic acid QTL on LGs A08 and B07 with the trait enhancing alleles from the HE parents (Table 4, Supplementary Figure 3). The remaining QTL that did not merge into C-QTL were identified as “population-specific” QTL (Supplementary Table 3).
The allelic effects of the seven parents used for the development of the mapping populations were calculated for the nine C-QTL. Comparative allelic strength of each of the seven parents for the nine C-QTL was determined following Chardon et al. (2014). This analysis revealed that the parent J8, possessing highest oil content (~49%), had the most positive alleles at five of the nine C-QTL, followed by TM-4 and Varuna possessing most positive alleles at three and one C-QTL, respectively (Table 5).
Table 5. Estimation of allelic effects (in percent) of nine Oil C-QTL involving five SE mapping populations and comparative allelic status of the seven parental lines.
QTL Mapping of Oil Content in ZE Populations
Of the three ZE populations, the QTL analysis of EJ8Z population identified a total of 10 QTL from three environments distributed over five LGs (Supplementary Table 5). The phenotypic variance explained by the detected Oil QTL ranged from 4.5 to 21.1%. Of these, three were detected as major QTL with R2 greater than 10.0%. Eight of the 10 QTL had the positive allele from J8, while EH-2 contributed positive alleles to the remaining two QTL.
The DEZ population detected a total of seven QTL from three year data. These QTL were distributed over four LGs (Supplementary Table 5). All the seven QTL were detected as major effect QTL with R2 varying from 10.8 to 24.2%. Six of the seven QTL had the trait enhancing allele from the high oil parent, Donskaja-IV.
The QTL re-analysis in VHZ population (Jagannath et al., 2011) with phenotyping data of one location identified a total of two QTL on LGs A02 and A03, of which the QTL on LG A02 was detected as a major QTL. Both the QTL had the trait enhancing allele from the low oil content parent Heera (Supplementary Table 5).
Oil QTL analysis from the three ZE populations detected a total of 19 QTL including 10 QTL from EJ8Z and seven QTL from DEZ populations which were phenotyped in three environments, and two QTL from VHZ population that was phenotyped only in one environment. However, in all the three ZE populations, no Oil QTL was detected in the erucic acid QTL regions of LGs A08 and B07 unlike the five SE populations that had most of the major effect QTL influencing oil content in that region. Of the 19 QTL detected in ZE populations, 11 were detected as major QTL with R2 more than 10.0%. Due to lack of sufficient number of component QTL from more than one mapping population in each of the 18 linkage groups, QTL meta-analysis of oil content in ZE populations was not undertaken.
Identification of Novel and Common QTL Between SE and ZE Populations
QTL meta-analysis between the Oil QTL of SE populations and ZE populations was performed to identify the “common” QTL between the two types of populations. The analysis revealed the merger of 17 SE Oil QTL with 13 ZE Oil QTL and resulted in identification of nine “common” Oil QTL between SE and ZE populations distributed over seven LGs (A02, A03, A08, B01, B04, B06, and B08) (Figure 1). The remaining six QTL discerned in the ZE populations that did not merge to form “common” QTL, were identified as “novel” QTL—specific to ZE populations (Table 6). Of these six “novel” QTL, J8 had the trait enhancing alleles at four QTL, while Donskaja-IV and EH-2 contributed the positive alleles at one “novel” QTL, each.
Figure 1. “Common” Oil QTL between SE and ZE populations identified by QTL meta-analysis. Meta-QTL regions are highlighted in red. Values in parentheses represent the phenotypic variance explained (R2) of the respective component QTL. Asterisk on LG A08 marks the map position of the FAE1.1 gene.
The allelic strength of the five parental lines at four of the nine “common” QTL merging the QTL from the same source populations (LG A08 in DE, B01 in EJ8, B04 in DE, and B06 in EJ8) identified— J8 and Donskaja-IV possessing the positive alleles for two “common” QTL, each. The allelic effects of these parental lines at corresponding “common” QTL were observed to have increased in ZE populations as compared to the SE populations (Table 7). At eight of the nine “common” QTL, a comparison of average R2 values of these nine “common” QTL between SE and ZE populations showed a general trend of increase in R2 values of ZE Oil QTL as compared to the SE Oil QTL (Figure 2). Identification of the “novel” QTL, and increased allelic effects and R2 values of SE QTL at “common” QTL therefore, revealed a significant change in the landscape of oil QTL in ZE background.
Table 7. Comparison of allelic effects (%) of Oil QTL underlying the “common” QTL between SE and ZE background in the mapping populations resourced from the same parental crosses.
Figure 2. Comparison of average phenotypic variance explained (R2) by “common” QTL identified between SE and ZE populations after QTL meta-analysis.
The 15 multi-population QTL identified in the present study (Table 8) were projected on to the reference genome of B. juncea (Yang et al., 2016), B. rapa (Wang et al., 2011), B. nigra (Yang et al., 2016), and B. napus (Chalhoub et al., 2014), using the sequences of the PCR primers flanking the identified QTL through “BLASTn” tool available in the BRAD database (http://brassicadb.org/brad/blastPage.php). Mapping of the QTL detected on A-subgenome of B. juncea identified 17 corresponding QTL (Liu et al., 2016; Raboanatahiry et al., 2017) distributed over three A-subgenome LGs of B. napus genome. In the absence of any QTL data from B. nigra and B-subgenome of B. carinata for oil content, no such comparative analysis could be undertaken for the B-subgenome Oil QTL of B. juncea.
Table 8. Summary of Oil QTL identified in B. juncea and their physical positions on the reference genomes of B. rapa, B. nigra, B. juncea, and B. napus.
Pleiotropic Effect of Erucic Acid On Oil Content
As the investigated traits—erucic acid and oil content were positively correlated (Supplementary Table 1) and were identified as major QTL for both the traits in overlapping intervals on LGs A08 and B07, Mt-CIM (Jiang and Zeng, 1995) was used to distinguish whether the correlation was an artifact due to two closely-linked QTL or a single QTL was affecting both the traits through pleiotropy. The analysis was undertaken for the five SE populations (EJ8A8B7, EPJA8B7, VHA8B7, DEB7, and TDA8) for LGs A08 and B07. In all the five SE populations, a joint analysis on LGs A08 and B07, revealed a similar pattern of LOD profile for the joint trait and for the two individual traits. These QTL were significant for both the individual as well as joint traits in all the environments (Figure 3), indicating existence of strong pleiotropy between the two traits. Hence the two C-QTL, Oil-C-A8-1 on LG A08 and Oil-C-B7-1 on LG B07 in SE populations, which were not detected in the corresponding ZE populations, were designated as pleiotropic QTL. The effect of pleiotropy was also evident from a significant change in QTL landscape in ZE background (Figure 2, Table 6).
Figure 3. Multi-trait mapping (Mt-CIM) of seed oil content and erucic acid content in five SE bi-parental mapping populations of B. juncea showing pleiotropism between erucic acid and oil content (A: LG A08 of EJ8A8B7 population, B: LG B07 of EJ8A8B7 population, C: LG A08 of EPJA8B7 population, D: LG B07 of EPJA8B7 population, E: LG A08 of VHA8B7 population, F: LG B07 of VHA8B7 population, G: LG B07 of DEB7 population, and H: LG A08 of TDA8 population). The horizontal lines indicate LOD significance threshold.
Breeding for higher seed oil content in both HE and ZE mustard is critical for enhancing the productivity of the mustard crop. To address this problem effectively, seven inbred parental lines containing significant phenotypic variation for the trait (ranging from 37% to 49%) in the B. juncea germplasm were used for the development of two sets of bi-parental mapping populations—one set of bi-parental mapping populations segregating for erucic acid content (SE populations) and the other set consisting of only zero-erucic DH lines (ZE populations). A comparative QTL analysis between the two types of mapping populations has provided a detailed insight into the QTL landscape for seed oil content trait in B. juncea.
Seven inbred parental lines from which five SE and three ZE mapping populations were sourced consisted of lines with very high (49%), intermediate (41 to 44%), and low (37%) oil content-almost representing the maximum phenotypic variation reported for the trait in the B. juncea germplasm. The SE mapping populations showed higher mean values and presence of transgressive segregation in both the directions as compared to the ZE populations for oil content in the study. This difference in the mean and the spectrum of transgressive segregation observed between the two types of populations could probably be due to the influence of erucic acid on oil content as these two traits are positively correlated (Supplementary Table 1). Gu et al. (2017) reported positive correlation (r = 0.483) between oil content and long chain fatty acid particularly erucic acid in a set of lines of B. napus consisting of ultrahigh (>55%), high (50–55%), medium (40–50%), and low (30–40%) oil groups.
QTL analysis in five SE populations revealed distribution of Oil QTL throughout the genome (in 16 of the 18 LGs) with a sizable number of minor QTL having R2 less than 10%. Our results in B. juncea, support the earlier findings in rapeseed that the Oil QTL are distributed throughout the genome (Wang et al., 2013; Jiang et al., 2014). Comparison of Oil QTL data between SE and ZE populations revealed that none of the 23 QTL identified in SE populations on LGs A08 and B07 that co-mapped to the erucic acid QTL regions were detected in ZE populations. These Oil loci constitute 76% of the major QTL detected in the five SE populations (19 of the 25 major QTL) and merge into two consensus QTL in SE populations (Oil-C-A8-1 on LG A08 and Oil-C-B7-1 on LG B07). The multi-trait analysis for erucic acid and oil content on these two LGs also indicated existence of a strong pleiotropic effect between these two traits. Hence, these two C-QTL were designated as “Pleiotropic” QTL. It indicates that these two consensus Oil QTL, although the key QTL in SE populations, cannot be used to improve the oil content in ZE B. juncea.
Pleiotropic effect of erucic acid on oil content has also been reported earlier in both B. napus (Ecke et al., 1995; Burns et al., 2003) and B. juncea (Jagannath et al., 2011). Ecke et al. (1995) interpreted this occurrence owing to increase in the molecular mass of the carbon atoms of fatty acids during the elongation of oleic acid to erucic acid and resulting in a corresponding increase in the oil content in B. napus. As erucic acid is the major fatty acid (≥45% of the fatty acid pool) in B. juncea as in case of non-Canola B. napus and shows high positive correlation with oil content, similar mechanisms should be operational in B. juncea as well. A supporting evidence to this fact was provided by Katavic et al. (2000) by over-expression of FAE1 (fatty acid elongase) gene in transgenic B. napus that showed an increase in erucic acid content up to 16.8% and a simultaneous increase in the oil content by 10.9% in non-Canola B. napus cv. Hero. The pleiotropic Oil QTL on the LG A08 (Oil-C-A8-1) identified in the present study in B. juncea also corresponds to the physical position in B. napus where maximum number of Oil QTL from multiple populations have been mapped (Table 8) (Liu et al., 2016; Raboanatahiry et al., 2017) indicating that these QTL in B. napus could also be “pleiotropic” QTL. However, Delourme et al. (2006) while comparing the QTL data of oil content from two bi-parental populations of B. napus interpreted that the Oil QTL detected in the erucic acid region of LG A08 could be due to linkage and not due to pleiotropy. This contention was based on the fact that a population (DY population) segregating for erucic acid detected no Oil QTL in the erucic acid region of LG A08 while another population (RNSL population) that was not segregating for erucic acid detected an Oil QTL in the corresponding region. Our results in the present study also detected an Oil QTL (Oil-A8-1) in the LG A08 in both SE and ZE populations of DE which does not segregate for the A08 erucic acid locus. The physical position of this Oil QTL is 5.98 Mb away from the pleiotropic QTL locus Oil-C-A8-1 (Table 8).
In addition to the two pleiotropic QTL (Oil-C-A8-1 and Oil-C-B7-1), this study has identified seven C-QTL (Table 5) and two major effect “population-specific” QTL–Oil-A3-1-TD and Oil-B3-1-TD (Supplementary Table 3) with significant contribution to the trait variability. The comparative allelic effects deciphered at the seven C-QTL identified J8 as the most desirable parental line for four loci (Oil-C-A7-1, Oil-C-B1-1, Oil-C-B3-1, and Oil-C-B6-1) and TM-4 for the remaining three loci (Oil-C-A10-1, Oil-C-B2-1, and Oil-C-B4-1)—these loci can be deployed in marker assisted breeding for improvement of oil content in HE mustard.
In the present study, we observed significant changes in the genetic landscape of Oil QTL when the influence of the erucic acid locus was absent. Eight out of the nine common Oil QTL between SE and ZE populations showed an upward trend in the phenotypic variance explained and also increase in the allelic effects of the respective parent in the ZE populations. It could probably be due to sharing of that part of the phenotypic variance owing to pleiotropic effect in SE populations that is reallocated among the Oil QTL in ZE populations when the pleiotropic effect of erucic acid is absent. Due to this reason, five QTL on LGs A03 (Oil-A3-1), A08 (Oil-A8-1), B01 (Oil-B1-1), B04 (Oil-B4-1), and B06 (Oil-B6-1) which were detected as minor QTL in SE populations were detected as major QTL in ZE populations (Figure 2). In addition, the appearance of six “novel” QTL in ZE populations (Table 6) could probably be due to exclusion of the effect of erucic acid on seed oil content. Jagannath et al. (2011) also reported a similar observation in B. juncea while comparing the SE and ZE QTL information from two mapping populations that were derived from one bi-parental cross. The present study, therefore, suggests that the genetic analysis of oil content in populations with high erucic acid content may not be relevant in breeding for high oil content in ZE mustard.
The present study incorporated multiple mapping populations derived from seven parental lines with significant variation for oil content, including the pureline J8 with an oil content of ~49%. At four of the “common” QTL and four of the “novel” QTL identified in this study, J8 holds stronger alleles as compared to the other parental lines, positioning it as an excellent source for improvement of oil content in ZE background. The nine “common” QTL (Figure 1) and the six “novel” QTL (Table 6) identified here can be highly useful in compensating for the reduction in oil content in “0” erucic lines. Fine mapping of the important loci identified here will allow identification of the candidate genes and their interactions involved in oil content variation in B. juncea.
KR and BY contributed equally and carried out the phenotyping and mapping work. SY helped in field experiments and data analysis. VG helped with genotyping. AM developed microspore derived F1DH mapping populations. DP and AP conceived and supervised the overall study and along with KR wrote the manuscript. All authors read and approved the final manuscript.
The study was funded by the Department of Biotechnology (DBT), Government of India through the award of two projects—(1) Centre of Excellence on Genome Mapping and Molecular Breeding of Brassicas (BT/01/COE/08/06/-II) and (2) DBT-UDSC Partnership Centre on Genetic Manipulation of Brassicas (BT/01/NDDB/UDSC/2016).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Support of the National Dairy Development Board (NDDB) for the research till September, 2016 is also acknowledged. KR and BY acknowledge the receipt of research fellowship from the NDDB and UGC, respectively. DP is supported with a J.C. Bose fellowship from the Department of Science & Technology (DST) and by the Council of Scientific and Industrial Research (CSIR) as a Distinguished Scientist.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01448/full#supplementary-material
Ackman, R. G., Eaton, C. A., Sipos, J. C., Loew, F. M., and Hancock, D. (1977). Comparison of fatty acids from high levels of erucic acid of RSO and partially hydrogenated fish oil in non-human primate species in a short-term exploratory study. Nutr. Diet 25, 170–185.
Arcade, A., Labourdette, A., Falque, M., Mangin, B., Chardon, F., Charcosset, A., et al. (2004). BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics 20, 2324–2326. doi: 10.1093/bioinformatics/bth230
Burns, M. J., Barnes, S. R., Bowman, J. G., Clarke, M. H. E., Werner, C. P., and Kearsey, M. J. (2003). QTL analysis of an inter-varietal set of substitution lines in Brassica napus: seed oil content and fatty acid composition. Heredity 90, 39–48. doi: 10.1038/sj.hdy.6800176
Cao, Z., Tian, F., Wang, N., Jiang, C., Lin, B., Xia, W., et al. (2010). Analysis of QTLs for erucic acid and oil content in seeds on A8 chromosome and the linkage drag between the alleles for the two traits in Brassica napus. J. Genet. Genomics 37, 231–240. doi: 10.1016/S1673-8527(09)60041-2
Chalhoub, B., Denoeud, F., Liu, S., Parkin, I. A., Tang, H., Wang, X., et al. (2014). Early allopolyploid evolution in the post-neolithic Brassica napus oilseed genome. Science 346, 950–953. doi: 10.1126/science.1253435
Chardon, F., Jasinski, S., Durandet, M., Lécureuil, A., Soulay, F., Bedu, M., et al. (2014). QTL meta-analysis in Arabidopsis reveals an interaction between leaf senescence and resource allocation to seeds. J. Exp. Bot. 65, 3949–3962. doi: 10.1093/jxb/eru125
Cheung, W. Y., Landry, B. S., Raney, P., and Rakow, G. F. W. (1998). Molecular mapping of seed quality traits in Brassica juncea L. czern. and coss. Acta Hortic. 459, 139–148. doi: 10.17660/ActaHortic.1998.459.15
Delourme, R., Falentin, C., Huteau, V., Clouet, V., Horvais, R., Gandon, B., et al. (2006). Genetic control of oil content in oilseed rape (Brassica napus L.). Theor. Appl. Genet. 113, 1331–1345. doi: 10.1007/s00122-006-0386-z
Dhaka, N., Mukhopadhyay, A., Paritosh, K., Gupta, V., Pental, D., and Pradhan, A. K. (2017). Identification of genic SSRs and construction of a SSR-based linkage map in Brassica juncea. Euphytica 213, 15. doi: 10.1007/s10681-016-1814-z
Ecke, W., Uzunova, M., and Weißleder, K. (1995). Mapping the genome of rapeseed (Brassica napus L.). II. Localization of genes controlling erucic acid synthesis and seed oil content. Theor. Appl. Genet. 91, 972–977. doi: 10.1007/BF00223908
Fu, Y., Zhang, D., Gleeson, M., Zhang, Y., Lin, B., Hua, S., et al. (2017). Analysis of QTL for seed oil content in Brassica napus by association mapping and QTL mapping. Euphytica 213, 17. doi: 10.1007/s10681-016-1817-9
Gu, J., Chao, H., Wang, H., Li, Y., Li, D., Xiang, J., et al. (2017). Identification of the relationship between oil body morphology and oil content by microstructure comparison combining with QTL analysis in Brassica napus. Front. Plant Sci.7:1989. doi: 10.3389/fpls.2016.01989
Gupta, V., Mukhopadhyay, A., Arumugam, N., Sodhi, Y. S., Pental, D., and Pradhan, A. K. (2004). Molecular tagging of erucic acid trait in oilseed mustard (Brassica juncea) by QTL mapping and single nucleotide polymorphisms in FAE1 gene. Theor. Appl. Genet. 108, 743–749. doi: 10.1007/s00122-003-1481-z
Jagannath, A., Sodhi, Y. S., Gupta, V., Mukhopadhyay, A., Arumugam, N., Singh, I., et al. (2011). Eliminating expression of erucic acid-encoding loci allows the identification of “hidden” QTL contributing to oil quality fractions and oil content in Brassica juncea (Indian mustard). Theor. Appl. Genet. 122, 1091–1103. doi: 10.1007/s00122-010-1515-2
Jiang, C., Shi, J., Li, R., Long, Y., Wang, H., Li, D., et al. (2014). Quantitative trait loci that control the oil content variation of rapeseed (Brassica napus L.). Theor. Appl. Genet. 127, 957–968. doi: 10.1007/s00122-014-2271-5
Katavic, V., Friesen, W., Barton, D. L., Gossen, K. K., Giblin, E. M., Luciw, T., et al. (2000). Utility of the Arabidopsis FAE1 and yeast SLC1-1 genes for improvements in erucic acid and oil content in rapeseed. Biochem. Soc. Trans. 28, 935–937. doi: 10.1042/bst0280935
Kim, H., Choi, S. R., Bae, J., Hong, C. P., Lee, S. Y., Hossain, M. J., et al. (2009). Sequenced BAC anchored reference genetic map that reconciles the ten individual chromosomes of Brassica rapa. BMC Genomics 10:43. doi: 10.1186/1471-2164-10-432
Liu, S., Fan, C., Li, J., Cai, G., Yang, Q., Wu, J., et al. (2016). A genome-wide association study reveals novel elite allelic variations in seed oil content of Brassica napus. 129, 1203–1215. doi: 10.1007/s00122-016-2697-z
Mahmood, T., Rahman, M. H., Stringam, G. R., Yeh, F., and Good, A. G. (2006). Identification of quantitative trait loci (QTL) for oil and protein contents and their relationships with other seed quality traits in Brassica juncea. Theor. Appl. Genet. 113, 1211–1220. doi: 10.1007/s00122-006-0376-1
Mika, V., Tillmann, P., Kopma, R., Nerusil, P., and Kucera, V. (2003). Fast prediction of quality parameters in whole seeds of oilseed rape (Brassica napus L.). Plant Soil Environ. 49, 141–145. doi: 10.17221/4104-PSE
Mukhopadhyay, A., Arumugam, N., Sodhi, Y. S., Gupta, V., Pradhan, A. K., and Pental, D. (2007). “High frequency production of microspore derived doubled haploids (DH) and its application for developing low glucosinolate lines in Indian Brassica juncea,” in 12th International Rapeseed Congress (Wuhan), 333–335.
Padmaja, K. L., Arumugam, N., Gupta, V., Mukhopadhyay, A., Sodhi, Y. S., Pental, D., et al. (2005). Mapping and tagging of seed coat color and the identification of microsatellite markers for marker-assisted manipulation of the trait in Brassica juncea. Theor. Appl. Genet. 111, 8–14. doi: 10.1007/s00122-005-1933-8
Panjabi, P., Jagannath, A., Bisht, N. C., Padmaja, K. L., Sharma, S., Gupta, V., et al. (2008). Comparative mapping of Brassica juncea and Arabidopsis thaliana using Intron Polymorphism (IP) markers: homoeologous relationships, diversification and evolution of the A, B and C Brassica genomes. BMC Genomics 9:113. doi: 10.1186/1471-2164-9-113
Paritosh, K., Gupta, V., Yadava, S. K., Singh, P., Pradhan, A. K., and Pental, D. (2014). RNA-seq based SNPs for mapping in Brassica juncea (AABB): Synteny analysis between the two constituent genomes A (from B. rapa) and B (from B. nigra) shows highly divergent gene block arrangement and unique block fragmentation patterns. BMC Genomics 15:396. doi: 10.1186/1471-2164-15-396
Pradhan, A. K., Gupta, V., Mukhopadhyay, A., Arumugam, N., Sodhi, Y. S., and Pental, D. (2003). A high density linkage map in Brassica juncea (Indian mustard) using AFLP and RFLP markers. Theor. Appl. Genet. 106, 607–614. doi: 10.1007/s00122-002-1083-1
Qiu, D., Morgan, C., Shi, J., Long, Y., Liu, J., Li, R., et al. (2006). A comparative linkage map of oilseed rape and its use for QTL analysis of seed oil and erucic acid content. Theor. Appl. Genet. 114, 67–80. doi: 10.1007/s00122-006-0411-2
Raboanatahiry, N., Chao, H., Guo, L., Gan, J., Xiang, J., Yan, M., et al. (2017). Synteny analysis of genes and distribution of loci controlling oil content and fatty acid profile based on QTL alignment map in Brassica napus. BMC Genomics 18:776. doi: 10.1186/s12864-017-4176-6
Ramchiary, N., Padmaja, K. L., Sharma, S., Gupta, V., Sodhi, Y. S., Mukhopadhyay, A., et al. (2007). Mapping of yield influencing QTL in Brassica juncea: implications for breeding of major oilseed crop of dryland areas. Theor. Appl. Genet. 115, 807–817. doi: 10.1007/s00122-007-0610-5
Rout, K., Sharma, M., Gupta, V., Mukhopadhyay, A., Sodhi, Y. S., Pental, D., et al. (2015). Deciphering allelic variations for seed glucosinolate traits in oilseed mustard (Brassica juncea) using two bi-parental mapping populations. Theor. Appl. Genet. 128, 657–666. doi: 10.1007/s00122-015-2461-9
Schranz, M. E., Lysak, M. A., and Mitchell-Olds, T. (2006). The ABC's of comparative genomics in the Brassicaceae, building blocks of crucifer genomes. Trends Plant Sci. 11, 535–542. doi: 10.1016/j.tplants.2006.09.002
Sun, M., Hua, W., Liu, J., Huang, S., Wang, X., Liu, G., et al. (2012). Design of New genome-and gene-sourced primers and identification of QTL for seed oil content in a specially high-oil Brassica napus cultivar. PLoS ONE 7:e47037. doi: 10.1371/journal.pone.0047037
Thies, W. (1971). Schnelle und einfache analysen der fettsaurzusammensetzung in einzelnen Raps-Kotyledonen I. Gaschromatographische and papier-chromatographischemethoden. Z. Pflanzenzuchtg. 65, 181–202.
Wang, S. C., Basten, C. J., and Zeng, Z. B. (2012). Windows QTL Cartographer 2.5. Department of Statistics; North Carolina State University, Raleigh, NC. Available online at: http://statgen.ncsu.edu/qtlcart/WQTLCart.htm
Wang, X., Wang, H., Long, Y., Li, D., Yin, Y., Tian, J., et al. (2013). Identification of QTLs associated with oil content in a high-oil Brassica napus cultivar and construction of a high-density consensus map for QTLs comparison in B. napus. PLoS ONE 8:e80569. doi: 10.1371/journal.pone.0080569
Xu, J., Qian, X., Wang, X., Li, R., Cheng, X., Yang, Y., et al. (2010). Construction of an integrated genetic linkage map for the A genome of Brassica napus using SSR markers derived from sequenced BACs in B. rapa. BMC Genomics 11:594. doi: 10.1186/1471-2164-11-594
Yadava, S. K., Arumugam, N., Mukhopadhyay, A., Sodhi, Y. S., Gupta, V., Pental, D., et al. (2012). QTL mapping of yield-associated traits in Brassica juncea: meta-analysis and epistatic interactions using two different crosses between east European and Indian gene pool lines. Theor. Appl. Genet. 125, 1553–1564. doi: 10.1007/s00122-012-1934-3
Yang, J., Liu, D., Wang, X., Ji, C., Cheng, F., Liu, B., et al. (2016). The genome sequence of allopolyploid Brassica juncea and analysis of differential homoeolog gene expression influencing selection. Nat. Genet. 48, 1225–1232. doi: 10.1038/ng.3657
Keywords: Brassica juncea, seed oil content, seed erucic acid content, QTL mapping, meta-analysis, pleiotropy
Citation: Rout K, Yadav BG, Yadava SK, Mukhopadhyay A, Gupta V, Pental D and Pradhan AK (2018) QTL Landscape for Oil Content in Brassica juncea: Analysis in Multiple Bi-Parental Populations in High and “0” Erucic Background. Front. Plant Sci. 9:1448. doi: 10.3389/fpls.2018.01448
Received: 08 June 2018; Accepted: 12 September 2018;
Published: 16 October 2018.
Edited by:Maoteng Li, Huazhong University of Science and Technology, China
Reviewed by:Leonardo Velasco, Instituto de Agricultura Sostenible (IAS), Spain
Harsh Raman, New South Wales Department of Primary Industries, Australia
Copyright © 2018 Rout, Yadav, Yadava, Mukhopadhyay, Gupta, Pental and Pradhan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Akshay K. Pradhan, firstname.lastname@example.org