International Journal of Genomics
Volume 2016 (2016), Article ID 2128010, 8 pages
http://dx.doi.org/10.1155/2016/2128010

Research Article

Mapping of Novel QTL Regulating Grain Shattering Using Doubled Haploid Population in Rice (Oryza sativa L.)

Gyu-Ho Lee,¹ In-Kyu Kang,² and Kyung-Min Kim¹

¹Division of Plant Biosciences, School of Applied Biosciences, College of Agriculture & Life Science, Kyungpook National University, Daegu 702-701, Republic of Korea
²Department of Horticultural Science, College of Agriculture and Life Sciences, Kyungpook National University, Daegu 702-701, Republic of Korea

Received 15 April 2016; Revised 19 May 2016; Accepted 26 May 2016

Academic Editor: Rupesh Deshmukh

Copyright © 2016 Gyu-Ho Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The critical evolutionary step during domestication of major cereals was elimination of seed shattering because the easy-to-shatter trait in wild relatives results in a severe reduction in yield. In this study, we analyzed the QTLs associated with shattering employing a high-density genetic map in doubled haploid (DH) population of rice (Oryza sativa L.). A genetic linkage map was generated with 217 microsatellite markers spanning 2082.4 cM and covering 12 rice chromosomes with an average interval of 9.6 cM between markers based on 120 DHLs derived from a cross between Cheongcheong indica type cultivar and Nagdong japonica type cultivar. In the QTL analysis, five QTLs pertaining to the breaking tensile strength (BTS) were detected in 2013 and 2015. Two regions of the QTLs related to BTS on chromosome 1 and chromosome 6 were detected. Several important genes are distributed in 1 Mbp region of the QTL on chromosome 6 and they are related to the formation of abscission layer. We decide to name this QTL qSh6 and the candidate genes in the qSh6 region can be employed usefully in further research for cloning.

1. Introduction

Major cereals were first cultivated approximately 10,000 years ago, when several plant characteristics, such as adaptation of flowering time to regional areas, increase in the number of seeds, variation in plant architecture and seed shape, response to increased fertility, and prevention of seed shattering, were altered by selection. These efforts of selection have generated a wide variety of crops [1, 2]. Rice, which is the largest food crop cultivated worldwide, meets the caloric needs of millions of people and is used as a staple food by half of the world’s population. There are two distinct types of domesticated rice, Oryza sativa or Asian rice and Oryza glaberrima or African rice; each has unique domestication histories [3]. The critical evolutionary step during domestication of major cereals was elimination of seed shattering because the easy-to-shatter trait in wild relatives results in a severe reduction in yield [1, 3] and wild rice also disperses seeds freely at maturity to guarantee the propagation, while cultivated rice retains seeds on the straw to make harvest easy and decrease the loss of production [4]. The degree of shattering can be distinguished into several types, easy-shattering, moderate shattering, and hard or nonshattering. The type of shattering is important to consider when selecting combine machines to prevent loss of yield during harvest. Specifically, small head-feeding combines are usually employed to harvest hard or nonshattering varieties, while large combine harvester-threshers are used for the moderate shattering varieties. Therefore, measurement of grain shattering is critical to enable the best harvesting method and reduce yield loss [5]. A digital force gauge to measure the degree of shattering is useful [6–10]. The device can measure the greatest amount of force against the direction of pulling during removal of seeds from a plant (Kg/pull). In the different methods, the shattering trait can be measured as the percentage of fallen spikelets among total spikelets of each panicle [11, 12]. The bagging panicles method is similar method to evaluate grain shattering for scoring the proportion of seeds naturally shed, shed by hand, and not shed by hand as 2, 1, and 0, respectively [13]. Also shattering habit was evaluated in shattering scale by hand gripping and recorded it on 1−9 scale [14]. However, there are several inevitable factors such as water content, growth stage, and man-made errors.

The shattering habit of rice is considered to be under relatively simple genetic control when compared with other characteristics related to domestication. Several recessive genes associated with the formation of an abscission layer, sh2, sh4, and sh-h on chromosomes 1, 3, and 7, have been reported [5, 15, 16]. Sh3 and SHA1 on chromosome 4 from wild rice were also reported [3, 17]. The grain shattering is considered to be caused by abscission layer [5]. The morphology of the abscission layer can differ in many rice varieties that show varying degrees of shattering [18–20]. Accordingly, it is important to elucidate the molecular mechanism to determine why some varieties do not have abscission layers and have an easy-shattering trait.

