1Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreports

epistatic interactions of major effect drought QTLs with genetic background loci determine grain yield of rice under drought stressShailesh Yadav1, Nitika sandhu1, Ratna Rani Majumder1, Shalabh Dixit1, Santosh Kumar2, s. p. singh3, N. p. Mandal4, s. p. Das5, Ram Baran Yadaw6, Vikas Kumar singh7, Pallavi sinha8, Rajeev K. Varshney 8 & Arvind Kumar1

Epistatic interactions of QTLs with the genetic background and QTL-QTL interaction plays an important role in the phenotypic performance of introgression lines developed through genomic-assisted breeding (GAB). In this context, NIL pairs developed with various drought QTL (qDTY) combinations in the genetic background of IR64, TDK1-Sub1 and Savitri backgrounds were utilized to study the interactions. Multi-season phenotyping of NIL pairs harboring similar qDTY combinations provided contrasting performance for grain yield under drought (RS) (classified as high and low yielding NILs) but nearly similar performance under non-stress(NS) conditions. Genome wide genotyping data revealed a total of 16, 5 and 6 digenic interactions were detected under RS conditions in low yielding NILs of IR64, TDK1-Sub1 and Savitri respectively while no significant interaction was found in high yielding NILs under RS and NS conditions in any of the genetic backgrounds used in this study. It is evident from this study that existence of epistatic interactions between QTLs with genetic background, QTL-QTL interaction and interactions among background markers loci itself on different chromosomes influences the expression of a complex trait such as grain yield under drought. The generated information will be useful in all the GAB program of across the crops for precise breeding.

Drought is one of the most severe climatic threat to rainfed rice production affecting on more than 23 million hectares of South and Southeast Asia1. A further increase in drought intensity has anticipated due to ongoing global climatic changes. Development of climate-adapted crop varieties is an urgent need to counter the effect of more than 50% yield losses due to abiotic stresses2,3. Widely grown high yielding rice varieties are highly sensitive to drought. The recent advances in molecular marker and genomics technologies have played a major role in selection of drought-tolerant traits. Rice drought breeding programme at IRRI has successfully identified and val-idated some of the major and consistent effect grain yield QTLs under reproductive stage drought (qDTYs) such as qDTY1.1

4,5, qDTY2.16, qDTY2.2

7, qDTY3.16, qDTY3.2

4, qDTY6.18 and qDTY12.1

9 using molecular markers. A consid-erable effort has been made to improve drought susceptible mega varieties of rice through marker-assisted intro-gression of qDTYs10–13. Assembling multiple alleles/QTLs/genes in a genetic background using marker-assisted introgression may be an appropriate approach to achieve the expected phenotypic variance in enhancing rice grain yield under drought stress7,12,14, submergence tolerance15, salt tolerance16, cold tolerance17, resistance to blast18–20, bacterial blight21–23, brown planthopper24,25. The introgression of major effect QTLs/genes and its com-binations do not always result in desirable improvement for the targeted trait(s). In some cases, the effect of identified genetic loci in a QTL/gene pyramiding program is not sufficient to fully explain the desired phenotypic variation due to various interactions occurring for a complex quantitative trait. The possible interactions between

1Rice Breeding Platform, International Rice Research Institute, DAPO Box 7777, Metro Manila, Philippines. 2icAR Research complex for eastern Region, Patna, Bihar, india. 3Bihar Agricultural University, Sabour, Bihar, india. 4central Rainfed Upland Rice Research station, national Rice Research institute, Hazaribagh, Jharkhand, india. 5icAR Research complex for neH Region, tripura centre, Lembucherra, tripura, india. 6Regional Agriculture Research Station, nARc, Bara, nepal. 7international Rice Research institute, South Asia Hub, icRiSAt, Patancheru, Hyderabad, india. 8International Crops Research Institute for the Semi-arid Tropics, Patancheru, Hyderabad, India. correspondence and requests for materials should be addressed to A.K. (email: [email protected])

Received: 22 August 2018

Accepted: 17 January 2019

Published: xx xx xxxx

opeN

2Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

the genetic make-up, epistasis, pleiotropy and linkage among/between the introgressed loci, with the genetic background and environment could have negative/positive effects on the expression of the traits that should never be underestimated26–28. QTL pyramiding approach has enormous potential to improve our understanding of interactions among QTLs and also help to develop new strategies to increase the efficiency of marker-assisted selection program. In this context, development of NILs using marker assisted selection (MAS) is a way forward to verify the effects of one or few major effect QTLs as well as increase the recovery of recurrent parent genome at a faster rate29. Rice NILs for complex traits such as drought7,13, blast19,30 and brown planthopper25,31 were devel-oped and evaluated in different genetic backgrounds of popular rice varieties.

NILs harboring QTLs/genes for abiotic/biotic stresses in the elite backgrounds with superior performances are useful material for immediate release to farmer’s field. NILs are also useful source to study the stable epistatic interactions between the introgressed donor segment and the rest of the recurrent parent genome which gives an important clue about differences in performance of the same set of NILs under variable conditions. The expres-sions of QTLs/genes for complex traits are strongly affected by the genetic backgrounds32,33. The estimation of genetic backgrounds or recurrent parent genome (RPG) recovery with different types of markers and number of markers also affect the estimation of RPG recovery as a good number of markers with good coverage estimates better information19. Better estimation of RPG recovery at early generation of selection can be useful to eliminate linkage drags which can cause unexpected traits in the MAS breeding products because of neighboring genes of the target QTLs. Also, the favorable allele at one locus in a particular background may become unfavorable in another genetic background due to the epistatic interaction34. Role of epistasis phenomena in plant breeding has been discussed for a long time, but the extent of the expression level of quantitative traits is due to this com-plex genetic phenomena are poorly understood35. Epistatic or non-allelic interaction is a genetic factor and may change the magnitude of phenotypic expression of QTLs by suppressing/enhancing the loci effect associated with a trait to the other loci of genetic background34,36. The epistatic effect is a deviation from the sum of independent effects of QTLs/genes therefore biased estimation of QTL effects is possible in genetic models assuming no inter-action in QTL mapping studies37,38. The epistatic interaction of QTLs and background or modifying loci has been recently reported for quantitative traits in various crops38–41.

Genomic interactions play an important role in deciding the performance of MAS lines for grain yield under drought and various significant QTL x QTL interactions in pyramided lines have been reported for drought-related traits12,28,42,43. Most of these studies highlighted capture of such interactions for increased per-formance of QTLs pyramided lines under drought. But none of these studies reported reasons for a lower yield of NILs with same combinations of pyramided QTLs in the same or different recipient backgrounds. Keeping in this view, the present study is an effort (i) to evaluate the performance of drought NILs of IR64, TDK1-Sub1 and Savitri backgrounds under reproductive stage moderate and severe drought stress and non-stress conditions (ii) to identify promising drought tolerant NILs with higher recovery of recipient genome (iii) to study epistatic interaction of major effect grain yield drought QTLs and with background marker loci.

ResultsPhenotypic performance of drought NILs possessing qDTY under RS and NS condi-tions. Significant phenotypic differences under RS and NS were observed among the parents and NILs in IR64, TDK1-Sub1 and Savitri backgrounds. The results of ANOVA, means, LSD, heritability of IR64 NILs with a qDTY and qDTY combinations under RS (SS and MS) and NS are summarized in Table 1. Among IR64 NILs, the average yield reduction was found 87% under SS and 74% under MS as compared to NS conditions. Trial mean of NILs in IR64 background for GY was 690 kg ha−1 under SS in 2015 DS while it was 1358 kg ha−1 under MS in 2017 DS. Broad-sense heritability (H) of GY under RS ranged from 0.78 to 0.89 in the years of 2015 and 2017 DS. The NIL pairs harboring same qDTY or qDTY combinations with on par background recovery, pro-vided significantly different grain yield (GY) under similar level of drought stress were classified into high and low yielding NILs (Table 1). For instance, nine IR64 NILs possessing qDTY12.1 QTL yielded 25 to 3375 kg ha−1 under SS situation, whereas no such difference was observed under NS situation (4490 to 6470 kg ha−1). Another four IR64 NILs possessing qDTY2.3 yielded 480 kg ha−1 to 2495 kg ha−1 under MS. In NS conditions, trial mean for GY was 5360 kg ha−1 and ranged from 4375 kg ha−1 to 6958 kg ha−1. Nine NILs with same combina-tion of qDTY1.2 + qDTY12.1 have shown a large variation for grain yield under drought and ranged from 89 to 2141 kg ha−1. The huge variation in grain yield under drought even having same QTL combinations was might be due to some hidden interactions and analysis of such interaction were explored and discussed in this study. Similarly three NILs with qDTY2.2 + qDTY2.3 yielded 572 to 1455 kg ha−1; four NILs with qDTY4.1 + qDTY12.1 yielded 369 to 936 kg ha−1, two NILs with qDTY1.1 + qDTY1.2 + qDTY12.1 yielded 289 to 611 kg ha−1 and two NILs with qDTY4.1 yielded 414 to 1031 kg ha−1

