Very low polymorphism in the germplasm typically used by breeding programs poses a significant bottleneck with regards to molecular breeding and the exploitation of breeding materials for quantitative trait analyses. California rice cultivars, derived from a very small base of temperate japonica germplasm and having a relatively brief breeding history, are a good example. In this study, we employed a reduced representation sequencing approach called Restriction Enzyme Site Comparative Analysis (RESCAN) to simultaneously identify and genotype single nucleotide polymorphisms (SNPs) in forty-five rice cultivars representing the majority of the 100 year-old breeding history in California. Over 20,000 putative SNPs were detected relative to the Nipponbare reference genome which enabled the identification and analysis of inheritance of pedigree haplotypes. Haplotype blocks distinguishing modern California cultivars from each other and from the ancestral short grain temperate japonica cultivars were easily identified. Reduced representation sequencing methods such as RESCAN are a valuable alternative to SNP chip genotyping and low coverage whole genome sequencing.
References
[1]
International Rice Genome Sequencing Project (2005) The map-based sequence of the rice genome. Nature 436: 793–800.
[2]
Yu J, Hu S, Wang J, Wong GK, Li S, et al. (2002) A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science 296: 79–92.
[3]
McNally KL, Childs KL, Bohnert R, Davidson RM, Zhao K, et al. (2010) Genome wide SNP variation reveals relationships among landraces and modern varieties of rice. Proc Natl Acad Sci USA 106: 12273–12278.
[4]
Yamamoto T, Nagasaki H, Yonemaru J, Ebana K, Nakajima M, et al. (2010) Fine definition of the pedigree haplotypes of closely related rice varieties by means of genome-wide discovery of single-nucleotide polymorphisms. BMC Genomics 11: 267.
[5]
Zhao K, Wright M, Kimball J, Eizenga G, McClung A, et al. (2010) Genomic diversity and introgression in O. sativa reveal the impact of domestication and breeding on the rice genome. PLoS One 5: e10780.
[6]
Zhao K, Tung CW, Eizenga GC, Wright MH, Ali ML, et al. (2011) Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat. Comm 2: 467.
[7]
Huang X, Wei X, Sang T, Zhao Q, Feng Q, et al. (2010) Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet 42: 961–967.
[8]
Huang X, Zhao Y, Wei X, Li C, Wang A, et al. (2011) Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet 44: 32–39.
[9]
Huang X, Kurata N, Wei X, Wang Z, Wang A, et al. (2012) A map of rice genome variation reveals the origin of cultivated rice. Nature 490: 497–501.
[10]
Davey JW, Hohenlohe PA, Etter PD, Boone JQ, Catchen JM, et al. (2011) Genome-wide genetic marker discovery and genotyping using next generation sequencing. Nat Rev Genet 12: 499–510.
[11]
Deschamps S, Llaca V, May GD (2012) Genotyping-by-sequencing in plants. Biology 1: 460–483.
[12]
Poland J, Rife TW (2012) Genotyping-by-sequencing for plant breeding and genetics. Plant Genome 5: 92–102.
[13]
Altshuler D, Pollara VJ, Cowles CR, Van Etten WJ, Baldwin J, et al. (2000) An SNP map of the human genome generated by reduced representation shotgun sequencing. Nature 407: 513–516.
[14]
Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, et al. (2008) Rapid SNP Discovery and Genetic Mapping Using Sequenced RAD Markers. PLoS One 3: e3376.
[15]
Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, et al. (2011) A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One 6: e19379.
[16]
van Orsouw NJ, Hogers RC, Janssen A, Yalcin F, Snoeijers S, et al. (2007) Complexity reduction of polymorphic sequences (CRoPS): a novel approach for large-scale polymorphism discovery in complex genomes. PLoS One 2: e1172.
[17]
Van Tassell CP, Smith TP, Matukumalli LK, Taylor JF, Schnabel RD, et al. (2008) SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries. Nat Methods 5: 247–252.
[18]
Andolfatto P, Davison D, Erezyilmaz D, Hu TT, Mast J, et al. (2011) Multiplexed shotgun genotyping for rapid and efficient genetic mapping. Genome Res 21: 610–617.
[19]
Monson-Miller J, Sanchez-Mendez DC, Fass J, Henry IM, Tai TH, et al. (2012) Reference genome-independent assessment of mutation density using restriction enzyme-phased sequencing. BMC Genomics 13: 72.
[20]
Dilday RH (1990) Contribution of ancestral lines in the development of new cultivars of rice. Crop Sci 30: 905–911.
[21]
McKenzie KS, Johnson CW, Tseng ST, Oster JJ, Brandon DM (1994) Breeding improved rice cultivars for temperate regions: a case study. Aust J Exp Agr 34: 897–905.
[22]
Rutger JN, Peterson ML, Hu CU (1977) Registration of 'Calrose 76' rice. Crop Sci 17: 987.
[23]
Carnahan HL, Johnson CW, Tseng ST, Mastenbroek JJ (1978) Registration of 'M9' rice. Crop Sci 18: 357–358.
[24]
Tseng ST, Carnahan HL, Johnson CW, Oster JJ, Hill JE, et al. (1984) Registration of 'L-202' rice. Crop Sci 24: 1213.
[25]
Seymour DK, Filiault DL, Henry IH, Monson-Miller J, Ravi M, et al. (2012) Arabidopsis doubled haploids-rapid homozygous lines for quantitative trait locus mapping. Proc Natl Acad Sci USA 109: 4227–32.
[26]
DeAngelis MM, Wang DG, Hawkins TL (1995) Solid-phase reversible immobilization for the isolation of PCR products. Nucleic Acids Res 23: 4742–4743.
[27]
Lennon NJ, Lintner RE, Anderson S, Alvarez P, Barry A, et al. (2010) A scalable, fully automated process for construction of sequence-ready barcoded libraries for 454. Genome Biol 11: R15.
Kim SI, McKenzie K, Tai TH (2009) A molecular survey of SD1 alleles used in U.S. rice cultivars. SABRAO J 41: 25–40.
[31]
Marchini J, Howie B, Myers S, McVean G, Donnelly P (2007) A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet 39: 906–913.
[32]
Marchini J, Howie B (2010) Genotype imputation for genome-wide association studies. Nat Reviews 11: 499–511.
[33]
Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25: 1754–60.
[34]
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, et al. (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25: 2078–2079.
[35]
Cavalli-Sforza LL, Edwards AW (1967) Phylogenetic analysis. Models and estimation procedures. Am J Hum Genet 3: 233–57.
[36]
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution 28: 2731–2739.
[37]
Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–59.
[38]
Falush D, Stephens M, Pritchard JK (2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164: 1567–87.