Various segregation patterns ranging from being monogenic to multigenic and depending on the crosses have been found, indicating that the shattering degree of cereals is a complex trait influenced by many genetic and environmental factors [21]. Many genetic studies for the seed shattering habits have been carried out using populations of natural variation from a cross between indica variety with shattering type and japonica variety with nonshattering type [2, 6, 11, 13, 14, 22]. A major QTL of grain shattering in rice, qSH1, was revealed to encode a BEL1-type homeobox gene [1]. A single-nucleotide polymorphism in the 5′ regulatory region of the qSH1 gene has also been shown to lead to the loss of seed shattering because of the absence of an abscission layer formation. In this study, we tried to detect QTLs associated with shattering traits with candidate genes by development of a genetic linkage map using 122 DHLs derived from a cross between the indica cultivar Cheongcheong and japonica cultivar Nagdong.

2. Materials and Methods

2.1. Plant Materials

The DH population for construction of the genetic map was developed by anther culture of the F₁ derived from a cross between Cheongcheong and Nagdong. Cheongcheong is a Tongil type rice variety that has resistance to the brown plant hopper, while Nagdong is the main variety cultivated in the regional area and has a partial abscission layer on the pedicel tissues [10]. The CNDH population was cultivated in a paddy field in 2013 and 2015 in the experimental fields of Kyungpook National University at Gunwie in Korea.

2.2. DNA Isolation and Genotyping

The genomic DNA of the parents cultivars and 120 DHLs were isolated from fresh leaves by CTAB method. In total, 217 microsatellite markers were used to determine genotype of the population for construction of a genetic linkage map. A total of 24 μL consisted of 2 μL of 10–20 ng/μL template DNA, 1 μL of 5–10 pM of each primer, 0.1 μL of Taq polymerase (Imclone Biotech Co., IN5001), 0.375 μL of dNTPs mixture, 2.4 μL of 10x Ex buffer, and 17.125 μL of nuclease-free water (QIAGEN, category number 129114). Analyses were conducted using GeneAmp PCR System 2700 (Applied Biosystems, USA). The reaction consisted of initial denaturation for 5 minutes at 96°C followed by 34 cycles of 30 seconds at 96°C, 30 seconds at 55°C, and 1 minute at 72°C, followed by final extension for 8 minutes at 72°C and then storage at 4°C PCR. Finally, the PCR products were confirmed by electrophoresis (QIAGEN, QIAxcel) for genotyping.

2.3. Construction of a Genetic Map

The data derived from PCR genotyping were input into Mapmaker 3.0 to calculate the genetic distance by applying the Kosambi function. A total of 217 microsatellite markers were used to generate a genetic linkage map spanning the 12 rice chromosomes using Mapmaker/EXP 3.0. The microsatellite markers were provided by the National Institute of Crop Science (NICS), RDA. The linkage among markers was determined by the group command at an LOD threshold of 3.0 and a recombination fraction of 0.30.

2.4. Phenotype Evaluation and QTL Analysis

The panicles samples were harvested at 40 days after heading (DAH) and then placed at room temperature for a week in 2013 while all the lines of the population were harvested at 40 DAH and placed at room temperature for three days in 2015. Five panicles per line were harvested for the measurement and shattering degree of spikelet was investigated as breaking tensile strength. To increase the precision and consistency, the shattering degree of spikelet was measured using a digital force gauge (Imada, Japan). The panicles were fixed to the force gauge to measure breaking tensile strength by pulling and bending grain spikelet using forceps (Figure 1). In total, 50 grains with 5 pedicles were removed by pulling and bending randomly. The Windows QTL cartographer 2.5 was used for analysis of QTLs associated with grain shattering. Composite interval mapping (CIM) was conducted for the entire genome by WinQTLCart 2.5 at a threshold of LOD 2.5 [23] after input of all required data and the heritability is set for calculation at 0.8 as a default value by the program. QTLs detection was attempted three times according to the date of harvest.

(a)

(b)

(c)

(d)

Figure 1: Measurement of seed shattering degree using a force gauge. (a) The device used to measure shattering degree was fixed to a stand reversely. (b) A holder was used to link tightly with a panicle stem. (c, d) The breaking tensile strength (BTS) was measured by bending and pulling a spikelet.

3. Results

3.1. Construction of the Genetic Map