. The mean plant height (PH) of IR64 NILs under RS was higher than the recurrent background IR64 while PH of NILs under NS was very similar to that of IR64 (Table 1). Days to flowering (DTF) were affected significantly with severity of stress and flowering delay of 10 days under SS and 7 days under MS was observed.

Mean performances of measured traits of TDK1-Sub1 NILs, heritability and LSD is presented in Table 2. Grain yield under severe stress (GYSS) ranged from 0 kg ha−1 to 479 kg ha−1 with trial mean of 327 kg ha−1. The recurrent parent TDK1-Sub1 yielded 14.4 kg ha−1 while the donor parent IR55419-04 performed well with GY of 881 kg ha−1 under SS situation. GY mean of MS trial (GYMS) was 889 kg ha−1, with range of 300 kg ha−1 to 1555 kg ha−1. The average yield reduction was more than 90% under SS while 84% less grain yield was recorded under MS among TDK1-Sub1 NILs. Heritability for GY was recorded 0.89 to 0.91 in RS and 0.58 in NS. The pyr-amided three NILs with qDTY6.1 + qDTY6.2 combination yielded 0 to 479 kg ha−1 under SS situation. However, the performance of both lines was 300 and 1555 kg ha−1 respectively in MS. Similarly, NIL pair consisting of qDTY3.1 + qDTY6.1 + qDTY6.2 combination yielded 86 kg ha−1 and 479 kg ha−1 under SS and 1085 kg ha−1 to 1555 kg ha−1 under MS. A reduction in PH was observed under SS and MS conditions compared to NS condition.

3Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

The mean PH under SS, MS and NS were 64 cm, 69 cm and 92 cm, respectively. Significant differences in DTF with trial mean of 94 days and flowering differences of 7 days between entries was observed under RS. Earliest line flowered in 80 days while most of the lines were of late maturing type and it flowered in 96 days. As observed in the case of IR64 NILs, TDK1-Sub1 NILs was also produced nearly similar grain yield under NS situation.

The same traits were also observed in the Savitri NILs (Table 2). A reduction in mean GY under SS and MS compared with mean GY under NS was 90% and 61%, respectively. The recurrent parent Savitri yielded 48 kg ha−1 while the donor parent IR77298-5-6-18 yielded 236 kg ha−1 under SS situation. NIL pair carrying single qDTY3.2 yielded 134 to 517 kg ha−1 under SS, whereas four NILs carrying two qDTY (qDTY3.2 + qDTY12.1) yielded 18 to 627 kg ha−1. The significant deviation observed in GY suggesting some negative interaction between background markers loci.

NILs QTL and QTLs combinations

GY (kg ha−1) DTF (days) PH (cm)

Classification

2015DS 2017DS 2017DS 2015DS 2017DS 2017DS 2015DS 2017DS 2017DS

SS MS NS SS MS NS SS MS NS

IR 102793:1-11-66-3-1-1 qDTY1.2 + qDTY12.1 2141 1455 4375 74 73 71 75 67 91 High yielding

IR 102796-14-140-1-1-1 qDTY1.2 + qDTY12.1 788 2740 4928 85 79 78 64 66 92 High yielding

IR 102793:1-11-192-1-1-4 qDTY1.2 + qDTY12.1 1386 1110 4748 84 82 79 64 64 89 High yielding

IR 102793:1-11-64-1-1-2 qDTY1.2 + qDTY12.1 2002 2750 6145 75 75 76 67 66 90 High yielding

IR 102784:2-42-17-2-1-1 qDTY1.2 + qDTY12.1 540 1810 6191 92 83 96 55 87 92 High yielding

IR 102796-14-132-1-1-3 qDTY1.2 + qDTY12.1 89 515 5855 85 84 81 55 58 86 Low yielding

IR 102796-14-65-2-1-1 qDTY1.2 + qDTY12.1 138 1050 5990 86 77 76 58 60 92 Low yielding

IR 102784:2-42-16-1-1-2 qDTY2.3 + qDTY3.2 1426 1860 5430 78 72 70 60 71 92 High yielding

IR 102784:2-62-481-2-1-1 qDTY2.3 + qDTY3.2 191 1255 5200 95 84 82 59 78 92 Low yielding

IR 102783:2-70-112-2-1-4 qDTY12.1 1580 1640 4615 84 80 75 63 64 93 High Yielding

IR 102783:2-70-112-3-1-1 qDTY12.1 3375 2210 5143 74 75 71 75 68 96 High yielding

IR 102784:2-42-88-2-1-2 qDTY12.1 1675 2280 6390 83 77 75 62 62 92 High yielding

IR 102783:2-70-135-4-1-1 qDTY12.1 996 1225 4561 77 76 71 61 56 80 High yielding

IR 102783:2-70-139-2-1-2 qDTY12.1 1194 2060 5848 77 76 73 67 63 97 High yielding

IR 102784:2-42-99-2-1-3 qDTY12.1 209 1330 4820 96 83 79 55 63 92 Low yielding

IR 102784:2-90-385-1-1-3 qDTY12.1 160 1485 6470 94 79 78 54 66 87 Low yielding

IR 102784:2-62-434-1-1-1 qDTY12.1 25 645 4490 99 84 83 55 61 89 Low yielding

IR 102784:2-90-385-3-1-1 qDTY12.1 39 810 4780 96 83 81 52 51 93 Low yielding

IR 102784:2-62-66-1-1-2 qDTY1.1 + qDTY1.2 + qDTY12.1 611 845 5255 88 81 79 100 87 115 High yielding

IR 102784:2-118-22-1-1-2 qDTY1.1 + qDTY1.2 + qDTY12.1 289 1595 6068 91 80 80 55 61 86 Low yielding

IR 102784:2-42-3-1-1-2 qDTY2.2 + qDTY2.3 1455 1960 5465 75 73 70 64 78 84 High yielding

IR 102796-14-124-1-1-3 qDTY2.2 + qDTY2.3 996 2235 5220 85 73 69 62 87 91 High yielding

IR 102784:2-42-136-1-1-3 qDTY2.2 + qDTY2.3 572 1255 4988 87 75 70 57 78 81 Low yielding

IR 102784:2-42-127-3-1-1 qDTY2.3 1692 730 4497 94 75 73 58 82 84 High yielding

IR 102793:1-11-69-1-1-3 qDTY2.3 1667 2495 5600 74 73 71 61 80 88 High yielding

IR 102784:2-42-136-2-1-1 qDTY2.3 33 480 6865 103 84 78 49 82 85 Low yielding

IR 102784:2-42-138-1-1-1 qDTY2.3 157 1650 6270 98 75 75 59 83 85 Low yielding

IR 102793:1-11-189-2-1-2 qDTY4.1 1031 1720 6290 89 83 73 63 91 98 High yielding

IR 102784:2-118-15-1-1-2 qDTY4.1 414 1625 6958 92 82 80 57 84 91 Low yielding

IR 102784:2-118-549-2-1-2 qDTY4.1 + qDTY12.1 852 2125 6378 88 82 79 57 84 90 High yielding

IR 102793:1-11-229-3-1-1 qDTY4.1 + qDTY12.1 656 1900 5618 87 81 78 59 83 90 High yielding

IR 102784:2-89-284-3-1-1 qDTY4.1 + qDTY12.1 936 1175 5502 87 83 80 59 83 90 High yielding

IR 102784:2-89-284-1-1-3 qDTY4.1 + qDTY12.1 369 1600 6145 97 83 82 52 82 90 Low yielding

IR64 — 642 715 5910 88 79 78 63 65 93 Background

IR 86918-B-315 qDTY1.1, qDTY1.2 654 800 4715 78 88 84 101 88 107 Donor

Way Rarem qDTY12.1 271 510 5050 98 80 96 86 80 122 Donor

IR 77298-14-1-2-17 qDTY2.2, qDTY4.1 490 1930 3868 87 82 81 93 79 99 Donor

