Analysis of genetic architecture and favorable allele usage of agronomic traits in a large collection of Chinese rice accessions

Genotyping and phenotyping large natural populations provide opportunities for population genomic analysis and genome-wide association studies (GWAS). Several rice populations have been re-sequenced in the past decade; however, many major Chinese rice cultivars were not included in these studies. Here, we report large-scale genomic and phenotypic datasets for a collection mainly comprised of 1,275 rice accessions of widely planted cultivars and parental hybrid rice lines from China. The population was divided into three indica/Xian and three japonica/Geng phylogenetic subgroups that correlate strongly with their geographic or breeding origins. We acquired a total of 146 phenotypic datasets for 29 agronomic traits under multi-environments for different subpopulations. With GWAS, we identified a total of 143 significant association loci, including three newly identified candidate genes or alleles that control heading date or amylose content. Our genotypic analysis of agronomically important genes in the population revealed that many favorable alleles are underused in elite accessions, suggesting they may be used to provide improvements in future breeding efforts. Our study provides useful resources for rice genetics research and breeding.

[1]  Wei Chen,et al.  Genome-wide association analyses provide genetic and biochemical insights into natural variation in rice metabolism , 2014, Nature Genetics.

[2]  Xiaomin Feng,et al.  Improving rice grain length through updating the GS3 locus of an elite variety Kongyu 131 , 2018, Rice.

[3]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[4]  Yan Li,et al.  Sequencing and de novo assembly of a near complete indica rice genome , 2017, Nature Communications.

[5]  Ming Zhou,et al.  Integrated analysis of phenome, genome, and transcriptome of hybrid rice uncovered multiple heterosis-related loci for yield increase , 2016, Proceedings of the National Academy of Sciences.

[6]  Shan Wang,et al.  Divergent Hd1, Ghd7, and DTH7 Alleles Control Heading Date and Yield Potential of Japonica Rice in Northeast China , 2018, Front. Plant Sci..

[7]  Xiaomin Feng,et al.  Improving the blast resistance of the elite rice variety Kongyu-131 by updating the pi21 locus , 2019, BMC Plant Biology.

[8]  M. Matsuoka,et al.  Loss-of-function mutations of the rice GAMYB gene impair alpha-amylase expression in aleurone and flower development. , 2004, The Plant cell.

[9]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[10]  Qifa Zhang,et al.  Genome-wide association studies of 14 agronomic traits in rice landraces , 2010, Nature Genetics.

[11]  Yuanli Song,et al.  Interaction between temperature and photoperiod in regulation of flowering time in rice , 2012, Science China Life Sciences.

[12]  David H. Alexander,et al.  Fast model-based estimation of ancestry in unrelated individuals. , 2009, Genome research.

[13]  Y. Ouyang,et al.  Genome-wide dissection of segregation distortion using multiple inter-subspecific crosses in rice , 2019, Science China Life Sciences.

[14]  Kenneth L. McNally,et al.  Genomic variation in 3,010 diverse accessions of Asian cultivated rice , 2018, Nature.

[15]  Kazuyuki Doi,et al.  Ehd1, a B-type response regulator in rice, confers short-day promotion of flowering and controls FT-like gene expression independently of Hd1. , 2004, Genes & development.

[16]  Q. Qian,et al.  Natural Variation in the Promoter of GSE5 Contributes to Grain Size Diversity in Rice. , 2017, Molecular plant.

[17]  Qian Qian,et al.  Rational design of high-yield and superior-quality rice , 2017, Nature Plants.

[18]  Lei Zhang,et al.  Genomic analysis of hybrid rice varieties reveals numerous superior alleles that contribute to heterosis , 2015, Nature Communications.

[19]  Qian Qian,et al.  Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm , 2011, Nature Genetics.

[20]  Chengzhi Liang,et al.  Dissecting the genetic basis of heavy panicle hybrid rice uncovered Gn1a and GS3 as key genes , 2018, Theoretical and Applied Genetics.

[21]  Mark H. Wright,et al.  Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa , 2011, Nature communications.

[22]  Haiyang Liu,et al.  Validation and characterization of Ghd7.1, a major quantitative trait locus with pleiotropic effects on spikelets per panicle, plant height, and heading date in rice (Oryza sativa L.). , 2013, Journal of integrative plant biology.

[23]  Updating the elite rice variety Kongyu 131 by improving the Gn1a locus , 2017, Rice.

[24]  Lijun Luo,et al.  Genome-Wide Association Studies of Image Traits Reveal Genetic Architecture of Drought Resistance in Rice. , 2018, Molecular plant.

