- Research
- Open access
- Published:
Dissection of QTLs conferring drought tolerance in B. carinata derived B. juncea introgression lines
BMC Plant Biology volume 23, Article number: 664 (2023)
Abstract
Background
Drought is one of the important abiotic stresses that can significantly reduce crop yields. In India, about 24% of Brassica juncea (Indian mustard) cultivation is taken up under rainfed conditions, leading to low yields due to moisture deficit stress. Hence, there is an urgent need to improve the productivity of mustard under drought conditions. In the present study, a set of 87 B. carinata-derived B. juncea introgression lines (ILs) was developed with the goal of creating drought-tolerant genotypes.
Method
The experiment followed the augmented randomized complete block design with four blocks and three checks. ILs were evaluated for seed yield and its contributing traits under both rainfed and irrigated conditions in three different environments created by manipulating locations and years. To identify novel genes and alleles imparting drought tolerance, Quantitative Trait Loci (QTL) analysis was carried out. Genotyping-by-Sequencing (GBS) approach was used to construct the linkage map.
Results
The linkage map consisted of 5,165 SNP markers distributed across 18 chromosomes and spanning a distance of 1,671.87 cM. On average, there was a 3.09 cM gap between adjoining markers. A total of 29 additive QTLs were identified for drought tolerance; among these, 17 (58.6% of total QTLs detected) were contributed by B. carinata (BC 4), suggesting a greater contribution of B. carinata towards improving drought tolerance in the ILs. Out of 17 QTLs, 11 (64.7%) were located on the B genome, indicating more introgression segments on the B genome of B. juncea. Eight QTL hotspots, containing two or more QTLs, governing seed yield contributing traits, water use efficiency, and drought tolerance under moisture deficit stress conditions were identified. Seventeen candidate genes related to biotic and abiotic stresses, viz., SOS2, SOS2 like, NPR1, FAE1-KCS, HOT5, DNAJA1, NIA1, BRI1, RF21, ycf2, WRKY33, PAL, SAMS2, orf147, MAPK3, WRR1 and SUS, were reported in the genomic regions of identified QTLs.
Conclusions
The significance of B. carinata in improving drought tolerance and WUE by introducing genomic segments in Indian mustard is well demonstrated. The findings also provide valuable insights into the genetic basis of drought tolerance in mustard and pave the way for the development of drought-tolerant varieties.
Background
Global warming and climate change present a severe threat to crop production, including the occurrence of various biotic and abiotic stresses [1, 2]. Drought, among the abiotic stresses, significantly impacts the growth and productivity of crop plants [3, 4]. Over the last 6–7 decades, India has experienced a consistent decline in summer monsoon rainfall, leading to an increased risk of droughts [5]. Between 1951 and 2016, droughts became more frequent and widespread, particularly in central India, the southwest coast, the southern peninsula, and the north-eastern parts of the country, which experienced an average of two droughts per decade [5]. Furthermore, the drought-affected area has expanded by 1.3% per decade during this period [5]. Rapeseed-mustard, an important group of edible oilseed crops, is cultivated worldwide and occupies 41.64 million hectares (Mha) with a production of approximately 87.30 million metric tons (MMT). It contributes 13.91% to global oilseed production (627.44 MMT). In India, Indian mustard (B. juncea) dominates the rapeseed-mustard group, covering over 90% of the total acreage (9.00 Mha) and accounting for 28.45% of the country's oilseed production (40.42 MMT) in 2022–23 [6].
Higher production of Indian mustard is needed to meet the edible oil demand of the ever-growing population, which can be well achieved by improving productivity and reducing the yield losses caused by different biotic and abiotic stresses. Efforts have been made in the past to improve seed and oil yields to achieve self-sufficiency, but a huge quantity of edible oil is still imported annually. According to the Indian Vegetable Oil Producers’ Association, India imported 14.38 million metric tons (MMT) of oil, costing over 18 billion dollars accounting for more than 70% of the country's total edible oil demand in the year 2022–23 (Source: https://www.indiastatagri.com/). The rise in per capita oil consumption is driven by factors like population growth, increasing income and changing dietary preferences, and is expected to further increase [7]. To meet the edible oil demand of growing population, it is estimated that approximately 16.4–20.5 MMT of rapeseed-mustard needs to be produced [8, 9], while the current production stands at only 11.5 MMT [6].
To achieve self-sufficiency in edible oil, the productivity of Indian mustard needs to be urgently improved. However, the narrow genetic base of this species poses a significant constraint to its improvement [10]. The crop's susceptibility to various pests, diseases, and environmental stresses contributes to inconsistent production patterns. Approximately 24% of Indian mustard cultivation is taken up in rainfed areas, resulting in substantial yield losses due to moisture deficit stress [11]. Most commercially released B. juncea varieties are sensitive to drought, leading to critical seed yield losses. Although a few drought-tolerant varieties have been developed for drought-prone regions, their poor yield potential limits adoption. Therefore, the search for new sources imparting drought tolerance is urgently needed to reduce yield losses, particularly in the drought-prone areas of the eastern and western parts of the country [12].
Brassica species and their close relatives possess beneficial traits such as tolerance to cold, salinity and drought, which could be incorporated into present-day cultivars [13,14,15,16,17,18]. Wide hybridization has been considered a novel approach for the successful transfer of desirable traits and generating selectable genetic variability in cultivated species [19,20,21]. Ethiopian mustard, B. carinata (BBCC; 2n = 4x = 34), possesses resistance/tolerance to various abiotic and biotic stresses, including drought, heat, aphid, white rust, Sclerotinia rot, Alternaria black spot, powdery mildew and blackleg [22, 23]. Exploiting B. carinata as a donor for drought tolerance can enhance the performance and stability of Indian mustard cultivars under diverse arid environments [23,24,25].
Drought tolerance is a complex, polygenic trait reported in several crop species, including mustard [18]. To improve drought tolerance, approaches like introgressions from wild relatives, quantitative trait loci (QTL) mapping and marker-assisted backcross breeding (MABB) are being employed [26]. QTL mapping requires a large number of molecular markers with genome-wide coverage for high resolution. Simple sequence repeats (SSRs) have been extensively used for QTL mapping [27], but their limited number and uneven distribution in the genome hinder resolution. This has been greatly achieved by the advent of next-generation sequencing (NGS) techniques.
Next generation sequencing platforms offer a cost-effective solution for genotyping by sequencing single nucleotide polymorphic (SNP) sites, offering thousands to millions of molecular markers. Whole genome resequencing is an effective method when a reference genome is available, but it can be costly for crop species with large genomes. Therefore, genotyping-by-sequencing (GBS), a reduced representation genotyping approach based on restriction site-associated DNA sequencing (RAD-seq), has emerged as a novel and alternative method [28]. It has been successfully employed for high-resolution linkage mapping, particularly in introgression lines derived from related species of Brassicas [28,29,30]. However, very few studies reported QTLs associated with drought tolerance in rapeseed-mustard [18, 31];) and the mapping of QTLs for drought tolerance or improved water use efficiency in B. carinata or its derived lines remain largely unexplored.
With a view to improve the ability of B. juncea to withstand moisture deficit stress conditions, efforts were made at the Indian Council of Agricultural Research-Indian Agricultural Research Institute (ICAR-IARI), New Delhi to broaden the genetic base of B. juncea through the development of introgression lines (ILs) carrying genomic segments from B. carinata. The present study also aimed to identify QTLs/genes responsible for imparting drought tolerance and/or improved water use efficiency (WUE). To the best of our knowledge, this is the first instance where ILs possessing useful genomic segment(s) from B. carinata were developed and deployed for genetic studies conducted under moisture deficit stress conditions. This diverse set of ILs, evaluated under moisture deficit stress conditions in different environments, was used to map the QTL(s) and decipher the underlying candidate genes associated with WUE, drought tolerant indices and yield contributing traits. This study will lay the foundation for future investigations into molecular mechanisms underlying drought tolerance/WUE in mustard.
Materials and methods
Development and cytological investigation on B. carinata derived B. juncea introgression lines
A set of 87 B. juncea introgression lines, in F6 generation, carrying the genomic segments of B. carinata was developed. These lines were derived from the cross between B. juncea cv. DRMRIJ 31 and B. carinata acc. BC 4, and following biparental mating among phenotypically selected plants within F2 population (Fig. 1). Phenotypic selection of desirable segregants having a close resemblance to B. juncea parent in the subsequent filial generations led to the development of B. carinata derived B. juncea introgression lines (ILs). These ILs were identified as homozygous and cytologically stable, as reported in our earlier study on the expression of heterosis by them [32]. Further, molecular analysis on these ILs has confirmed the presence of introgression segments from B. carinata [33].
Phenotypic evaluation of ILs under moisture deficit stress conditions
A set of ILs along with their parents were evaluated at ICAR-IARI, New Delhi (L-1) and ICAR-Directorate of Rapeseed and Mustard Research (ICAR-DRMR), Bharatpur, Rajasthan (L-2). During 2018–19 crop season, the experiment was conducted at L1 location under rainfed (RE1) and irrigated (IE1) conditions, whereas during 2020–21 crop season it was done at both L1 and L2 locations under rainfed (RE2 and RE3) and irrigated (IE2 and IE3) conditions, thus consisting of a total of six environments. An augmented Randomised Complete Block Design (RCBD) with four blocks and three checks (Pusa Agrani, DRMRIJ 31 and Pusa Mustard 30) was used for the morphological characterization of ILs differing in metric traits and response to moisture deficit stress. Checks were replicated in each block under both rainfed and irrigated conditions. Each IL was raised in a paired row plot of four meters in length. Within each plot, row-row and plant-plant spacing were kept at 30 cm and 15 cm, respectively, and plots were separated by 60 cm. No irrigation was applied in the rainfed plots, while two irrigations of 50 mm were applied in the irrigated plots at 45 and 90 days after sowing (DAS). All recommended agronomic practices were followed for raising the crop.
Data were recorded on 15 quantitative traits, viz., plant height (cm), primary branches/plant, secondary branches/plant, main shoot length (cm), total siliquae on main shoot, siliqua length (cm), seeds/siliqua, seed yield/plant (g), biological yield/plant (g), seed yield/plot (g), harvest index, oil content (%), 1,000-seed weight (g), days to 50% flowering, and days to maturity. The data were recorded on five randomly selected competitive plants from each plot in each replication, except for seed yield per plot (g), days to 50% flowering and days to maturity, where observations were recorded on a plot basis. The biological yield/plant (g) was recorded at maturity after completely drying the harvested plants.
Statistical analysis
An analysis of variance for augmented randomized complete block design (Augmented-RCBD) was performed for all the studied traits following the statistical model,
where, Yijk is the mean values of a trait in ith environment of jth genotype in kth block; μ is the population mean; G, C and Ei is genotype, check and ith environment, respectively; G × Ei and C × Ei is the effect of the genotype by environment interaction and check by environment interaction, respectively; Bk (Ei) is the effect of the kth block within ith environment; and eilk is the residual error.
Best linear unbiased predictors (BLUP) values for all studied traits were estimated to perform combined analysis under both rainfed and irrigated conditions across the environments using Augmented Complete Block Design in R (ACBD-R) software [34, 35]. The descriptive statistics such as mean, range and broad-sense heritability ((h2); [36]) were estimated using Plant Breeding Tools software [37].
Assessment of drought tolerance
Adjusted mean values of all the studied traits under both rainfed and irrigated conditions were used for the estimation of drought susceptibility index (DSI), drought tolerance index (DTI), tolerance index (TOL) and mean relative performance (MRP) using the following formulae:
-
Drought Susceptibility Index (DSI) = (1–Ys / Yns) / D; [38]
-
Drought tolerance index (DTI) = (Ys × Yns) / (Ynsm).2; [39]
-
Tolerance index (TOL) = Yns–Ys; [40]
-
Mean relative performance: MRP = [Ys/Ysm + Yns/Ynsm]
Where, D is stress intensity and calculated as, D = 1– (Ysm / Ynsm); Ys = seed yield of genotype under moisture deficit stress conditions (drought stress); Yns = seed yield of genotype under irrigated or non-stress conditions; Ysm = mean seed yield of all genotypes under moisture deficit stress conditions; Ynsm = mean seed yield of all genotypes under irrigated or non-stress conditions.
Water-use efficiency (WUE) for ILs was estimated from the following formulae:
The minimum and maximum temperatures (°C) and rainfall (mm) were recorded during the crop growth seasons during 2018–19 and 2020–21 by the meteorological observatories located within 500 m distance from the trial sites at ICAR-IARI, New Delhi and ICAR-DRMR, Bharatpur (Fig. 2). The effective rainfall during crop season from November to March was calculated for respective season using Cropwat software (Version 8.0) by USDA SCS method. The water use efficiency in ILs was then estimated using effective rainfall observed.
Genotyping by sequencing
A total of 87 introgression lines (ILs) were used for genotyping by sequencing [GBS; [41]]. The genomic DNA was extracted from each genotype following the modified Cetyl Trimethyl Ammonium Bromide (CTAB) method suggested by Saghai-Maroof et al. [42]. Genomic services were outsourced from NxGenBio Life Sciences Private Limited. An optimized GBS library was prepared using high-quality DNA following the protocol given by Elshire et al. [41], and raw sequence data were generated using Illumina True Seq sequencing. Sequence reads were processed and aligned to the reference genome of B. juncea cv. Varuna UDSC Var 1.1 [43] using bwa version 0.7.17 [44]. Polymorphic SNPs were identified using the TASSEL-GBS pipeline [45], which further filtered for minor allele frequency of > 0.05 as well as missing genotypic frequency of < 0.01. The linkage map was constructed using polymorphic non-redundant binned SNPs with the least missing data using IciMapping 4.2 tool ([http://www.isbreeding.net; [46]).
Mapping of QTLs and candidate gene analysis
Agro-morphological traits recorded on ILs across the environments and molecular data were used to map Quantitative Trait Loci (QTL) in different environments using BIP and MET function in IciMapping 4.2 software [46]. To increase the authenticity and reliability of the detected QTLs, Logarithm of the Odds (LOD) was kept at ≥ 3. QTLs were named as per the previously described method [47, 48]. The name of the identified QTL begins with the initial "q" followed by the short form of the trait and the linkage group (A1 to A10 and B1 to B8). When any linkage group contains more than one QTL, the detected QTLs are numbered in order of their physical locations. Candidate genes within the region of identified QTLs were predicted following the method suggested by Wang et al. [49]. The identification of candidate genes was based on the annotation of the Arabidopsis thaliana genome and the physical positions of SNPs in the B. juncea reference genome [43].
Results
Phenotypic evaluation of introgression lines
Wide range of phenotypic variations were observed for all the 15 agro-morphological traits in the ILs evaluated under both rainfed and irrigated conditions. All the traits were normally distributed (Fig. 3), indicating the suitability of ILs for QTL mapping. The notched boxplots obtained between rainfed and irrigated conditions for each trait were compared following the Wilcoxon test, with the corresponding level of significance depicted through P values (Fig. 3). Pooled analysis over the environments indicated a significant mean difference between rainfed and irrigated conditions for studied traits, viz., plant height (cm), primary branches/plant, secondary branches/plant, total siliquae on main shoot, biological yield per plant (g), harvest index, seed yield per plot (g), 1,000-seed weight (g), days to 50% flowering and days to maturity. The mean values of harvest index and 1,000 seed weight (g) were found to be significantly higher under rainfed than in irrigated conditions (Fig. 3).
The data recorded on all the traits under three different environments were pooled, and analysis revealed a significant effect of genotypes (P ≤ 0.01) under both rainfed and irrigated conditions (Table 1). This indicates the presence of sufficient genetic variability in the ILs for all the studied traits across environments. Significant G × E interactions were also observed for all the traits studied (Table 1). A wide range of variation, along with transgressive segregants, was recovered for almost all the traits recorded under rainfed and irrigated conditions (Table 1). The broad-sense heritability across the environments ranged from 0.19 to 0.95 and from 0.28 to 0.91 under rainfed and irrigated conditions, respectively. Under rainfed conditions, all the traits exhibited high heritability except for secondary branches per plant, seed yield per plant and biological yield per plant, which exhibited low to moderate heritability (Table 1).
Water use efficiency (WUE) and drought tolerance indices
Water use efficiency (WUE) and drought tolerance indices, viz., drought susceptibility index (DSI), drought tolerance index (DTI), tolerance index (TOL) and mean relative performance (MRP) for seed yield (kg/ha) were estimated under three different environments. WUE among the parental lines was observed to vary from 1.23 to 2.88 kg m−3, 1.41 to 5.23 kg m−3 and 5.06 to 5.70 kg m−3 in RE1, RE2 and RE3 environments, respectively, under a moisture deficit stress condition. On the other hand, WUE among ILs ranged from 0.72 to 4.09 kg m−3, 1.83 to 6.19 kg m−3 and 1.77 to 10.56 kg m−3 in RE1, RE2 and RE3 environments, respectively. The highest and lowest values for this trait were recorded for IL94 and IL154 in RE1, IL160 and IL135 in RE2 and IL160 and IL125 in RE3 environments, respectively (Supplementary Table). In all the environments, ILs had a higher WUE than their B. juncea parent (Fig. 4). Results also revealed that in all environments, the mean WUE of parents and ILs were higher under rainfed conditions than under irrigated conditions (Fig. 4).
Important seed yield contributing traits such as siliqua length (cm), seed yield per plant (g), biological yield per plant (g), harvest index, seed yield per plot (g), 1,000 seed weight (g) and oil content (%) exhibited higher DTI and MRP for ILs than their parents (Fig. 5). DTI for seed yield varied from 0.17 to 2.73, 0.18 to 2.21 and 0.18 to 2.12 in ILs, whereas it varied from 0.12 to 1.46, 0.10 to 1.40 and 0.35 to 1.25 in their parents in E1, E2 and E3 environments, respectively (Supplementary Table). Similarly, MRP varied from 0.82 to 2.68, 0.63 to 2.37 and 1.31 to 2.36 in parents, whereas it ranges from 0.92 to 3.68, 0.85 to 2.98 and 0.96 to 3.09 in ILs in E1, E2 and E3 environments, respectively (Supplementary Table). ILs, viz., IL104, IL106, IL119, IL121, IL126, IL128, IL131 and IL135 have a DTI value of > 1 and a MRP value of > 2, indicating their superiority under both rainfed and irrigated conditions across all the environments (Supplementary Table).
Construction of linkage map
The linkage map was constructed using 5,165 genome-wide SNPs obtained from GBS and grouped into 838 unique bins (Table 2). The linkage map spanned over a length of 1,671.87 cM with an average marker interval of 2.00 cM and an average marker density of 3.09 cM. The number of genetic bins on each chromosome varied from 25 (10A) to 65 (9A) and markers from 159 (10A) to 438 (9A) per chromosome. The shortest chromosome was 4A, which carried 204 markers with a genetic length of 5.21 cM, an average marker interval of 0.12 cM, and a marker density of 39.16 per marker. The 3B chromosome was the longest, with 277 markers, 245.27 cM genetic length and an average marker interval of 6.45 cM. Marker density was found to be lowest in chromosome 3B (1.13 cM per marker) and highest in 4A (39.16 cM per marker), followed by 5B (10.25 cM per marker). The A and B genomes contain 2932 and 2233 markers with a total length of 736.17 cM and 935.70 cM, respectively. Among the two genomes of B. juncea, the A genome had a higher marker density (3.98 cM per marker and 56.77% of total markers) than the B genome (2.39 cM per marker and 43.23% of total markers) (Table 2).
Mapping of QTLs for agro-morphological traits under drought conditions
In three different environments, namely 2018–19 (Delhi; E1), 2020–21 (Delhi; E2), and 2020–21 crop season (Bharatpur; E3), created under moisture deficiency stress conditions, a total of 24 additive QTLs linked with seed yield and its contributing attributes were discovered (Table 3). These QTLs were distributed over eleven chromosomes (3A, 7A, 8A, 9A, 10A, 1B, 2B, 3B, 5B, 6B and 8B) and explained about 5.49 to 20.65% of the total phenotypic variation (PV). Among 24 QTLs, 14 (58%) loci carried alleles from the B. carinata parent (BC 4) for increasing phenotypic values and to withstand moisture deficit stress. While the remaining 10 (42%) loci carried alleles from the B. juncea parent (Table 3).
Three additive effect QTLs were identified, one in each environment for plant height, explaining PV from 11.55% to 16.38% and distributed over three chromosomes (9A, 2B and 3B). For secondary branches per plant, three QTLs were detected on three different chromosomes (7A, 10A and 3B) and explained 5.49 to 15.66% PV. A QTL regulating main shoot length under moisture deficit stress was identified to exhibit high PV (16.17%). For siliqua length two major QTLs, viz., qSL.2B.1 on chromosome 2B and qSL.8A.1 on chromosome 8A were detected to explain PV of 12.96 and 20.65%, respectively. The QTL qSL.2B.1, flanked by markers 2B_65115276 and 2B_64577224, was found to regulate siliqua length in two of the three environments. All the QTLs associated with siliqua length observed positive additive effects, indicating the contribution of B. carinata parent BC 4 to increase siliqua length under moisture deficit stress conditions (Table 3).
For seeds per siliqua, a major QTL (qSS.9A.1), which is flanked by markers 9A_38513927 and 9A_5962669, was detected in two different environments and explained PV up to 18.21%. Two QTLs on chromosome 3B were found to be associated with seed yield per plant under moisture deficit stress conditions, and these explained PV ranging from 14.15 to 18.34% (Table 3). Three QTLs were discovered to be associated with biological yield per plant; they were present on chromosomes 3A and 3B and explained PV ranging from 5.53 to 11.91% under moisture deficit stress conditions. A QTL flanked by markers 3B_9101086 and 3B_22690875 (qBY.3B.2) was detected to explain high PV in two different environments, whereas, another QTL on chromosome 8A explained 11.49% variation for oil content. Three QTLs, viz., qTSW.3B.1, qTSW.3B.2 and qTSW.6B.1, were found to regulate 1,000-seed weight (g). Two of these were found to be located on chromosome 3B and the other on chromosome 6B, explaining PV from 11.78 to 14.58%, respectively. Another QTL for 1000-seed weight (qTSW.3B.2) flanked by 3B_16490510 and 3B_15005559 markers was found to explain 14.40 to 14.58% PV in two environments. For days to 50% flowering, three QTLs were detected, one in each environment, explaining PV from 7.42 to 17.84% and were located on three different chromosomes (1B, 3B and 5B). QTLs viz., qDM.8B.1 and qDM.8B.1 were reported to be associated with shorter maturity duration under moisture deficit stress in E2 and E3 environments, explaining 12.41 and 16.35% of PVs, respectively (Table 3).
Several QTLs with a wide range of QTL × environment interactions were detected across the environments under moisture deficit stress conditions (Table 4). For plant height, QTL detected on chromosome 2B (qPH.2B.1) had a high Logarithm of Odds (LOD) score for additive effect [LOD (A)] and a small LOD score for additive × environment effect [LOD (A × E)], which indicated weak QTL-environment interactions. QTL, qPH.3B.1 exhibited a high LOD (A × E) and low LOD (A) values, implying significant interactions between the QTL and the environment. QTLs for primary and secondary branches per plant showed very high QTL × environment interactions, which is evident by the small value of LOD (A) and high LOD (A × E). In the case of main shoot length, out of three QTLs detected, one on 3B (qMSL.3B.2) exhibited very strong QTL × environment interactions. All the QTLs identified for siliqua length, seeds per siliqua, 1000-seed weight and oil content displayed high values for LOD (A) and small LOD (A × E) values. The QTLs for seed yield per plant and biological yield per plant exhibited medium-to-high QTL × environment interactions. For days to 50% flowering, all identified QTLs showed small QTL × environment interactions. Only one QTL present on chromosome 8B (qDM.8B.1), explaining variation for days to maturity under moisture deficit stress conditions, had strong QTL × environment interactions, while the rest were responsible for controlling variation for this trait only in one or the other environments (Table 4).
Mapping of QTLs for water use efficiency and drought tolerance indices
A total of five additive QTLs associated with different WUE and drought tolerance indices were identified in three environments under moisture deficit stress conditions (Table 5). These were located on four different chromosomes (7A, 2B, 7B and 8B) and explained PV ranging from 5.42 to 24.18%. Out of these five QTLs, three had a positive additive effect, which indicated that these loci carry alleles from the B. carinata parent (BC 4) and are responsible for improving drought tolerance in the ILs. The rest of the two QTLs identified in these ILs were inherited from B. juncea cultivar DRMRIJ 31.
A QTL (qWUE.7B.1) with flanking markers 7B_30548045 and 7B_870706 was found to improve WUE under moisture deficit stress conditions and explained up-to 13.40% PV. Two QTLs, viz., qDTI.2B.1 and qDTI.7A.1, were reported to control the DTI and explained high PV (12.51 and 24.18%, respectively). Similarly, two QTLs (qMRP.2B.1 and qMRP.8B.1) were reported to explain 12.63 and 5.42% PV for MRP, respectively (Table 5). Another QTL (qDTI.2B.1) located on the 2B chromosome, flanked by markers 2B_3291748 and 2B_1636089, was identified to explain more than 10% PV for DTI and MRP traits (Table 5).
Significant QTL × environment interactions were also observed for QTLs explaining WUE and drought tolerance indices across the environments (Table 6). A total of 16 QTLs were reported to influence WUE under moisture deficit stress and drought tolerance indices. Six QTLs responsible for improving WUE in this study were identified on three different chromosomes, and three of them (qWUE.7A.1, qWUE.7A.2 and qWUE.3B.1) had higher LOD (A) than LOD (A × E), indicating small QTL × environment interactions. The other three QTLs (qWUE.7B.1, qWUE.7B.2 and qWUE.7B.3), on the other hand, have significant QTL × environment interactions. All five QTLs responsible for improved DTI exhibited higher LOD (A) than LOD (A × E), suggesting smaller QTL × environment interactions. Four QTLs associated with higher MRP, namely qMRP.7A.1, qMRP.1B.1, qMRP.2B.1 and qMRP.3B.1, had higher LOD (A) than LOD (A × E), and had small QTL × environment interactions. Higher QTL × environment interaction was also observed for qMRP.8B.1 QTL (Table 6).
Identification of co-localized QTLs and candidate genes in the QTL regions
Analysis was performed for the identification of two or more co-localized QTLs, at the same position in the genome, governing agro-morphological traits, WUE and drought tolerance indices under moisture deficit stress conditions (Table 7; Fig. 6). Four QTLs, viz., qSS.7A.1, qSL.7A.1, qWUE.7A.2 and qDTI.7A.1 were found to be located together in the small confidence interval on chromosome 7A. Two QTLs (qSL.9A.1 and qSS.9A.1) controlling siliqua length and seeds per siliqua, respectively, were found co-localized on the 9A chromosome. Two QTLs (qDTI.1B.1 and qMRP.1B.1) on 1B and two (qDF.2B.1 and qDM.2B.1) on the 2B chromosome were detected at the same position in the genome (Fig. 6). The hotspot carrying eight QTLs, viz., qBY.3B.2, qMSL.3B.1, qSY.3B.1, qSL.3B.1, qDF.3B.1, qPH.3B.1, qDM.3B.1 and qWUE.3B.1 was identified on chromosome 3B in a small confidence interval. Furthermore, QTLs qDF.7B.1 and qDM.7B.1 were also reported to be located at the same position on chromosome 7B (Table 7; Fig. 6). To make any B. juncea line tolerant to moisture deficit stress, QTL hotspots identified in the present study can be further transferred from ILs through MAS.
The search for candidate genes in the different QTL regions through in silico analysis on B. juncea chromosomes identified a total of seventeen genes in Indian mustard with their orthologs in Arabidopsis thaliana (Table 7). SOS2like and NPR1 genes were identified in the qMSL.6A.1 QTL region, while qSS.7A.1 and qSL.7A.1 carried the FAE1 KCS gene. Three candidate genes, viz., HOT5, DNAJA1 and NIA1, were reported in the qWUE.7A.2 and qDTI.7A.1 QTL regions. However, QTL qMRP.7A.1 identified to govern MRP under moisture deficit stress conditions, harbored the BRI1 (brassinosteroid insensitive 1) gene in the same chromosomal region. The genes RF21 and ycf2 were identified to be co-localized in the genomic region of QTLs qSL.9A.1 and qSS.9A.1, whereas genes PAL and SAMS2 were detected in the qSB.3B.1, qBY.3B.2 and qMSL.3B.1 QTL regions. Furthermore, QTL qSB.3B.1 harbors the WRKY-33 gene on the 3B chromosome (Table 7).
Similarly, the QTL responsible for higher DTI (qDTI.6B.1) carried candidate genes SOS2like, ERFI, BRI and WRR1 in the same genomic region, whereas the MAPK3 gene was co-localized with QTLs qSL.3B.2 and qTSW.3B.1. The gene SUS (sucrose synthase 09) was found to be located within qWUE.7B.1, a QTL identified to improve WUE under moisture deficit stress conditions (Table 7). The association among identified QTLs and already reported candidate genes that are directly or indirectly involved in the pathways associated with drought tolerance in plants was established in this study. Further, the QTLs carrying putative genes in the same genomic region may even be considered as candidate loci in the future studies.
Discussion
Interspecific hybridization is a prominent method for generating selectable genetic variation in cultivated Brassica species [67]. B. carinata being tolerant to various abiotic stresses, including drought, is most suitable for broadening the genetic base of other Brassicas [23, 33]. Introgression lines, carrying genomic segments from an alien species into the genetic background of cultivated genotype(s), are highly desirable for the identification and mapping of novel genes/QTLs. Simultaneously, such efforts provide valuable insights into the genetic basis of drought tolerance. The present study was, thus, aimed to develop B. carinata-derived B. juncea ILs possessing drought tolerance/improved WUE under water stress conditions and elucidating responsive QTLs/genes.
The existence of significant differences (P < 0.001) among the ILs for all agro-morphological traits under both rainfed and irrigated conditions indicated the wide range of variability created for these traits, as reported in different environments (Table 1). Furthermore, the pooled analysis revealed significant genotype-by-environment (G × E) interactions for all the traits studied, suggesting the influence of the environment on the expression of these traits. Notably, under both rainfed and irrigated conditions, a wide range of variation along with the recovery of transgressive segregants was observed for nearly all the agro-morphological traits (Table 1), indicating the distribution of desirable alleles in both parents.
Water use efficiency and drought tolerance indices
In general, higher values of WUE were observed under rainfed conditions than in irrigated conditions and ILs had higher WUE than their parents in both test environments (Fig. 4). Associations among the drought tolerance indices suggested the usefulness of DTI and MRP in the indirect selection of genotypes suitable for drought tolerance in B. carinata-derived B. juncea lines [32]. Therefore, in the present study, DTI and MRP were calculated from the average values of ILs and parents. Several ILs exhibited high DTI (> 1) and MRP (> 2) values, indicating their superiority in both rainfed and irrigated conditions across all environments (Supplementary Table). Furthermore, seed yield contributing traits such as siliqua length, seed yield per plant, biological yield per plant, harvest index, seed yield per plot, 1000-seed weight and oil content observed higher DTI and MRP for ILs than their respective parents (Fig. 5), implying that improved WUE in these ILs might be due to the complementation of beneficial alleles [32]. These ILs have the potential to become cultivars after large-scale testing and/or be involved in hybridization to develop varieties for drought-prone areas [32, 68].
Construction of linkage map and mapping of QTLs governing drought tolerance
The linkage map is a prerequisite for the precise mapping of QTLs associated with drought tolerance [69]. In the present study, it was constructed using 5,165 GBS-based SNP markers widely distributed on 18 chromosomes. It covered a total length of 1671.87 cM with an average interval of 2.0 cM between marker loci (Table 2), indicating the availability of a large number of genome-wide SNP markers for the identification of genomic regions governing the traits under study. The GBS approach has earlier been successfully used for the construction of linkage maps [43, 70,71,72] and mapping of QTLs in B. juncea and its derived lines [29, 30, 73, 74].
A total of 29 QTLs were identified in the present study. Out of 29, 5 were responsible for WUE and drought tolerance, while the remaining 24 were responsible for additive effect explaining variation in various agro-morphological traits (Table 3 & 5). Results revealed that B. carinata has contributed about 17 QTLs (58.6% of total QTLs detected) which regulates 10 of the studied traits, suggesting that alleles from the B. carinata largely contributed to improving the drought tolerance in these ILs. DRMRIJ 31 contributed 12 (41.4%) of the 29 QTLs identified in different environments. This implies that positive loci governing agro-morphological traits under moisture deficit stress conditions were majorly contributed by B. carinata parent in the ILs.
Out of 17 QTLs identified in ILs, 10 (58.8% of total QTLs detected) were contributed by the B genome of B. carinata, while the rest seven (41.2%) were contributed by the A genome. It is likely that pairing and recombination between the B genomes of both species, as well as the higher frequency of recombination between the A and C genomes, resulted in the introgression of the genomic regions between B. carinata and B. juncea [75, 76]. The present study reports a larger number of QTLs for drought tolerance and improved WUE in the B genome. As also revealed in the previous study, a larger number of introgressed segments were observed in the B genome than in the A genome of the ILs [33]. This might be due to the conservation of novel alleles for drought tolerance in B. carinata [77], and their elimination in B. juncea during the course of evolution. It is expected as a consequence of nucleo-cytoplasmic interactions leading to considerable changes in the B genome of B. juncea, whereas it remains intact in B. carinata [78, 79].
In the present study, 23 major QTLs expressing high PV (> 10%) for agro-morphological traits, WUE and drought tolerance indices under moisture deficit stress conditions were identified (Table 3 & 5). A major QTL expressing PV upto 20.65% and flanked by markers 2B_65115276 and 2B_64577224 was identified to regulate the siliqua length (qSL.2B.1) trait in two different environments. A QTL governing seed per siliqua (qSS.9A.1), flanked by markers 9A_38513927 and 9A_5962669, explained PV upto 18.21% was detected in two different environments (Table 3). Both QTLs exhibited a positive additive effect, indicating that QTLs derived from the B. carinata parent (BC 4) are responsible for increasing siliqua length and seeds per siliqua in ILs under moisture deficit stress. Most of the QTLs reported to explain seed yield contributing traits showing positive additive effects, in the present study, were derived from B. carinata (BC 4), thus highlighting the usefulness of B. carinata in stabilizing productivity traits under drought stress conditions (Table 3). A few QTLs identified in current study for seed yield-contributing traits were found on the same chromosomes as in prior studies by Dhaka et al. [80], Rout et al., [72], and Aakanksha et al. [81]. Notably, these QTLs were associated with imparting drought tolerance and improving WUE under moisture deficit stress conditions in B. carinata derived introgression lines.
QTL × environment interactions for morpho-physiological traits in ILs
The ability of genotypes to produce different phenotypes in a wide range of environments is mainly due to phenotypic plasticity arising from the interaction of QTLs with environments [82]. Therefore, understanding QTLs × environment interaction will help to select stable genotypes across environments, which further improves crop productivity [83, 84]. Previous studies also reported significant QTLs × environment interaction for yield-contributing traits in several crops, including mustard, and demonstrated the varied range of QTL expression with a change in environment [69, 80].
The present study also demonstrated the QTLs × environment interactions for morpho-physiological traits in B. carinata derived ILs under moisture deficit stress conditions (Table 4 & 6). Fifty QTLs were identified to govern agro-morphological traits (Table 4) and 16 QTLs for WUE/drought tolerance indices (Table 6) collectively demonstrated substantial QTL × environment interactions, indicating the existence of selectable genetic variations for phenotypic plasticity among these ILs. In addition, all the QTLs governing plant height, secondary branches per plant, main shoot length, seed yield/plant, oil content, days to 50% flowering and days to maturity, as well as a few QTLs for siliqua length, seeds/siliqua, biological yield/plant and 1000-seed weight were detected in specific environments, indicating their varied range of expression in different environments [69, 85]. The QTLs detected for siliqua length, seeds per siliqua, 1,000-seed weight and oil content observed small LOD (A x E) and high values for LOD (A), indicating their stable response to different environments (Table 4). This finding is consistent with previous studies of Singh et al. [86]; Singh et al. [87]; Singh et al. [88]; Binod et al. [89], which demonstrate the preponderance of additive gene action governing these traits.
The majority of QTLs for WUE and drought tolerance indices had low QTL × environment interactions (Table 6), revealing their stability across the environments. However, QTLs associated with WUE and MRP, namely qWUE.7B.1, qWUE.7B.3, qWUE.7B.2 and qMRP.8B.1, had the highest QTL × environment interactions, implying differential expression pattern in response to different environments. The effect of environment on QTL expression demonstrates that QTL × environment interactions is a key component of genetic variation which can play an important role in defining future mustard breeding programs. The QTLs identified in the present study can be used in a wide range of environments or any specific environment based on their degree of QTL × environment interactions.
Identification of co-localized QTLs and candidate genes in the QTL regions
The genomic regions containing multiple QTLs for different traits, also called QTL hotspots, enable simultaneous selection and accelerating the breeding progress through MAS [90]. The present study reported four co-localized QTLs, viz., qSS.7A.1, qSL.7A.1, qWUE.7A.2 and qDTI.7A.1 on chromosome 7A (Table 7; Fig. 6). QTLs for siliqua length (qSL.9A.1) and seeds per siliqua (qSS.9A.1), on the other hand, were reported at the same position with the same set of flanking markers, indicating a pleiotropic effect. QTLs, viz., qBY.3B.2, qMSL.3B.1, qSY.3B.1, qSL.3B.1, qDF.3B.1, qPH.3B.1, qDM.3B.1 and qWUE.3B.1 were observed to be co-localized on chromosome 3B, revealing a QTL hotspot that can be further exploited through MAS for improving trait values under moisture deficit stress conditions (Table 7).
Seventeen candidate genes, known to regulate various pathways related to biotic and abiotic stress tolerance, were identified within QTL regions through in silico analysis in this study (Table 7). The genomic region of qMSL.6A.1 contains SOS2like and NPR1 genes associated with salt tolerance and salicylic acid-mediated systemic acquired resistance (SAR) pathways, respectively [50, 51]. QTLs qSS.7A.1 and qSL.7A.1 carry the FAE1 KCS gene involved in erucic acid production [91, 92]. Three candidate genes, HOT5, DNAJA1, and NIA1, were identified in the qWUE.7A.2 and qDTI.7A.1 QTL regions, known for their role in inducing abiotic stress tolerance [53,54,55]. The qMRP.7A.1 QTL, associated with moisture deficit stress, harbours the BRI1 gene linked to drought stress tolerance [56]. RF21 and ycf2 genes were co-localized in the qSL.9A.1 and qSS.9A.1 QTL regions, respectively while WRKY-33 gene was found in the qSB.3B.1 QTL, associated with abiotic stress responses [58]. PAL and SAMS2 genes were detected in qBY.3B.2 and qMSL.3B.1 QTL regions, respectively, known for their involvement in drought and salt tolerance [59,60,61,62, 93, 94]. Lastly, the SOS2like gene involved in plant response to abiotic stresses found to be co-localized with qSY.3B.2 and qSL.3B.1 QTLs for seed yield and siliqua length respectively [50].
Similarly, the genomic region of co-localized QTLs qSL.3B.2 and qTSW.3B.1 encompasses the MAPK3 gene, which is reported to impart drought tolerance in tomato [63]. The drought tolerance index QTL qDTI.6B.1 contained four genes: SOS2 [50]; ERF [64]; BRI1 [56]; and WRR1 [65]. Three of these four genes, viz., SOS2, ERF and BRI1, were reported to be involved in plant responses to abiotic stresses [50, 56, 64]. The SUS gene encoding sucrose synthase 09 [66] was identified within the genomic region of QTL qWUE.7B.1, and identified to be responsible for improving water use efficiency under moisture deficit stress conditions (Table 7). Thus, QTLs detected in multiple environments and carrying already known genes associated with drought tolerance can also be used to identify drought-responsive key candidate genes and their markers following the transcriptomic approach by mapping the transcripts over the reference sequence as reported in crops including Brassicas [95, 96]. Furthermore, the co-localized QTLs conferring drought tolerance and/or improved WUE can be further subjected to fine mapping and validation for their wider applications in mustard breeding through marker-assisted selection. The material and information generated from this study are expected to have long-term implications for the development of drought-tolerant mustard varieties.
Conclusion
The study demonstrated the potential of B. carinata-derived B. juncea ILs for improving drought tolerance and WUE in B. juncea cultivar DRMRIJ 31 and provided insights into the genetic basis of yield-contributing traits under moisture deficit stress conditions. Significant differences for yield-contributing agro-morphological traits were observed among the ILs, evaluated under both rainfed and irrigated conditions, indicating the creation of sufficient genetic variability through interspecific hybridization. These ILs exhibited higher WUE and drought tolerance indices than their parents; moreover, some of these lines also showed higher phenotypic values under both rainfed and irrigated conditions. The genotyping by sequencing (GBS)-based linkage map helped identify the genomic regions associated with drought tolerance and improved WUE. Twenty-nine QTLs for seed yield-related traits were discovered, with a significant contribution from the B genome of B. carinata. Significant QTL-environment interactions were also observed, indicating the influence of the environment on the expression of the studied traits. QTL hotspots for various traits were identified in the present study, which will further provide opportunities for marker-assisted selection. In silico analysis identified 17 candidate genes involved in stress tolerance pathways within the QTL regions. The material and information generated from the present study have demonstrated the usefulness of interspecific hybridization among Brassicas in creating selectable novel genetic variability and paved the way for the development of high-yielding varieties with better WUE and improved productivity in water-scarce regions.
Availability of data and materials
All data generated or analyzed during this study are included in the paper or supplementary information.
References
Kogan F, Guo W, Yang W. Drought and food security prediction from NOAA new generation of operational satellites. Geomatics Nat Hazards Risk. 2019;10(1):651–66.
Hassan T, Hamid M, Wani SA, Malik AH, Waza SA, Khuroo AA. Substantial shifts in flowering phenology of Sternbergia vernalis in the Himalaya: Supplementing decadal field records with historical and experimental evidences. Sci Total Environ. 2021;795: 148811.
Sehgal A, Sita K, Bhandari K, Kumar S, Kumar J, Vara Prasad PV, Nayyar H. Influence of drought and heat stress, applied independently or in combination during seed development, on qualitative and quantitative aspects of seeds of lentil (Lens culinaris Medikus) genotypes, differing in drought sensitivity. Plant Cell Environ. 2019;42(1):198–211.
Rasheed A, Hassan M, Aamer M, Batool M, Sheng F, Ziming W, Huijie L. A critical review on the improvement of drought stress tolerance in rice (Oryza sativa L.). Not Bot Horti Agrobot Cluj-Napoca. 2020;48(4):1756.
Krishnan R, Sanjay J, Gnanaseelan C, Mujumdar M, Kulkarni A, Chakraborty S. Assessment of climate change over the Indian region: a report of the ministry of earth sciences (MOES), government of India. Singapore: Springer Nature; 2020. https://doi.org/10.1007/978-981-15-4327-2.
Anonymous. USDA foreign agricultural service, Global agricultural information network Circular Series. Global Market Analysis International Production Assessment Division Washington. 2023. WAP 5–23. https://ipad.fas.usda.gov/.
Chandra K, Pandey A, Mishra SB. Interrelationship and cause-effect of morphological traits with grain yield and oil content among Indian mustard (Brassica juncea L. Czern. & Coss.) Genotypes under non-irrigated and irrigated condition. J Exp Agri Int. 2018;28:1–14.
DRMR. VISION 2030p30. Directorate of Rapeseed-Mustard Research, Bharatpur 321303 Rajasthan. 2011. Available at https://www.drmr.res.in/pub/.
RR Valiyaveettil GK Jha B Kathayat 2023 Pushing for self-sufficiency in edible oils in India in the aftermath of recent global events Natl Acad Sci Lett https://doi.org/10.1007/s40009-023-01294-z.
Chauhan JS, Singh KH, Singh VV, Kumar S. Hundred years of rapeseed-mustard breeding in India: accomplishments and future strategies. Indian J Agric Sci. 2011;81(12):1093–109.
Chauhan JS, Tyagi MK, Kumar A, Nashaat NI, Singh M, Singh NB, et al. Drought effects on yield and its components in Indian mustard (Brassica juncea L). Plant Breed. 2007;126(4):399–402.
Sharma P, Sharma HO, Rai PK. Strategies and technologies for enhancing rapeseed-mustard production and farmer income. Indian Farming. 2018;68(01):44–8.
Malik RS. Prospects for Brassica carinata as an oilseed crop in India. Exp Agri. 1990;26(1):125–9.
Warwick SI. Wild species in the tribe Brassicaceae (Cruciferae) as sources of agronomic traits. In: Guide to the wild germplasm of Brassica and allied crops. Technical Bulletin 17E, 1–19. Center for Land and Biological Resources Research Branch, Agriculture Canada, Ottawa, Ontario, Canada. 1993.
Kasuga M, Liu Q, Miura S, Yamaguchi-Shinozaki K, Shinozaki K. Improving plant drought, salt, and freezing tolerance by gene transfer of a single stress-inducible transcription factor. Nat Biotechnol. 1999;17(3):287–91.
Dalal M, Tayal D, Chinnusamy V, Bansal KC. Abiotic stress and ABA-inducible Group 4 LEA from Brassica napus plays a key role in salt and drought tolerance. J Biotechnol. 2009;139(2):137–45.
Seo YJ, Park JB, Cho YJ, Jung C, Seo HS, Park SK, Nahm BH, Song JT. Overexpression of the ethylene-responsive factor gene BrERF4 from Brassica rapa increases tolerance to salt and drought in Arabidopsis plants. Mol Cells. 2010;30(3):271–7.
Fletcher RS, Mullen JL, Heiliger A, McKay JK. QTL analysis of root morphology, flowering time, and yield reveals trade-offs in response to drought in Brassica napus. J Exp Bot. 2015;66(1):245–56.
Prakash S. Non-homologous meiotic pairing in the A and B genomes of Brassica: its breeding significance in the production of variable amphidiploids. Genet Res Cambridge. 1973;21:133–7.
Inomata N. Wide hybridization and meiotic pairing. In: Kalia HR, Gupta SK, editors. Recent Advances in Oilseed Brassicas Kalyani Publisher, Ludhiana, India. 1997. p. 53–76.
Zhang X, Lu G, Long W, Zou X, Li F, Nishio T. Recent progress in drought and salt tolerance studies in Brassica crops. Breed Sci. 2014;64(1):60–73.
Raman R, Qiu Y, Coombes N, Song J, Kilian A, Raman H. Molecular diversity analysis and genetic mapping of pod shatter resistance loci in Brassica carinata L. Front Plant Sci. 2017;30:1765.
Thakur AK, Singh KH, Sharma D, Parmar N, Nanjundan J. Breeding and genomics interventions in Ethiopian mustard (Brassica carinata A Braun) improvement-A mini review. S Afri J Bot. 2019;125:457–65.
Lohani N, Golicz AA, Singh MB, Bhalla PL. Genome-wide analysis of the Hsf gene family in Brassica oleracea and a comparative analysis of the Hsf gene family in B. oleracea, B. rapa and B. napus. Funct Integr Genom. 2019;19:515–31.
Mafakheri, M, Kordrostami M. Newly Revealed Promising Gene Pools of Neglected Brassica Species to Improve Stress-Tolerant Crops. In: Hasanuzzaman M, editor. The Plant Family Brassicaceae. Singapore: Springer; 2020. https://doi.org/10.1007/978-981-15-6345-4_4.
Mwadzingeni L, Shimelis H, Dube E, Laing MD, Tsilo TJ. Breeding wheat for drought tolerance: progress and technologies. J Integr Agric. 2014;15(5):935–43.
Vieira MLC, Santini L, Diniz AL, Munhoz CDF. Microsatellite markers: what they mean and why they are so useful. Genet Mol Biol. 2016;39:312–28.
Rana K, Atri C, Akhatar J, Kaur R, Goyal A, Singh MP, Banga SS. Detection of first marker trait associations for resistance against Sclerotinia sclerotiorum in Brassica juncea-Erucastrum cardaminoides introgression lines. Front Plant Sci. 2019;10:1015.
Rana K, Atri C, Gupta M, Akhatar J, Sandhu PS, Kumar N, Jaswal R, et al. Mapping resistance responses to Sclerotinia infestation in introgression lines of Brassica juncea carrying genomic segments from wild Brassicaceae Bfruticulosa. Sci Rep. 2017;7(1):1–12.
Atri C, Akhatar J, Gupta M, Gupta N, et al. Molecular and genetic analysis of defensive responses of Brassica juncea-B. fruticulosa introgression lines to Sclerotinia infection. Sci Rep. 2019;9(1):1–12.
Li F, Chen B, Xu K, Wu J, Song W, Bancroft I, Wu X. Genome-wide association study dissects the genetic architecture of seed weight and seed quality in rapeseed (Brassica napus L). DNA Res. 2014;21(4):355–67.
Limbalkar OM, Singh R, Kumar P, Nanjundan J, Parihar CM, Vasisth P, Yadava DK, Chinnusamy V, Singh N. Deployment of Brassica carinata A. Braun derived Brassica juncea (L.) Czern lines for improving heterosis and water use efficiency under water deficit stress conditions. Front Plant Sci. 2021;12:765645.
Vasisth P, Singh N, Limbalkar OM, Sharma M, Dhanasekaran G, Meena ML, Singh R, et al. Introgression of heterotic genomic segments from Brassica carinata into Brassica juncea for enhancing productivity. Plants. 2023. https://doi.org/10.3390/plants12081677.
Rodriquez F, Alvarado G, Pacheco A, Burgueno J. ACBD-R, Augmented Complete Block Design with R for Windows Version 40. CIMMYT Res Data Softw Repos Netw. 2018. Corpus ID: 65413032.
Lozada DN, Carter AH. Accuracy of single and multi-trait genomic prediction models for grain yield in US Pacific Northwest winter wheat. Crop Breed Genet Genom. 2019;1:e190012. https://doi.org/10.20900/cbgg20190012.
Lush JL. Intra-sire correlations or regressions of offspring on dam as a method of estimating heritability of characteristics. J Anim Sci. 1940;1940(1):293–301. https://doi.org/10.2527/jas1940.19401293x.
IRRI. Plant Breeding Tools (PBTools) Version 1.4. Biometrics and Breeding Informatics, PBGB Division, International Rice Research Institute, Los Banos, Philippines. 2014.
Fischer RA, Maurer R. Drought resistance in spring wheat cultivars. I. Grain yield responses. Aus J Agric Re. 1978;29(5):897–912.
Fernandez GC. Effective selection criteria for assessing plant stress tolerance. In: Proceeding of the International Symposium on Adaptation of Vegetables and other Food Crops in Temperature and Water Stress Aug13–16 Shanhua Taiwan. 1992;257–70. https://doi.org/10.22001/wvc.72511.
Rosielle AA, Hamblin J. Theoretical aspects of selection for yield in stress and non-stress environment. Crop Sci. 1981;21(6):943–6.
Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, Mitchell SE. A robust, simple Genotyping-by-Sequencing (GBS) approach for high diversity species. PLoS ONE. 2011;6(5):e19379.
Saghai-Maroof MA, Soliman KM, Jorgensen RA, Allard R. Ribosomal DNA spacer-length polymorphisms in barley: mendelian inheritance, chromosomal location, and population dynamics. Proc Natl Acad Sci. 1984;81(24):8014–8.
Paritosh K, Yadava SK, Singh P, Bhayana L, Mukhopadhyay A, Gupta V, Bisht NC, et al. A chromosome-scale assembly of allotetraploid Brassica juncea (AABB) elucidates comparative architecture of the A and B genomes. Plant Biotechnol J. 2021;19(3):602–14.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinform. 2009;25(14):1754–60.
Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, Buckler ES. TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline. PLoS ONE. 2014;9(2): e90346.
Meng L, Li H, Zhang L, Wang J. QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015;3(3):269–83.
McCouch SR. Report on QTL nomenclature. Rice Genet News. 1997;14:11–3.
Zhao X, Yu K, Pang C, Wu X, Shi R, Sun C, Wang X, et al. QTL analysis of five silique-related traits in Brassica napus L across multiple environments. Front Plant Sci. 2021;12:766271–766271.
Wang H, Zaman QU, Huang W, Mei D, Liu J, Wang W, Hu Q. QTL and candidate gene identification for silique length based on high-dense genetic map in Brassica napus L. Front Plant Sci. 2019;10:1579.
Liu J, Ishitani M, Halfter U, Kim CS, Zhu JK. The Arabidopsis thaliana SOS2 gene encodes a protein kinase that is required for salt tolerance. Proc Natl Acad Sci. 2000;97(7):3730–4. https://doi.org/10.1073/pnas.97.7.3730.
Backer R, Naidoo S, Van den Berg N. The Nonexpressor of Pathogenesis-Related Genes 1 (NPR1) and related family: mechanistic insights in plant disease resistance. Front Plant Sci. 2019;10:102.
Ghanevati M, Jaworski JG. Active-site residues of a plant membrane-bound fatty acid elongase β-ketoacyl-CoA synthase, FAE1 KCS. Biochim Biophys Acta Mol Cell Biol Lipids. 2001;1530(1):77–85.
Leterrier M, Chaki M, Airaki M, Valderrama R, Palma JM, Barroso JB, Corpas FJ. Function of S-nitrosoglutathione reductase (GSNOR) in plant development and under biotic/abiotic stress. Plant Signal Behav. 2011;6(6):789–93.
Kong F, Deng Y, Zhou B, Wang G, Wang Y, Meng Q. A chloroplast-targeted DnaJ protein contributes to maintenance of photosystem II under chilling stress. J Exp Bot. 2014;65(1):143–58.
Zhao C, Cai S, Wang Y, Chen ZH. Loss of nitrate reductases NIA1 and NIA2 impairs stomatal closure by altering genes of core ABA signalling components in Arabidopsis. Plant Signal Behav. 2016;11(6):1456–69.
Fabregas N, Lozano-Elena F, Blasco-Escámez D, Tohge T, Martínez-Andújar C, Albacete A, et al. Overexpression of the vascular brassinosteroid receptor BRL3 confers drought resistance without penalizing plant growth. Nat Commun. 2018;9(1):1–13.
Nakai M. Reply: the revised model for chloroplast protein import. Plant Cell. 2020;32(3):543–6.
Li W, Pang S, Lu Z, Jin B. Function and mechanism of WRKY transcription factors in abiotic stress responses of plants. Plants. 2020;9(11):1515.
Rasool F, Uzair M, Naeem MK, Rehman N, Afroz A, Shah H, Khan MR. Phenylalanine ammonia-lyase (PAL) genes family in wheat (Triticum aestivum L): genome-wide characterization and expression profiling. Agronomy. 2021;11(12):2511.
Zhang X, Bao Z, Gong B, Shi Q. S-adenosylmethionine synthetase 1 confers drought and salt tolerance in transgenic tomato. Environ Exp Bot. 2020;179: 104226.
Hua Y, Zhang B, Cai HY, Bai X, Ji W, Zhu Y. Stress-inducible expression of GsSAMS2 enhances salt tolerance in transgenic Medicago sativa. African J Biotechnol. 2012;11(17):4030–8.
Sun F, Ma J, Wang P, Yang Y. Genome-wide identification of the SAMS gene family in upland cotton (Gossypium hirsutum L) and expression analysis in drought stress treatments. Genes. 2022;13(5):860.
Muhammad T, Zhang J, Ma Y, Li Y, Zhang F, Zhang Y, Liang Y. Overexpression of a mitogen-activated protein kinase SlMAPK3 positively regulates tomato tolerance to cadmium and drought stress. Mol. 2019;24(3):556.
Muller M, Munné-Bosch S. Ethylene response factors: a key regulatory hub in hormone and stress signaling. Plant Physiol. 2015;169(1):32–41.
Heena A, Padmaja KL, Kumar P, Nitika M, Tewari AK, Arundhati M, Deepak P. BjuWRR1, a CC-NB-LRR gene identified in Brassica juncea, confers resistance to white rust caused by Albugo candida. Theor Appl Genet. 2019;132(8):2223–36.
Li M, He Q, Huang Y, Luo Y, Zhang Y, Chen Q, Wang Y, et al. Sucrose synthase gene family in Brassica juncea: genomic organization, evolutionary comparisons, and expression regulation. Peer J. 2021;9: e10878.
Katche E, Quezada-Martinez D, Katche EI, Vasquez-Teuber P, Mason AS. Interspecific hybridization for Brassica crop improvement. Crop Breed Genet Genom. 2019;1(1).
Yohannes T, Abraha T, Kiambi D, Folkertsma R, Hash CT, Ngugi K, Mutitu E, et al. Marker-assisted introgression improves Striga resistance in an Eritrean farmer-preferred sorghum variety. Field Crops Res. 2015;173:22–9.
Hussain W, Baenziger PS, Belamkar V, Guttieri MJ, Venegas JP, Easterly A, Poland J. Genotyping-by-sequencing derived high-density linkage map and its application to QTL mapping of flag leaf traits in bread wheat. Sci Rep. 2017;7(1):1–15.
Yadava SK, Arumugam N, Mukhopadhyay A, Sodhi YS, Gupta V, Pental D, Pradhan AK. 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. 2012;125(7):1553–64.
Yang Y, Shen Y, Li S, Ge X, Li Z. High density linkage map construction and QTL detection for three silique-related traits in Orychophragmus violaceus derived Brassica napus population. Front Plant Sci. 2017;8:1512.
Rout K, Yadav BG, Yadava SK, Mukhopadhyay A, Gupta V, Pental D, Pradhan AK. QTL landscape for oil content in Brassica juncea: analysis in multiple bi-parental populations in high and “0” erucic background. Front Plant Sci. 2018;9:1448.
Wu J, Chen P, Zhao Q, Cai G, Hu Y, Xiang Y, Zhou Y. Co-location of QTL for Sclerotinia stem rot resistance and flowering time in Brassica napus. Crop J. 2019;7(2):227–37.
Song J, Li B, Cui Y, Zhuo C, Gu Y, Hu K, Wen J, Yi B, Shen J, Ma C, Fu T. QTL mapping and diurnal transcriptome analysis identify candidate genes regulating Brassica napus flowering time. Int J Mol Sci. 2021;22(14):7559.
Mason AS, Huteau V, Eber F, Coriton O, Yan G, Nelson MN, Chèvre AM. Genome structure affects the rate of autosyndesis and allosyndesis in AABC, BBAC and CCAB Brassica interspecific hybrids. Chrom Res. 2010;18:655–66.
Katche E, Gaebelein R, Idris Z, Vasquez-Teuber P, Lo YT, Nugent D, Mason AS. Stable, fertile lines produced by hybridization between allotetraploids Brassica juncea (AABB) and Brassica carinata (BBCC) have merged the A and C genomes. New Phytol. 2021;230(3):1242–57.
Jesske T, Olberg B, Schierholt A, Becker HC. Resynthesized lines from domesticated and wild Brassica taxa and their hybrids with B. napus L.: genetic diversity and hybrid yield. Theor Appl Genet. 2013;126:1053–65.
Song KM, Osborn TC, Williams PH. Brassica taxonomy based on nuclear restriction fragment length polymorphisms (RFLPs) 1. Genome evolution of diploid and amphidiploid species. Theor Appl Genet. 1988;75:784–94.
Song K, Lu PI, Tang K, Osborn TC. Rapid genome change in synthetic polyploids of Brassica and its implications for polyploid evolution. Proc Natl Acad Sci. 1995;92(17):7719–23.
Dhaka N, Mukhopadhyay A, Paritosh K, Gupta V, Pental D, Pradhan AK. Identification of genic SSRs and construction of a SSR-based linkage map in Brassica juncea. Euphytica. 2017;213(1):1–13.
Aakanksha, Yadava SK, Yadav BG, Gupta V, Mukhopadhyay A, Pental D, Pradhan AK. Genetic analysis of heterosis for yield influencing traits in Brassica juncea using a doubled haploid population and its backcross progenies. Front Plant Sci. 2021;12:721631.
El-Soda M, Malosetti M, Zwaan BJ, Koornneef M, Aarts MG. Genotype× environment interaction QTL mapping in plants: lessons from Arabidopsis. Trends Plant Sci. 2014;19(6):390–8.
Zhang Y, Li YX, Yang W, Liu ZZ, Cheng L, Bo P, Yu L. Stability of QTL across environments and QTL-by-environment interactions for plant and ear height in maize. Agric Sci China. 2010;9(10):1400–12.
Zhao F, Xu S. Genotype by environment interaction of quantitative traits: a case study in barley. G3: Genes Genomes, Genet. 2012;2(7):779–88.
Yang DL, Jing RL, Chang XP, Li W. Identification of quantitative trait loci and environmental interactions for accumulation and remobilization of water-soluble carbohydrates in wheat (Triticum aestivum L) stems. Genetics. 2007;176(1):571–84.
Singh KH, Gupta MC, Srivastava KK, Kumar PR. Combining ability and heterosis in Indian mustard. J Oilseeds Res. 2003;20:35–9.
Singh M, Singh L, Srivastava SBL. Combining ability analysis in Indian mustard (Brassica juncea L. Czern & Coss). J Oilseed Brassica. 2010;1(1):23–7.
Singh N, Yadava DK, Sujata V, Singh R, Giri SC, Dass B, Barun S, Prabhu KV. Combining ability and heterobeltiosis for yield and yield contributing traits in high quality oil Indian mustard (Brassica juncea) genotypes. Indian J Agric Sci. 2015;85(4):498–503.
Binod K, Anil P, Shivramakrishnan R, Vinoth R, Singh SK, Mishra SB. Combining ability for yield and oil quality traits in Indian mustard [Brassica juncea (L.) Czern & Coss] using line× tester analysis. Adv Res J Crop Imp. 2016;7(1):22–31. https://doi.org/10.15740/HAS/ARJCI/7.1/22-31.
Varshney RK, Gaur PM, Chamarthi SK, Krishnamurthy L, Tripathi S, Kashiwagi J, Samineni S, et al. Fast-track introgression of “QTL-hotspot” for root traits and other drought tolerance traits in JG 11, an elite and leading variety of chickpea. Plant Genome. 2013;6:plantgenome2013-07. https://doi.org/10.3835/plantgenome2013.07.0022.
Han J, Lühs W, Sonntag K, Zähringer U, Borchardt DS, Wolter FP, Frentzen M. Functional characterization of [beta]-ketoacyl-CoA synthase genes from Brassica napus L. Plant Mol Biol. 2001;46(2):229.
Yan G, Li D, Cai M, Gao G, Chen B, Xu K, Wu X. Characterization of FAE1 in the zero erucic acid germplasm of Brassica rapa L. Breed Sci. 2015;65(3):257–64.
Khakdan F, Alizadeh H, Ranjbar M. Molecular cloning, functional characterization and expression of a drought inducible phenylalanine ammonia-lyase gene (ObPAL) from Ocimum basilicum L. Plant Physiol Biochem. 2018;130:464–72.
Liu F, Xie L, Yao Z, Zhou Y, Zhou W, Wang J, Gong C. Caragana korshinskii phenylalanine ammonialyase is up-regulated in the phenylpropanoid biosynthesis pathway in response to drought stress. Biotechnol Biotechnol Equip. 2019;33(1):842–54.
Iquebal MA, Sharma P, Jasrotia RS, Jaiswal S, Kaur A, Saroha M, Angadi UB, Sheoran S, Singh R, Singh GP, Rai A. RNAseq analysis reveals drought-responsive molecular pathways with candidate genes and putative molecular markers in root tissue of wheat. Sci Rep. 2019;9(1):13917.
Bergmann T, Menkhaus J, Ye W, Schemmel M, Hasler M, Rietz S, Leckband G, Cai D. QTL mapping and transcriptome analysis identify novel QTLs and candidate genes in Brassica villosa for quantitative resistance against Sclerotinia sclerotiorum. Theor Appl Genet. 2023;136(4):86.
Acknowledgements
OML is grateful to the University Grants Commission of India and Post-Graduate School, ICAR-Indian Agricultural Research Institute (IARI), New Delhi for granting a Senior Research Fellowship during the pursuit of a doctoral degree. The authors are thankful to the Division of Genetics, ICAR-IARI, New Delhi for providing the necessary facilities during the evaluation of ILs. The authors are also thankful to the ICAR-DRMR, Bharatpur, Rajasthan for extending help and resources for evaluating ILs.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
Oml conducted the field trials, analysed the data and prepared the draft manuscript. PV recorded observations and analysed data. PJ helped in the bioinformatics analyses. RS maintained the ILs, conducted the trial and recorded observations. MS, GS, MK and MLM helped in recording the observations in multilocation field trials. MAI and SJ guided in the bioinformatics analyses. MR supervised the cytological analysis. AW revised the manuscript. KHS conducted field trials at Bharatpur. RB and DK guided the research. NS developed the ILs, conceptualized & supervised the research, analysed the data and wrote the manuscript. All authors read the manuscript and agreed with its content.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Field as well as laboratory studies were conducted in full compliance with local regulations. This article did not involve any studies with human, animals, or endangered species. The collection of plant materials and the experimental procedures conducted in this study adhered to institutional, national and international guidelines and legal requirements.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Supplementary Table. Seed yield, water use efficiency (WUE), drought tolerance index (DTI), and mean relative performance (MRP) of introgression lines (ILs) along with their parents evaluated under rainfed (RE) and irrigated (IE) conditions in three environments.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Limbalkar, O.M., Vasisth, P., Singh, G. et al. Dissection of QTLs conferring drought tolerance in B. carinata derived B. juncea introgression lines. BMC Plant Biol 23, 664 (2023). https://doi.org/10.1186/s12870-023-04614-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-023-04614-z