Identification of Major Effect QTLs for Agronomic Traits and CSSLs in Rice from Swarna/Oryza nivara Derived Backcross Inbred Lines

Backcross inbred lines (BILs) derived from elite x wild crosses are very useful for basic studies and breeding. The aim of this study was to map quantitative trait loci (QTLs) associated with yield and related traits and to identify chromosomal segment substitution lines (CSSLs) from unselected BC2F8 BILs of Swarna/Oryza nivara IRGC81848. In all, 94 BILs were field evaluated in 2 years (wet seasons, 2014 and 2015) for nine traits; days to 50% flowering, days to maturity (DM), plant height (PH), number of tillers, number of productive tillers, panicle weight, yield per plant, bulk yield, and biomass. BILs were genotyped using 111 polymorphic simple sequence repeats distributed across the genome. Fifteen QTLs including 10 novel QTLs were identified using composite interval mapping, Inclusive composite interval mapping and multiple interval mapping (MIM). O. nivara alleles were trait-enhancing in 26% of QTLs. Only 3 of 15 QTLs were also reported previously in BC2F2 of the same cross. These three included the two major effect QTLs for DM and PH detected in both years with 13 and 20% phenotypic variance. Further, a set of 74 CSSLs was identified using CSSL Finder and 22 of these showed significantly higher values than Swarna for five yield traits. CSSLs, 220S for panicle weight and 10-2S with consistent high yield in both years are worthy of large scale field evaluation. The major QTLs and 22 significantly different CSSLs are a useful resource for rice improvement and dissecting yield related traits.

Backcross inbred lines (BILs) derived from elite x wild crosses are very useful for basic studies and breeding. The aim of this study was to map quantitative trait loci (QTLs) associated with yield and related traits and to identify chromosomal segment substitution lines (CSSLs) from unselected BC 2 F 8 BILs of Swarna/Oryza nivara IRGC81848. In all, 94 BILs were field evaluated in 2 years (wet seasons, 2014 and 2015) for nine traits; days to 50% flowering, days to maturity (DM), plant height (PH), number of tillers, number of productive tillers, panicle weight, yield per plant, bulk yield, and biomass. BILs were genotyped using 111 polymorphic simple sequence repeats distributed across the genome. Fifteen QTLs including 10 novel QTLs were identified using composite interval mapping, Inclusive composite interval mapping and multiple interval mapping (MIM). O. nivara alleles were trait-enhancing in 26% of QTLs. Only 3 of 15 QTLs were also reported previously in BC 2 F 2 of the same cross. These three included the two major effect QTLs for DM and PH detected in both years with 13 and 20% phenotypic variance. Further, a set of 74 CSSLs was identified using CSSL Finder and 22 of these showed significantly higher values than Swarna for five yield traits. CSSLs, 220S for panicle weight and 10-2S with consistent high yield in both years are worthy of large scale field evaluation. The major QTLs and 22 significantly different CSSLs are a useful resource for rice improvement and dissecting yield related traits.