Vandana qDTY2.3, qDTY2.3 2074 1890 4760 75 73 72 89 76 92 Donor

Trial Mean 690 1358 5360 90 83 80 69 82 93

LSD0.05 848.44 998 1416 9.52 8.52 5.9 10.92 8.92 10.16

H 0.78 0.89 0.58 0.82 0.85 0.92 0.91 0.80 0.88

Table 1. Mean performance for grain yield (GY), days to flowering (DTF) and plant height (PH) of IR64 NILs under severe stress (SS), moderate stress (MS) and non-stress (NS) conditions. Note: LSD0.05: least significant difference at the 5% confidence level; H: Broad sense heritability.

4Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

Mean PH of Savitri NILs was significantly reduced under drought stress conditions ranged from 61 to 85 cm in SS, 58 to 75 cm in MS and from 78 to 102 cm in NS conditions. PH was reduced by 17 to 20 cm in drought stress. The mean days to flowering (DTF) in Savitri NILs ranged from 72 to 81 days in SS, 76 to 82 days in MS and from 74 to 81 days in NS conditions. DTF was severely affected under drought for recipient parent Savitri varied from 103 days in NS while 113 and 96 days was recorded in SS and MS conditions respectively. The heritability (H) was high for DTF (0.84,0.85,0.92) and medium to high for PH (0.72,0.80,0.88) and GY (0.74,0.89,0.58) under SS, MS and NS conditions for Savitri NILs (Table 2).

Promising drought NIL lines with the maximum genetic background coverage of recipient parents. Recurrent parent genome recovery (RPG) were estimated using 111, 107 and 89 polymorphic SSR markers among NILs carrying qDTYs in IR64, TDK-Sub1 and Savitri backgrounds, respectively. The background recovery among the NILs of multiple genetic backgrounds ranged from 85 to 93% in the case of IR64 NILs, 86 to 90% in the case of TDK-Sub1 NILs and 80 to 89% in the case of Savitri NILs. Most of the IR64 NIL lines showed more than 90% RPG recovery followed by TDK-Sub1 and least recovered lines belong to Savitri back-ground. Graphical representation of the selected pyramided NILs achieving higher genome recovery of respective backgrounds is depicted in Supplementary Fig. S1 (IR64 NILs with qDTY2.3, qDTY3.2), Supplementary Fig. S2 (TDK1-Sub1 NILs with qDTY6.1, qDTY6.2) and Supplementary Fig. S3 (Savitri NILs with qDTY3.2, qDTY12.1). The selected promising NILs capturing higher recurrent genome background and yielded well under SS, MS and NS conditions is presented in Table 3. These superior NILs with excellent drought tolerance could be used a drought tolerant line/variety to be released after multi-location evaluation in national and provincial coordinated trials in target countries.

QTL–QTL and QTL-background interaction in drought NILs. A total of 16 digenic interactions between marker loci were detected in low yielding NILs of IR64 background under SS and MS situation. Out of 16 digenic interactions, we identified 6 interactions between drought grain yield QTLs (qDTY2.2, qDTY2.3 and qDTY3.2) and background markers loci while 9 interactions were found among background markers loci of different chromosome. A Significant negative QTL-QTL interaction was also found between introgressed qDTY1.2 (RM543-RM212) located on chromosome 1 with qDTY2.3 (RM263-RM573) identified on chromosome 2. In context of QTL and background marker loci interactions, qDTY2.2 (RM154-RM236) interacted negatively with two genetic loci (RM7-RM332) on chromosome 3 and chromosome 4 (RM142-RM119), while qDTY2.3 (RM263-RM573) was interacted negatively with genetic loci (RM108-RM215) on chromosome 9 and chromo-some 10 (RM271-RM269) under SS situation (Table 4, Fig. 1a). Similarly, qDTY3.2 interacted negatively with

NILS QTL and QTLs combinations

GY (kg ha−1) DTF (days) PH (cm)

Background Classification

2015DS 2017DS 2017DS 2015DS 2017DS 2017DS 2015DS 2017DS 2017DS

SS MS NS SS MS NS SS MS NS

IR 102777-5-83-1-2-7 qDTY6.1 + qDTY6.2 316 465 5590 95 78 78 93 96 126 TDK-Sub1 High yielding

IR 102774-26-8-3-2-5 qDTY6.1 + qDTY6.2 0 300 4915 96 99 85 66 65 108 TDK-Sub1 Low yielding

IR 102777-6-86-2-2-7 qDTY6.1 + qDTY6.2 163 540 5695 89 83 84 67 60 104 TDK-Sub1 Low yielding

IR 102777-6-86-2-2-11 qDTY3.1 + qDTY6.1 + qDTY6.2 479 1555 4966 80 77 74 71 66 93 TDK-Sub1 High yielding

IR 102774-15-32-3-1-2 qDTY3.1 + qDTY6.1 + qDTY6.2 86 1085 4958 84 76 76 93 84 129 TDK-Sub1 Low yielding

TDK1-Sub1 — 14.4 527 5412 70 80 82 61 71 110 Recipient —

IR55419–04 qDTY3.1, qDTY6.1, qDTY6.2 881.5 652 5361 72 89 102 68 70 92 Donor —

Trial Mean 327 889 5847 94 93 88 64 69 92

LSD0.05 305 929 1371 6.8 2.7 1.15 17.8 5.08 12.8

H 0.91 0.89 0.58 0.93 0.85 0.92 0.89 0.80 0.88

IR 106523-21-28-1-2-B qDTY3.2 517 2565 7098 72 78 78 68 61 90 Savitri High yielding

IR 106531-4-22-3-3-B qDTY3.2 134 2185 6478 78 82 78 74 68 90 Savitri Low yielding

IR 106522-41-8-3-B qDTY3.2 + qDTY12.1 627 1650 6035 72 76 74 76 65 92 Savitri High yielding

IR 106523-3-9-3-2-B qDTY3.2 + qDTY12.1 545 2165 6058 75 77 81 63 58 79 Savitri High yielding

IR 106531-13-31-2-3-B qDTY3.2 + qDTY12.1 145 1425 5725 76 78 77 61 60 78 Savitri Low yielding

IR 106529-15-33-1-2-B qDTY3.2 + qDTY12.1 18 1090 5035 81 80 79 85 75 102 Savitri Low yielding

Savitri — 48 903 6641 113 96 103 48 66 104 Recipient —

IR 77298-5-6-18 qDTY3.2 236 1767 6393 96 84 88 57 57 84 Donor —

IR74371-46-1-1 qDTY12.1 999 1715 4238 76 73 70 74 73 98 Donor

Trial Mean 563 2319 6069 88 80 83 61 64 87

LSD0.05 540 522 808 8.8 3.21 2.12 6.46 4.33 7.07

H 0.74 0.89 0.58 0.84 0.85 0.92 0.72 0.80 0.88

Table 2. Mean performances for grain yield (GY), days to flowering (DTF) and plant height (PH) of TDK-Sub1 and Savitri NILs under severe stress (SS), moderate stress (MS) and non- stress (NS). Note: LSD0.05: least significant difference at the 5% confidence level; H: Broad sense heritability.

5Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

background markers loci at RM307-RM537 on chromosome 4 and at RM206-RM254 located on chromosome 11 under MS condition (Table 4, Fig. 1b). Out of 9 digenic background markers interaction detected in present study, two negative interactions with marker interval (RM7-RM232; chromosome 3) was found at single loci position at 65 cM with background markers located on chromosomes (9 and10). Similarly, two negative interactions with marker interval (RM246-RM543; chromosome 1) were found at single loci position at 141 cM with background markers located on chromosomes (9 and 11). The PVE% of epistatic QTLs ranged from 3.47 to 25.46% under SS to MS. Surprisingly, none of high yielding NILs expressed any interactions in SS and MS conditions (Fig. 1c,d). Also, no significant interactions were also observed under NS conditions for IR64 NILs including both high yielding as well as low yielding NILs (Fig. 1e). A heat map according to LOD score along the chromosomes reflecting all the significant interactions among low yielding IR64 NILs under SS and MS situations is presented in Supplementary Fig S4.