[25]  Qin He,et al.  Combinations of the Ghd7, Ghd8 and Hd1 genes largely define the ecogeographical adaptation and yield potential of cultivated rice. , 2015, The New phytologist.

[26]  Lei Wang,et al.  Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice , 2008, Nature Genetics.

[27]  Tao Huang,et al.  Genomic architecture of heterosis for yield traits in rice , 2016, Nature.

[28]  Heng Li,et al.  A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data , 2011, Bioinform..

[29]  F. Wang,et al.  The expression of Orysa;CycB1;1 is essential for endosperm formation and causes embryo enlargement in rice , 2009, Planta.

[30]  W. Zhai,et al.  TAC1, a major quantitative trait locus controlling tiller angle in rice. , 2007, The Plant journal : for cell and molecular biology.

[31]  M. Yano,et al.  Hd1, a Major Photoperiod Sensitivity Quantitative Trait Locus in Rice, Is Closely Related to the Arabidopsis Flowering Time Gene CONSTANS , 2000, Plant Cell.

[32]  Weibo Xie,et al.  Genome-wide association studies reveal that members of bHLH subfamily 16 share a conserved function in regulating flag leaf angle in rice (Oryza sativa) , 2018, PLoS genetics.

[33]  M. Matsuoka,et al.  Loss-of-Function Mutations of the Rice GAMYB Gene Impair α-Amylase Expression in Aleurone and Flower Development Article, publication date, and citation information can be found at www.plantcell.org/cgi/doi/10.1105/tpc.017327. , 2004, The Plant Cell Online.

[34]  Haiyang Wang,et al.  Days to heading 7, a major quantitative locus determining photoperiod sensitivity and regional adaptation in rice , 2014, Proceedings of the National Academy of Sciences.

[35]  Quan Xu,et al.  The DENSE AND ERECT PANICLE 1 (DEP1) gene offering the potential in the breeding of high-yielding rice , 2016, Breeding science.

[36]  Yu Li,et al.  Stability of QTL Across Environments and QTL-by-Environment Interactions for Plant and Ear Height in Maize , 2010 .

[37]  HuiHui Yu,et al.  Expression patterns of photoperiod and temperature regulated heading date genes in Oryza sativa , 2013, Comput. Biol. Chem..

[38]  Xuehui Huang,et al.  Ef-cd locus shortens rice maturity duration without yield penalty , 2019, Proceedings of the National Academy of Sciences.

[39]  L. Xiong,et al.  Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice , 2014, Nature Communications.

[40]  A. Fujiyama,et al.  A map of rice genome variation reveals the origin of cultivated rice , 2012, Nature.

[41]  Weibo Xie,et al.  Genetic Architecture of Natural Variation in Rice Chlorophyll Content Revealed by a Genome-Wide Association Study. , 2015, Molecular plant.

[42]  H. Kang,et al.  Variance component model to account for sample structure in genome-wide association studies , 2010, Nature Genetics.

[43]  Shaoqing Li,et al.  Artificial Selection of Gn1a Plays an Important role in Improving Rice Yields Across Different Ecological Regions , 2015, Rice.

[44]  Kosuke M. Teshima,et al.  Variations in Hd1 proteins, Hd3a promoters, and Ehd1 expression levels contribute to diversity of flowering time in cultivated rice , 2009, Proceedings of the National Academy of Sciences.

[45]  Benjamin E. Gaddy,et al.  Surfactant-enabled epitaxy through control of growth mode with chemical boundary conditions. , 2011, Nature communications.

[46]  Qifa Zhang,et al.  Genetic and molecular bases of rice yield. , 2010, Annual review of plant biology.

[47]  D. Reich,et al.  Principal components analysis corrects for stratification in genome-wide association studies , 2006, Nature Genetics.

[48]  L. Xiong,et al.  Breeding signatures of rice improvement revealed by a genomic variation map from a large germplasm collection , 2015, Proceedings of the National Academy of Sciences.

[49]  Susan McCouch,et al.  Genome-wide association and high-resolution phenotyping link Oryza sativa panicle traits to numerous trait-specific QTL clusters , 2016, Nature Communications.

[50]  F. Teng,et al.  The evidence for non‐additive effect as the main genetic component of plant height and ear height in maize using introgression line populations , 2009 .

[51]  Qun Xu,et al.  Pan-genome analysis highlights the extent of genomic variation in cultivated and wild rice , 2018, Nature Genetics.

[52]  Xiaomin Feng,et al.  Updating the Genome of the Elite Rice Variety Kongyu131 to Expand Its Ecological Adaptation Region , 2019, Front. Plant Sci..

[53]  N. Kurata,et al.  OsHAP3 genes regulate chloroplast biogenesis in rice. , 2003, The Plant journal : for cell and molecular biology.