Total markers consisting of 217 microsatellites were used to construct a genetic linkage map of the CNDH population derived from a cross between the indica cultivar, Cheongcheong, and japonica cultivar, Nagdong. The genetic linkage map of 12 chromosomes is generated by Mapmaker 3.0 covering 2082.4 cM in the entire genome with an average of 9.6 cM between markers (Supplemental Figure  S1) (see Supplementary Material available online at http://dx.doi.org/10.1155/2016/2128010).

3.2. Phenotype Evaluation

The bending and pulling strength of the DH population at 40 days after heading (DAH) in 2013 and 2015 were measured by unit of gravity force (gf) using a digital force gauge. The determination of shattering degree in 2013 and 2015 has shown that the shattering value of Cheongcheong is of easy-shattering type and Nagdong is of moderately difficult-shattering type (Figure 1). The two different breaking tensile strengths (BTSs) by pulling and bending have a correlation over 87% with significance at the level of 0.01 (Table 2, Supplemental Figure  S3). In addition, the BTS has shown the other grain characteristics including grain width, thickness, and thousand grain weight are related to the BTS below 30% with 0.05 of significance level (Table 2). The pulling strength in 2015 only was similar to normal distribution.

3.3. Analysis of QTLs

A total of 3 QTLs for pulling strength and 4 QTLs for bending were confirmed on chromosomes 1, 4, 6, 9, and 10 in 2013 while 3 QTLs for pulling strength on chromosomes 1, 2, 6, and 4 QTLs for bending strength were detected during 2015 (Table 2, Figure 3). The qPS1 was represented as 5.81 as the highest LOD value with 11% of phenotypic variation in 2013 while 6.93 of LOD for qPS6 with 13% of phenotypic variation was the highest in 2015. The qPS1 and qPS10 were detected in the same region of qBS1 on RM11966–RM11849 and qBS10 on RM25219–RM25036 in 2013. The QTLs confirmed on chromosomes 1 (qPS1) and 6 (qPS6) also were placed on regions, RM11966–RM11849, and RM20632–RM439 in 2015, respectively. The position of qPS6 in 2015 was equal to qPS6 in 2013 and the QTLs for pulling and bending strength with around 3.7 of LOD and 0.08% of phenotypic variation on chromosome 10 were detected on the same region, RM25219–RM25181, in 2013.

4. Discussion

Grain shattering is considered to be an important trait because it can affect patterns of harvest using combine machine and yield [5]. In this study, 217 microsatellite markers were selected from 1402 microsatellite markers by polymorphism test for construction of a genetic linkage map and it was employed to analyze QTLs associated with grain shattering. The grain shattering can be measured in different methods. The bagging method is to obtain the ratio of grains naturally shed because the pedicels of rice are sensitive to shaking in harvest period. Therefore, the shattering was measured in the ratio of fallen spikelets in the entire spikelet of each panicle [11, 12]. It was also evaluated using two bagged panicles per plant and the proportion of seeds naturally shed, shed by hand gripping, and nonshed by hand gripping was recorded by scores 2, 1, and 0, respectively [13]. The shattering degree was measured by hand gripping and recorded on 1–9 scale [14]. However, many factors such as growth environment, growth stage, water content, and human errors can influence the shattering evaluations. Recently, the grain shattering has also been measured using a digital force gauge that can provide more precise data than other methods by pulling or bending grain pedicels [1, 5, 10]. In the shattering investigation, the parental cultivars of the DH population, Cheongcheong and Nadong, have shown that they are represented as of easy-shattering and moderately difficult shattering type, respectively (Table 1). However, the shattering degree of the parents is represented to be increased in 2015 (Figure 2). It can be assumed that it is due to the different drying time after harvest. There is a gap for three days between 2013 and 2015.

Table 1: Evaluation for shattering degree by force gauge.

Trial
2013
2015

PS (gf)
BS (gf)
PS (gf)
BS (gf)

Cheongcheong
98.5 ± 31.0
10.9 ± 3.5
157.4 ± 19.1
6.45 ± 4.05

Nagdong
219.7 ± 24.6
38.1 ± 2.9.
261.6 ± 16.5
34.7 ± 10.21

DH population
126.7 ± 79.2
18.9 ± 16.0
161.2 ± 74.6
17.0 ± 14.8

Gravity force; mean ± standard deviation.

Table 2: Analysis of the correlation coefficients for BTS and grain characteristics.

Year
Trait
PS
BS
GL
GW
GT
RLW
TGW

2013
PS
1
0.882
−0.077
0.245
0.226
0.157
0.221

BS

1
−0.146
0.285
0.228
−0.221
0.119

GL

1
−0.444
−0.339
−0.861
0.368

GW

1
0.784
−0.827
0.554

GT

1
−0.633
0.600

RLW

1
−0.076

GW

1

2015
PS
1
0.872
−0.164
0.292
0.306
−0.231
0.203

BS

−0.204
0.250
0.236
−0.245
0.116

GL

−0.444
−0.339
0.861
0.368

GW

1
0.784
−0.827
0.554

GT

1
−0.633
0.600

RLW

1
−0.0765

TGW

1

Significant at the level of 0.05 and significant at the level of 0.01. PS: pulling strength; BS: bending strength; GL: grain length; GW: grain width; GT: grain thickness; RLW: ratio of length/width; TGW: thousand grain weight.

(a)

(b)

Figure 2: The frequency distribution for BTS in DH population. BTS was measured as pulling strength (a) and bending strength (b) for two years. The black arrows are represented as Cheongcheong while the white arrow indicates Nagdong.

Figure 3: The chromosomal locations of the QTLs related to pulling strength (PS) and bending strength (BS) in 2013 and 2015. The rectangles and diamonds are the QTL associated with pulling and bending strength.

Many QTLs associated with shattering trait have been reported from the previous studies (Table 4). In this study, five QTLs were detected on chromosomes 1, 4, 6, 9, and 10 in 2013 and also five QTLs were identified on chromosomes 1, 2, 6, and 8 in 2015. The QTLs associated with pulling strength and bending strength on chromosome 1 were named qPS1 and qBS1 according to the QTL nomenclature [4, 24]. However, the qPS1 and qBS1 detected in 2013 and 2015 have shared the region of qSH1 on chromosome 1 with 68.6% of phenotypic variation [1]. The SHA1 on chromosome 4 can be considered to be a major gene related to grain shattering and also the QTL for bending strength on chromosome 4 (qBS4) has been detected on the opposite region of chromosome 4 compared with SHA1. However, there was only one QTL related to BTS showing relatively low LOD value among the detected QTLs in the present study. Therefore, we tried to focus on other QTLs for grain shattering newly found on chromosome 6 (qPS6) with the highest LOD value in 2015 (Table 3). Although two QTLs, sh6 and qSH6-1, related to grain shattering on chromosome 6 were found in other researches [9, 14], the QTLs on chromosome 6 in current study have been detected newly at different position. The qPS6 was detected from the QTL analysis in 2013 and 2015 and the regions of qPS6 and qBS6 were same location flanked by RM20632–RM439. The correlation analysis has shown that pulling and bending strength have over 80% of highly significant correlation at the level of 0.01 (Table 2). Thus, the two QTLs, the qPS6 and qBS6, were changed in the equal name of qSh6 and we tried to find candidate genes for qSH6 distributed on the 1 Mbp region.

Table 3: The QTL analysis related to pulling and bending strength in 2013 and 2015.

Year
Chr.
Locus
LOD
Add.

Marker interval
Allele

2013
1
qPS1
5.81
29.57
0.11
RM11966-RM11849
Cheongcheong

6
qPS6
4.25
24.64
0.08
RM20632-RM439

10
qPS10
3.73
22.80
0.07
RM25219-RM25181

1
qBS1
3.43
4.12
0.06
RM11966-RM11849
Cheongcheong

4
qBS4
3.27
4.79
0.08
RM1205

9
qBS9
2.78
3.75
0.05
RM566-RM24288

10
qBS10
3.76
4.77
0.08
RM25219-RM25181

2015
1
qPS1
5.14
26.15
0.11
RM11966-RM11849
Cheongcheong

2
qPS2
4.68
23.29
0.09
RM12915-RM12925

6
qPS6
6.93
28.70
0.13
RM20632-RM439

1
qBS1-1
4.95
6.39
0.14
RM10458-RM11194
Cheongcheong

1
qBS1-2
4.26
5.16
0.11
RM11966-RM11849

6
qBS6
3.13
4.33
0.07
RM20632-RM439

8
qBS8
3.03
4.69
0.08
RM149-RM23191

Chromosome, additive effect, and phenotypic variation.

Table 4: Comparison of QTL identified from previous studies.

Traits
Chr. number
Population type (size)
Marker (number)
Parental cultivars
QTL name
References

PNS
1, 4, 5, 8
F₂ (172)
RFLP, SSR, AFLP (348)
Aijiao Nante (indica), P16 (O. rufipogon)
sh1, sh4, sh5, sh8
Xiong et al. (1999) [11]

PNS
1, 4, 8, 11
RILs (125)
RFLP, isozyme morphological (147)
Pei-Kuh (indica), W1944 (O. rufipogon)
qSHT-1, qSHT-4, qSHT-8, qSHT-11
Cai and Morishima (2000) [13]

PNS
1
DH (151)
SSR (68)
Miara (japonica), C6′ (japonica)
Sh-2
Bres-Patry et al. (2001) [12]

PNS
1, 3, 6
RIL (120)
SSR (124)
Hwayeongbyeo (japonica), W1944 (O. rufipogon)
sh1, sh3, sh6
Lee et al. (2005) [14]

PS
1, 2, 5, 11, 12
F₂ (156)
RFLP, RAPD (609)
Kasalath (indica), Nipponbare (japonica)
qSH1
Konish et al. (2006) [1]

PS
3, 4, 8
F₂ (304)
SSR
CL16 (indica), IRGC 80470 (O. nivara)
sh3, sh4, sh8
Li et al. (2006) [7]

PS, BS
7
F₂ (240)
SSR (77)
Hsh (japonica), Blue & Gundil (japonica)
sh-h
Ji et al. (2006) [5]

SR
1, 3, 4, 11, 12
F₂ (376)
SSR
Teqing (indica), YJCWR (O. rufipogon)
SHA1
Lin et al. (2007) [3]

PS
1, 6
CSSL (103)
SSR (132)
93-11 (indica), Nipponbare (japonica)
qSH1-1, qSH6-1
Zhu et al. (2009) [9]

SS, AL, PS, BS
1, 3, 4, 5, 8, 9
DH (120)
SSR, STS (182)
Samgang (indica), Nagdong (japonica)
Qss1, Qss3, Qss4, Qss5-5, Qal1, Qal5-1, Qps8, Qbs9
Qin et al. (2010) [10]

PS: pulling strength, BS: bending strength, PNS: percentage of natural shed, SR: shattering rate, SS: shattering scale, and AL: abscission layer.

The pulling and bending strength of pedicel in shattering trait have a considerable relationship with the existence of abscission layer [1, 5]. The abscission layer of pedicel is essential for seed propagation and it is well known that it is affected by many factors such as ethylene [25] and injury due to insect or fungal attack [26]. Therefore, Os06g6991100 encoding AP2 superfamily and describing ethylene responsive factor 121 (ERF121) and APETALA2/ethylene-responsive element binding protein 61 (AP2/EREBP61) is considered as one of the important candidate genes. The SHAT1 gene, which encodes APETALA2 transcription factor on chromosome 4, is required for seed shattering through specifying abscission zone development in rice [27]. Also the existence of many protein kinase-related and pathogen-related genes indicates that the shattering habit is affected by signaling, which can change in hormone, from many environmental factors [27]. Os06g069700 and Os06g0700700 encode heavy metal ATPase 1 (HMA1) and heavy metal ATPase 2 (HMA2). They are related to transport heavy metal, cadmium, zinc, and copper. The high accumulation of cadmium rice for phytoremediation indicated that cultivars can be developed with resistance in shattering and lodging by marker-assisted selection [28]. Os06g0701600 and Os06g0701700 that encode high-affinity K⁺ transporter 2;4 (HKT2;4) and HKT1 are related to control transportation of cNa⁺, K⁺, and Ca²⁺. The HKT rice gene, named OsHKT2;4, can function as both ion carriers and channels [29]. Calcium is one of the components of cell wall and known to maintain both cell wall and membrane in the abscission process. In particular the calcium pectate functions as the cementing substance between cells in the bean leaf [30, 31]. Furthermore, calcium can affect inhibiting the formation of abscission layer [32]. These studies indicate that OsHKT2;4 can affect the formation of abscission layer indirectly. The Os06g0701900 and Os06g0702000 encode small auxin-up RNA (SAUR). The exact function of SAUR is unknown but some of them play a role in auxin signaling with calcium and calmodulin [33]. Therefore, we can conclude that some of candidate genes distributed on qSh6 are involved in formation of abscission layer and recognition of ethylene directly or indirectly by many complex signalings. Among the 165 candidate genes that can be searched in the qSh6 region with 1 Mpb, some genes are considered as important in shattering habit. The present study will be useful in the separation and introduction of the genes to the parental cultivars for confirmation of gene function. We can presume the influence of the genes for grain shattering habit by previous researches indirectly. Thus, the comparison between expression of the genes and phenotype via gene introduction to parental cultivar after cloning the genes will be helpful to understand their function.

5. Conclusion

The shattering habit is really complex phenomenon for propagation of seed. We could find newly detected QTLs related to breaking tensile strength. Among the seven QTLs in 2015, the qSh6 that explains 13% of the phenotypic variation with 6.93 of LOD will be exploited for cloning in further research. The qSh6 was detected in QTL for pulling and bending strength region on chromosome 6 equally in 2015. About seven candidate genes describing ethylene response factor (ERF), heavy metal ATPase (HMA), high-affinity K+ transporter (HKT), and small auxin-up (SAUR) were found among the 165 gene loci in 1 Mbp region. They will be employed in further cloning research for introduction of the gene to the parental plant and for rice breeding.

Competing Interests

The authors declare that they have no competing interests regarding the publication of this paper.

Acknowledgments

This work was supported by a grant from the Next-Generation BioGreen 21 Program (no. PJ011257012016), Rural Development Administration, Republic of Korea.

References

S. Konishi, T. Izawa, S. Y. Lin et al., “An SNP caused loss of seed shattering during rice domestication,” Science, vol. 312, no. 5778, pp. 1392–1396, 2006.

T. Sang and S. Ge, “The puzzle of rice domestication,” Journal of Integrative Plant Biology, vol. 49, no. 6, pp. 760–768, 2007.

Z. Lin, M. E. Griffith, X. Li et al., “Origin of seed shattering in rice (Oryza sativa L.),” Planta, vol. 226, no. 1, pp. 11–20, 2007.

S. R. McCouch, Y. G. Cho, M. Yano et al., “Report on QTL nomenclature,” Rice Genetics Newsletter, vol. 14, pp. 11–12, 1997.

H.-S. Ji, S.-H. Chu, W. Jiang et al., “Characterization and mapping of a shattering mutant in rice that corresponds to a block of domestication genes,” Genetics, vol. 173, no. 2, pp. 995–1005, 2006.

W. C. Kennard, R. L. Phillips, and R. A. Porter, “Genetic dissection of seed shattering, agronomic, and color traits in American wildrice (Zizania palustris var. interior L.) with a comparative map,” Theoretical and Applied Genetics, vol. 105, no. 6-7, pp. 1075–1086, 2002.

C. Li, A. Zhou, and T. Sang, “Genetic analysis of rice domestication syndrome with the wild annual species, Oryza nivara,” New Phytologist, vol. 170, no. 1, pp. 185–194, 2006.

C. Li, A. Zhou, and T. Sang, “Rice domestication by reducing shattering,” Science, vol. 311, no. 5769, pp. 1936–1939, 2006.

W. Zhu, J. Lin, D. Yang et al., “Development of chromosome segment substitution lines derived from backcross between two sequenced rice cultivars, Indica recipient 93-11 and Japonica donor Nipponbare,” Plant Molecular Biology Reporter, vol. 27, no. 2, pp. 126–131, 2009.

Y. Qin, S.-M. Kim, X. Zhao et al., “Identification for quantitative trait loci controlling grain shattering in rice,” Genes & Genomics, vol. 32, no. 2, pp. 173–180, 2010.

L. Z. Xiong, K. D. Liu, X. K. Dai, C. G. Xu, and Q. F. Zhang, “Identification of genetic factors controlling domestication-related traits of rice using an F₂ population of a cross between Oryza sativa and O. rufipogon,” Theoretical and Applied Genetics, vol. 98, no. 2, pp. 243–251, 1999.

C. Bres-Patry, M. Lorieux, G. Clément, M. Bangratz, and A. Ghesquière, “Heredity and genetic mapping of domestication-related traits in a temperate japonica weedy rice,” Theoretical and Applied Genetics, vol. 102, no. 1, pp. 118–126, 2001.

H.-W. Cai and H. Morishima, “Genomic regions affecting seed shattering and seed dormancy in rice,” Theoretical and Applied Genetics, vol. 100, no. 6, pp. 840–846, 2000.

S.-J. Lee, C.-S. Oh, J.-P. Suh, S. R. McCouch, and S.-N. Ahn, “Identification of QTLs for domestication-related and agronomic traits in an Oryza sativa × O. rufipogon BC₁F₇ population,” Plant Breeding, vol. 124, no. 3, pp. 209–219, 2005.

S. Oba, N. Sumi, F. Fujimoto, and T. Yasue, “Association between grain shattering habit and formation of abscission layer controlled by grain shattering gene sh-2 in rice (Oryza sativa L.),” Japanese Journal of Crop Science, vol. 64, no. 3, pp. 607–615, 1995.

Y. Fukuta and T. Yagi, “Mapping of a shattering resistance gene in a mutant line SR-5 induced from an indica rice variety, Nan-jing11,” Breeding Science, vol. 48, no. 4, pp. 345–348, 1998.

M. Eiguchi and Y. Sano, “A gene complex responsible for seed shattering and panicle spreading found in common wild rices,” Rice Genetics Newsletter, vol. 7, pp. 105–107, 1990.

I. D. Jin, H. Terao, and J. Inouye, “On the cracking of abscission layer in Asian rice cultivar (Oryza sative L.),” Japanese Journal of Crop Science, vol. 51, no. 4, pp. 542–545, 1982.

I. Jln, J. Inouye, and N. Q. Ng, “Histological peculiarities of the abscission layers of African rice, Oryza glaberrima Steud. and its relation with degree of grain shedding,” Japanese journal of crop science, vol. 59, no. 3, pp. 475–480, 1990.

I. D. Jin, Y. H. Bae, and J. Inouye, “Formation and development of abscission layer between pedicel and rachilla, and changes in grain shedding during ripening in African rice, Oryza glaberrima steud,” Korean Journal of Crop Science, vol. 40, no. 1, pp. 103–112, 1995.

L. H. Tang and H. Morishima, “Variation and inheritance of seed shedding in weed rice,” Rice Genetics Newsletter, vol. 6, pp. 72–73, 1989.

M. J. Thomson, T. H. Tai, A. M. McClung et al., “Mapping quantitative trait loci for yield, yield components and morphological traits in an advanced backcross population between Oryza rufipogon and the Oryza sativa cultivar Jefferson,” Theoretical and Applied Genetics, vol. 107, no. 3, pp. 479–493, 2003.

Z.-B. Zeng, “Precision mapping of quantitative trait loci,” Genetics, vol. 136, no. 4, pp. 1457–1468, 1994.

S. R. McCouch and CGSNL, “Gene nomenclature system for rice,” Rice, vol. 1, no. 1, pp. 72–84, 2008.

T. H. Goodspeed and P. Davidson, “Controlled pollination in Nicotiana,” University of California Publications in Botany, vol. 5, pp. 429–434, 1918.

J. E. Coit and R. W. Hodgson, “An investigation of the abnormal shedding of young fruits of Washington navel orange,” University of California Publications in Agricultural Science, vol. 3, pp. 283–368, 1919.

Y. Zhou, D. Lu, C. Li et al., “Genetic control of seed shattering in rice by the APETALA2 transcription factor SHATTERING ABORTION1,” The Plant Cell, vol. 24, no. 3, pp. 1034–1048, 2012.

T. Abe, F. Taguchi-Shiobara, Y. Kojima et al., “Detection of a QTL for accumulating Cd in rice that enables efficient Cd phytoextraction from soil,” Breeding Science, vol. 61, no. 1, pp. 43–51, 2011.

W.-Z. Lan, W. Wang, S.-M. Wang et al., “A rice high-affinity potassium transporter (HKT) conceals a calcium-permeable cation channel,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 15, pp. 7089–7094, 2010.

H. C. Sampson, “Chemical changes accompanying abscission in Coleus blumei,” Botanical Gazette, vol. 66, no. 1, pp. 32–53, 1918.

F. T. Addicott, “Environmental factors in the physiology of abscission,” Plant Physiology, vol. 43, no. 9, pp. 1471–1470, 1968.

B. W. Poovaiah and A. C. Leopold, “Inhibition of abscission by calcium,” Plant Physiology, vol. 51, no. 5, pp. 848–851, 1973.

T. Yang and B. W. Poovaiah, “Molecular and biochemical evidence for the involvement of calcium/calmodulin in auxin action,” The Journal of Biological Chemistry, vol. 275, no. 5, pp. 3137–3143, 2000.