Low yielding NILs of TDK1-Sub1 background had shown 5 significant epistatic interactions and all these interactions were detected between introgressed QTL qDTY6.1 (RM587-510) and background markers loci located across chromosomes 2, 3, 5, 6 and 11. Two negative interactions were found between QTL qDTY6.1 (RM587-510) and marker pairs RM115-RM275 and RM473-RM206 located on chromosomes 6 and 11 respectively while 3 interactions were found between background markers loci (RM300-RM475 on chromosome 2; RM411-RM319 on chromosome 3; RM592-RM437 on chromosome 5) with a common loci (RM587-510) associated with qDTY6.1 (Table 4, Fig. 2a). The PVE% of interacting QTLs ranged from 12.04 to 19.95%. Low yielding NILs under MS (Fig. 2b) and high yielding NILs under SS (Fig. 2c) and MS (Fig. 2d) have not shown any significant interac-tions. Also, no significant interactions were observed under well-watered (NS) conditions for TDK1-Sub1 NILs including both high yielding and low yielding NILs (Fig. 2e). A heat map according to LOD score along the chromosomes reflecting the significant digenic interactions found among low yielding TDK1-Sub1 NILs under SS situations is presented in Supplementary Fig S5. In total, 6 epistatic interactions were identified among back-ground markers loci under SS and MS in low yielding lines of Savitri background (Table 4, Fig. 3a,b). The inter-action between introgressed QTLs (qDTY3.2, qDTY12.1) and background markers loci was not identified in this background. Despite, away from position of introgressed QTL (qDTY3.2; 10 cM) flanking markers (RM7-RM251, 67 cM) interacted negatively with background marker pair RM55-570 on chromosome 3 under SS situation. Under MS, on chromosome 4 (RM537- RM335) and chromosome 10 (RM244- RM239) were interacted with a common genetic locus (RM136- RM275) on chromosome 6 with explaining PVE upto 12.76%.No interactions were found in high yielding NILs under SS and MS conditions (Fig. 3c,d). Also, no significant interactions were found under NS condition for Savitri NILs including both high yielding as well as low yielding NILs (Fig. 3e). A heat map according to LOD score along the chromosomes reflecting the significant digenic interactions found among low yielding Savitri NILs under SS and MS situations is presented in Supplementary Fig S6.

DiscussionEpistasis is a genetic phenomenon of interaction which may enhance or reduce the expression (depending on degree and direction) of interacting loci underlying QTLs associated with the complex trait44,45. QTL mapping and interaction studies in recent years in different crops suggested the possible presence of epistasis such as (1) interactions between introgressed QTLs (2) interactions between introgressed QTLs and ‘background’ (modify-ing) loci and (3) interactions between complementary background loci46 affecting overall expressions of trait. To further understand this interaction (QTL-QTL and QTL to background), we have selected range of rice drought NILs developed in different genetic backgrounds.

Several large and consistent effect QTLs for grain yield under drought (qDTY1.1, qDTY2.2, qDTY3.1, qDTY12.1) were identified and validated in various rice genetic backgrounds using the QTL associated markers4,6,7,9. Exploiting the available molecular markers for introgression of consistent drought grain yield QTLs (qDTYs)

Background Promising NILs Parentage Foreground

RPG recovery %

GY (kg ha−1)

Severestress

Moderate stress

Non-stress

IR64 IR 102784:2-42-16-1-1-2 IR 99621-181/IR 99620-158//IRRI 149 qDTY2.3 + qDTY3.2 93 1426 1860 5430

IR64 IR 102793:1-11-229-3-1-1 IR 99622-312/IR 99620-158//IRRI 149 qDTY4.1 + qDTY12.1 92 656 1900 5618

IR64 IR 102793:1-11-192-1-1-4 IR 99621-14/IR 99619-361//IRRI 149 qDTY1.2 + qDTY12.1 91 1386 1110 4748

IR64 IR 102784:2-62-66-1-1-2 IR 99621-181/IR 99620-158//IRRI 149 qDTY1.1 + qDTY1.2 + qDTY12.1 91 611 845 5255

IR64 IR 102796-14-140-1-1-1 IR 99622-306/IR 99620-228//IRRI 149 qDTY1.2 + qDTY12.1 90 788 2740 4928

IR64 IR 102796-14-124-1-1-3 IR 99622-306/IR 99620-228//IRRI 149 qDTY2.2 + qDTY2.3 90 996 2235 5220

IR64 IR 102783:2-70-112-3-1-1 IR 99621-14/IR 99619-361//IRRI 149 qDTY12.1 90 3375 2210 5143

IR64 IR 102793:1-11-64-1-1-2 IR 99621-181/IR 99620-158//IRRI 149 qDTY1.2 + qDTY12.1 89 2002 2750 6145

IR64 IR 102793:1-11-66-3-1-1 IR 99622-312/IR 99620-158//IRRI 149 qDTY1.2 + qDTY12.1 89 2141 1455 4375

TDK1-Sub1 IR 102777-5-83-1-2-7 IR 90266-B-542-1/IR07F289 qDTY6.1 + qDTY6.2 90 316 465 5590

TDK1-Sub1 IR 102777-6-86-2-2-11 IR 90266-B-542-1/IR07F289 qDTY3.1 + qDTY6.1 + qDTY6.2 89 479 1555 4966

Savitri IR 106523-21-28-1-2-B IR 102773-1061//IR 90252-B-548-2-B/IR 90250-B-475-1-B qDTY3.2 89 517 2565 7098

Table 3. Selected promising drought NILs identified with maximum recurrent parent genome recovery. Note: RPG: Recurrent parent genome recovery.

6Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

through MAS is most desirable step for generation of drought-tolerant lines/varieties. In this study, contrasting performance of a pair of NILs possessed similar qDTY and qDTY combinations were analyzed to identify the possible epistatic interactions of introgressed qDTYs QTLs and background marker loci under different level of stress and non -stress conditions.

Genetic backgrounds IR64, TDK1-Sub1 and Savitri, used in this study to develop NILs which produces higher grain yield under drought stress condition, widely cultivated and accepted by farmers in India, Lao PDR and Nepal, respectively. Foreground selection using the flanking markers associated with qDTYs QTLs ensured successful incorporation of the target QTLs in NILs of respective backgrounds (Supplementary Table S1). Marker-assisted background selection using uniformly distributed polymorphic markers have been extensively used to track recovery of background genome19,47,48. Rice drought breeding has witnessed some of the success-ful MAS introgression of qDTYs QTLs in elite backgrounds7,8,12,13 using SSR and SNP markers for background selection. The Broad-sense heritability estimates in Tables 1, 2 were recorded moderate to high for grain yield in all the experiments varied from 0.58 (58%) to 0.89 (89%). Many of our previous studies have also reported

NILs Background Trait Name

ChromoSome1* Position 1†

FlankingMarkers (position1)

ChromoSome 2‡ Position 2⁑

FlankingMarkers (position 2)

LODepi§

PVEepi(%)¶

LODadd§§

PVE add(%)¶¶ Add1⁑⁑ Add2‡‡

Add1 by Add2††

IR64 GYMS 3 10 RM231-RM22(qDTY3.2)

4 0 RM307- RM537 7.67 25.46 3.17 4.89 124.65 458.93 −298.92

IR64 GYMS 3 10 RM231-RM22(qDTY3.2)

11 110 RM206- RM254 7.32 15.64 3.025 3.00 214.96 231.23 −235.02

IR64 GYSS 1 146 RM543-RM212(qDTY1.2)

2 124.7 RM263- RM573(qDTY2.3)

3.20 4.98 2.60 2.57 625.07 −248.3 −972.50

IR64 GYSS 1 141 RM246-RM543 9 86.7 RM108-RM215 3.33 4.75 2.50 2.60 603.78 −269.83 −947.59

IR64 GYSS 1 141 RM246-RM543 11 0 RM286-RM332 3.60 4.66 2.68 2.47 −363.62 579.67 −974.07

IR64 GYSS 2 124.7 RM263- RM573(qDTY2.3)

9 86.7 RM108-RM215 3.82 4.43 2.79 2.77 −396.48 507.01 −990.95

IR64 GYSS 2 9.8 RM154-RM236(qDTY2.2)

3 65 RM7- RM232 3.26 3.92 2.94 2.65 −426.60 481.96 −966.28

IR64 GYSS 2 9.8 RM154-RM236(qDTY2.2)

4 75 RM142-RM119 3.20 4.66 2.68 2.47 −363.62 579.67 −974.07

IR64 GYSS 2 84.8 RM341- RM475 10 21.3 RM222-RM311 3.67 3.54 2.84 2.58 434.01 −424.80 −962.57

IR64 GYSS 2 84.8 RM341-RM475 9 86.7 RM108-RM215 3.75 4.63 2.99 2.60 553.25 −393.57 −992.71

