https://orcid.org/0000-0002-9880-1257
Figures
Abstract
Grain shape strongly influences the economic value and grain yield of rice. Thus, identifying quantitative trait loci (QTLs) for grain shape has been a longstanding goal in rice genetic research and breeding programs. Single nucleotide polymorphism (SNP) markers are ubiquitous in the rice genome and are more abundant and evenly distributed on the 12 rice chromosomes than traditional markers. An F2 population was genotyped using the RICE6K SNP array to elucidate the mechanisms governing grain shape. Thirty-five QTLs for grain shape were detected on 11 of 12 chromosomes over 2 years. The major QTL cluster qGS7 was detected in both years and displayed strong genetic effects on grain length and width, showing consistency with GL7/GW7. Some minor QTLs were also detected, and the effects of four QTLs on seed size were then validated using BC1F6 populations with residual heterozygous lines in each QTL region. Our findings provide insights into the molecular basis of grain shape as well as additional resources and approaches for producing hybrid high-yield rice varieties.
Citation: Chen J, Liu K, Zha W, Zhou L, Li M, Xu H, et al. (2021) Identification and verification of grain shape QTLs by SNP array in rice. PLoS ONE 16(11): e0260133. https://doi.org/10.1371/journal.pone.0260133
Editor: Zhukuan Cheng, The Institute of Genetics and Developmental Biology (IGDB) of the Chinese Academy of Sciences (CAS), China, CHINA
Received: July 18, 2021; Accepted: November 4, 2021; Published: November 22, 2021
Copyright: © 2021 Chen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: This work was supported by grants from the National Natural Science Foundation (32001529), the Hubei Academy of Agricultural Science Foundation (2021NKYJJ02), Key R&D program of Hubei Province (2020BBA031), Natural Science Foundation of Guangdong Province (2020A1515011390) and the Hubei Key Laboratory of Food Crop Germplasm and Genetic Improvement Foundation (2019lzjj02).
Competing interests: The authors have declared that no competing interests exist.
Introduction
It is critical to increase crop productivity through efficient breeding to meet the challenge of feeding a rapidly growing world population in a smaller area [1]. Rice (Oryza sativa) is a staple food crop worldwide and a key model cereal crop. Grain shape in rice is controlled by a combination of grain length, grain width, and grain thickness and is a major determinant of grain yield and grain appearance quality [2]. Therefore, it is essential to explore the genetic basis of grain shape and identify new quantitative trait loci (QTLs) for the improvement of grain yield and quality.
Most QTLs influence grain shape by regulating cell division or proliferation, and several major QTLs have been cloned and functionally characterized. GS3, the first cloned gene for grain shape, encodes a putative transmembrane protein and function as a negative regulator of grain length [3, 4]. The QTL qGL3/qGL3.1 encodes a protein phosphatase with a Kelch-like repeat domain and negatively controls grain length [5]. For grain width, GW2 encodes a novel RING-type protein with E3 ubiquitin ligase activity and regulates grain width [6]. GW5/GSE5 codes for a plasma membrane-associated protein and governs grain width; a deletion in the promoter region causes decreased expression and increases grain width [7, 8]. GS5 encodes a putative serine carboxypeptidase and is a positive regulatory factor of grain size [9], and the GS5 allele with a higher expression level resulted in increased grain size. Similarly, GW8, which encodes the transcription factor OsSPL16 with an SBP domain, is also a positive regulator of grain size [10]. GL7/GW7, contains a 17.1-kb tandem duplication and encodes a protein homologous to Arabidopsis thaliana LONGIFOLIA proteins, increases grain length and reduces grain width [11, 12]. The cloning and functional characterizations of these genes have greatly enriched our knowledge of the molecular mechanisms determining grain shape and improve genomic breeding in rice [13, 14]. However, additional genes with minor effects on grain shape remain largely unknown.
QTL mapping is frequently used in genetic studies to identify new genes and QTLs. Compared with traditional molecular markers, single nucleotide polymorphisms (SNPs) are more abundant and uniformly distributed on the 12 rice chromosomes and are regarded as the most desirable molecular marker [15]. Next-generation sequencing technologies have identified millions of DNA sequences, improving the discovery of SNPs [16, 17]. Whole-genome sequencing and DNA microarrays [18] are high-throughput genotyping platforms for constructing high-density genetic linkage maps with SNPs. Both methods have been applied in gene discovery [19–21], and microarrays have been widely used in studies of humans, Arabidopsis, and rice to produce high-resolution linkage maps [22–24].
This study constructed a high-resolution genetic linkage map using a high-density SNP array. Linkage analyses were performed in F2 and its derived F3 populations to determine the genetic basis of grain size in rice. Our results provide insights into the construction of linkage maps with gene arrays, and the identified QTLs can potentially improve rice grain size.
Materials and methods
Plant materials and field planting
An F2 population was constructed by crossing the long-grain variety, O. sativa cv. Pusa, with the medium-grain variety, O. sativa cv. H2613S. The male parent (Pusa) has the longest grain in a collection of 529 core germplasms. The female parent H2613S is a well-known male sterile line characterized by excellent plant architecture and several advantages regarding quality and resistance (Fig 1A). Five hundred F2 plants were obtained from the cross between H2613S and Pusa. Two markers closely linked to the male sterile gene tms5 (TMS1F: AAGTTGCCACCCTCTTTCAG, TMS1R: TGTGTGAAGGGGTGCTACAG; TMS2F: GCAAAAGCTCAAGCCAGAGT, TMS2R: TCTCAGGCACCGTCAATGTA) were used to detect the fertility of these F2 lines. One hundred fertile plants with a homozygous genotype at the tms5 locus were used as the mapping population. The F2 population was planted in 2017, and the F3 population derived from the F2 population was planted in 2018, during normal growing seasons at the experimental field in Wuhan, China. H2613S was used as the recurrent parent to construct the backcross BC1F1 populations. Then, four respectively derived plants in the BC1F1 were self-pollinated to obtain BC1F2 populations. BC1F6 residual heterozygous lines were generated using seed descent and marker-assisted selection and were planted in 2020 under the same experimental conditions to confirm the effects of QTLs. All accessions were sown in the seedling nursery, and 25-day-old seedlings were transplanted into two-row plots (12 plants, with 26.5 cm of space between the rows and 16.5 cm of space between plants within a row).
Trait measurements and statistical analyses
Harvested rice grains were air-dried and stored at room temperature for at least 3 months, and all fully-filled rice grains were collected from each plant to investigate grain shape traits, including grain length, grain width, grain length-to-width ratio, and thousand-grain weight. These traits were measured at a high-throughput rice phenotyping facility [25]. P-values for phenotypic coefficients were calculated using a two-sided t-test and the cor.test function in R software [26].
Genotyping and linkage map construction
The F2 population was genotyped using the RICE6K array, which contained 5,102 SNP and InDel markers. There were 1,596 unique variants between Pusa and H2613S (S1 Fig). The genetic linkage map of the population was calculated using the MapMaker 3.0 program. Analyses of main-effect QTLs were conducted in the mapping population by composite interval mapping using WinQTLCart version 2.5. In the analyses, the likelihood ratio and t-test were combined to test the significance of the single-locus QTL additive effects. The likelihood ratio and t values corresponding to P = 0.001, equivalent to LOD (log likelihood value) ≥ 3.0, were used as the threshold for selecting the putative main-effect QTLs. The peak points of the likelihood ratio in the linkage map were considered the putative positions of the QTLs. The relative contribution of a genetic component (%) was calculated as the proportion of phenotypic variance explained by that component in the selected model [27].
Results
Grain shape variation in the F2 and its derived F3 populations
The phenotypic values of grain length and width were both higher in 2018 than in 2017 (Fig 1B–1D). The grain length and thousand-grain weight were normally distributed and showed a continuous variation in both years (Fig 1B and 1E), and grain width also showed a continuous variation and approximately followed a normal distribution (Fig 1C). The distribution of the length-to-width ratio was bimodal in both years, and the two peaks were close to each other (Fig 1D).
Data on grain length, width, length-to-width ratio, and thousand-grain weight are shown in Tables 1 and 2. These phenotypes were strongly correlated with each other, except thousand-grain weight and grain length-to-width ratio. Grain width was significantly and negatively correlated with grain length and thousand-grain weight was significantly and positively associated with grain length and grain width in both years.
Mapping of QTLs for grain shape
LOD thresholds for selecting QTLs for grain shape traits were determined by 1000 permutations at P = 0.05. Thirty-five QTLs for grain shape traits were distributed on all chromosomes, except chromosome 11, in the two years (Table 3 and Fig 2). The phenotypic variance explained by each QTL ranged from 0.12% to 74.26%.
GL, grain length; GW, grain width; LWR, grain length-to-width ratio; TGW, thousand-grain weight.
For grain length, a major QTL was detected on chromosome 7. The QTL (qGL7) was flanked by R0721009 and R0726296 and was stabled in both years; it explained 59.20% of the phenotypic variation in 2017 and 52.84% in 2018. In addition, eight minor QTLs (qGL2.2, qGL6, qGL8.1, qGL1, qGL2.1, qGL4, qGL8.2, and qGL10) were detected on chromosomes 1, 2, 4, 6, 8, and 10. Of these, qGL2.2, qGL6, and qGL8.1 were detected in both years and explained 0.12% to 7.23% of the phenotypic variation, and qGL1, qGL2.1, qGL4, qGL8.2, and qGL10 were detected in one year and explained 0.23% to 9.24% of the phenotypic variation.
For grain width, a major QTL was detected on chromosome 7 in the two years. The QTL (qGW7) was flanked by R0721009 and R0726296 and explained 46.81% of the phenotypic variation in 2017 and 30.97% in 2018. Moreover, eight minor QTLs (qGW1, qGW3, qGW9 qGW12, qGW2.1, qGW2.2, qGW6.1, and qGW6.2) were detected on chromosomes 1, 2, 3, 6, 9, and 12. Of these, qGW1, qGW3, qGW9, and qGW12 were detected in both years and explained 1.58% to 17.29% of the phenotypic variation, and qGW2.1, qGW2.2, qGW6.1, and qGW6.2 were detected in one year and explained 2.88% to 26.70% of the phenotypic variation.
Marker-trait linkage analyses revealed that eight QTLs for the grain length-to-width ratio were detected on chromosomes 1, 2, 4, 7, 8, 9, and 12. Among them, the QTL flanked by F0721726 and R0726296 on chromosome 7 (qLWR7) was found in the two years, and explained 74.26% of the phenotypic variation in 2017 and 44.96% in 2018. qLWR12 was flanked by R1203509 and R1207724 on chromosome 12 and was detected in both years. This QTL explained 8.36% of the phenotypic variation in 2017 and 9.63% in 2018. Moreover, one and five QTLs were detected in both years and one year, respectively; these explained 0.99% to 13.54% of the phenotypic variation.
Nine QTLs for thousand-grain weight were detected on chromosomes 2, 3, 4, 5, 7, 8, 10, and 12. Three minor QTLs were detected in both years and explained 2.68% to 14.47% of the phenotypic variation. The remaining QTLs were detected in one year and explained 2.57% to 23.64% of the phenotypic variation.
High precision of QTL mapping via high-density SNP array
The region flanked by R0710169 and R0726296 on chromosome 7 contained four QTLs, qGL7, qGW7, qLWR7, and qTGW7.2, which were collectively termed qGS7. These QTLs controlled grain length, grain width, grain length-to-width ratio, and thousand-grain weight, respectively. The top points of qGL7, qGW7, and qLWR7 were all located between R0724564 and R0724700, corresponding to 136 kb (Fig 3). According to the Nipponbare 6.0 reference genome, 15 candidate genes are located in this 136-kb region, including a major grain size gene, GL7/GW7. Some minor candidate QTLs possessing gene linkage or pleiotropy were observed as well. The QTL qGW12 was identified in both years, and the top point was located between R1220728 and R1221279, corresponding to 551 kb (S2 Fig). Exploring and utilizing the linkage or pleiotropy of these QTLs underlying grain traits would be beneficial for genetic research and grain yield improvement. These results support the high precision of QTL mapping via high-density SNP arrays.
Red line indicates the position of the peak SNP.
Validation of the genetic effects of major and minor QTLs
The BC1F6 residual heterozygous line populations of qGW1, qGW3, qGS7, and qGW9 were analyzed to validate the genetic effects of these QTLs (Fig 4). The qGW1 locus from female parents increased grain width by 0.05 mm. The qGW3 locus from female parents increased grain width by 0.04 mm and had little effect on grain length. The qGS7 locus from male parents increased grain length by 0.20 mm and decreased grain width by 0.15 mm. All these QTLs had significant effects on grain length-to-width ratio. These findings suggest that these QTLs have great potential for grain size improvement.
The BC1F6 residual heterozygous line populations of each QTL were analyzed to validate the genetic effects of these QTLs. Blue bars represent alleles from pusa, Green bars represent alleles from H2613s. Error bars are based on standard deviation of each genotype. P-values based on two-tailed t-test.
Discussion
Rice is a staple food, feeding more than half of the world’s population [1]. Rice grain shape is one of the most important factors contributing to rice yield. Therefore, the identification of major or minor QTLs for grain shape is a longstanding goal in rice genetic research and breeding programs. In this study, a RICE6K SNP array with high-density SNP markers was used in genotyping to identify QTLs for grain shape traits.
Grain size and grain weight are agronomic traits with high heritability in rice. However, the phenotypic values of grain length, grain width and grain weight were all lower in 2017 than in 2018 (Fig 1) because these quantitative traits are largely affected by the environment. Ambient temperature was higher in the summer of 2017 than in the summer of 2018, potentially affecting grain filling. Genetic factors may also contribute to phenotypic differences between the two study years because the dominance effect was higher in F2 than F3.
In the present study, 35 QTLs for grain shape were detected in two years. Among them, a QTL cluster on chromosome 7 contributed to grain length and grain width and displayed strong genetic effects, and the qGS7 allele tended to increase grain length and decrease grain width. This region also contained the gene GL7/GW7, cloned previously [11, 12]. GL7/GW7 encodes a TONNEAU1-recruiting motif protein and acts as a semi-dominant locus. Upregulation of GL7/GW7 expression leads to a decrease in cell width and an increase in cell length of epidermal cells of the inner and outer glumes. In our study, the genomic DNA fragments from the candidate gene in Pusa and H2613S were sequenced and analyzed using primers NGSP11F and 210QCF, and a 1,471 bp fragment with a functional allele associated with longer grain length was amplified from the male parent, consistent with previous results (S3 Fig). Thus GL7/GW7 might be the candidate gene for qGS7.
Several minor QTLs for grain shape appeared as novel QTLs that have not been reported previously. In both years, the qGL2.2 region was close to the SNP marker KASP0224 and showed moderate effects on grain length and grain width. The qGW9 region was close to the SNP marker KASP0915 and was responsible for grain width and length-to-width ratio in both years. Similarly, qGW1 was identified in both years and affected on grain width. These regions contained several minor QTLs for grain width, consistent with our previous study in other populations. GS3 and GW5 are two major genes respectively for grain length and grain width, and strongly contribute to grain length and width variance in rice landraces [28]. However, these genes were not identified in QTL mapping in our samples, probably because the genotypes of GS3 and GW5 in the two parental lines were identical.
Since the sequencing of two rice genome references (93–11 and Nipponbare), millions of SNPs are currently available on the NCBI website and are used as genetic markers. Compared with traditional markers, such as RFLP and SSR, SNPs are ubiquitous in the rice genome and are highly abundant and evenly distributed on the 12 rice chromosomes, improving the development of linkage maps using RFLP and SSR markers. In addition, SNP genotyping using a SNP array is cheaper than whole-genome sequencing and can efficiently obtain the genotypes of genome-wide SNPs for several individuals from an experimental population [29]. Nonetheless, the development of high-quality linkage maps is limited by the population size and number of markers. Thus, the resolution of QTL mapping is positively correlated with the density of the genetic map and population size [30], and the number of polymorphic markers evaluated in QTL analyses increases as the use of SNP arrays expands. In our study, the top SNPs (qGL7, qGW7, and qLWR7) were located between R0724564 and R0724700, a span of 136 kb, and this region is highly consistent with GL7/GW7. These studies suggest that SNP arrays are an efficient and reliable genotyping method to identify new genes or QTLs.
Most rice cultivars approved in China show medium grain size, good quality, and high yield; nonetheless, these characteristics can be improved in many cultivars [31]. Therefore, breeding rice cultivars with a combination of desirable traits is challenging. Although some genes have been cloned and many QTLs have been identified in rice cultivars, the genetic basis and molecular mechanism underlying grain shape and quality are incompletely understood [32]. A previous study showed that mature grains of the near-isogenic line NIL-gs9 were slender than those of Nipponbare, and the transparency of mature rice was greatly improved in NIL-gs9, which showed no or very little chalkiness [33]. The authors designed QTL pyramiding based on combinations of alleles qgw8 and qgs3 by molecular marker-assisted selection and developed a slender rice variety, substantially improving grain quality [10]. These results suggest that grain width may be correlated with the degree of chalkiness. The major QTL qGS7 can potentially produce slender and less chalky grains. The QTL qGW3 affected grain width but not grain length, potentially influencing chalkiness. qGW9 had a strong effect on grain length but a slight effect on grain width. Additional characterization and validation studies combined with fine mapping are necessary to identify novel genes and alleles and improve our understanding of the relationship between these candidate genes and the observed phenotypes, providing researchers with effective genetic tools to select beneficial alleles.
Supporting information
S1 Fig. Haplotype maps of example lines/plants detected by RICE6K array.
(a) One Pusa/H2613S F2 line. (b) The whole Pusa/H2613S F2 lines. Each short line at the chromosomes indicates the position of a single nucleotide polymorphism (SNP); AA, female parental homozygous genotype; BB, male parental homozygous genotype; and AB, heterozygous genotype.
https://doi.org/10.1371/journal.pone.0260133.s001
(TIF)
S2 Fig. Local qGW12 mapping result of grain width.
Red line indicates the position of the peak SNP.
https://doi.org/10.1371/journal.pone.0260133.s002
(TIF)
S3 Fig. The gel image of the PCR amplification with the 210QCF and NGSP11F primer pair using DNA samples from Pusa, H2613s and C815s.
https://doi.org/10.1371/journal.pone.0260133.s003
(TIF)
S1 Table. List of primers used in this study.
https://doi.org/10.1371/journal.pone.0260133.s004
(DOCX)
S2 Table. The genotyping results of qGW1, qGW3, qGS7 and qGW9 in the BC1F6 residual heterozygous lines.
https://doi.org/10.1371/journal.pone.0260133.s005
(XLSX)
Acknowledgments
We gratefully thank Mei Yang for logistic work, and two anonymous reviewers for critical evaluation of the study.
References
- 1. Khush GS. What it will take to feed 5.0 billion rice consumers in 2030. Plant Mol Biol. 2005; 59(1):1–6. pmid:16217597
- 2. Huang RY, Jiang LR, Zheng JS, Wang TS, Wang HC, Huang YM, et al. Genetic bases of rice grain shape: so many genes, so little known. Trends Plant Sci. 2012; 18:218–226. pmid:23218902
- 3. Fan CC, Xing YZ, Mao HL, Lu TT, Han B Xu CG, et al. GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein. Theor Appl Genet. 2006; 112(6):1164–1171. pmid:16453132
- 4. Takano-Kai N, Jiang H, Kubo T, Sweeney M, Matsumoto T, Kanamori H, et al. Evolutionary history of GS3, a gene conferring grain length in rice. Genetics. 2009; 182(4):1323–1334. pmid:19506305
- 5. Qi P, Lin YS, Song XJ, Shen JB, Huang W, Shan JX, et al. The novel quantitative trait locus GL3.1 controls rice grain size and yield by regulating Cyclin-T1;3. Cell Res. 2012; 22(12):1666–1680. pmid:23147796
- 6. Song XJ, Huang W, Shi M, Zhu MZ, and Lin HX. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat Genet. 2007; 39(5):623–630. pmid:17417637
- 7. Shomura A, Izawa T, Ebana K, Ebitani T, Kanegae H, Konishi S, et al. Deletion in a gene associated with grain size increased yields during rice domestication. Nat Genet. 2008; 40(8):1023–1028. pmid:18604208
- 8. Weng JF, Gu SH, Wan XY, Gao H, Guo T, Su N, et al. Isolation and initial characterization of GW5, a major QTL associated with rice grain width and weight. Cell Res. 2008; 18(12):1199–209. pmid:19015668
- 9. Li YB, Fan CC, Xing YZ, Jiang YH, Luo LJ, Sun L, et al. Natural variation in GS5 plays an important role in regulating grain size and yield in rice. Nat Genet. 2011; 43:1266. pmid:22019783
- 10. Wang SK, Wu K, Yuan QB, Liu XY, Liu ZB, Lin XY, et al. Control of grain size, shape and quality by OsSPL16 in rice. Nat Genet. 2012;44(8):950–954. pmid:22729225
- 11. Wang SK, Li S, Liu Q, Wu K, Zhang JQ, Wang SS, et al. The OsSPL16-GW7 regulatory module determines grain shape and simultaneously improves rice yield and grain quality. Nat Genet. 2015; 47(8):949–954. pmid:26147620
- 12. Wang YX, Xiong GS, Hu J, Jiang L, Yu H, Xu J, et al. Copy number variation at the GL7 locus contributes to grain size diversity in rice. Nat Genet. 2015; 47(8):944–948. pmid:26147619
- 13. Hu J, Wang YY, Fang YX, Zeng LJ, Xu J, Yu HP, et al. A Rare Allele of GS2 Enhances Grain Size and Grain Yield in Rice. Mol Plant. 2015; 8:1455–1465. pmid:26187814
- 14. Che RH, Tong HN, Shi BH, Liu YQ, Fang SR, Liu DP, et al. Control of grain size and rice yield by GL2-mediated brassinosteroid responses. Nat Plants. 2015; (2):15195. pmid:27250747
- 15. Gupta PK, Rustgi S, and Mir RR. Array-based high-throughput DNA markers for crop improvement. Heredity. 2008; 101(1):5–18. pmid:18461083
- 16. Ansorge WJ. Next-generation DNA sequencing techniques. N Biotechnol. 2009; 25:195–203. pmid:19429539
- 17. Huang XH, Feng Q, Qian Q, Zhao Q, Wang L, Wang AR, et al. High-throughput genotyping by whole-genome resequencing. Genome Res. 2009; 19(6):1068–1076. pmid:19420380
- 18. Yu HH, Xie WB, Li J, Zhou FS, and Zhang QF. A whole-genome SNP array (RICE6K) for genomic breeding in rice. Plant Biotechnolo J. 2013; 12:8–37. pmid:24034357
- 19. Huang XH, Wei XH, Sang T, Zhao Q, Feng Q, Zhao Y, et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet. 2010; 42(11):961. pmid:20972439
- 20. Huang XH, Zhao Y, Wei XH, Li CY, Wang AH, Zhao Q, et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2011; 44:32–39. pmid:22138690
- 21. Yu H, Xie W, Wang J, Xing Y, Xu C, X Li, et al. Gains in QTL Detection Using an Ultra-High Density SNP Map Based on Population Sequencing Relative to Traditional RFLP/SSR Markers. PLoS One. 2011; 6(3):e17595. pmid:21390234
- 22. Sham A, Moustafa K, Al-Shamisi S, Alyan S, Iratni R, and AbuQamar S. Microarray analysis of Arabidopsis WRKY33 mutants in response to the necrotrophic fungus Botrytis cinerea. PLoS One. 2017; 12(2):e0172343. pmid:28207847
- 23. Tan C, Han ZM, Yu HH, Zhan W, Xie WB, Chen X, et al. QTL scanning for rice yield using a whole genome SNP array. J Genet Genomics. 2013; 40(12):629–638. pmid:24377869
- 24. Pasquini G, Barba M, Hadidi A, Faggioli Fet, Negri R, Sobol I, et al. Oligonucleotide microarray-based detection and genotyping of Plum pox virus. J Virol Methods. 2008; 147(1):118–126. pmid:17920703
- 25. Yang WN, Guo ZL, Huang CL, Duan LF, Chen GX, Jiang N, et al. Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice. Nat Commun. 2014; 5:5087. pmid:25295980
- 26. Yndgaard F, and Solberg SO. Analysing genebank collections using "R": making trait information widely available to users. Gene Conserve. 2016; 15:01–52.
- 27. Xia D, Zhou H, Qiu L, Jiang HC, Zhang QL, Gao GJ, et al. Mapping and verification of grain shape QTLs based on an advanced backcross population in rice. PLoS One. 2017; 12(11):e0187553. pmid:29145412
- 28. Zhou H, Li PB, Xie WB, Hussain S, Li YB, Xia D, et al. Genome-wide Association Analyses Reveal the Genetic Basis of Stigma Exsertion in Rice. Mol Plant. 2017; (10):634–644. pmid:28110091
- 29. Thomson MJ, Singh N, Dwiyanti MS, Wang DR, Wright MH, Perez FA, et al. Large-scale deployment of a rice 6 K SNP array for genetics and breeding applications. Rice. 2017; 10(1):40. pmid:28856618
- 30. Chen MS, Presting G, Barbazuk WB, Goicoechea JL, Blackmon B, Fang GC, et al. An Integrated Physical and Genetic Map of the Rice Genome. The Plant Cell. 2002; 14(3):537–545. pmid:11910002
- 31. Zhang QF. Strategies for developing Green Super Rice. Proc Natl Acad Sci USA. 2007; 104(42):16402–16409. pmid:17923667
- 32. Chen JX, Zhou H, Gu Y, Xia D, Wu B, Gao GJ, et al. Mapping and verification of grain shape QTLs based on high-throughput SNP markers in rice. Mol Breeding. 2019; 39(3):42. https://doi.org/10.1007/s11032-019-0955-x
- 33. Zhao DS, Li QF, Zhang CQ, Zhang C, Yang QQ, Pan LX, et al. GS9 acts as a transcriptional activator to regulate rice grain shape and appearance quality. Nat Communications. 2018; 9(1):1240. pmid:29588443.