INTRODUCTION
Rice (Oryza sativa L.) is the most important cereal food crop cultivated worldwide. The rapid growth in the world's population, which is expected to reach 9.1 billion by 2050 demands global rice production to be doubled (Ray et al., 2013). Over the years, intensive breeding to develop high yielding varieties quickly using few elite parents has resulted in loss of genetic variability in the cultivars. On the other hand, land races and wild species of rice are sources of abundant genetic variation that can be tapped to increase rice yield sustainably. The wild species of Oryza are an important source of genetic variability for tolerance to biotic and abiotic stresses and for improvement of yield as well (Brar and Khush, 1997;Brar and Singh, 2011;Singh K. et al., 2016). Tanksley and Nelson (1996) proposed advanced backcross quantitative trait loci (AB-QTL) strategy to simultaneously identify and transfer valuable alleles from unadapted germplasm into the elite cultivars and it has been followed in several crops. The annual wild species, O. nivara adapted to seasonally dry habitats is the closest progenitor of O. sativa (Sharma and Shastry, 1965). It is also a potential source of favorable alleles for agriculturally important traits. One O. nivara accession IRGC101508 from Uttar Pradesh, India was identified as the only accession resistant to grassy stunt virus after screening 5000 accessions and 1000 breeding lines (Khush, 1977). O. nivara also contributed resistance to bacterial blight (Cheema et al., 2008) and blast (Eizenga et al., 2013). Li et al. (2006) mapped QTLs for domestication traits in F 2 population derived from CL16/O. nivara IRGC80470. QTLs for seedling vigor, yield and quality traits were reported from an advanced backcross population derived from M-202/O. nivara IRGC100195 (Eizenga et al., 2015). More recently, a set of 131 ILs were developed from a cross between 93-11 and O. nivara W2014 and 65 QTLs were identified for 13 agronomic traits using whole genome resequencing (Ma et al., 2016).
At Indian Institute of Rice Research, two accessions of O. nivara IRGC81832 and IRGC81848 which were genetically distinct from 22 other accessions were extensively used for QTL mapping of yield in BC 2 F 2 and quality traits in BC 2 F 3 seed (Sarla et al., 2003). Kaladhar et al. (2008) identified 17 major effect QTLs for different yield traits in BC 2 F 2 population derived from Swarna/O. nivara IRGC81832 including yldp8.1 with LOD score of 8.76 which increased yield by 5.8 g per plant and grain number by 426 grains per plant. Significant yield enhancing QTLs qyldp2.1, qyldp3.1, qyldp8.1, qyldp9.1, qyldp11.1 were reported from the other BC 2 F 2 population derived from Swarna/O. nivara IRGC81848 (Swamy, 2009;Swamy et al., 2014). QTLs for stem diameter qSD7.2, qSD8.1, qSD9.1, rachis diameter qRD9.1 and number of secondary branches qNSB1.1 were identified as good targets for use in MAS . QTLs qmp1.2, qkw3.1, qkw6.1, qklac12.1 for grain quality traits such as milling percentage, kernel width and kernel length after cooking were also identified in these two populations (Swamy et al., 2012). One IL IET21542 (RPBio4918-248) derived from Swarna/O. nivara IRGC81848 gave mean yield of 5.5 t/ha for three consecutive years across several locations and a maximum yield of 10.6 t/ha in Coimbatore in wet season 2011 (Annual Progress Report, 2012. It was released as DRR Dhan 40 for three important rice growing states of India -West Bengal, Maharashtra and Tamil Nadu in 2013 (Sarla, 2014). This variety has two QTLs yld9.1 for yield and nfg9.1 for number of filled grains from O. nivara. Also, six ILs 212S, 215S, 221S, 224S, 228S, and 230S from the same cross were BPH resistant (Lakshmi et al., 2010) and 228S in particular was resistant to multiple pests when evaluated in 10 green house conditions and 46 field tests against 11 pests in multiple resistance screening trials (MRST) (Annual Progress Report, 2012. Two ILs 166S and 75S were tolerant to drought and salinity (Rai et al., 2010). Selected ILs from these two populations were analyzed for G × E interaction and two stable ILs 166S and 14S were identified (Divya et al., 2016). The usefulness of O. nivara in improving varieties for yield and other traits is thus quite evident.
Chromosome segment substitution lines (CSSLs), which carry a single or few chromosome segments from the donor in the genetic background of the recurrent parent with the whole donor genome, are ideal for QTL mapping/cloning and also as a genetic resource for pyramiding target segments and breeding (Ali et al., 2010). Several CSSLs have been developed in rice and characterized for yield related traits (Ando et al., 2008;Ookawa et al., 2016). Most of the CSSLs used in these reports were derived from crosses within O. sativa. The other cultivated species of rice O. glaberrima can also be used to improve O. sativa (Sarla and Swamy, 2005). A set of 34 CSSLs were developed from O. glaberrima and 105 QTLs for 10 yield traits related to grain yield, plant stature and maturity were identified (Shim et al., 2010).
Chromosomal segment substitution lines have also been developed using wild rice species. CSSLs using O. rufipogon helped identify QTLs for yield and related traits (Qiao et al., 2016). A total of 40 CSSLs were developed from O. longistaminata in the background of Taichung 65 and evaluated for yield traits (Ramos et al., 2016). Twenty six CSSLs harboring O. nivara genomic segments in the genetic background of Koshihikari were developed and evaluated for agriculturally important traits . A set of 131 ILs carrying, a total of 767 chromosomal segments from O. nivara (W2014) in the genetic background of 93-11 were recently reported (Ma et al., 2016). The O. nivara allele conferred positive effects at 37% of yield associated QTLs. A new gene-specific InDel marker LQ30 for a gene LOC_Os03g14850 for improved stigma length was developed using a single segment substitution line SSSL14 (Liu et al., 2015). Thus CSSLs can be easily used for gene discovery. We also wanted to know if the major QTLs detected in BC 2 F 2 remain major effect when the BC 2 F 2 are advanced to BC 2 F 8 or do new QTLs become major QTLs with the change in background from BC 2 F 2 to BC 2 F 8 of the same cross. Keeping in view, the importance of ILs as a prebreeding material harboring loci for yield enhancing QTLs/genes and CSSLs as a genetic resource for crop improvement, the present study was aimed to (i) detect major QTLs for yield traits in BC 2 F 8 population derived from Swarna/Oryza nivara (ii) compare them with those identified in BC 2 F 2 from the same cross and (iii) identify a set of complete chromosomal segment substitution lines from BC 2 F 8 BILs.

Plant Material
A set of BC 2 F 5 BILs, derived from a cross between an elite rainfed lowland cultivar Swarna (O. sativa) also known as MTU 7029 as a recurrent parent and a wild accession O. nivara IRGC81848 as a donor parent was developed by Swamy (2009) and these BILs were self-pollinated in consecutive generations to obtain BC 2 F 8 families by single panicle selection (Supplementary Table 1).

Phenotypic Evaluation
Ninety four BC 2 F 8 lines along with recurrent parent Swarna were grown in wet seasons (

Trait Association
Correlation was calculated based on two sample t-test with equal variances. CORREL function of the analysis tool pack was used and then imported in excel to find the correlation coefficient between two variables. Correlation among traits was computed at P < 0.05 and P < 0.01, respectively.

Genotyping
Total genomic DNA was isolated from fresh leaf samples of 94 BC 2 F 8 plants and the parents following CTAB method (Doyle and Doyle, 1987). In all, 324 markers were tested for parental polymorphism. One hundred and twenty four SSRs were polymorphic of which 111 gave clear bands and also segregated in BC 2 F 8 (Supplementary Table 2). PCR was carried out in thermal cycler (G-STORM, United States) with a final 1 http://bbi.irri.org/products reaction volume of 10 µl containing 15 ng of genomic DNA, 1X assay buffer, 200 µM of dNTPs, 1.5 mM MgCl 2 , 10 pmol of forward and reverse primer and 1 unit of Taq DNA polymerase (Thermo Scientific). PCR cycles were programmed as follows: initial denaturation at 94 • C for 5 min followed by 35 cycles of 94 • C for 30 s, 55 • C for 30 s, 72 • C for 1 min and a final extension of 10 min at 72 • C. Amplified products were resolved in 3% agarose gel prepared in 0.5 × TBE buffer and electrophoresed at 120 V for 2 h. Gels were stained with ethidium bromide and documented using gel documentation system (Alpha Imager, United States).

Linkage Mapping and QTL Analysis
Linkage map was constructed based on genotypic data of 94 BILs using 111 polymorphic SSR markers on all chromosomes using MAP function (BC 2 RIL) of QTL IciMapping v4.1 2 using the Kosambi mapping function (Kosambi, 1944 (Joehanes and Nelson, 2008). Analysis was undertaken using automatic parameter setting and controlling marker forward stepwise. The threshold of LOD for declaring the presence of significant QTL for each trait was determined using 1000 permutations and α = 0.05 in QGene 4.4.0 software. At α = 0.05, the LOD threshold values ranged from 2.8 to 8.8.

Identification of CSSL Set
Chromosomal segment substitution lines were identified using genotypic data of 111 polymorphic loci in 94 BILs in the background of the recurrent parent Swarna using the software CSSL Finder 3 . Statistical analysis was carried out using PB tools (Version 1.4 4 ) for the test of significance.

Phenotypic Evaluation
The mean values of parents and BILs in each year for nine traits and number of BILs showing significant trait increase over recurrent parent Swarna is given in

Trait Correlations
Significant correlations were observed among the traits in both the years at P < 0.05 and P < 0.01 (Supplementary Table 3). In both years, significant positive correlation was observed between DFF and DM; PH with PW, YLDP and BM; NT and NPT; PW and YLDP; YLDP with BY and BM; and BM and BY. Significant negative correlation was observed between PH and NPT in both years. In 2014, DFF was positively correlated with NT and YLDP, but negatively with BY in 2015. NT and NPT showed significant positive correlation with YLDP, BY, and BM only in 2014. In 2015, positive significant correlation was observed between PH and BY; PW and BY, but were negatively correlated in 2014. Negative correlation was observed for DM with PW and BY; PH with NT in both years, but significant only in 2015. Correlation analysis of both years mean data showed that YLDP has significant positive correlation with PH, NT, NPT, PW, BM, and BY. Significantly correlated traits DFF and DM were negatively correlation with both PW and BY. PH had significant positive correlation with PW, YLDP, and BM but negative correlation with NT and NPT. BM also showed significant positive correlation with PH, NT, NPT, and PW.

Molecular Characterization of BILs
Considering all the 111 loci and 94 BILs, a total of 7488 alleles were detected covering all chromosomes. Swarna homozygous alleles accounted for 84%, O. nivara homozygous alleles 11% and heterozygous alleles 4%. The highest number of O. nivara alleles were detected on chromosome 1 (128 alleles), chromosome 2 (120 alleles) and chromosome 8 (108 alleles) and lowest number of O. nivara alleles were detected on chromosome 10 (10 alleles). The O. nivara introgressions among BILs ranged from 1.15 to 27.0%, with a mean of 10.1%. Number of heterozygotes at any locus ranged from 0 to 40 (42%). RM125 on chromosome 7 followed by RM519 on chromosome 12 showed the highest number of heterozygotes. BIL 138S had was highest number of 24 heterozygotic loci (21% of all loci).

Construction of Linkage Map and QTL Mapping
Linkage map was constructed based on genotypic data of 94 BILs using 111 polymorphic SSR markers on all chromosomes. In all, 15 QTLs were identified for all traits except for BM, considering data of both years using CIM, ICIM, and MIM (Figure 1). Five QTLs were identified in 2014 and 10 QTLs in 2015 (Tables 2, 3). Two QTLs were common in both years for DM and PH considering CIM, ICIM and MIM. PV explained by these 15 QTLs ranged from 12 to 56%. O. nivara alleles were trait enhancing in 26% of QTLs for 3 traits.

Common QTLs Identified in 2014 and 2015
In all, 15 QTLs were identified, of which only 2 QTLs were common in both years and their PV ranged from 13 to 27%. These were qDM2.1 and qPH1.1. The common QTL qDM2.1 was identified at chromosomal region RM8080-RM6318 and explained PV of 13% with LOD score 2.9. A major QTL qPH1.1 at marker interval RM226-RM431 was identified in both years with average LOD of 4.5 explaining PV of 20%. O. nivara allele was trait-enhancing in qPH1.1 and Swarna allele at qDM2.1 in both years.

QTLs Identified in 2014
A total of 5 QTLs were identified, of which 2 were also identified in 2015. One QTL qPH1.2 was identified at RM128-RM226 region with PV of 13%. One QTL each for NPT (qNPT1.1) and PW (qPW9.1) were identified with PV 13% each and LOD score of 3 and 3.5, respectively. The increasing effect of these three QTLs qPH1.2, qNPT1.1 and qPW9.1 was from O. nivara.

QTLs Identified in 2015
In all, 10 QTLs were identified, of which 2 were also identified in 2014. The PV explained by these 10 QTLs ranged from 12 to 56%. Three QTLs each were identified for YLDP and BY. Three QTLs qYLDP1.1, qYLDP5.1, and qYLDP12.1 showed PV ranging from 36 to 56% and trait enhancing alleles from Swarna. Of these three QTLs, one major effect QTL qYLDP12.1 was identified at chromosomal region RM519-RM235 with highest LOD of 16.7 and 56% PV. Three major QTLs qBY2.1, qBY6.1, and qBY8.1 were identified with PV ranging from 27 to 33%. One QTL each for DFF (qDFF2.1) and NT (qNT1.1) were identified with PV 12 and 22.5%. Two common QTLs qDM2.1 and qPH1.1 were identified in both years.

Identification of CSSLs
Chromosomal segment substitution line Finder output showed that 74 CSSLs had homozygous chromosome segments from O. nivara substituting Swarna segments (Figure 2). These 74  CSSLs together showed 89% coverage of the wild genome based on the 111 SSRs marker data. Small regions on chromosomes 4, 6, 7, 9, 11, and 12 were not represented in the 74 CSSLs. The average number of substituted segments per chromosome in these was 6.5. The number of substituted segments in each CSSL ranged from 2 to 18, with an average of 9.5% (Supplementary Figure 3). Chromosome 1 was represented by 8 CSSLs, chromosome 2 by 13 CSSLs, chromosome 3 by 9 CSSLs and chromosome 10 by 4 CSSLs. Chromosomal segment substitution line 142S had only three substituted segments and they were on chromosome 10 and showed significantly higher BM than Swarna. Two lines 220S and 166-23-1S showed significantly higher PW than Swarna in both the years.

DISCUSSION
In the present study, transgressive segregants were obtained with about 15% improvement over Swarna for many yield related traits which indicates that alleles from O. nivara were favorable in the genetic background of Swarna in BC 2 F 8 generation also. It is significant that 10 new QTLs were identified in BC 2 F 8 which were not detected previously in BC 2 F 2 population   2014). Of 15 QTLs identified in the study, three QTLs for PH, DM, and YLDP on chromosomes 1, 2 and 12, respectively, were identified in BC 2 F 2 also. Most of the QTL regions identified in our study were associated with two or more traits. For example, chromosomal region RM226-RM431 on chromosome 1 had QTLs for PH and NT. Likewise, the chromosomal region RM1220 on chromosome 1 had QTLs for NPT and YLDP; RM8080-RM6318 on chromosome 2 had QTLs for DFF and DM.
Correspondingly, there was also significant correlation between DFF and DM in both years; NPT and YLDP in 2015. This may be due to the pleiotropic effects of same QTL alleles/genes controlling PH, NPT and YLDP or presence of two or more adjacent QTL alleles at same locus each controlling a different trait. These QTLs can be used for marker assisted improvement of different traits simultaneously because of their strong and consistent linkage with yield. However, it may be noted that though QTLs for correlated traits were colocalised but the effect of QTLs depends on the genetic context and the direction of their association (Dufey et al., 2015).
Common QTLs in BC 2 F 2 and BC 2 F 8 Yield enhancing QTLs have been reported previously from O. nivara in BC 2 F 2 families (Kaladhar et al., 2008;Swamy et al., , 2014. A total of 28 QTLs for yield traits were reported in BC 2 F 2 population derived from the same cross of Swarna/O. nivara IRGC81848 and 78% of the loci from O. nivara were trait enhancing (Swamy et al., 2014). Three common QTLs for traits PH (qPH1.1-RM431) on chromosome 1, DM (qDM2.1-RM3874) on chromosome 2, and YLDP (qYLDP12.1-RM519) on chromosome 12 were identified in both BC 2 F 2 (previous study) and BC 2 F 8 (present study) and these three are thus significant major effect QTLs. The first common QTL qPH1.1 is close to the well known semi-dwarf locus sd-1, the green revolution gene and explained PV of 22 and 17.5% in 2014 and 2015, respectively. Thus, it shows the robustness of our study as these major QTLs were identified in BC 2 F 2 and BC 2 F 8 and in both years. The major effect QTL qPH1.1 was identified in both years and in both generations (BC 2 F 8 and BC 2 F 2 ) indicating that it is also a stable QTL. RM431 the peak marker is within the region of QTL qDTY1.1 reported previously for yield under drought and being transferred into submergence tolerant versions of three high yielding mega rice varieties Swarna-Sub1, Samba Mahsuri-Sub1, and IR 64-Sub1 using MAS (Vikram et al., 2011;Singh R. et al., 2016). In our study, O. nivara allele of qPH1.1 increases PH. The semi dwarfing gene sd-1, has been widely utilized in rice production, however, attention has to be paid to breed taller plants with strong culm to enhance plant yield by increasing biomass (Han et al., 2017). They identified QTLs for PH in a RIL population derived from Zhenshan 97 and Xizang 2. Two QTLs qph1 and qph7.1 detected in 3 years, explained 13% PV with large additive effect of 12 cm from Xizang 2 allele. In our study, O. nivara alleles increased PH by 11 cm in both the years. RM431 flanks PH QTLs reported previously from elite/ wild crosses using accessions of O. nivara, O. rufipogon and from landraces also (Wickneswari et al., 2012;Eizenga et al., 2013;Mohammadi et al., 2013). In addition to height, RM431 was linked with QTLs for grain number/panicle, harvest index (Li et al., 2012). This indicates that alleles at the vicinity of RM431 are associated with different yield related traits in different genetic backgrounds.
The second common QTL qDM2.1 at chromosomal region RM3874 falls within a metaQTL MQTL2.3 reported for PW . The third common QTL qYLDP12.1 at RM519 region also harbors QTLs qrd12.1 for rachis diameter, qnpt12.1 for number of productive tillers, qnsp12.1 for number of spikelets and qnfg12.1 for number of filled grains (Swamy, 2009;Swamy et al., 2014). There are two previous reports on identification of common QTLs in BC 2 F 2 and later generations. Rangel et al. (2008) identified one common QTL for grain yield on chromosome 1 at RM1 locus from BG90-2/RS-16 (O. glumaepatula) in both BC 2 F 2 and BC 2 F 8 . Similarly, Wickneswari and Bhuiyan (2014) identified two common QTLs qSPL-1-1 and qSPL-8 for spikelets per plant in BC 2 F 2 and BC 2 F 5 from MR219/Oryza rufipogon Griff. IRGC105491. QTL Regions in BC 2 F 8 Which Were Linked with Different Agronomic Traits in BC 2 F 2 Comparing QTLs in BC 2 F 8 with BC 2 F 2 we found there were few other chromosomal regions common in BC 2 F 8 and BC 2 F 2 but linked with different traits. In BC 2 F 8, RM434 was linked with qPW9.1 where as in BC 2 F 2 RM434 was linked with qyldp9.1. Likewise, RM128 was linked with qPH1.2 in BC 2 F 8 and qnpt1.1 in BC 2 F 2 ; RM223 for qBY8.1 in BC 2 F 8 and qnfg8.1and qyldp8.1 in BC 2 F 2 ; RM519 for qYLDP12.1 in BC 2 F 8 and qnpt12.1, qnsp12.1 and qnfg12.1 in BC 2 F 2 . There could be several reasons for this observation. Earlier, 100 SSRs were used to genotype 227 BC 2 F 2 mapping population whereas in this study 111 different SSRs were used to genotype 94 BC 2 F 8 BILs and only 37 markers were common with the previous study. The detection of novel QTLs in BC 2 F 8 might be due to changes in genetic background from BC 2 F 2 to BC 2 F 8 or differences in the size of mapping population and the distribution of SSR markers used in the two studies. Also environment conditions were different in the two generations and this might have contributed to differences in phenotype and QTLs.
Considering that the available population size was small, the phenotypic variation explained by QTLs is quite likely to be overestimated. However, since many of the major effect QTLs detected in our study were also identified for the same or different trait previously, it indicates the robustness and importance of these QTLs in regulating the phenotype of yield associated traits. The significant QTLs identified in this study can be further evaluated for use in marker assisted transfer to adapted varieties to improve not only yield but several related traits as QTLs for drought and salinity are being transferred into high yielding mega varieties (Singh R. et al., 2016).

Identification of CSSLs
Chromosomal segment substitution line are an important genetic resource in rice to discover novel genes by focusing on small chromosomal regions (Ali et al., 2010;Subudhi et al., 2015;Ramos et al., 2016). A set of 74 CSSLs with substituted chromosomal segments of O. nivara in the genetic background of Swarna were identified. Graphical genotypes of CSSL library revealed 89% coverage of O. nivara genome. A few chromosomal regions were not represented by substituted segments on chromosomes 4, 7, 9, and 11. This could be due to less number of BILs used or the presence of lethal alleles, hybrid sterility and gametophyte lethal genes or even low recombination in these chromosomal locations. Introgression of O. nivara chromosomal segments in the genetic background of 93-11 was reported by Ma et al. (2016). The coverage of the O. nivara genome by the ILs was 94.96%. Several sets of CSSLs have been developed using other wild species such as O. rufipogon (Furuta et al., 2014;Qiao et al., 2016;Ogawa et al., 2016), O. minuta (Guo et al., 2013), O. meridionalis (Arbelaez et al., 2015), O. longistaminata (Ramos et al., 2016). A set of 198 CSSLs was developed from a cross between 93-11 and O. rufipogon, and introgressed segments covered 84.9% of the wild rice genome (Qiao et al., 2016). Subudhi et al. (2015) developed 74 CSSLs covering 99% of the weedy rice donor PSRR-1 genome. In their study, donor segments per line ranged from 1 to 3 with 64% of CSSLs with single homozygous donor segments. In our study number of donor segments per line ranged from 3 to 18 as CSSLs were identified from a set of available BILs at BC 2 F 8 generation without MAS. A few backcrosses can be made along with MAS of target regions to recover more background genome of this set or selected CSSLs.
In our previous studies, two stable lines 166S and 14S were reported to be efficient in compartmentalization of Na+ in leaf tissue and grain yield of 166S was least affected by salt stress (Divya et al., 2016;Pushpalatha et al., 2016). Also, three lines 24S, 70S, 14-3S were identified as heat tolerant lines for spikelet fertility and YLDP in both wet and dry seasons (Prasanth et al., 2016). In present study, these five lines 14S, 24S, 70S, 166S and 14-3S were part of the set of CSSLs and they were found to have 11.5-18.6% chromosomal segments from O. nivara. Thus, several elite BILs such as these can be used in prebreeding programs for identification of candidate genes for different yield related traits. Among the CSSLs, 220S showed significantly higher PW over the parent Swarna in both years. 10-2S showed highest YLDP in both years over Swarna. These two lines, 220S and 10-2S showed significantly higher yield than Swarna, so they can be used to dissect target QTL regions. In addition, another 20 CSSLs showed different yield traits as significantly higher than Swarna and these marker defined CSSLs can be further utilized for fine mapping.

CONCLUSION
Our results provide evidence that O. nivara has novel, stable, major effect QTL alleles for PH, DM not only in BC 2 F 2 but even in BC 2 F 8 . ILs with O. nivara alleles for increasing PH or decreasing DM are potential donors for transfer into other popular lines. The set of 74 CSSLs is being phenotyped for other agronomic traits including disease and pest resistance and is an important genetic resource to discover novel alleles for several traits in rice.

ETHICS STATEMENT
The authors declare that the experiments comply with the current laws of the country in which they were performed and in compliance with ethical standards.