IR64 GYSS 2 124.7 RM263- RM573(qDTY2.3)

10 66.3 RM271-RM269 3.52 4.62 2.74 2.58 −400.88 545.46 −999.2

IR64 GYSS 3 65 RM7-RM232 9 86.7 RM108-RM215 3.81 4.63 2.94 2.59 554.08 −392.93 −991.94

IR64 GYSS 3 65 RM7-RM232 10 36.3 RM311-RM271 3.52 4.43 2.76 2.58 504.44 −394.71 −991.78

IR64 GYSS 4 75 RM142-RM119 9 86.7 RM108-RM215 3.60 4.63 2.92 2.59 553.86 −393.27 −992.10

IR64 GYSS 8 90 RM515-RM419 9 91.7 RM108-RM215 3.27 4.64 2.92 2.58 555.68 −391.84 −992.32

IR64 GYSS 10 16.3 RM222-RM311 10 71.3 RM269-RM228 3.39 3.47 2.68 2.50 −501.44 328.77 −1027.18

TDK1-Sub1 GYSS 2 91.9 RM300- RM475 6 10 RM587- RM510(qDTY6.1)

31.06 14.09 13.302 4.09 −73.18 72.46 −72.46

TDK1-Sub1 GYSS 3 150 RM411- RM319 6 10 RM587- RM510(qDTY6.1)

30.38 12.40 11.162 5.80 −22.97 50.06 −26.05

TDK1-Sub1 GYSS 5 36.7 RM592- RM437 6 10 RM587- RM510(qDTY6.1)

30.40 17.90 13.89 3.93 −122.05 113.11 −113.11

TDK1-Sub1 GYSS 6 10 RM587- RM510(qDTY6.1)

6 115 RM275- RM30 28.0 14.88 10.87 4.98 45.55 −42.51 −44.08

TDK1-Sub1 GYSS 6 10 RM587- RM510(qDTY6.1)

11 100 RM473-RM206 30.14 19.95 14.30 4.76 40.54 −35.12 −35.17

Savitri GYMS 4 8.5 RM537-RM335 6 65 RM136-RM275 3.43 10.56 2.41 3.90 2.62 −211.62 −382.62

Savitri GYMS 10 15 RM244-RM239 6 65 RM136-RM275 3.43 12.76 2.58 4.22 −211.62 2.62 −382.62

Savitri GYSS 2 49.8 RM71- RM290 3 182 RM55- RM570 4.31 10.78 2.80 5.01 10.66 −58.66 −238.33

Savitri GYSS 3 67 RM7- RM251 3 187 RM55- RM570 4.10 12.76 2.61 7.94 −58.66 10.66 238.33

Savitri GYSS 3 167 RM319- RM55 7 73.2 RM10- RM47 6.33 14.31 2.44 2.12 −69.38 144.08 −94.60

Savitri GYSS 3 182 RM55- RM570 11 0 RM4B- RM332 4.21 11.78 2.88 4.01 −58.66 10.66 −238.33

Table 4. Epistatic interaction of major effect loci of grain yield under drought with background loci of IR64, TDK1-Sub1 and Savitri NILs. Note: GYMS: Grain yield under moderate stress; GYSS: Grain yield under severe stress; *Chromosome ID at the first scanning position; †Scanning position in cM of the first flanking marker pair; ‡Chromosome ID at the second scanning position; ⁑Scanning position in cM of the second flanking marker pair; §LOD score caused by epistatic effects; ¶PVE(%): Phenotypic variation explained by epistatic effects; §§LOD score caused by additive effects; ¶¶PVE(%): Phenotypic variation by the additive effects; ⁑⁑Estimated additive effect of position 1; ‡‡Estimated additive effect of position 2; ††Additive by additive epistatic effect at the two scanning positions.

7Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

moderate to high heritability for grain yield under drought stress for QTL identification4,8. In present study, 600 genome-wide SSR markers were used for estimating the similarity between backgrounds and their derived NILs. Promising NILs have captured 89–93% recovery of respective backgrounds (Table 3) while 44 NILs used in this study possessed 85–93% genome recovery. Phenotypic selection combined with marker-assisted background selection approach after qDTY introgression was followed, that could be one of the reasons not to achieve the expected recovery in some of the NILs. Similar concern was discussed earlier by Dixit et al.49 in his study of marker-assisted breeding to combine submergence and drought tolerance in rice. However, in future drought NILs advancement, donor segment can be further minimized by using the recombinant selection markers and generation of whole genome profiling data through SNP-chips19 of respective DTY QTLs. A similar approach was also explained by Jena et al.25 in context to rice NILs developed for brown planthopper (BPH) resistance having low recovery (82–89.5%) of recurrent parent genome IR24 even after three back crosses.

Role of digenic interactions between DTY QTLs affecting the performance of one QTL over the other drought QTL has been well documented in various QTL mapping and marker-assisted drought breeding pro-gramme12,42,43. Positive interaction of qDTY12.1 with qDTY2.3 and qDTY3.2 have been reported by Dixit et al.42 which increases grain yield under drought significantly while Shamsudin et al.12 identified positive digenic inter-actions of qDTY2.2 and qDTY3.1 with qDTY12.1 enhancing overall expression of drought-related traits. Interaction of qDTY3.2 with qDTY1.1 and qDTY12.1 reduces the flowering duration and subsequently increases GY under stress43. A positive epistatic interaction of qDTY4.1 and qDTY9.1 loci with qDTY7.1 have been reported in the pyr-amided line of Samba Mahsuri enhancing the grain yield under drought28. Apart from these known qDTYs QTL x QTL interactions, there are many unknown genomic interactions exist with background loci called as back-ground noise which should be identified for effective deployment of drought QTLs in MAS breeding program.

Interaction between the target QTL/gene with the genetic background has been discussed earlier for vari-ous complex traits in rice18,50,51. However, there have been limited attempts to identify particular background markers loci interacting with introgressed QTLs and influencing the overall expression of target trait. There is an urgent need to identify such epistatic gene interactions which complicate the genotype-phenotype relationship of complex traits such as drought. In present study, we have detected negative interactions between drought QTLs

Figure 1. Cyclic illustrations of epistatic QTLs of IR64 NILs for grain yield under various levels of drought stress. The dotted lines indicate marker pairs interacting significantly on same or different chromosomes with their corresponding LOD value due to epistatic effect (a) Grain yield under severe stress (GYSS): chromosome 1 (qDTY1.2 RM543-RM212) showed epistatic negative interaction with qDTY2.3 (RM263-RM573) on chromosome 2 while qDTY2.2 (RM154-RM236) and qDTY2.3 (RM263-RM573) showed epistatic interaction with background markers loci on chromosomes 3, 4, 9 and 10 (b) Grain yield under moderate stress (GYMS): chromosome 3 (qDTY3.2 RM232-RM22, 10 cM) showed epistatic interaction with background loci of chromosomes 4 and 11 (c,d) GYMS and GYSS of high yielding NILs: no significant interactions was found (e) Grain yield under non- stress (GYNS): no interaction.

8Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

(qDTY2.2, qDTY2.3 and qDTY3.2) with IR64 background markers loci located on chromosomes (2, 3, 4, 10 and 11), among background markers loci located on different chromosomes and between drought grain yield QTL (qDTY1.2) with qDTY2.3, which might be one of the reason of poor performance of low yielding NILs carrying DTY QTLs (Table 1). The influence of epistatic drought QTL interactions explained PVE% ranging from 3.47% to 25.46% which could be very crucial in determining the grain yield under drought. Five digenic interactions were found in low yielding NILs of TDK1-Sub1 background, where qDTY6.1 (RM587-RM510) had shown negative interaction with background marker located on five different chromosomes. Savitri background NILs had shown 6 epistatic interactions between background markers loci in low yielding lines. Many of such epistatic interactions are background specific. However, in our study, the high effect of many of the qDTY QTLs interaction with the genetic background marker loci indicate that these are biological interactions and not just statistical. The QTLs qDTY2.2, qDTY2.3, qDTY3.2 and qDTY6.1 have shown interactions with the background markers on the same chro-mosome in three other genetic backgrounds (data unpublished).Some earlier studies32,52 highlighted the impor-tance of epistatic QTL effect over the additive QTLs are influencing the expression of introgressed trait in crops like Arabidopsis and rice. Thomson et al.53 reported a loss of Saltol QTL effect in developed NILs and suggested an interaction between Saltol and other background loci. Similarly, Babu et al.54 observed the difference in level of seedling stage salt tolerance among Saltol introgressed advanced lines and assumed the existence of a possible interaction between Saltol QTL and genetic background but these reports did not report on the particular back-ground markers loci interacting negatively with Saltol QTL, resulting in variable expression of Saltol introgressed lines. Revealing such epistatic interactions in present study, between introgressed QTL and background markers loci sheds more light in understanding the differences in phenotypic expression of MAS introgressed/pyramided lines for quantitative traits that can be useful in future MAS programs.

