Bangale Umakanth1†
Balija Vishalakshi1†
P. Sathish Kumar1
S. J. S. Rama Devi1
Vijay Pal Bhadana2
P. Senguttuvel3
Sudhir Kumar4
Susheel Kumar Sharma5
Pawan Kumar Sharma5
M. S. Prasad6
Maganti S. Madhav1*- 1Biotechnology Division, ICAR-Indian Institute of Rice Research, Hyderabad, India
- 2Plant Breeding, ICAR-Indian Institute of Rice Research, Hyderabad, India
- 3Hybrid Rice Division, ICAR-Indian Institute of Rice Research, Hyderabad, India
- 4Plant Breeding Section, ICAR Research Complex for NEH Region, Manipur Centre, Imphal, India
- 5Plant Pathology Section, ICAR Research Complex for NEH Region, Manipur Centre, Imphal, India
- 6Plant Pathology Division, ICAR-Indian Institute of Rice Research, Hyderabad, India
North-East (NE) India, the probable origin of rice has diverse genetic resources. Many rice landraces of NE India were not yet characterized for blast resistance. A set of 232 landraces of NE India, were screened for field resistance at two different hotspots of rice blast, viz., IIRR-UBN, Hyderabad and ICAR-NEH, Manipur in two consecutive seasons. The phenotypic evaluation as well as gene profiling for 12 major blast resistance genes (Pitp, Pi33, Pi54, Pib, Pi20, Pi38, Pita2, Pi1, Piz, Pi9, Pizt, and Pi40) with linked as well as gene-specific markers, identified 84 resistant landraces possessing different gene(s) either in singly or in combinations and also identified seven resistant landraces which do not have the tested genes, indicating the valuable genetic resources for blast resistance. To understand the molecular diversity existing in the population, distance and model based analysis were performed using 120 SSR markers. Results of both analyses are highly correlated by forming two distinct subgroups and the existence of high diversity (24.9% among the subgroups; 75.1% among individuals of each subgroup) was observed. To practically utilize the diversity in the breeding program, a robust core set having an efficiency index of 0.82 which consists of 33 landraces were identified through data of molecular, blast phenotyping, and important agro-morphological traits. The association of eight novel SSR markers for important agronomic traits which includes leaf and neck blast resistance was determined using genome-wide association analysis. The current study focuses on identifying novel resources having field resistance to blast as well as markers which can be explored in rice improvement programs. It also entails the development of a core set which can aid in representing the entire diversity for efficiently harnessing its properties to broaden the gene pool of rice.
Introduction
Rice is a major staple crop in the world and demand for rice is increasing every year (Ray et al., 2013). Biotic stresses continue to be the constraint in rice production and becoming severe in the climate change regime. One of the most identifiable major biotic stresses is the blast disease caused by Magnaporthe oryzae. It infects rice leaves, nodes, color, and panicles at different stages of crop growth thus, decreasing the overall yield. In India, neck blast is becoming severe in many agro-ecological zones and causing more threats to the rice production (Laha, 2017). Although 100 major blast resistance genes (R-genes) have been identified, mapped and their tightly linked DNA markers are available (Miah et al., 2013), only one major gene (Pb1) has been reported for neck blast (Hayashi et al., 2010). In addition to this, rapid evolution and high variability of Magnaporthe with the rise in temperatures relates to its high pathogenicity. Hence identification of durable resistance sources and resistance genes is a continuous process in rice improvement programs.
Landraces are genetically dynamic and display equilibrium with both the environment and pathogens (Harlan, 1975). Many potential landraces of rice are being replaced by high yielding varieties to meet the food requirements. Despite being less productive they are known to have a high genetic variance for several biotic stresses (Hanamaratti et al., 2008), so they can be explored for rice improvement. North-East (NE) India, known to have highly diverse rice germplasm and the landraces belonging to this region are expected to have high genetic variability for many biotic stresses including the rice blast. Many rice landraces of NE India were not yet characterized for blast resistance.
Genetic diversity among the populations and their genetic relationships aids in conservation and parental selection in the improvement programs. Identification of populations with a high level of genetic variation will become a valuable resource for broadening the genetic base since it enables the identification of superior alleles for several traits including blast resistance. Association mapping is the best approach to identify the chromosomes segments harboring genes/QTLs (quantitative trait loci) controlling key traits and the genotype–phenotype associations in germplasm. For the successful application of association mapping, understanding of population structure is essential to avoid type I and type II errors between molecular markers and traits of interest (Zhang P. et al., 2011). Determination of population structure also helps in categorizing populations into subsections which are genetically related. The rapid advances in marker technology along with morphological traits can precisely estimate the genetic diversity among the populations (Kumbhar et al., 2015). For estimation of genetic diversity, population structure, linkage disequilibrium, and mapping various traits, molecular markers have been successfully used by many researchers (Ried et al., 2011; Nachimuthu et al., 2015; Anandan et al., 2016). SSRs aid in accurately estimating the genetic diversity among germplasm and are reported to be more efficient than single-nucleotide polymorphic markers (SNPs; Das et al., 2013; Singh et al., 2013; Nachimuthu et al., 2015). Development of a core set representing the total variation of the population with novel and superior traits is essential for successful utilization of the population. Core and mini core sets of rice were earlier made using morphological traits and molecular markers (Zhang H. et al., 2011; Liu et al., 2015). However, no core set is available exclusively for the NE landraces. Evaluation of genetic diversity among NE Indian germplasm was reported earlier (Das et al., 2013; Roy et al., 2014, 2015) but very few studies exist on the understanding the diversity in relation to the blast resistance (Ghaley et al., 2012; Mahender et al., 2012).
In the current study, we have systematically performed phenotyping for blast resistance in a span of two years. Phenotyping along with gene profiling for major blast genes identified novel genetic resources which might contain the novel resistance genes. For clear understanding, the genetic diversity of NE landraces were analyzed with distance as well as model based diversity approaches. To apply the existing diversity in breeding programs, a core set was derived. Using 232 landraces, the marker–trait association was identified for blast resistance as well as yield-related traits. The information generated in the present study has a lot of value in the current rice improvement programs.
Materials and Methods
A total of 232 landraces (Supplementary Table S1) were collected from the NE region of India (Figure 1) for the current study. Along with the landraces, three varieties of Basmati, five varieties of Indica, seven varieties of Temperate Japonica, and five varieties of Tropical Japonica were also selected for the comparison of landrace alleles. The line harboring particular blast resistance gene was considered as a positive control, and the susceptible cultivars namely BPT, Swarna, and CO39 which do not have any blast resistant gene were considered as negative controls for gene profiling.
FIGURE 1. Geographical distribution of diverse rice landraces (232) across North Eastern states of India.
Evaluation of Blast Disease at Two Hotspot Locations
Approximately 40–50 seeds of each landrace were sown in uniform blast nursery (UBN) for phenotyping of leaf blast disease at ICAR-IIRR, Hyderabad, India. The entire nursery block, as well as each entry, was surrounded by HR12 (susceptible cultivar) on four sides for spreading the disease evenly. A mixture of M. oryzae isolates collected from different blast hotspot regions of Andhra Pradesh and Telangana (Surapu et al., 2017) were used to screen landraces for the leaf blast disease. The plants were spread with inoculum at the two-leaf stage and the score was taken once the susceptible check gets totally infected with blast disease. The score was taken based on 0–9 scale (Standard evaluation system, SES, IRRI, 2002; Prasad et al., 2011). Scores 0–3 were considered resistant (R), 4–5 as moderately resistant (MR), and 6–9 as susceptible (S). Field screening of blast disease was done at ICAR Research Complex for NEH Region, Manipur under natural disease conditions and the disease pressure was increased by inoculating with the mixture of Magnaporthe isolates that were collected from different blast hotspots of NE India. For augmenting the disease pressure for Neck blast, sowings were delayed than usual season and also each plant was injected with a spore suspension of a mixture of isolates collected from NE India using a syringe (Aglawe et al., 2017). Disease score for each plant was recorded for at least three times and the highest score out of these was considered. Neck blast was evaluated based on the percentage of infections on the neck of the panicle at physiological maturity of the plant (Madhav et al., 2013). The highest disease score observed in both the seasons was considered. The screening was done in the same way for two seasons (Kharif 2012 and 2013) in two replications.
Agro-morphological Evaluation of Landraces
The landraces were grown in a block prepared by uniform distribution of soil enriched with nutrients required for the plant growth. A spacing of 20 cm between each row and 15 cm between plants was adopted at the time of planting. Twenty plants of each landrace were sown in rows and the data was recorded for five selected plants. Data was taken on various parameters like plant height, presence and absence of awns, seed color, tiller number, and spikelet number. At the time of harvesting, biomass and yield per plant were also recorded.
DNA Isolation
DNA was isolated from the leaf tissue by following Murray’s (Murray and Thompson, 1980) protocol. The quantity of DNA was estimated in nanodrop (Thermo Fisher Scientific, United States) as well as in agarose gel.
Gene Profiling
PCR analysis was done for landraces together with the positive controls and susceptible checks to identify the presence of 12 major known blast resistance genes namely Pitp, Pi33, Pi54, Pib, Pi20, Pi38, Pita2, Pi1, Piz, Pi9, Pizt, and Pi40 using molecular markers as described in Devi et al. (2015). Gene-based (SNP and STS) markers were used for six genes, viz., Pi54, Pi40, Pita2, Pi9, Piz, and Pizt and the remaining genes were profiled with tightly linked SSR markers. The markers used for the gene profiling and positive controls for respective genes are listed in Supplementary Table S2. Genetic diversity analysis was carried out with 120 SSR markers which are uniformly spread across the 12 chromosomes (Supplementary Table S3).
Diversity Analysis
For diversity analysis, all the alleles derived from each SSR marker were scored as “1” for the presence of amplicons and “0” for its absence in the respective allele position. The size of the amplicons was also determined based on their migration relative to standard molecular weight marker (50 bp DNA ladder from Fermentas). Allele number, allele frequency, heterozygosity, and polymorphic information content (PIC) values for all the markers were calculated using Power marker V3.25 software (Liu and Muse, 2005).
Clustering of Genotypes
For clustering of genotypes, a pairwise distance matrix was computed using DARwin V6 software. An unweighted neighbor-joining tree was constructed using dissimilarity index method. The genetic distance was calculated by Jaccard’s coefficient with bootstrap analysis of 1000 iterations using the SSR marker data (Perrier and Jacquemoud-Collet, 2006).
Population Structure Determination
The structure of the population was determined by a Bayesian-based approach using Structure V2.3.4 software (Pritchard et al., 2000). The number of subpopulations was classified based on the true K value. For accurate determination of subpopulations, admixture model was followed. In this model, the burnin period length was set to as high as 150,000 with Markov Chain Monte Carlo (MCMC) 150,000 replications. Each K value was run for 10 times with K value varying from 1 to 25. To obtain clear peak at the ΔK (ad hoc quantity) to K value, the output generated from structure software was loaded in structure harvester (http://taylor0.biology.ucla.edu/structureHarvester/).
Molecular Variance
The presence of molecular variance among the population and individuals estimated by structure was analyzed through analysis of molecular variance (AMOVA). The estimation of AMOVA based on Nei’s distance matrix (Nei, 1973) along with the pairwise FST and the principal component analysis (PCoA) to characterize the subgroups of the germplasm set was computed using GenAlEx 6.5 software (Peakall and Smouse, 2012).
Establishment of Core Set
A mini-core collection was developed using Powercore V1.0 software (Kim et al., 2007). For both core and mini-core collections, Nei genetic diversity index and Shannon–Weaver diversity index was calculated. The significant difference percentage between core set and the entire collection was analyzed for evaluating mean difference percentage (MD%) and variance difference percentage (VD%) of traits. To evaluate the properties of a core set against the entire collection, the coincidence rate percentage (CR%) and variable range percentage (VR%) were evaluated.
Association Analysis
The association tests between the marker and the agronomic traits in the population were run based on two models, general linear model (GLM) and mixed linear model (MLM) using TASSEL V3.0 (Bradbury et al., 2007). GLM was run by considering the population structure of accessions as described by Kang et al. (2008). The MLM was performed using the Kinship file generated from genotypic data and Q matrix generated by running population structure at optimized K value in order to minimize the false positive associations (Yu et al., 2006; Gupta et al., 2014). The significant marker–trait association was determined by P ≤ 0.05 and the magnitude of the QTL effects by marker R2.
Results
Leaf and Neck Blast Phenotyping
The stringent blast screening of 232 landraces at two locations in the year 2012 revealed 166 landraces showing resistance to leaf blast and 148 to neck blast. Whereas in 2013, 172 landraces were resistant to leaf blast and 154 to neck blast. On average, 153 landraces showed resistance to leaf blast and 135 landraces resistance to neck blast and interestingly, 91 showed resistance to both leaf and neck blast. The mean scores of resistant landraces ranged between 0 and 3 while the susceptible check HR12 gave the score of 9.
Gene Profiling
The gene profiling of landraces revealed that Pi9 was found in the maximum number of landraces, i.e., 102 (43.9%), while the Pi20 gene is present in the least, i.e., 11 landraces (4.74%). One unique landrace Apaghi june found to be a rich source for blast genes as it has seven genes, viz., Pi9, Piz, Pi33, Pib, Pi38, Pi40, and Pi54. Another landrace, i.e., Chalhtssia has six genes, viz., Pi1, Pi9, Pi54, Pib, Pitp, and Piz. Seven landraces have the combinations of five genes whereas 14 landraces have four gene combinations. Similarly, the landraces containing three gene combinations, two gene combinations and single genes were also identified; these details are given in Supplementary Table S1 (Figure 2). Among the 91 blast resistant landraces, seven landraces, viz., Meghalaya lakang, Chingphourel, Manuikhamei, Kemenya kepeyu, Wainem, Thekrulha, and Koyajang do not have any of the tested major blast resistant genes (Figure 3).
FIGURE 2. (A) Distribution of 12 major blast resistance genes in a set of 232 rice landraces collected from North-East region of India (∗gene having a gene-based marker). (B) Landrace harboring different number of blast resistance gene.
FIGURE 3. Phenotypic reaction of unique seven blast resistant landraces at IIRR’s uniform blast nursery with HR12 as a control for susceptibility.
Genetic Diversity
Diversity analysis with the allelic data generated from 120 SSR markers detected a total of 360 alleles, as the alleles per loci ranged from 2 to 8 with an average of 3 alleles per locus. Among the total alleles, 15% rare alleles (having allele frequency <5%) were found. The SSR marker RM22837 (Chr-8) recorded the highest number of alleles (8), which indicates its high polymorphic nature, while many markers showed two alleles. The major allele frequency ranged from 0.31 to 0.98 with an average of 0.59 per locus, and the average observed heterozygosity was 0.09. Where as expected heterozygosity ranged from 0.05 to 0.80 with an average of 0.51 per locus. The relative informativeness of each marker has been identified by PIC value, highest PIC value (0.77) was found for RM22837 while least (0.05) for RM21052 with an average PIC value of 0.44 for markers (Supplementary Table S3).
The neighbor-joining tree showed the presence of two distinct major clusters, of which 59 landraces grouped in cluster-1 and 173 landraces in cluster-2. The cluster-1 is again divided into four subclusters—1A, 1B, 1C, and 1D, whereas cluster-2 divided into five subclusters, i.e., 2A, 2B, 2C, 2D, and 2E. The landraces present in the subcluster 1B shares alleles with Japonica and Tropical Japonica varieties while landraces present in 1C shares with Japonica and Basmati accessions. The landraces present in the subcluster 2D shares with Indica, Japonica, and Tropical Japonica varieties while landraces of 2E shares with Japonica and Basmati accessions. The sub-clusters 1A, 1D, 2A, 2B, and 2C formed exclusively with landraces (Figure 4).
FIGURE 4. Unweighted neighbor-joining tree of landraces grouped into two major clusters representing the genetic relationship among the landraces.
The population structure analysis of landraces revealed the log likelihood value (ΔK) maximized to the highest value of 2260.4 at K = 2 (Figure 5A), showing a clear peak indicating the classification of entire population into two distinct subgroups SG1 (subgroup1) and SG2 (subgroup2) (Supplementary Table S1). The landraces sharing the membership fractions with a probability of greater than 80% were allotted to corresponding subgroups (Figure 5B). Among the 232 landraces, 33(14.2%) landraces were categorized under SG1 whereas the majority of the landraces that is 158 (68.1%) were grouped in SG2, and the remaining 41(17.7%) landraces were grouped as admixture. Subgroup2 which has a high number of landraces (158) was further subjected to structure analysis. This analysis revealed the classification into four sub-subgroups (SSG) (Supplementary Figure S1) as the ΔK value showed a clear peak at K = 4. The SSG1 consists of 24 landraces, 13 in SSG2, 10 in SSG3, 12 SSG4, and 99 landraces in admixture.
FIGURE 5. Representation of population structure dividing the landraces in two subgroups based on K value (A,B). Landraces with >80% probability are assigned to corresponding subgroups and exceptions are shown as an admixture. (A) Population structure inferring the landraces divided into two subgroups based on K value. (B) Population structure inferring the landraces divided into two subgroups based on K value. Landraces with probability >80% were assigned to corresponding subgroups and remaining as an admixture.
The distribution of molecular variance exposed significant genetic differentiation between and within populations among the landraces. Molecular variation of 75.1% was observed among the individuals within the subgroup (population) and 24.9% variation was observed among the subgroups. The pairwise FST value 0.278 indicated the significant difference between two subgroups (Table 1).
TABLE 1. AMOVA between structured landraces and pairwise comparison using FST values (GenAlEX).
The PCoA was performed for characterizing the subgroups of landraces. The first two-dimensional scatter plots described the PCoA axes accounting for 12.48 and 4.04% of the genetic variation among subpopulations formed from 252 accessions (Supplementary Figure S2).
Powercore Analysis
To make the core set representing the entire diversity of landraces, 11 morphological characters and genotypic data of 120 SSR markers was used for the power core analysis using non-heuristic approach. This analysis identified 33 landraces (14.2%) that represent the core set (Figure 6 and Supplementary Table S4A). The core set has 4.42% MD%, and 100% CR%. The core collection with VD% 29.17% and VR% (variable rate of coefficient of variation) 118.18% was observed. The efficiency index of the core was 0.82 and the PIC was 0.99 estimated by phenotypic characters and SSR markers (Supplementary Table S4B). Among phenotypic traits, biomass contributed to the highest Nei’s diversity index of 0.74 while the presence of awn contributed the lowest (0.21). The SSR markers showed an average Nei’s diversity of 0.6 in the core set.
FIGURE 6. Depiction of (A) core set land races and (B) distribution of its subgroups.
Association Analysis
To identify the marker associated with the blast resistance as well as with other traits, we used 232 landraces for association study. Association analysis of landraces based on GLM and MLM using Tassel software revealed 20 marker–trait associations with a significant P-value of less than 0.05. One marker (RM480) present on chromosome-5 showed association with leaf blast with a significant r2 value of 0.21 and at this location, blast gene Pi26 was reported. Similarly, RM5401 located on chromosome 5 having an r2 value of 0.21 is associated with the neck blast resistance where no single QTL for blast resistance was reported. Two markers (RM21936 and RM3787) located on chromosome 7 and chromosome 9, respectively were associated with both leaf and neck blast resistance. About 15 markers were associated with yield-related traits namely number of tillers, a number of spikelets, plant height, yield per plant, and color. The marker–trait associations and their significance were listed in Table 2. The relative kinship estimates based on the 120 SSR data resulted in 83.3% of the pairwise kinship estimates ranging from 0.25 to 1.0 and the remaining estimates ranged from 1 to 2, with a low number of pairs filling in the higher estimate categories (Supplementary Figure S3). The kinship analysis revealed that the many of the accessions had a weak relationship with the other accessions in the landrace population. The linkage disequilibrium analysis among the locus pairs of 120 SSR markers revealed the D′ values ranging from 0 to 0.52 with an average of 0.33. About 26% of the marker pairs were having a significant LD (P < 0.01).
TABLE 2. Association analysis between agronomic traits and SSR markers.
Discussion
Successful rice breeding programs require the presence of genetic diversity in the existing germplasm. The current study focuses on ascertaining this phenomenon in NE Indian landraces. We used a completely unique and unexplored germplasm for identifying the best donors and markers associated with blast resistance. We demarcate the diversity into a manageable core set.
Systematic screening for blast disease at two different locations for two consecutive seasons led to the identification of 91 landraces resistant to both leaf and neck blast. Blast screening at different locations increases the consistency of phenotyping as different virulent blast isolates have diverse pathogenic behavior (Madan et al., 2012). Previously a similar strategy was used for identification of best introgression lines for blast resistance (Devi et al., 2015). A mixture of blast isolates from geographically distinct locations (Hyderabad and Manipur) led to the identification of stable QTLs in one of the NE landrace (Aglawe et al., 2017). Neck blast is becoming more problematic than the leaf blast at several locations in India (Aglawe et al., 2017; Laha, 2017) so, it has become imperative in breeding programs to choose the accession having resistance for two phases of blast resistance (Abhijeet et al., 2013). So far only one neck blast resistance gene (Pb1) was identified (Hayashi et al., 2010) and many of leaf blast resistant accessions may not hold the resistance for the neck blast, hence screening for neck blast was given a priority. Moreover, the NE hilly areas have compatible climatic conditions for a natural infestation of neck blast.
In addition to the phenotyping, gene profiling using SNP, SSR, STS, and indel markers was also performed so as to identify the presence of valuable genes to the blast resistance (Cho et al., 2008) (Supplementary Table S2). These results indicated that Pi9 was major gene followed by Piz and Pizt existing in landraces. The Pi9 is predicted to be persistent in landraces from a long time as it was originated from the wild species (Oryza minuta) of rice. Anupam et al. (2017) identified the frequency of seven blast genes ranging from 0 to 80% in 74 local landraces collected from Tripura by profiling with gene-based STS and SNP markers. Prevalence of Piz among 32 NE Indian landraces by gene profiling with STS and SNP markers was reported (Imam et al., 2014). Shikari et al. (2013) also evaluated the 54 landraces collected from different parts of India for the presence of Pita and Pita2 genes with two gene-based and three linked markers and found 11 landraces containing Pita and nine landraces containing Pita2 genes. Gene-specific markers of eight blast resistance genes (Pib, Piz, Pizt, Pi9, Pi40, Pi5, Pia, and Pita) were also explored by Mahender et al. (2012) for characterizing the various accessions of Manipur (another state of NE India) and found landraces containing two to seven blast resistance genes, which again indicates the richness of blast diversity in landraces of NE India. In the current study, the landraces used were entirely different from the previous work (Mahender et al., 2012; Imam et al., 2014; Roy et al., 2015), we also checked their presence in the national gene bank (National Bureau of Plant Genetic Resources, NBPGR, New Delhi) and found 189 landraces that does not exist in it. Around 85% of landraces (197) harbors blast genes ranging from 1 to 7, indicating the rich diversity for blast resistance and these landraces serves as good genetic resources for blast resistant genes. Interestingly, seven resistant landraces did not show the presence of any of the tested genes, may contain the novel genes (Supplementary Table S1 and Figures 2, 3). Identification of new genes having durability and its transfer into a cultivable variety may increase its resistance toward the disease since most of the identified genes are losing their effectiveness to blast disease (Das et al., 2012).
Genetic diversity analysis was done for determining the diverse nature of the landraces. To precisely differentiate closely related individuals, 120 SSR markers which are well distributed in the genome were used in the current study. The extent of polymorphism was detected by calculating PIC values and for evaluation of genetic diversity across the chromosomes; an allelic number was also calculated. Average PIC value of markers was 0.44 with an average of three alleles per locus (Supplementary Table S3). An average PIC value of 0.47 was reported in 74 Tripura rice landraces by using 30 SSR markers (Anupam et al., 2017). Roy et al. (2016) found a PIC value of 0.62 with the NE Himalayan landraces. Yang et al. (1994) reported average PIC value of 0.58 in 238 landraces collected from India, China, Japan, and IRRI. Giarrocco et al. (2007) also found average PIC value of 0.69 in 68 accessions of Argentina. A significant amount of rare alleles were found (15%) which might contribute to the overall genetic diversity of the population. Nachimuthu et al. (2015) reported the presence of 5% rare alleles in 192 rice germplasm upon screening with 61 SSR loci. Anupam et al. (2017) reported the existence of 11 rare alleles in 74 germplasm lines upon analysis with 30 SSR loci.
In the present study, a gene diversity of 0.51 was observed among the total population (232 landraces), which is comparable to most of the diversity panel consisting of global accessions (gene diversity of 0.5–0.7) (Ni et al., 2002; Garris et al., 2005; Ali et al., 2011; Roy et al., 2015, 2016). Gene diversity of NE landraces is higher than overall gene diversity of rice core collection (0.544) from China, Philippines (Zhang H. et al., 2011) and US accession panel (average gene diversity of 0.43; Agrama and Eizenga, 2008). These results clearly indicating that the NE landraces of India are showing the high genetic diversity in rice that exists in India (Supplementary Table S3).
The distance based clustering analysis grouped the landraces into two major clusters (Figure 4). Choudhury et al. (2013) used the distance based clustering method to group the 100 cultivable rice varieties released in different states into their respective decade periods. Similarly, rice genotypes collected from Tamil Nadu were classified into indica and japonica types based on the distance based clustering (Vanniarajan et al., 2012). In addition to the distance based clustering of the population, the genetic architecture of diverse material was also done by the model based structure analysis, as it offers better architecture of population with molecular markers (Powell et al., 1996; Horst and Wenzel, 2007; Varshney et al., 2007). Several researchers have analyzed the population structure by implementing model-based approach using Structure for numerous rice accessions (Garris et al., 2005; Agrama et al., 2007; Zhang et al., 2007, Zhang P. et al., 2011; Agrama and Eizenga, 2008; Jin et al., 2010; Ali et al., 2011; Chakhonkaen et al., 2012; Courtois et al., 2012). The present study classifies the landraces into two subgroups SG1 and SG2 using the threshold value >80%. Similar results were observed in grouping 3024 Chinese rice landraces (Zhang et al., 2009), 64 landraces of NE Himalayan region of India (Roy et al., 2016), and 192 rice accessions of varied origin (Nachimuthu et al., 2015). But, several studies were reported on the grouping of rice accessions which varied from two to eight groups (Garris et al., 2005; Agrama et al., 2007; Zhang et al., 2007, 2009; Zhang P. et al., 2011; Ali et al., 2011; Chakhonkaen et al., 2012). Roy et al. (2015) divided 107 NE landraces into three different subgroups by structure analysis using 40 SSRs. Das et al. (2013) used 23 SSR markers and distinguished 91 rice lines (which include 83 NE landraces and eight checks) into four subpopulations. Thirty-seven chakhao landraces of NE were also divided into six subgroups based on genetic structure analysis (Roy et al., 2014). The current study grouped 232 accessions based on the maximum membership probability of >80%, into 33 accessions under SG1, 158 accessions under SG2 and 41 accessions as an admixture (Figure 5). Further population structure analysis in SG2 classified the landraces into four subgroups indicating that independent structure analysis defines the precise genetic structure of the accessions. The results of the unweighted neighbor-joining tree constructed on the basis of UPGMA and PCoA harmonized with the structure analysis of the landraces. The PCoA cited the resistant and susceptible genotypes into four different quadrants at different positions.
AMOVA is useful for studying molecular variance within the species. It is used to study the patterns and degree of relatedness revealed by multidimensional scaling and the clustering dendrogram. It also summarizes the population structure, while remaining flexible enough to accommodate different types of assumptions about the evolution of the genetic system. AMOVA of the population showed greater variability among the individuals within the populations. The variation observed among the individuals within populations was higher (75.1%) than the variation present between the populations (24.9%). The FST value of 0.27 indicated that the higher divergence was observed between subgroups. Nachimuthu et al. (2015) upon analyzing 192 rice accessions (landraces and varieties collected from different parts of the world) and grouped them in two subgroups and reported 14% variations among the groups and 86% within the group. Similarly, Zhang et al. (2009) also grouped the 3024 rice accessions, into two subgroups and found the variation among the groups which was little higher (36.65%) and variation among individuals within populations was 40.80%.
The development of core set is very much essential for better management of germplasm lines in crop improvement strategies. It is helpful to precisely characterize, explore and conserve resources, monitor the genetic drift during preservation and identify gaps in genetic diversity (Frankel and Brown, 1984). Choudhury et al. (2014) developed a core set of 701 accessions (10%) among 6984 germplasm lines collected from NE India using molecular markers and agronomical traits. Zhang H. et al. (2011) established a core set of 189 accessions from 4310 Chinese accessions which include landraces. In the present study, a core set of 14% (33 landraces) was developed which represents the existing entire genetic diversity of the total landrace collection (Supplementary Table S4A and Figure 6). The core set developed indicates the homogeneous distribution ranges of the phenotypic traits as MD% is less than 20% and CR% is larger than 80% (Kim et al., 2007). The VD% and VR% present within the range represents the diversity as the original collection (Hu et al., 2000). This shows that the core collection represents the entire collection. At present, there is no universally accepted method for the development of core set as various factors such as sampling percentage, data type, number of traits observed, genetic diversity of germplasm, grouping method, and sampling method will play in demarcating the core group (Upadhyaya et al., 2006; Diez et al., 2012; Rao et al., 2012).
Association mapping is a novel approach which helps in identifying marker–trait relation based on linkage disequilibrium (Zondervan and Cardon, 2004). Based on GLM, we identified 20 marker associations for blast disease and other agronomic traits. Of which, eight markers loci are novel since, in their location, QTLs for those traits were not reported. These newly identified trait-associated SSR markers can be used for the discovery of novel genes/QTLs. While with MLM analysis, only 17 marker–trait associations were found and all these associations were reported in the GLM. Associations detected in GLM were not found in MLM, since in later analysis effect of kinship may play a role. The markers RM5499 and RM190 which are previously reported for grain yield also showed association in the present study and in the same way two markers (RM286 and RM6051) for plant height. Few markers showed association with more than one trait which is of great interest in breeding programs. Two markers RM21936 and RM3787 were newly identified to be associated significantly (P < 0.05) with leaf and neck blast disease (Table 2). These markers may be helpful for mapping genes for neck blast with durable resistance which does not exist in the currently available germplasm. Exploring novel gene pool for neck blast resistance and identification of novel genes which can offer resistance to leaf and neck blast resistance has become quite essential. Anandan et al. (2016) reported 16 marker associations in 629 rice genotypes using 39 SSR markers for identification of genes associated with early seedling vigor using GLM. In an association analysis of yield and grain quality traits in Embrapa Rice Core Collection (ERiCC) of 242 accessions (Inbred lines, cultivars, upland, and lowland accessions) from Brazil with 86 SSR markers and field data, eight markers were identified which were associated with four different traits (Borba et al., 2010). Twenty marker associations were identified for 29-grain metabolites in rice using 218 SSR markers in a set of 48 rice germplasm lines from the Chinese core collection (Lou et al., 2011). Twenty-six SSR markers associated with blast resistance (P < 0.01) in a set of 276 indica landraces from China and few from different parts of the world were reported (Wu et al., 2016). The identified markers associated with specific traits can be validated for their effectiveness in various genetic backgrounds and those can be helpful to the breeders in the pyramiding of QTLs from different sources through marker-assisted selection. Hence association mapping became an unleashing tool for identification of markers linked to different traits.
Conclusion
Host and pathogen arms race is a continuous phenomenon in evolution. The premise of collecting, analyzing and characterizing unexplored NE landraces is for identification of donors containing multiple blast resistance genes which can be readily used in the breeding programs. The unique seven landraces (Meghalaya lakang, Chingphourel, Manuikhamei, Kemenya kepeyu, Wainem, Thekrulha, and Koyajang) having resistance to leaf and neck blast can be explored for the discovery of novel genes for broadening the gene pool to combat the pathogen competition. Distance and model based diversity analysis are highly correlated and genetic architecture of diverse material was effectively determined. The rich diversity existing in the landraces were divided into manageable subgroups and identified the molecular variance present among the populations and individuals which can be harnessed in the selection of distant individuals. To utilize the diverse landrace germplasm, an efficient core set which represents the entire diversity was made. Further, association analysis using a core set identified novel eight markers associated with the blast and other agronomically important traits, these markers are the candidates for the validation in the larger populations. The information emanating this work will help in conservation and utilization of valuable rice germplasm of NE India.
Author Contributions
MSM designed all experiments, developed the structure and arguments for the paper, made critical revisions, and approved the final version. BU and BV carried out the major work and prepared manuscript and contributed to the refinement of the manuscript. PSK helped in analyzing the data and SRD helped in data collection of agro-morphological traits. PS reviewed and approved the final manuscript. VPB, SK, SKS, and PKS helped in the collection and phenotypic screening of rice landraces for neck blast resistance at the North-East region. MSP helped in the screening of leaf blast resistance at IIRR. All authors have read and approved the final manuscript.
Funding
The present research work manuscript was supported by Department of Biotechnology, Government of India with the sanction number BT/188/NE/TBP/2011 dated 26/4/2012.
Conflict of Interest Statement
The authors declare that the research was conducted by the support of DBT, Government of India.
Acknowledgments
The authors acknowledge the Department of Biotechnology (DBT), New Delhi, India for providing funds for carrying out the research work. BU and BV thank Dr. P. Sudhakar and Dr. A. Krishna Satya, Department of Biotechnology, Acharya Nagarjuna University (ANU), Guntur, India for the Ph.D. registration.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.01500/full#supplementary-material
References
Abhijeet, G., Willocquet, L., Savary, S., and Kumar, J. (2013). Variability in aggressiveness of rice blast (Magnaporthe oryzae) isolates originating from rice leaves and necks: a case of pathogen specialization? PLoS ONE 8:e66180. doi: 10.1371/journal.pone.0066180
Aglawe, S. B., Bangale, U., Ramadevi, S. J. S., Balija, V., Pal, B. V., Kumar, S. S., et al. (2017). Identification of novel QTLs conferring field resistance for rice leaf and neck blast from a unique landrace of India. Gene Rep. 7, 35–42. doi: 10.1016/j.genrep.2017.01.007
Agrama, H., and Eizenga, G. (2008). Molecular diversity and genome-wide linkage disequilibrium patterns in a worldwide collection of Oryza sativa and its wild relatives. Euphytica 160, 339–355. doi: 10.1007/s10681-007-9535-y
Agrama, H., Eizenga, G., and Yan, W. (2007). Association mapping of yield and its components in rice cultivars. Mol. Breed. 19, 341–356. doi: 10.1007/s11032-006-9066-6
Ali, M., McClung, A. M., Jia, M. H., Kimball, J. A., McCouch, S. R., and Georgia, C. E. (2011). A rice diversity panel evaluated for genetic and agro-morphological diversity between subpopulations and its geographic distribution. Crop. Sci. 51, 2021–2035. doi: 10.2135/cropsci2010.11.0641
Anandan, A., Anumalla, M., Pradhan, S. K., and Ali, J. (2016). Population structure, diversity and trait association analysis in rice (Oryza sativa L.) germplasm for early seedling vigor (ESV) using trait linked SSR markers. PLoS ONE 11:e0152406. doi: 10.1371/journal.pone.0152406
Anupam, A., Imam, J., Quatadah, S. M., Siddaiah, A., Das, S. P., Variar, M., et al. (2017). Genetic diversity analysis of rice germplasm in Tripura State of Northeast India using drought and blast linked markers. Rice Sci. 24, 10–20. doi: 10.1016/j.rsci.2016.07.003
Barman, S. R., Gowda, M., Venu, R. C., and Chatto B. B. (2004). Identification of a major blast resistance gene in the rice cultivar ‘Tetep’. Plant Breed. 123, 300–302. doi: 10.1111/j.1439-0523.2004.00982.x
Berruyer, R., Adreit, H., Milazzo, J., Gaillard, S., Berger, A., Dioh, W., et al. (2003). Identification and fine mapping of Pi33, the rice resistance gene corresponding to the Magnaporthe grisea avirulence gene ACE1. Theor. Appl. Genet. 107, 1139–1147. doi: 10.1007/s00122-003-1349-2
Borba, T. C. O., Brondani, R. P. V., Breseghello, F., Coelho, A. S. G., Mendonca, J. A., Rangel, P. H. N., et al. (2010). Association mapping for yield and grain quality traits in rice (Oryza sativa L.). Genet. Mol. Biol. 33, 515–524. doi: 10.1590/S1415-47572010005000065
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Bucker, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Chakhonkaen, S., Pitnjam, K., Saisuk, W., Ukoskit, K., and Muangprom, A. (2012). Genetic structure of Thai rice and rice accessions obtained from the international rice research institute. Rice 5:19. doi: 10.1186/1939-8433-5-19
Cho, Y. C., Jeung, J. U., Park, H. J., Yang, C. I., Choi, Y. H., Choi, I. B., et al. (2008). Haplotype diversity and durability of resistance genes to blast in Korean Japonica rice varieties. J. Crop Sci. Biotech. 11, 205–214.
Choudhury, B., Khan, M. L., and Dayanandan, S. (2013). Genetic structure and diversity of Indigenous Rice varieties (Oryza sativa) in Eastern Himalayan region of Northeast India. Springer Plus 2, 228–237. doi: 10.1186/2193-1801-2-228
Choudhury, D. R., Singh, N., Singh, A. K., Kumar, S., Srinivasan, K., Tyagi, R. K., et al. (2014). Analysis of genetic diversity and population structure of rice germplasm from north-eastern region of India and development of a core germplasm set. PLoS ONE 9:e113094. doi: 10.1371/journal.pone.0113094
Courtois, B., Frouin, J., Greco, R., Bruschi, G., Droc, G., Hamelin, C., et al. (2012). Genetic diversity and population structure in a European collection of rice. Crop Sci. 52, 1663–1675. doi: 10.2135/cropsci2011.11.0588
Das, A., Soubam, D., Singh, P. K., Thakur, S., Singh, N. K., and Sharma, T. R. (2012). A novel blast resistance gene, Pi54rh cloned from wild species of rice, Oryza rhizomatis confers broad spectrum resistance to Magnaporthe oryzae. Funct. Integr. Genomics 12, 215–228. doi: 10.1007/s10142-012-0284-1
Das, B., Sengupta, S., Parida, S. K., Roy, B., Ghosh, M., Prasad, M., et al. (2013). Genetic diversity and population structure of rice landraces from Eastern and North Eastern States of India. BMC Genet. 14:71. doi: 10.1186/1471-2156-14-71
Devi, S. J. S. R., Kuldeep, S., Umakanth, B., Vishalakshi, B., Renuka, P., Sudhakar, K. V., et al. (2015). Development and identification of novel rice blast resistant sources and their characterization using molecular markers. Rice Sci. 22, k300–308. doi: 10.1016/j.rsci.2015.11.002
Diez, C. M., Imperato, A., Rallo, L., Barranco, D., and Trujillo, I. (2012). Worldwide core collection of olive cultivars based on simple sequence repeat and morphological markers. Crop Sci. 52, 211–221. doi: 10.2135/cropsci2011.02.0110
Fjellstrom, R., McClung, A. M., Conaway-Bormans, C. A., Anna, M. M., Marchetti, M. A., Shank, A. R., et al. (2004). Development of DNA markers suitable for marker assisted selection of three Pi genes conferring resistance to multiple Pyricularia grisea pathotypes. Crop Sci. 44, 1790–1798. doi: 10.2135/cropsci2004.1790
Fjellstrom, R., McClung, A. M., and Shank, A. R. (2006). SSR markers closely linked to the Pi-z locus are useful for selection of blast resistance in a broad array of rice germplasm. Mol. Breed. 17, 149–157. doi: 10.1007/s11032-005-4735-4
Frankel, O. H., and Brown, A. H. D. (1984). “Plant genetic resources today: a critical appraisal,” in Crop Genetic Resources: Conservation and Evaluation, eds J. H. W. Holden and J. T. Williams (London: George Allen and Unwin Ltd), 249–257.
Garris, A. J., Tai, T. H., Coburn, J., Kresovich, S., and McCouch, S. (2005). Genetic structure and diversity in Oryza sativa L. Genetics 169, 1631–1638. doi: 10.1534/genetics.104.035642
Ghaley, B. B., Christiansen, J. L., and Andersen, S. B. (2012). Genetic diversity for blast resistance of Bhutan rice landraces. Euphytica 184, 119–130. doi: 10.1007/s10681-011-0582-z
Giarrocco, L. E., Marassi, M. A., and Salerno, G. L. (2007). Assessment of the genetic diversity in Argentine rice cultivars with SSR markers. Crop Sci. 47, 853–858. doi: 10.2135/cropsci2005.07.0198
Gowda, M., Barman, S. R., and Chattoo, B. B. (2006). Molecular mapping of a novel blast resistance gene Pi38 in rice using SSLP and AFLP markers. Plant Breed. 125, 596–599. doi: 10.1111/j.1439-0523.2006.01248.x
Gupta, P. K., Kulwal, P. L., and Jaiswal, V. (2014). Association mapping in crop plants: opportunities and challenges. Adv. Genet. 85, 109–147. doi: 10.1016/B978-0-12-800271-1.00002-0
Hanamaratti, N. G., Prashanthi, S. K., Salimath, P. M., Hanchinal, R. R., Mohankumar, H. D., Parameshwarappa, K. G., et al. (2008). Traditional landraces of rice in Karnataka: reservoirs of valuable traits. Curr. Sci. 94, 242–247.
Harlan, J. R. (1975). Crops and Man. Madison, WI: American Society of Agronomy and Crop Science Society of America, 295.
Hayashi, K., Yoshida, H., and Ashikawa, I. (2006). Development of PCR-based allele-specific and InDel marker sets for nine rice blast resistance genes. Theor. Appl. Genet. 113, 251–260. doi: 10.1007/s00122-006-0290-6
Hayashi, N., Inoue, H., Kato, T., Funao, T., Shirota, M., Shi mizu, T., et al. (2010). Durable panicle blast-resistance gene Pb1 encodes an atypical CC-NBS-LRR protein and was generated by acquiring a promoter through local genome duplication. Plant J. 64, 498–510. doi: 10.1111/j.1365-313X.2010.04348.x
Hittalmani, S., Parco, A., Mew, T. W., Zeigler, R. S., and Huang, N. (2000). Fine mapping and DNA marker-assisted pyramiding of the three major genes for blast resistance in rice. Theor. Appl. Genet. 100, 1121–1128. doi: 10.1007/s001220051395
Horst, L., and Wenzel, G. (2007). Molecular Marker Systems in Plant Breeding and Crop Improvement: Biotechnology in Agriculture and Forestry, Vol. 55. Berlin: Springer.
Hu, J., Zhu, J., and Xu, H. M. (2000). Methods of constructing core collections by stepwise clustering with three sampling strategies based on the genotypic values of crops. Theor. Appl. Genet. 101, 264–268. doi: 10.1007/s001220051478
Imam, S., Alam, N. P., Mandal, M., Variar, P., and Shukla, P. (2014). Molecular screening for identification of blast resistance genes in north east and eastern Indian rice germplasm (Oryza sativa L.) with PCR based markers. Euphytica 196, 199–211. doi: 10.1007/s10681-013-1024-x
Jia, Y., Redus, M., Wang, Z., and Rutger, J. N. (2004). Development of a SNLP marker from the Pi-ta blast resistance gene by tri-primer PCR. Euphytica 138, 97–105. doi: 10.1023/B:EUPH.0000047079.42768.4d
Jia, Y., Wang, Z., and Singh, P. (2002). Development of dominant rice blast Pi-ta resistance gene markers. Crop Sci. 42, 2145–2149.
Jin, L., Lu, Y., Xiao, P., Sun, M., Corke, H., and Bao, J. (2010). Genetic diversity and population structure of a diverse set of rice germplasm for association mapping. Theor. Appl. Genet. 121, 475–487. doi: 10.1007/s00122-010-1324-7
Kang, H. M., Zaitlen, N. A., Wae, C. M., Kirby, A., Heckerman, D., Daly, M. J., et al. (2008). Effective control of population structure in model organism association mapping. Genetics 178, 1709–1723. doi: 10.1534/genetics.107.080101
Kim, K. W., Chung, H. K., Cho, G. T., Ma, K. H., Chandrabalan, D., Gwag, J. G., et al. (2007). PowerCore: a program applying the advanced M strategy with a heuristic search for establishing core sets. Bioinformatics 23, 2155–2162. doi: 10.1534/genetics.107.080101
Kumbhar, S. D., Kulwal, P. L., Patil, J. V., Sarawate, C. D., Gaikwad, A. P., and Jadhav, A. S. (2015). Genetic diversity and population structure in landraces and improved rice varieties from India. Rice Sci. 22, 99–107. doi: 10.1016/j.rsci.2015.05.013
Laha, G. S. (2017). Draft Proceedings of 52nd Annual Rice Research Group Meetings, All India Coordinated Rice Improvement Project (AICRIP). New Delhi: Indian Council of Agricultural Research.
Li, W., Lei, C., Cheng, Z., Jia, Y., Huang, D., Wang, J., et al. (2008). Identification of SSR markers for a broad-spectrum blast resistance gene Pi20(t) for marker-assisted breeding. Mol. Breed. 22, 141–149. doi: 10.1007/s11032-008-9163-9
Liu, K., and Muse, S. V. (2005). PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21, 2128–2129. doi: 10.1093/bioinformatics/bti282
Liu, W., Shahid, M. Q., Bai, L., Lu, Z., Chen, Y., Jiang, L., et al. (2015). Evaluation of genetic diversity and development of a core collection of wild rice (Oryza rufipogon Griff.) populations in China. PLoS ONE 10:e0145990. doi: 10.1371/journal.pone.0145990
Lou, Q., Ma, C., Wen, W., Zhou, J., Chen, L., Feng, F., et al. (2011). Profiling and association mapping of grain metabolites in a subset of the core collection of Chinese rice germplasm (Oryza sativa L.). J. Agric. Food. Chem. 59, 9257–9264. doi: 10.1021/jf201602g
Madan, K. M., Madhav, M. S., Prasad, M. S., RamaDevi, S. J. S., RamKumar, G., and Viraktamath, B. C. (2012). Analysis of population structure of Magnaporthe grisea using genome specific microsatellite markers. Curr. Trends Biotechnol. Pharm. 6, 173–182.
Madhav, M. S., Laha, G. S., Padmakumari, A. P., Somasekhar, N., Mangrauthia, S. K., and Viraktamath, B. C. (2013). “Phenotyping Rice for Molecular Plant Breeding,” in Phenotyping for Plant Breeding, eds S. Panguluri and A. Kumar (New York, NY: Springer), doi: 10.1007/978-1-4614-8320-5_1
Mahender, A., Swain, D. M., Gitishree, D., Subudhi, H. N., and Rao, G. J. N. (2012). Molecular analysis of native Manipur rice accessions for resistance against blast. Afr. J. Biotechnol. 11, 1321–1329. doi: 10.5897/AJB11.2178
Miah, G., Rafii, M. Y., Ismail, M. R., Puteh, A. B., Rahim, H. A., Islam, K. N., et al. (2013). A review of microsatellite markers and their applications in rice breeding programs to improve blast disease resistance. Int. J. Mol. Sci. 14, 22499–22528. doi: 10.3390/ijms141122499
Murray, M. G., and Thompson, W. F. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 8, 4321–4325. doi: 10.1093/nar/8.19.4321
Nachimuthu, V. V., Muthurajan, R., Duraialaguraja, S., Sivakami, R., Pandian, B. A., Ponniah, G., et al. (2015). Analysis of population structure and genetic diversity in rice germplasm using SSR markers: an initiative towards association mapping of agronomic traits in Oryza sativa. Rice 8:30. doi: 10.1186/s12284-015-0062-5
Nei, M. (1973). Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. U.S.A. 70, 3321–3323. doi: 10.1073/pnas.70.12.3321
Ni, J., Colowit, P. M., and Mackill, D. J. (2002). Evaluation of genetic diversity in rice subspecies using microsatellite markers. Crop Sci. 42, 601–607. doi: 10.2135/cropsci2002.6010
Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research—an update. Bioinformatics 28, 2537–2539. doi: 10.1093/bioinformatics/bts460
Perrier, X., and Jacquemoud-Collet, J. (2006). DARwin Software. Available at: http://darwin.cirad.fr/darwin
Powell, W., Morgante, M., Andre, C., Hanafey, M., Vogel, J., Tingey, S., et al. (1996). The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis. Mol. Breed. 2, 225–238. doi: 10.1007/BF00564200
Prasad, M. S., Madhav, M. S., Laha, G. S., Ladhalakshmi, D., Krishnaveni, D., Satendrakumar, M., et al. (2011). Rice Blast Disease and Its Management. Hyderabad: Directorate of Rice Research.
Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959.
Qu, S., Liu, G., Zhou, B., Bellizzi, M., Zeng, L., Dai, L., et al. (2006). The broad-spectrum blast resistance gene Pi9 encodes a nucleotide-binding site-leucine-rich repeat protein and is a member of a multigene family in rice. Genetics 172, 1901–1914. doi: 10.1534/genetics.105.044891
Ramkumar, G., Srinivasarao, K., Mohan, K. M., Sudarshan, I., Sivaranjani, A. K. P., Gopalakrishna, K., et al. (2011). Development and validation of functional marker targeting an InDel in the major rice blast disease resistance gene Pi54 (Pikh). Mol. Breed. 27, 129–135. doi: 10.1007/s11032-010-9538-6
Rao, E. S., Kadirvel, P., Symonds, R. C., Geethanjali, S., and Ebert, A. W. (2012). Using SSR markers to map genetic diversity and population structure of Solanum pimpinellifolium for development of a core collection. Plant Genet. Resour. 10, 38–48. doi: 10.1017/S1479262111000955
Ray, D. K., Mueller, N. D., West, P. C., and Foley, J. A. (2013). Yield trends are insufficient to double global crop production by 2050. PLoS ONE 8:e66428. doi: 10.1371/journal.pone.0066428
Ried, A., Hof, L., Felix, G., Rucker, B., Tams, S., Milczynska, E., et al. (2011). Construction of an integrated microsatellite and key morphological characteristic database of potato varieties on EU common catalogue. Euphytica 182, 239–249. doi: 10.1007/s10681-011-0462-6
Roy, S., Banerjee, A., Mawkhlieng, B., Misra, A. K., Pattanayak, A., Harish, G. D., et al. (2015). Genetic diversity and population structure in aromatic and quality rice (Oryza sativa L.) landraces from North-Eastern India. PLoS ONE 10:e0141405. doi: 10.1371/journal.pone.0129607
Roy, S., Banerjee, A., Pattanayak, A., Roy, S. S., Rathi, R. S., Misra, A. K., et al. (2014). Chakhao (delicious) rice landraces (Oryza sativa L.) of North-east India: collection, conservation and characterization of genetic diversity. Plant Genet. Resour. 12, 264–272. doi: 10.1017/S1479262113000580
Roy, S., Marndi, B. C., Mawkhlieng, B., Banerjee, A., Yadav, R. M., Misra, A. K., et al. (2016). Genetic diversity and structure in hill rice (Oryza sativa L.) Landraces collected from the North-Eastern Himalayas of India. BMC Genet. 17:107. doi: 10.1186/s12863-016-0414-1
Sharma, T. R., Chauhan, R. S., Singh, B. M., Paul, R., Sagar, V., and Rathore, R. (2002). RAPD and pathotype analysis of Magnaporthe grisea population from north-western Himalayan region of India. J. Phytopathol. 150, 649–656. doi: 10.1046/j.1439-0434.2002.00812.x
Shikari, A. B., Khanna, A., Krishnan, S. G., Singh, U. D., Rathore, R., Tonapi, V., et al. (2013). Molecular analysis and phenotypic validation of blast resistance genes Pita and Pita 2 in landraces of rice (Oryza sativa L.). Indian J. Genet. 73, 131–141. doi: 10.5958/j.0975-6906.73.2.020
Singh, N., Choudhury, D. R., Singh, A. K., Kumar, S., Srinivasan, K., Tyagi, R. K., et al. (2013). Comparison of SSR and SNP markers in the estimation of genetic diversity and population structure of Indian Rice varieties. PLoS ONE 8:e84136. doi: 10.1371/journal.pone.0084136
Surapu, R. B., Srinivas, P., Aruna, J., Vijay, S., Durga Rani, P. Ch, Reddy, N., et al. (2017). Survey of Magnaporthe grisea isolates around Andhra Pradesh and Telangana states, India. Int. J. Curr. Microbiol. App. Sci. 6, 61–70. doi: 10.20546/ijcmas.2017.605.008
Upadhyaya, H. D., Gowda, C. L. L., Pundir, R. P. S., Reddy, V. G., and Singh, S. (2006). Development of a core subset of finger Millet Germplasm using geographical origin and data on 14 quantitative traits. Genet. Resour. Crop. Evol. 53, 679–685. doi: 10.1007/s10722-004-3228-3
Vanniarajan, C., Vinod, K. K., and Pereira, A. (2012). Molecular evaluation of genetic diversity and association studies in rice (Oryza sativa L.). J. Genet. 91, 9–19. doi: 10.1007/s12041-012-0146-6
Varshney, R. K., Chabane, K., Hendre, P. S., Aggarwal, R. K., and Graner, A. (2007). Comparative assessment of EST-SSR, EST-SNP and AFLP markers for evaluation of genetic diversity and conservation of genetic resources using wild, cultivated and elite barleys. Plant Sci. 173, 638–649. doi: 10.1016/j.plantsci.2007.08.010
Wu, Y. Y., He, J. B., Li, A. H., Fang, N. Y., He, W. W., Dang, L. L., et al. (2016). Population structure analysis and association mapping of blast resistance in indica rice (Oryza sativa L.) landraces. Genet. Mol. Res. 15, 1–11. doi: 10.4238/gmr.15038254
Yang, G. P., Maroof, M. A. S., Xu, C. G., Zang, Q., and Baiyashev, R. M. (1994). Comparative analysis of microsatellite DNA polymorphism in landraces and cultivars of rice. Mol. Gen. Genet. 245, 187–194. doi: 10.1007/BF00283266
Yu, J. M., Pressoir, G., Briggs, W. H., Bi, I. V., Yamasaki, M., Doebley, J. F., et al. (2006). A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 38, 203–208. doi: 10.1038/nq1702
Zhang, D., Zhang, H., Wang, M., Sun, J., Qi, Y., Wang, F., et al. (2009). Genetic structure and differentiation of Oryza sativa L. in China revealed by microsatellites. Theor. Appl. Genet. 119, 1105–1117. doi: 10.1007/s00122-009-1112-4
Zhang, D., Zhang, H., Wei, X., Qi, Y., Wang, M., Sun, J., et al. (2007). Genetic structure and diversity of Oryza sativa L. in Guizhou, China. Chin. Sci. Bull. 52, 343–351. doi: 10.1007/s11434-007-0063-x
Zhang, H., Zhang, D., Wang, M., Sun, J., Qi, Y., Li, J., et al. (2011). A core collection and mini core collection of Oryza sativa L. in China. Theor. Appl. Genet. 122, 49–61. doi: 10.1007/s00122-010-1421-7
Zhang, P., Li, J., Li, X., Liu, X., Zhao, X., and Lu, Y. (2011). Population structure and genetic diversity in a rice core collection (Oryza sativa L.) investigated with SSR markers. PLoS ONE 6:e27565. doi: 10.1371/journal.pone.0027565
Zhou, B., Qu, S., Liu, G., Dolan, M., Sakai, H., Lu, G., et al. (2006). The eight amino-acid differences within three leucine-rich repeats between Pi2 and Piz-t resistance proteins determine the resistance specificity to Magnaporthe grisea. Mol. Plant Microbe Interact. 19, 1216–1228. doi: 10.1094/MPMI-19-1216
Keywords: rice landrace, blast resistance, genetic diversity, population structure, core set, association mapping
Citation: Umakanth B, Vishalakshi B, Sathish Kumar P, Rama Devi SJS, Bhadana VP, Senguttuvel P, Kumar S, Sharma SK, Sharma PK, Prasad MS and Madhav MS (2017) Diverse Rice Landraces of North-East India Enables the Identification of Novel Genetic Resources for Magnaporthe Resistance. Front. Plant Sci. 8:1500. doi: 10.3389/fpls.2017.01500
Received: 20 May 2017; Accepted: 14 August 2017;
Published: 29 August 2017.
Edited by:
Manoj Prasad, National Institute of Plant Genome Research, IndiaReviewed by:
Rakesh Singh, National Bureau of Plant Genetic Resources (ICAR), IndiaRomesh Kumar Salgotra, Sher-e-Kashmir University of Agricultural Sciences and Technology of Jammu, India
Copyright © 2017 Umakanth, Vishalakshi, Sathish Kumar, Rama Devi, Bhadana, Senguttuvel, Kumar, Sharma, Sharma, Prasad and Madhav. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maganti S. Madhav, c2hlc2h1MjRAZ21haWwuY29t
†These authors have contributed equally to this work.