To develop and characterize introgression lines for leaf and neck blast resistance, 326 introgression lines were developed using various accessions of six different AA genome wild species in the genetic background of elite Indian varieties like PR114 and Pusa 44 and were screened for blast resistance. Stringent phenotyping coupled with genotyping using gene based markers led to the identification of four resistant introgression lines, which showed promising resistance and do not possess any of the tested genes. Furthermore, multi-location screening confirmed the field resistance of the four introgression lines to both leaf and neck blast. Molecular characterization of these introgression lines using genome-wide simple sequence repeat markers revealed the presence of small percentage of wild Oryza genome introgrssion. So these lines can be used for mapping and identification of novel leaf and neck blast resistance genes. Thus, these four introgression lines can be considered as new genetic resources for blast resistance.
Rice is the staple food for more than half of the world’ s population, and global rice demand is estimated to rise from 6.76 × 108 t in 2010 to 8.52 × 108 t in 2035 (Khush, 2013). To produce 1.76 × 108 t additional rice, it is need to increase the yield and also minimize the yield loss caused by various diseases and insect pests. Rice blast disease, caused by Magnaporthe oryzae, is a major constraint for sustainable rice production. To date, around 100 blast resistance genes have been identified, and many of them have been cloned and characterized, such as Pb1, Pia, Pib, Pid2, Pid3, Pik, Pik-h/Pi54, Pik-m, Pik-p, Pish, Pit, Pita, Piz-t, Pi1, Pi2/Piz-5, Pi5, Pi9, Pi21, Pi25, Pi36, Pi37, Pi35, Pi64, Pi56, Pi63 and PiCO39 (Devanna et al, 2014). Though many resistant varieties to M. oryzae have been developed, the resistance is not long lasting, because the high pathogen plasticity in the fields makes single resistance gene break down after three to five years of the cultivar release (Bonman et al, 1986; Lang et al, 2009). Hence, development of broad spectrum and durable blast resistant varieties is essential for combating this disease, which requires continuous efforts of breeders and pathologists. Wild species of Oryzacan be exploited to widen the gene pool of rice for biotic and abiotic stress (Ram et al, 2013). Various studies have demonstrated that wild species as a reservoir of hidden gene(s) can be used for crop improvement like grassy stunt virus resistance gene fromO. nivara (Brar and Khush, 1997), QTLs for bacterial leaf blight and brown planthopper resistance (Zhang and Xie, 2014), and yield enhancing QTLs from various wild species (Linh et al, 2008; Rangel et al, 2008; Chen et al, 2009). Though wild species have been used extensively for other important agronomic traits of rice, they have been rarely used in blast resistance breeding programmes, since only two resistance genes, Pi9and Pi40, were reported from wild species (Amante et al, 1992; Jeung et al, 2007). Thus, it is worth to develop introgression lines (ILs) from various wild species and characterize for blast resistance. Further, stable ILs will facilitate the use of wild species derived genes/alleles in breeding programmes (Lei et al, 2013; Zhou et al, 2014). The objective of this study was to identify novel genetic resources for blast disease resistance. Here we reported development and systematic screening of ILs across India for field resistance to rice blast. Since, field resistance allows effective control of a pathogen under natural field conditions and it is considered to be durable when exposed to new races of blast (Fukuoka and Okuno, 2001). We also made an attempt to identify the reported blast resistance (R) genes among the resistant ILs using gene based markers which are reported earlier (Hayashi et al, 2006; Lmam et al, 2014). In addition to the identification of blast R genes, we tried to identify the introgressed regions among the identified ILs using genome-wide simple sequence repeat (SSR) markers. SSR markers are most popular genetic markers with high polymorphism rate, high abundance and broad distribution throughout the genome, which are inherited in a Mendelian fashion as codominant markers (Miah et al, 2013). Also, the availability of many SSR markers in rice helps the breeders for crop improvement through molecular marker assisted breeding (McCouch et al, 2002). The strategic phenotyping for blast disease (field) resistance coupled with genotyping using markers has led to the identification of resistant ILs which are valuable genetic resources not only for blast resistance but also for novel gene identification leading to cultivar development (Ashikari and Matsuoka, 2006; Rahman et al, 2011).
The stable ILs was developed at Punjab Agricultural University, Ludhiana, India, with PR114 and Pusa 44 as recurrent parents and six wild species of Oryza (AA genome) as donors. These lines were developed through repeated backcrossing coupled with phenotypic selection for yield parameters. The donor species included 30 accessions of O. glaberrima, O. nivara, O. rufipogon, O. longistaminata, O. glumaepatulaandO. barthii(Supplemental Table 1). PR114 and Pusa 44 are the elite rice varieties with high yield and good cooking quality (Gaikwad et al, 2014).
From the 326 rice ILs (BC2-3F9) of PR114 and Pusa 44, single plant was selected based on the maximum features of its recurrent parent, and the seeds of single plant were used for further studies. Six monogenic lines carrying single blast R gene and one variety harbouring the rice blast R genes were used as resistant checks for identification of known gene allele among the four ILs (Table 1). The varieties viz, LTH, BPT5204, HR12 and Co39 were used as susceptible checks for gene profiling.
The seeds of ILs were sown on uniform blast nursery beds along with local susceptible checks HR12 (tall, low tillering variety with susceptibility to leaf blast), Co39 (universal blast susceptible variety), and BPT5204 (popularly known as Samba Mahsuri, mega variety with high yield and good quality, but susceptible to rice blast). Among the three checks, HR12 was used as spreader row. After 15 d of germination, the plants were inoculated with most virulent isolates of India i.e. NLR-1 maintained as DRR (Directorate of Rice Research) isolate (Padmavathi et al, 2005; Prasad et al, 2009; Ratna et al, 2011). The inoculum contained 50 000 conidia/mL with 0.5% glycerol and the nursery beds were maintained with 90% relative humidity. The ILs was extensively screened for blast resistance as follows: (i) Screened for resistance at the Indian Institute of Rice Research in three seasons (January to April in 2010 and 2011, and June to September in 2011), with three replications in each season; (ii) The selected ILs was screened at seven blast endemic areas to know their leaf blast resistance in the two seasons (June to September 2011 and 2012); (iii) The short listed ILs were further screened at hilly areas (which are hotspots for natural neck blast infestation in India) to check their resistance to neck blast (June to September in 2011 and 2012) (Supplemental Table 2); and (iv) To check the field resistance, the selected ILs were screened across different parts of the country through All India Co-ordinated Rice Improvement Programme (AICRIP) in 2012 and 2013 (http://www.drr icar.org/aicrip.htm). In AICRIP screening, the ILs were screened at 27 different centres for leaf blast and 8 centres for neck blast. The blast screenings of (ii), (iii) and (iv) were conducted in two replications per season and the mean scores were obtained. All the scores were obtained according to the standard evaluation system (Supplemental Table 3) (Panguluri et al, 2013). The selected four ILs were evaluated for the essential agronomic traits in the fields. The plants were grown at Indian Institute of Rice Research during Kharif in 2012 and 2013. Individual plants from the ILs were evaluated and the data for plant height, number of tillers per plant, number of panicles per plant, yield per plant and biomass were documented from 10 plants for each entry and the mean values were determined.
Genomic DNA was isolated using modified cetyltrimethyl ammonium bromide (CTAB) method (Saghai et al, 1984) and the quality of isolated DNA was checked using 0.8% agarose gel electrophoretically and quantified using the nanodrop (Thermo Fisher, USA). To identify the presence of blast resistance gene alleles among the ILs, gene profiling was carried out according to Lmam et al (2014). For this, seven effective blast resistance genes (Pi54, Piz-t, Pi1, Pi9, Pi2, Pita/Pita2and Pi40) were analyzed among blast resistant ILs using gene based markers (for Pi54, Piz-t, Pi9, Pita/Pita2, Pi40andPi2) as well as two flanking markers (for Pi1) (Table 1). PCR was performed using 1 U of Taq DNA polymerase and 10 × PCR buffer (Genei, India) in 15-µ L reaction volume with a thermal profile as described in Lmam et al (2014). The PCR products were resolved on agarose gel (4% for the products less than 500 bp and 2% for the products more than or equal to 1 kb) in 0.5 × TBE buffer, stained with ethidium bromide (0.5 pg/mL) and documented under UV light. Scoring was done based on presence (1) or absence (0) of the gene in comparison with the resistant and susceptible checks.
The four ILs (IL-1 to IL-4) were forwarded further for precise identification of the donor genome introgression. A set of 499 SSR markers (McCouch et al, 2002) which were spread uniformly across the twelve linkage groups were used. Scoring of the ILs was done by comparing the IL alleles with the recurrent parent allele and the donor allele.
The introgressed regions in the four ILs were analyzed based on the existing literature for the presence of meta QTLs (Ballini et al, 2008). Further, in silico analysis of the introgressed region was done by using RAP-DB (http://rapdb.dna.affrc.go.jp/). For this, the list of annotated genes between the introgressed regions were downloaded from the database and categorized them based on their molecular functions. Later the genes associated with disease resistance were listed.
The ILs were developed through backcrossing followed by recurrent selection. In brief, two or three backcrosses were carried out between donor and recurrent parents (PR114 and Pusa 44). After each backcross, the plants were selected based on phenotypic characters as well as yield related traits which were close to the recurrent parent. This is followed by continuous selfing so as to develop BC2-3F9 generation, at this stage, almost all the ILs were closely resemble to their respective recurrent parents. The process of developing the ILs is summarized in Supplemental Fig. 1.
The 326 ILs were screened for leaf blast resistance at DRR in three seasons (2010 to 2011), 50 ILs showed resistant reaction with mean score at 0 to 3 (Fig. 1 and Supplemental Table 4) and 276 were susceptible with mean score at 4 to 9), while both the recipient parents showed extreme susceptibility. Among the 50 resistant ILs, surprisingly majority of the lines (41) were under the PR114 background and only 9 lines were under the Pusa 44 background. Among all the wild species screened for blast, O. glaberrimacontributed major share (34%) of resistance, followed by O. longistaminata(32%), O. nivara, O. rufipogon, O. glumaepatula (each species for 10%) while the least fromO. barthii(4%). Among the 50 resistant ILs, top four ILs (designated as IL-1 to IL-4, Table 2) were selected (based on the phenotypic data ranging from 1.0 to 1.5) for screening at endemic regions of leaf and neck blast to check their field resistance across different parts of India (Supplemental Table 2). All the ILs showed high level of resistance, and their mean resistance score were in the range of 1.3 to 2.6. IL-1 showed the higher resistance (mean score of 1.3) than the others (Table 2). Furthermore, the screening at hilly areas for neck blast resistance indicated that all the lines were resistant to neck blast with the mean scores ranging from 1.6 to 2.3 (Table 2). These ILs screened for multiple disease resistance under AICRIP during 2012 and 2013 showed leaf and neck blast resistance with susceptibility index ranging from 3.3 to 3.8 and 2.0 to 3.0, respectively. The comparison of blast scores were also made with C101LAC (monogenic line for Pi1), C101A51 (monogenic line for Pi2) and HR12 (local susceptible check). The results showed that the ILs showed resistant reaction better than those of lines containing Pi1 and Pi2 genes, the resistant checks (Table 2). The agronomic trait evaluation of the top four ILs has revealed that the ILs performance was at par with that of the recurrent parents (Supplemental Table 5 and Supplemental Fig. 2).
To know the allelic status of gene in the ILs, genotyping was done using gene based markers associated with blast R genes. The ILs which showed marker alleles same to positive check were considered as positive for the gene (score as 1), similarly, the ILs which showed marker allele correlating to the susceptible checks were considered as negative for the gene (score as 0). A total of seven important blast resistance genes of India were genotyped in the present study, and none of these four ILs showed the presence of tested blast R gene alleles (Supplemental Table 6).
To identify the introgressed genomic regions from the wild Oryza species, four ILs were screened using 499 SSR markers (selected based on at least one marker per Mb) (Fig. 2). The results revealed the presence of 0.6% to 4.6% donor genome introgression (Supplemental Fig. 3). Most introgressions were observed at one region of particular chromosome except in IL-1 and IL-4. These two ILs showed the introgression at multiple chromosomes. The physical distances of introgressed regions among ILs varied from 2.3 to 9.0 Mb across different chromosomes (Table 3 and Supplemental Fig. 4). IL-2 had the minimum introgression (0.6%) of donor genome at chromosome 3 from 2.4 to 4.3 Mb region whereas IL-4 had the maximum introgression (4.6%) (Table 4). This information will facilitate in precise mapping and identification of gene(s) in these ILs.
The meta QTL analysis revealed the presence of meta QTLs in all four ILs where the donor genome introgression regions were found. In IL-3 and IL-4, we determined the presence of blast R genes. Through in silico, we observed the presence of nucleotide binding site-leucine rich repeat, nucleotide-binding adaptor shared by ARC domain containing protein, and serine/ threonine kinase genes in the target regions in all the ILs, which may be the candidate genes involved in resistance (Table 3).
Wild species are rich resources for a variety of important genes which have agronomical importance (Cheema et al, 2008; Nataraj et al, 2011) including disease resistance (Jena et al, 2006; Ramkumar et al, 2011). However, till today only two genes (Pi9 and Pi40) have been identified from wild species (Amante et al, 1992; Jeung et al, 2007). ILs have proven to be potential genetic resources for detecting rice blast resistance (Tsunematsu et al, 2000; Rahman et al, 2011) and also for exploring the agronomic traits (Fujita et al, 2009; Gu et al, 2012). Since we used different AA genome for the development of the ILs, it is expected to identify novel genetic resources for rice blast disease. Hence it is worth screening ILs for blast resistance, upon stringent screening at DRR with most virulent race in artificial screening conditions, 50 ILs were resistant and quite numbers of ILs were susceptible (Supplemental Table 1), since the advancement of generations of these lines were done earlier based on the agro-morphological and yield-related traits which resembles the recurrent parent. Majority of the resistant ILs were derived from O. glaberrima and O. longistaminata. Recently, Sié et al (2012) and He et al (2014) have demonstrated the utility of O. longistaminataand O. glaberrima for many valuable traits like rhizomatousness, disease resistance and drought tolerance that can be used to improve cultivated rice through transcriptomic, proteomic and metabolomic studies, and they exploited the several accessions of these wild species in breeding programs. Out of 50, only 4 lines were further advanced to check their field resistance at multiple locations, which included hotspots for leaf and neck blast. These hotspots, known to have diverse pathogen populations, generally offer field resistance, which is much better than testing with the single isolates of the pathogen. Further, these lines were also tested with artificial and natural infestations with predominant blast isolates of various test sites of AICRIP programme. The tested four ILs showed consistently resistant reactions in all the tested locations and seasons. Though there are many reports on the leaf blast resistance genes, less information is known towards neck blast, as only one gene Pb1 has been reported so far (Hayashi et al, 2010). In this context, the four ILs can be used as potential donors for neck blast resistance breeding programmes in India. The top four ILs were derived from O. glumaepatula (IL-1), O. glaberrima (IL-2), O. nivara (IL-3) and O. barthii(IL-4), respectively. Very interestingly, till date no major gene or major QTL has been identified from any of the four wild species from which the ILs have derived. The rest 46 resistant ILs were unexplored but the material is worth enough to look into further.
Apart from stringent phenotypic screening, the gene based markers were used to identify the known blast R alleles (Hayashi et al, 2006; Lmam et al, 2014). For this study, seven blast R genes which are being used more often in blast resistance breeding programs in India were selected (Prasad et al, 2011). Moreover, it is difficult and time-consuming process to conduct allelism tests for all the known blast R genes. Present study revealed that the selected four ILs does not have alleles of the seven tested blast R genes. Similar approach was followed earlier for identification of alleles or better performing alleles from the diverse germplasm through PCR based approach (Mahender et al, 2012; Ingole Kishor et al, 2014; Khan et al, 2014; Lmam et al, 2014).
Availability of large number of reliable SSR markers in rice and the availability of donor accessions, we can precisely locate the donor genome introgression of four ILs. Interestingly, the donor genome introgression was small in most of the ILs and these genomic regions were not reported earlier for presence of major blast resistance genes (at least in two ILs), indicating the presence of novel genes. Increasing the marker density along the chromosomes will enhance the precise identification of introgression regions. In the present study, 499 SSR markers (one marker per Mb) spread uniformly across the genome (except on chromosomes 8 and 9) were selected for this assay. Further increase in density of markers would certainly aid in the narrowing down the introgression regions. This observation was further strengthened by in silico analysis, which showed the presence of various disease related genes in the introgressed regions among the ILs (Table 3). In addition, we observed the presence of meta QTLs in the introgressed regions of all ILs (Ballini et al, 2008), indicating the importance of these genomic regions in the point of field resistance. This information will be useful in mapping of genes associated with resistance. All the ILs used in this study were stable lines, since they were in BC2-3F9 generation which can directly enter into the cultivar development programme. These lines can serve as good resources for the identification of novel genes for leaf and neck blast resistance. Two of these lines (IL-1 and IL-2) were registered (INGR15001 and INGR15002) by the National Bureau of Plant Genetic Resources (NBPGR), India, which is a responsible organization for repository and proper distribution of the ILs as open materials for further studies. Further, we identified that 46 ILs which are also resistant to blast and multiple disease, can be further evaluated for other biotic and abiotic stress. To conclude, we used the wild species in resistance breeding programmes by developing ILs, and we followed systemic approach to identify the ILs which showed the stable blast resistance. The scope of this study lies in the identification of four valuable and novel genetic resources which can be exploited further for blast resistance breeding programmes.
S. J. S. Rama Devi thanks the Council for Scientific and Industrial Research (CSIR), New Delhi, India and Department of Genetics, Osmania University, Hyderbad, Inida for providing Senior Research Fellowship.
The following materials are available in the online version of this article at http://www.sciencedirect.com/science/ journal/16726308; http://www.ricescience.org.
Supplemental Table 1. List of introgression lines used in this study.
Supplemental Table 2. List of leaf and neck blast hot spot areas in India.
Supplemental Table 3. Disease severity score and corresponding symptoms on plants.
Supplemental Table 4. Resistant introgression lines identified at Indian Institute of Rice Research.
Supplemental Table 5. Agronomic evaluation of the introgression lines.
Supplemental Table 6. Identification of rice blast resistance genes using gene specific markers.
Supplemental Fig. 1. Development and blast disease resistance evaluation of introgression lines.
Supplemental Fig. 2. Agronomic trait evaluation of the four introgression lines.
Supplemental Fig. 3. Representative gel pictures for identification of donor genome in introgression lines IL-1, IL-2 and IL-3.
Supplemental Fig. 4. Precise identification of donor genome introgression in selected introgression lines using genome wide simple sequence repeat markers.
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|
[36] |
|
[37] |
|
[38] |
|
[39] |
|
[40] |
|
[41] |
|
[42] |
|
[43] |
|
[44] |
|
[45] |
|
[46] |
|
[47] |
|
[48] |
|
[49] |
|
[50] |
|
[51] |
|
[52] |
|
[53] |
|
[54] |
|