ConclusionsPerformance of NILs carrying same major effect drought grain yield QTLs/QTLs combinations are strongly affected by interactions with the genetic background loci. This is the first report on in-depth analysis of loss of introgressed QTL effects in SS, MS due to negative interactions between drought QTLs and background markers loci, drought QTL-QTL interaction and among interacting background markers loci on various chromosomes in SS, MS in low yielding NILs but not under NS. For both low as well as high yielding NILs, no interaction was detected under NS. Absence of such negative interactions in high yielding NILs under MS, SS indicate future MAS program to carefully select against such interactions to increase grain yield under drought. Promising NILs in background of popular varieties (IR64, TDK1-Sub1 and Savitri) with absence of negative interactions and

Figure 2. Cyclic illustrations of epistatic QTLs of TDK1-Sub1 NILs for grain yield under various levels of drought stress (a) Grain yield under severe stress(GYSS) Chromosome 6 (qDTY6.1 RM587-RM510,10 cM) showed epistatic interaction with background loci of chromosomes 2, 3, 5 and 11 (b) Low yielding GYMS: no interaction, (c,d) GYMS and GYSS of high yielding NILs: no significant interactions was found (e) GYNS: no interaction.

9Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

performing well in drought stress and non-stress conditions might be a very useful lines for use in future MAS/variety to be released for drought-prone ecosystems.

MethodsPlant materials. A total of 44 drought NILs derived from three backgrounds (33 NILs of IR64, 5 NILs of TDK1-Sub1 and 6 NILs of Savitri) were used in this study (Supplementary Table S2). We initiated the process of NILs development with a total of 130 lines (70 lines from IR64, 32 lines from TDK1-Sub1 and 28 lines from Savitri) developed through marker-assisted breeding for the three genetic backgrounds for the present study. However, before including all the lines we attempted foreground selection with QTLs linked markers as well as background marker selection of all the lines with 600 genome wide background SSR markers and those lines which follows the minimum criteria of RPG recovery (>80%) were taken for the selection of NILs for further evaluation. As a result, finally we have selected a total of 44 drought NILs derived from three backgrounds (33 NILs of IR64, 5 NILs of TDK1-Sub1 and 6 NILs of Savitri) were used in this study. A total of eight sets of IR64 NILs (20 high yielding and 13 low yielding) possessing different QTL and QTL combinations was developed through crossing of IR64 with different drought tolerant donors contributing different grain yield QTLs. Four different drought tolerant donors were utilized for development of the selected NILs, for instance IR77298-14-1-2-17 (qDTY2.2 and qDTY4.1), IR 86918-B-315 (qDTY1.1 and qDTY1.2), Vandana (qDTY2.3 and qDTY3.2) and Way Rarem (qDTY12.1) were crossed with IR64 and through marker assisted backcross breeding the above mentioned NILs were developed7. The pyramided version of NILs was developed through intercrossing of NILs of different gene combinations7. Two sets of TDK1-Sub1 NILs (2 high yielding and 3 low yielding NILs) possessing differ-ent QTL combinations were utilized in the present study. The pyramided version of NILs were developed from the cross between IR55419-04/2*TDK1. After confirming the presence of three QTLs (qDTY3.1, qDTY6.1 and qDTY6.2) in pyramided line, it was utilized for backcrossing with submergence tolerant parent TDK1-Sub1 to develop NILs in TDK1-Sub1 background8,49. Similarly, two sets of Savitri NILs (3 high yielding and 3 low yield-ing NILs) were developed from intercrossing of two BC1 mapping population (IR 77298-5-6-18/2*Savitri IR 74371-46-1-1/2*Savitri) to develop lines carrying qDTY3.2 and qDTY12.1 QTLs in Savitri background10,11,13. All of the NILs used in the present study in all the three background were re-confirmed with earlier reported markers for the presence of QTL and QTL combinations.

Figure 3. Cyclic illustrations of epistatic QTLs of Savitri NILs for grain yield under various levels of drought stress (a) GYSS: chromosome 2 (RM71-RM290,50 cM) and chromosome 3 (RM7-RM251 67 cM, RM319-RM55 167 cM and RM55-RM570,182 cM) have shown epistatic interaction with background loci of chromosomes 3,7 and 11 (b) GYMS: chromosome 4(RM537-RM335) and chromosome 10(RM244-RM239) was interacted negatively with a common genetic loci (RM136-RM275) on chromosome 6 (c,d) GYMS and GYSS of high yielding NILs: no significant interactions was found (e) GYNS: no interaction.

1 0Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

Evaluation of NILs under reproductive-stage drought stress and non-stress conditions. To study the grain yield variability of NILs pyramided with same QTL and QTL combinations, a total of 33 NILs derived from IR64, 5 NILs from TDK1-Sub1 and 6 NILs from Savitri backgrounds were evaluated during the dry seasons (DS) of 2015 and 2017 in a transplanted lowland ecosystem under reproductive-stage drought stress (RS) and non-stress (NS) conditions at the Ziegler Experiment Station of IRRI, Los Baños, Laguna, Philippines (14 °30′N longitude 121 °15′E latitude). Seedlings were raised in a wet-bed nursery for 21 days before transplanting into the main experimental field. The field lay out was generated using the statistical software Plant Breeding Tools (PBTools v1.4) developed at IRRI. The experimental design was alpha lattice in 2 replications in two-row plots of 5 m length with row to row and plant to plant spacing of 20 × 20 cm. Continuous supply of irrigation at 5 cm standing water was provided till the crop maturity in non-stress experiment. Fertilizer nitrogen, phos-phorus, and potassium (NPK) were applied at recommended dose of 120:30:30 kg ha−1. In RS trials, two doses of nitrogen fertilizer (basal and first split) were applied before initiating stress and a third dose was applied in adjustment with life-saving irrigation.

Reproductive stage drought stress experiment was conducted as described previously4,6. In brief for RS exper-iment normal irrigation was maintained in the field up to 30 days after transplanting, after which the water was drained out to initiate stress and continued till the harvesting. Water table depth was measured using PVC pipe with a plugged hole at the bottom and inserted in soil to 1 m depth and 15 cm above the soil surface in exper-imental stress field at regular interval. Life-saving irrigation was provided at level of severe stress when all the susceptible checks showed severe leaf rolling with minimum probability to recover upon watering and water table remained below 1 meter. Life-saving irrigation provided through flash flooding was drained out after 24 hrs to initiate further cycle of stress.

The stress imposed was classified from severe to moderate based on the relative yield reduction to non-stress (NS) experiments explained earlier by42,55. Experiments with yield reduction of <70% were characterized as severe stress (SS) and experiments with yield reduction of 31–70% were characterized as moderate stress (MS) as earlier reported by Kumar et al.55.

Observations recorded. Data on days to 50% flowering (DTF), plant height (PH), and grain yield (GY) were recorded in RS and NS experiments. DTF was observed as the number of days from seeding to the 50% the plants in a plot flower. The PH of three random plants from each NILs was measured before maturity from the soil surface to the tip of the main tiller and then averaged for mean analysis. Grains from individual plot were harvested at physiological maturity, oven-dried to a moisture content of 14% before weighing6,9. The grain yield in grams was further converted in kg ha−1 before entry means analysis in PBTools.

NILs genotyping. The presence of introgressed drought QTLs and QTLs combinations in IR64 (qDTY1.1, qDTY1.2, qDTY2.2, qDTY2.3, qDTY3.2, qDTY4.1 and qDTY12.1) TDK1-Sub1 (qDTY3.1, qDTY6.1, and qDTY6.2) and Savitri backgrounds (qDTY3.2 and qDTY12.1) were confirmed in respective NILs using the earlier reported peak and flanking markers of respective grain yield drought QTLs4,6–9. Fresh and young leaves from 25 days old trans-planted seedlings were collected in plot-wise bulk from each NIL- and their respective parents. Genomic DNA was extracted from collected leaves using a modified CTAB protocol developed by Murray and Thompson56, dissolved in 200 μl of TE (Tris-EDTA) buffer and stored at −20 °C. PCR amplification with the molecu-lar markers was performed on a thermal cycler (G-Storm GS1, UK) at Genotyping service laboratory (GSL), IRRI, Philippines. Each PCR reaction mixture of 15 μl consisted of 10 ng of rice genomic DNA, 1 × PCR buffer, 100 μM dNTPs, 100 μM oligonucleotide primers and 1 unit of Taq polymerase. PCR products were resolved cast-ing high-resolution 8% (v/v) polyacrylamide gel electrophoresis (PAGE) (CBS scientific, model MGV-202–33) and resolved in a 1x TBE buffer at 90 volts for 1–2 h, depending on the PCR product sizes. The separated DNA fragments after electrophoresis were stained with SYBER Safe™ and visualized under UV trans-illuminator (AlphaImager™ System).

To study the genomic reconstitution of NILs with the similar QTL or QTLs combinations but showing vari-ation in grain yield under drought, background genotyping was carried out on the selected NILs and parents. A total of 600 SSR markers equally distributed on all 12 chromosomes were selected from the Gramene database (http://www.gramene.org/) for identification of polymorphic markers between parents, which were further uti-lized for background selection in order to estimate recurrent parent genome (RPG) recovery.

Statistical analysis, recurrent parent genome recovery and epistatic interaction analysis. The phenotypic data collected from drought stress and non-stress experiments were statistically analyzed for the computation of trial means and standard error of difference (SED) using PBTools v1.4. Least significant difference (LSD) “at the 5% and 1% significant levels” were used to compare the means of the test entries and to infer the significant differences of the traits studied between parents and each NILs. Linear mixed model for analysis of variance was calculated with following equation:

= + + + +µYijk Gi Rj BK Rj eijk( )

where, Yijk is measurement recorded in plot, µ is overall mean, Gi is effect of ith genotype, Rj is the effect of the jth replicate, BK (Rj) is the block effect of jth replicate and eijk is the error. The genotypes were considered as fixed and the replicates and block effects were random for estimating the entry means. The recurrent parent genome (RPG) recovery percentage of drought NILs of various background were estimated with the following formula57

=+

×N

Recurrent parent genome (RPG) recoery % 2(B) (H)2

100

1 1Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

Where, B=SSR marker loci homozygous for genetic background; H = no. of marker loci still in heterozygous; N = total no. of polymorphic SSR markers used for background estimation.

Graphical representation of background genome was constructed using software Graphical Genotype (GGT v 2.0)58 (https://www.plantbreeding.wur.nl/).

Epistatic (digenic) interactions between markers loci of all the NILs including high and low yielding NILs background wise were analyzed considering all the conditions (SS, MS and NS) in model using QTL IciMapping59 ver. 4.0.1. In epistasis mapping method, a two-stage stepwise regression strategy was adopted to identify the most significant markers and marker-pair followed by two-dimensional scanning to identify significant digenic epistasis using adjusted phenotypic values based on best fitted multiple regression model60. In the first stage, the significant markers and marker pairs explaining epistatic and additive variations were selected in model and then in next stage the stepwise regression with stricter probability levels were applied to the residuals from the first stage to select significant marker pairs and estimate their effects in model50. The stricter probability level was applied in the second stage to avoid over-fitting due to the large number of regression variables. Extensive simulations have shown epistatic QTLs and QTL epistatic networks can be identified efficiently by ICIM mapping even in case two interacting QTL have some additive effects60. The threshold LOD value was determined by a permutation test involving 1000 runs at a significance level of P = 0.01 to detect significant digenic interactions between marker loci. In mapping parameters, window size and walk speed used for the genome scan was 10 cM and 1 cM respectively.

Data AvailabilityThe data sets supporting the results of this article are included within the article.

References 1. Serraj, R. et al. Drought resistance improvement in rice: An integrated genetic and resource management strategy. Plant Prod. Sci.

14(1), 1–14 (2011). 2. Kumari, S. et al. Transcriptome map for seedling stage specific salinity stress response indicates a specific set of genes as candidate

for saline tolerance in Oryza sativa L. Funct. Integr. Genomic 9, 109–123 (2009). 3. Khan, S., Javed, M. A., Jahan, N. & Manan, F. A. A short review on the development of salt tolerant cultivars in rice. Int J Public

Health Sci (IJPHS) 5(2), 201–212 (2016). 4. Vikram, P. et al. qDTY1.1, a major QTL for Rice grain yield under reproductive-stage drought stress with a consistent effect in

multiple elite genetic backgrounds. BMC Genet. 12, 89 (2011). 5. Ghimire, K. H. Identification and mapping of QTL (qDTY1.1) with a consistent effect on GY under RS. Field Crops Res. 131, 88–96

(2012). 6. Venuprasad, R. et al. Identification and characterization of large-effect quantitative trait loci for grain yield under lowland drought

stress in rice using bulk-segregant analysis. Theor. Appl. Genet. 120, 177–190 (2009). 7. Swamy et al. Genetic, physiological, and gene expression analyses reveal multiple QTL enhance the yield of rice mega-variety IR64

under drought. PLoS One 8, e62795 (2013). 8. Dixit, S. et al. Multiple major QTL lead to stable yield performance of rice cultivars across varying drought intensities. BMC Genet.

15, 16 (2014). 9. Bernier, J., Kumar, A., Venuprasad, R., Spaner, D. & Atlin, G. N. A large-effect QTL for grain yield under reproductive-stage drought

stress in upland rice. Crop Sci. 47, 507–516 (2007). 10. Yadaw, R. B. et al. A QTL for high grain yield under lowland drought in the background of popular rice variety Sabitri from Nepal.

Field Crops Res. 144, 281–287 (2013). 11. Mishra, K. K. et al. qDTY 12.1: a locus with a consistent effect on grain yield under drought in rice. BMC Genet. 14, 12 (2013). 12. Shamsudin, N. A. A. et al. Pyramiding of drought yield QTLs into a high-quality Malaysian rice cultivar MRQ74 improves yield

under reproductive stage drought. Rice 9, 21 (2016). 13. Dixit, S., Yadaw, R. B., Mishra, K. K. & Kumar, A. Marker-assisted breeding to develop the drought-tolerant version of Sabitri, a

popular variety from Nepal. Euphytica 213, 184 (2017a). 14. Kumar, A. et al. Marker-assisted selection strategy to pyramid two or more QTLs for quantitative trait-grain yield under drought.

Rice 11, 35 (2018). 15. Septiningsih, E. M. et al. Development of Submergence-tolerant rice cultivars: the Sub1 locus and beyond. Ann. Bot. 103, 151–60

(2009). 16. Pang, Y. et al. Simultaneous improvement and genetic dissection of salt tolerance of rice (Oryza sativa L.) by designed QTL

pyramiding. Front. Plant Sci. 8, 1275 (2017). 17. Yang, T. et al. Identification and pyramiding of QTLs for cold tolerance at the bud bursting and the seedling stages by use of single

segment substitution lines in rice (Oryza sativa L.). Mol. Breed. 36, 96 (2016). 18. Singh, V. K. et al. Marker-assisted simultaneous but stepwise backcross breeding for pyramiding blast resistance genes Pi2 and Pi54

into an elite Basmati rice restorer line PRR78. Plant Breed. 132(5), 486–495 (2013). 19. Khanna, A. et al. Development and evaluation of near-isogenic lines for major blast resistance gene(s) in Basmati rice. Theor. Appl.

Genet. 128(7), 1243–1259 (2015). 20. Fukuoka, S. et al. Gene pyramiding enhances durable blast disease resistance in rice. Sci. Rep. 5, 7773 (2015). 21. Suh, J. P. et al. Development of breeding lines with three pyramided resistance genes that confer broad-spectrum bacterial blight

resistance and their molecular analysis in rice. Rice 6, 5 (2013). 22. Pradhan, S. K. et al. Pyramiding of three bacterial blight resistance genes for broad-spectrum resistance in deepwater rice variety,

Jalmagna. Rice 8(1), 51 (2015). 23. Das, G., Rao, G. N. J., Varier, M., Prakash, A. & Prasad, D. Improved Tapaswini having four BB resistance genes pyramided with six

genes/QTLs, resistance/tolerance to biotic and abiotic stresses in rice. Sci. Rep. 8, 2413 (2018). 24. Wang, Y. et al. Marker assisted pyramiding of Bph6 and Bph9 into elite restorer line 93–11 and development of functional marker for

Bph9. Rice 10, 51 (2017). 25. Jena, K. K., Hechanova, S. L., Verdeprado, H., Prahalada, G. D. & Kim, S. R. Development of 25 near-isogenic lines (NILs) with ten

BPH resistance genes in rice (Oryza sativa L.): Production, resistance spectrum, and molecular analysis. Theor. Appl. Genet. 1, 16 (2017).

26. Xu, Y. & Crouch, J. H. Marker-assisted selection in plant breeding: from publications to practice. Crop Sci. 48, 391–407 (2008). 27. Bovill, W. D. et al. Pyramiding QTL increases seedling resistance to crown rot (Fusarium pseudograminearum) of wheat (Triticum

aestivum). Theor. Appl. Genet. 121(1), 127–136 (2010). 28. Sandhu, N. et al. Positive interactions of major-effect QTLs with genetic background that enhances rice yield under drought. Sci.

Rep. 8, 1626 (2018).

1 2Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

29. Wijerathna, Y. M. Marker assisted selection: biotechnology tool for rice molecular breeding. Adv. Crop Sci. Technol. 3, 187 (2015). 30. Yanoria, M. J. et al. Development of near-isogenic lines of Japonica-type rice variety Lijiangxintuanheigu as differentials for blast

resistance. Breed. Sci. 60, 629–638 (2010). 31. Xiao, C. et al. Development and evaluation of near-isogenic lines for brown planthopper resistance in rice cv. 9311. Sci. Rep. 6, 38159

(2016). 32. Mei, H. W. et al. QTLs influencing panicle size detected in two reciprocal introgressive line (IL) populations in rice (Oryza sativa L.).

Theor. Appl. Genet. 112, 648–656 (2006). 33. Zheng, T. Q. et al. Genetic effects of background-independent loci for grain weight and shape identified using advanced reciprocal

introgression lines from Lemont/Teqing in rice (Oryza sativa L.). Crop Sci. 51, 2025 (2011). 34. Li, Z. K. et al. Epistatic for three grain yield components in rice (Oryza sativa L.). Genet. 145, 453–465 (1997). 35. Holland, J. B. Epistasis and plant breeding. Plant Breed. Rev. 21, 27–92 (2001). 36. Mackay, I. J. et al. An eight-parent multiparent advanced generation inter-cross population for winter-sown wheat: creation,

properties, and validation. G3 (Bethesda) 4, 1603–1610 (2014). 37. Falconer, D. S. Introduction to Quantitative Genetics. John Wiley and Sons, New York, 340 pp (1989). 38. Bocianowski, J. Using weighted multiple linear regressions in estimation of QTL-by-QTL epistasis effects. Genet. Mol. Biol 35,

802–809 (2012). 39. Eshed, Y. & Zamir, D. Less-than-additive epistatic interactions of quantitative trait loci in tomato. Genet. 143, 1807–1817 (1996). 40. Mao, D., Liu, T., Xu, C., Li, X. & Xing, Y. Epistasis and complementary gene action adequately account for the genetic bases of

transgressive segregation of kilo-grain weight in rice. Euphytica 180, 261–271 (2011). 41. Gowda, S. J. M. et al. Mapping of QTLs governing agronomic and field traits in chickpea. J. Appl.Genet. 52, 9–21 (2011). 42. Dixit, S. et al. Increased drought tolerance and wider adaptability of qDTY12.1 conferred by its interaction with qDTY2.3 and qDTY3.2.

Mol. Breed. 30, 1767–1779 (2012). 43. Vikram, P. et al. Linkages and interaction analysis of major effect drought grain yield QTLs in rice. PLoS One 11(3), e0151532 (2016). 44. Caicedo, A. L., Stinchcombe, J. R., Olsen, K. M., Schmitt, J. & Purugganan, M. D. Epistatic interaction between Arabidopsis FRI and

FLC flowering time genes generates a latitudinal cline in a life history trait. Proc. Natl. Acad. Sci. 101, 15670–15675 (2004). 45. Rahman, M. A. et al. Mapping QTLs using a novel source of salinity tolerance from Hasawi and their interaction with environments

in rice. Rice 10, 47 (2017). 46. Liao, Y. C., Wu, P., Hu, B. & Yi, K. K. Effects of genetic background and environment on QTLs and epistasis for rice (Oryza sativa L.)

panicle number. Theor. Appl. Genet. 103, 104–111 (2001). 47. Basavaraj, S. H. et al. Marker assisted improvement of bacterial blight resistance in parental lines of Pusa RH10, a superfine grain

aromatic rice hybrid. Mol. Breed. 26, 293–305 (2010). 48. Ahmed, F. et al. Recurrent parent genome recovery in different populations with the introgression of Sub1 gene from a cross between

MR219 and Swarna-Sub1. Euphytica 207, 605–618 (2016). 49. Dixit, S. et al. Combining drought and submergence tolerance in rice: marker-assisted breeding and QTL combination effects. Mol.

Breed. 37, 143 (2017b). 50. Koide, Y. et al. A set of near-isogenic lines for blast resistance genes with an Indica-type rainfed lowland elite rice (Oryza sativa L.)

genetic background. Field Crop Res. 123, 19–27 (2011). 51. Singh, V. K. et al. Incorporation of blast resistance into ‘PRR78’, an elite Basmati rice restorer line, through marker assisted backcross

breeding. Field Crops Res. 128, 8–16 (2012). 52. Wang, Z. & Zhang, J. Impact of gene expression noise on organismal fitness and the efficacy of natural selection. Proc. Natl. Acad.

Sci. 108, E67–E76 (2011). 53. Thomson, M. J. et al. Characterizing the Saltol quantitative trait locus for salinity tolerance in rice. Rice 3, 148–160 (2010). 54. Babu, N. N. et al. Marker aided incorporation of Saltol, a major QTL associated with seedling stage salt tolerance, into Oryza sativa

‘Pusa basmati 1121’. Front. Plant Sci. 8, 41 (2017). 55. Kumar, A. et al. Yield and yield-attributing traits of rice (Oryza sativa L.) under lowland drought and suitability of early vigor as a

selection criterion. Field Crops Res. 114, 99–107 (2009). 56. Murray, M. G. & Thompson, W. F. Carnegie Inst. Wash. Year Book 76, 255–259 (1980). 57. Sundaram, R. M. et al. Marker assisted introgression of bacterial blight resistance in Samba Mahsuri, an elite indica rice variety.

Euphytica 160, 411–422 (2008). 58. Van Berloo, R. GGT 2.0: versatile software for visualization and analysis of genetic data. J. Hered. 99(2), 232–236 (2008). 59. Wang, J., Li, H., Zhang, L., & Meng, L. Users’ manual of QTL IciMapping. The Quantitative Genetics Group, Institute of Crop

Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100081, China, and Genetic Resources Program, International Maize and Wheat Improvement Center (CIMMYT), Apdo. Postal 6–641, 06600 Mexico, D.F., Mexico (2014).

60. Li, H., Ribaut, J. M., Li, Z. & Wang, J. Inclusive composite interval mapping (ICIM) for digenic epistasis of quantitative traits in biparental population. Theor. Appl. Genet. 116, 243–260 (2008).

AcknowledgementsThe authors are thankful for the support received from genotyping service laboratory (GSL)-IRRI during the process of genotyping and technical support from M.T. Sta Cruz and P. Maturan in conducting the field experiments. This study was supported by the Bill and Melinda Gates Foundation through STRASA (Stress Tolerant Rice for Africa and South Asia) project. Authors thank BMGF for the financial support to the study.

Author ContributionsS.Y. was involved in genotyping, observations, experimental data analysis and drafting the manuscript; N.S. was involved interpretation of data and revising the manuscript; R.R.M. helped in genotyping; S.D., S.K., S.P.S., N.P.M., S.P.D., R.B.Y., V.K.S., P.S., R.K.V. were involved in drafting and correction of the manuscript R.K.V. and V.K.S. was also involved with conceptualization of the result from present study on other crops. A.K. was involved in the design of the experiment and in the critical revision of the manuscript. All authors approved the final version of the manuscript.

Additional InformationSupplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-39084-7.Competing Interests: The authors declare no competing interests.Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

13Scientific RepoRts | (2019) 9:2616 | https://doi.org/10.1038/s41598-019-39084-7

www.nature.com/scientificreportswww.nature.com/scientificreports/

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 Cre-ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per-mitted 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 license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2019