Identification of QTL TGW12 responsible for grain weight in rice based on recombinant inbred line population crossed by wild rice (Oryza minuta) introgression line K1561 and indica rice G1025
BMC Genetics volume 21, Article number: 10 (2020)
Limited genetic resource in the cultivated rice may hinder further yield improvement. Some valuable genes that contribute to rice yield may be lost or lacked in the cultivated rice. Identification of the quantitative trait locus (QTL) for yield-related traits such as thousand-grain weight (TGW) from wild rice speices is desired for rice yield improvement.
In this study, sixteen TGW QTL were identified from a recombinant inbred line (RIL) population derived from the cross between the introgression line K1561 of Oryza minuta and the rice cultivar G1025. TGW12, One of most effective QTL was mapped to the region of 204.12 kb between the marker 2,768,345 and marker 2,853,491 of the specific locus amplified fragment (SLAF). The origin of TGW12 was tested using three markers nearby or within the TGW12 region, but not clarified yet. Our data indicated thirty-two open reading fragments (ORFs) were present in the region. RT-PCR analysis and sequence alignment showed that the coding domain sequences of ORF12, one MADS-box gene, in G1025 and K1561 were different due to alternative slicing, which caused premature transcription termination. The MADS-box gene was considered as a candidate of TGW12.
The effective QTL, TGW12, was mapped to a segment of 204.12 kb using RILs population and a MADS-box gene was identified among several candidate genes in the segment. The region of TGW12 should be further narrowed and creation of transgenic lines will reveal the gene function. TGW12 could be applied for improvement of TGW in breeding program.
Rice (Oryza sativa L.) is the world’s most important cereal crop as a staple food . High rice yield is needed to meet the requirement of rapidly increasing population. Grain weight is important to rice yield, and it is usually represented by thousand-grain weight. Many QTLs for rice grain weight have been mapped to rice chromosomes [2,3,4,5,6,7]. A dozen of QTLs/genes affecting grain weight have be cloned and functionally charaterized, such as the cytokinin oxidase/dehydrogenase (CKX) OsCKX2(Gn1a) , the transmembrane protein GS3 [9, 10 ], and its homolog DEP1 , the RING-type E3 ubiquitin ligase GW2 (grain width and weight 2) , the arginine-rich domain nuclear protein qSW5/GW5 , the serine carboxypeptidase GS5/qTGW5a , the Kelch-like domain qGL3/ OsPPKL1 , the SBP domain transcription factor GW8 (OsSPL16) , the IAA (indole-3-acetic acid)- glucose hydrolase protein TGW6 , the GNAT-like protein GW6a (Grain weight of chromosome 6)/OsglHAT1 , the ABC1-like kinase OsAGSW1 , the AP2 transcript factor qHD5 and OsSNB [20, 21], the cytochrome P450 protein GNS4 , the otubain-like protease WTG1 (WIDE AND THICK GRAIN 1) , the 16-kDa α-amylase/trypsin inhibitor RAG2 , the GSK3/SHAGGY-Like Kinase qTGW3 . Some genes such as GW2, WTG1, OsCKX2, TGW3, TGW6, qSW5/GW5qSW5/GW5, GS3, DEP1, qGL3/ OsPPKL1,and OsSNB negatively regulate grain weight, while others like GW6a/OsglHAT1, GNS4, GS5/qTGW5a, OsAGSW1, RAG2, and GW8 (OsSPL16) function as positive regulator of grain weight.
To date, QTLs/genes associated with TGW have been mostly cloned from the cultivated rice. However, it is known that the genetic resource of cultivated rice turned quite limited during the process of wild rice domestication, which may hinder the yield improvement of the cultivated rice. Wild rice species should contain many valuable genes that can be used for genetic improvements of cultivated rice . Thus, the genetic resource of wild rice species should be explored and used for rice high-yield breeding. It would be an effective way to widen the genetic basis of cultivated rice by introduction and application of favorable wild rice genes.
Oryza minuta (O. minuta) is a tetraploid wild rice that possesses a number of favorable yield related genes . Rahman et al. mapped 22 novel yield-related QTLs for 16 agronomic traits using a set of introgression line (IL) of O. minuta, and demonstrated that 57% of the QTLs were derived from O. minuta . In a previous study, we also detected 28 QTLs for yield-related traits using ILs derived from the backcross of IR24 (O. sativa L) and O. minuta, and found that 46.4% of notable QTLs were from O. minuta .
To further identify the favorable yield-related genes from O. minuta, a recombinant inbred line has been developed by crossing of K1561 to indica rice G1025 . K1561 is one out of 192 ILs derived from backcross progenies (BC4F2) of IR24 and O. minuta. It was produced from one time cross of IR24 and O. minuta, then four times backcross with IR24 as recurrent parent, and four times self-cross. K1561 shows excellent agronomic traits such as long panicles and high TGW. G1025 is an excellent restorer line that is widely used in Guangxi Province of China with dense grains but light TGW. In this study, QTL mapping was conducted on the advanced RILs population by SSR and SLAF markers. Thirteen QTLs or TGWs responsible for TGW were detected under fiver environments, and the most effective QTL TGW12 was mapped to the segment of 204.12 kb based on the high-density genome map constructed with SLAF. The candidate genes of TGW12 were preliminarily concluded, and one gene encoding MADS-box protein was considered as the putative candidate based on sequence alignment. The TGW12 allele for increasing TGW might originated from O. minuta and likely used for rice yield improvement.
The two parents G1025 and K1561 showed highly significant differences in TGW under five environments including two locations in Nan-Ning (NN), and Wu-Han (WH), in China with an average of 16.01 g and 32.07 g (Table 1), respectively. TGW values of the 201 individuals were mostly distributed between 20 and 30 g with an average of 24.60 g, 25.69 g, 22.80 g, 24.97 g, and 25.69 g in 2013NN (2013 in Nan-Ning), 2014NN, 2015NN, 2016NN, 2016WH (in 2016 in Wu-Han)(Table 1, Fig. 1, Additional file 1: Table S1), respectively. TGW values of 104, 106, 104, 107, and 107 out of the 201 individuals were smaller than the average in 2013NN, 2014NN, 2015NN, 2016NN, 2016WH, and TGW of the remaining individuals was either equal to or larger than the average (Additional file 1: Table S1).
QTL mapping of TGW by simple sequence repeats (SSR)
TGW QTLs were preliminarily detected by 300 SSR markers evenly distributed on the 12 chromosomes. The population were F6, F7 RILs derived from the cross of G1025 and K1561 planted in NN in 2013, 2014. Four QTLs TGW3, TGW7, TGW9.2, and TGW12 were stably detected on the chromosomes 3, 7, 9 and 12 in the two environments (Table 2). TGW12 had the greatest effect, which located on the region of RM247 and RM7003 (Table 2), so it was selected for further analysis. There were other 166 SSR markers (Additional file 2: Table S2) in the region based on the genome sequencing data of Nipponbare . The polymorphism of the 166 SSRs was firstly detected between the parental lines G1025 and K1561. As a result, nine SSRs showed polymorphism but only five displayed clear bands. The five SSRs were further used to detect F6, F7 RILs population in 2013 and 2014. Finally, TGW12 was mapped to the 5.1 cM region between RM27638 and RM27748 (Fig. 2).
QTL mapping of TGW by SLAF markers
We have developed 5521 SLAF markers by SLAF sequencing . To further map TGW QTLs, those SLAF markers were used to screen F8 RILs in NN in 2015 and F9 RILs in NN and WH in 2016, respectively. A total of eight QTLs were detected in three environments, namely, TGW7, TGW9.1, and TGW12 in 2015NN; TGW7, TGW9.2, and TGW12 in 2016NN; TGW7 and TGW12 in 2016WH (Table 2, Fig. 3). TGW7 and TGW12 were both detected in three environments, and TGW9.1 and TGW9.2 were each detected once. TGW7 explained the phenotypes for 8.01, 10.76, and 10.43% inheritance with an LOD of 6.94, 7.69, and 7.48 in 2015NN, 2016NN, and 2016WH, respectively, whereas TGW12 showed 22.36, 17.48, and 17.95% inheritance explaining for the phenotypes with an LOD of 15.42, 11.99, and 11.96 in the three environments, respectively (Table 2). TGW12 had a greater effect than TGW7, which was consistent with the results analyzed by SSR mapping (mentioned above). Further analysis for TGW12 was conducted by comparing the linkage map constructed by SSR and SLAF markers. Consequentially, nine SLAF markers fell into the region of RM27638-RM27748, and TGW12 was further narrowed to 204.12 kb region between SLAFs Marker 2,768,345 and Marker 2,853,491 (Fig. 2).
Evaluation of TGW12 phenotype and genetic origin identification of TGW12 segment
In order to evaluate whether the phenotypes were determined by TGW12, 16 out of the 201 RILs containing the TGW12 region were identified by means of the markers nearby the region (Fig. 4). Then, the phenotypes and genotypes of the 16 RILs were compared. All the 16 RILs with one or two segments of K1561 showed TGW increase than the recurrent parent G1025, suggesting TGW12 control TGW (Fig. 4). To clarify whether the increasing effect of TGW12 was originated from O. minuta, the genotypes of G1025, K1561, IR24, and O. minuta were examined using markers nearby or within TGW12. The genotype of K1561 was the same as that of IR24 but different from that of G1025 and O. minuta on the sites of RM27638 and RM27748, which are nearby TGW12 (Figs. 2 and 5). However, the genotype of K1561 was the same as those of IR24 and O. minuta, but it was different from that of G1025 at the site of Marker 2,768,345, which is linked with TGW12. It was hard to draw a conclusion whether TGW12 originated from IR24 or O. minuta based on the above results. It has been suggested that translocation through centric break-fusion occurred more frequently than recombination in the introgression lines with interspecific cross, which didn’t always result in an O. minuta chromosome arm onto a complete or incomplete O. sativa chromosome [30, 31]. Thus, TGW12 origination remains to be determined in the near future. It is feasible to compare sequence of TGW12 candidate among O. minuta, IR24, and K1561 once it was fine mapped.
Preliminary prediction of candidate genes for TGW12
Analysis of annotated genes indicated that 32 ORFs located in the 204.12 kb region based on Nipponbare genome annotation (http://rice.plantbiology.msu.edu) (Table 3). Among them, 13 ORFs encoded functional proteins and 19 ORFs were annotated as transposon/retrotransposon proteins, hypothetical proteins, or expressed protein. It is worth noting that there were four transcription factors (TFs) among the functional proteins: two MADS-box proteins (ORF12, ORF14), one ZF-HD protein (ORF24), and one B-box zinc finger protein (ORF27). Because TFs play crucial roles in regulation of plant growth and development , the four TFs were considered preferentially as putative candidate genes of TGW12. Reverse transcript (RT)-PCR were conducted to amplify the CDS (coding domain sequence) of the four transcription factors from the parents G1025 and K1561. Sequence comparison indicated that the amplified sequence of ORF12 in K1561 was 56 bp shorter than that of G1025, which resulted in premature transcription termination, consequently leading to a peptide with 45 amino acid residues only in K1561, whereas ORF12 in G1025 encoded a protein with 202 amino acid residues (Additional file 3: Figure S1). Further analysis revealed that lack of 56 bp of ORF12 in K1561 was due to alternative splicing (AS) in the first extron (Fig. 6b, c), which causes premature termination of ORF12 translation (Additional file 3: Figure S1). There were no difference in the CDS of ORF14, ORF24, and ORF27 between K1561 and G1025 (data not shown). Thus, the MADS-box (ORF12) was likely one putative candidate of TGW12. However, possibility of other nine functional proteins and the hypothetical proteins, or expressed protein independently or collectively affecting TGW could not be excluded. Further investigation is required.
Rice is one of the most important staple food widely consumed by one-half world’s peopulation and more product is needed as population increases in future. However, further yield improvement of rice is constrained by the narrow genetic basis of cultivated rice varieties. Wild rice species are good candidates to explore the genetic resource for valuably genes to further enhance rice productivity .
O. minuta possesses a number of outstanding genes associated with resistance and yield . Lots of QTLs for yield related traits have been identified using the introgression lines (IL) consisting of O. minuta segments [27, 28]. In this study, sixteen TGWs were detected using the advanced RILs population under five environments (Table 2). Among them, the most effective QTL, TGW12, was mapped on chromosome 12. Eight QTLs for TGW (AQAG040, AQCF014, AQAG053, AQGP079, AQDR045, AQDR047, AQCY020, CQAS153) had been previously mapped to Chromosome 12 (http://archive.gramene.org/). Location comparison indicated that TGW12 was partially overlapped with AQDR045. TGW12 and AQDR045 were respectively located in the regions of 4,037,811–6,150,143 and 1,589,200–5,829,185 on Chromosome 12 in the physical map of Nipponbare genome , suggesting that there might be a major QTL controlling grain weigh in this region. However, AQDR045 was mapped by using the popu1ation derived from two cultivated rice Lemont and Teqing , whereas TGW12 was mapped using the population of one O. minuta introgression line and one cultivated line. Although it was uncertain whether TGW12 allele originated from O. minuta or IR24, it showed a great increasing effect for TGW and could be directly applied in the breeding program.
There were 32 annotated ORFs in the TGW12 region, among which four ORFs were identified as TFs based on sequence alignment and considered as TGW12 candidates due to their regulatory roles in plant growth and development. Sequence analysis indicated that ORF12 in K1561 was truncated, with 45 amino acid residues only, due to an AS event while the counterpart in G1025 possessed a full length of the protein. No sequence difference in the other three TFs (ORF14, ORF24, and ORF27) was found between K1561 and G1025. ORF12 encoded a MADS-box protein, which possesses a highly conserved DNA-binding MADS domain and is involved predominantly in developmental processes . In Arabidopsis, there are 107 genes encoding MADS-box proteins , and almost all of them are involved in the process of flower and seed development . In rice, 75 MADS-box genes were identified, and more than 20 were transcribed during the stages of panicle and seed development . In addition, alternative splicing of one MADS-box transcription factor OsMADS1 encoded by OsLG3b (Os03g0215400) controls grain length and yield in japonica rice . Our results suggested that ORF12 could be critical for the function of TGW12, even functioned as TGW12. However, further studies are required to examine the function of other ORFs located in the region.
In this study, an effective QTL TGW12 related to the trait of thousand-grain weight in rice was mapped to a segment with 204.12 kb using RILs population derived from the cross progenies of one O. minuta introgression line and one cultivated rice. Out of 32 ORFs located in the region of TGW12, ORF12 encoded a MADS-box protein could be crucial for the TGW12 function. Further investigation is required to validate this speculation.
Plant materials and field trials
O. minuta (Acc. No. 101133) and IR24 were kindly provided by the International Rice Germplasm Centre of the International Rice Research Institute. The parental line G1025 was kindly provide by Rice Research Institute of Guangxi Academy of Agricultural Science. The parental lines K1561 was developed by our lab in Rice Research Institute of Guangxi Academy of Agricultural Science. It is produced by one time cross of IR24 and O. minuta, then four times backcross with IR24 as recurrent parent, and four times self-cross. The parental lines G1025 and K1561 along with 201 F6, F7 RILs were planted in Nanning (NN) from February to July in 2013 and 2014, respectively. The parental lines along with 201 F8 RILs were planted in NN (February to July) in 2015, and parents and F9 RILs were planted in NN (February to July) and Wuhan (WH) from May to October in 2016, respectively. The phenotypes of parents and RILs were collected to map TGW based on SSR or SLAFs. Grain weight was calculated based on 200 grains and converted to TGW after harvesting and sun-drying. The mean values of ten plants were used as input data to identify QTLs (Additional file 1: Table S1).
SSR, linkage, and QTL analysis
DNA was extracted from fresh leaves following the CTAB procedure . SSR markers were used to analyze a polymorphism between the parents (Additional file 2: Table S2). SSR were synthesized according to published sequences . Polymerase chain reaction (PCR) was conducted in a 15 μL volume as follow: 50 ng of template DNA, 0.3 μL of 10 mM each dNTPs, 0.5 units of Taq DNA polymerase, 1.5 μL of 10 × PCR buffer with Mg2+, and 0.5 μL of 10 μM forward and reverse primers. The reaction conditions was carried out as an initial denaturation at 94 °C for 5 min, followed by 35 cycles of 30 s at 94 °C, 30 s at 56 °C, and 30 s at 72 °C, with a final extension at 72 °C for 10 min. PCR products were separated on 6% polyacrylamide denaturing gels, and the bands were revealed by the silver-staining protocol .
Single nucleotide polymorphism (SNP) genotyping, linkage map construction and QTL analysis
Genomic DNA was extracted from fresh leaves of the parents and RILs by CTAB . Quantified DNA was used for SLAF sequencing by an Illumina HiseqTM 2500 . SLAF markers, developed in previous work, were used for genotyping, linkage map construction and QTL analysis for TGW in this study as described by Zhu et al. .
Derived cleaved amplified polymorphic sequences (dCAPS) marker development
dCAPS marker was developed for SLAF Marker2768345. Primers were designed according to dCAPS Finder 2.0 (http://helix.wustl.edu/dcaps/dcaps.html). PCR were conducted in a 20 μL volume as follow: 100 ng of template DNA, 0.5 μL of 10 mM each dNTPs, 1 units of Taq DNA polymerase, 2.0 μL of 10 × PCR buffer with Mg2+, and 0.5 μL of 10 μM forward and reverse primers. The reaction conditions was carried out as an initial denaturation at 94 °C for 3 min, followed by 35 cycles of 30 s at 94 °C, 30 s at 60 °C, and 30 s at 72 °C, with a final extension at 72 °C for 5 min. PCR products were digested with EcoR V (Takara, China) for 4 h at 37 °C, then were resolved in 2% agrose gel to genotype.
Candidate gene prediction and RT-PCR analysis
The predicted genes in the target region of QTL were analyzed according to the annotation of Nipponbare reference genome2. RT-PCR was conducted as described by Sha et al. . In brief, total RNA was extracted from roots, leaves, stems, shoot apical meristem, young panicles (1–5 cm), mature panicles (15–20 cm) of K1561 and G1025 with Trizol (Invitrogen, Carlsbad, CA, USA). RNAs were digested with DNase I (Promega, USA) to eliminate genomic DNA contamination before cDNA synthesis. For cDNA synthesis, 5.0 μg total RNA was used for reverse transcribed by M-MLV Reverse Transcriptase (Promega, USA) using oligo-d (T) according to the user manual. For RT-PCR analysis, 2 μl of the first cDNA strand was used with gene-specific primer pairs in 20 μL reaction volumes, which contained 1 μM forward and reverse primer each, 200 μM dNTPs, and 1 U Taq enzyme (Takara, China). A DNA Engine peltier ThermalCycler (Bio-RAD) was used to perform the reaction. The PCR procedure was as follows: 94 °C for 2 min; 35 cycles of 94 °C for 30 s, 53 °C for 30 s, 72 °C for 1 min per 1 kb; 5 min for a final elongation at 72 °C. The PCR products were purified from agarose gel after electrophoresis and cloned into pMD-18 T vector (Takara, China). Clones were sequenced by Tianyihuiyuan Company (Wuhan, China). Sequences were aligned by ClustalW (https://www.genome.jp/tools-bin/clustalw).
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Derived cleaved amplified polymorphic sequences
Grain width and weight 2
Grain weight of chromosome 6
Open reading fragments
Polymerase chain reaction
Quantitative trait locus
Recombinant inbred line
Reverse transcript (RT)-PCR
Specific locus amplified fragment
Single nucleotide polymorphism
Simple sequence repeats
Wide and thick grain 1
Khush G. Green revolution: the way forward. Nat Rev Genet. 2001;2:815–22.
Guo L, Ma L, Jiang H, Zeng D, Hu J, Wu L, et al. Genetic analysis and fine mapping of two genes for grain shape and weight in rice. J Integr Plant Biol. 2009;51:45–51.
Liu T, Shao D, Kovi M, Xing Y. Mapping and validation of quantitative trait loci for spikelets per panicle and 1,000-grain weight in rice (Oryza sativa L.). Theor Appl Genet. 2010;120:933–42.
Bai X, Luo L, Yan W, Kovi M, Xing Y. Quantitative trait loci for rice yield-related traits using recombinant inbred lines derived from two diverse cultivars. J Genet. 2011;90:209–15.
Tang S, Shao G, Wei X, Chen M, Sheng Z, Luo J, et al. QTL mapping of grain weight in rice and the validation of the QTL qTGW3.2. Gene. 2013;27:201–6.
Xu F, Sun X, Chen Y, Huang Y, Tong C, Bao J. Rapid identification of major QTLs associated with rice grain weight and their utilization. PLoS One. 2015;10(3):e0122206.
Zhang H, Fan Y, Zhu Y, Chen J, Yu S, Zhuang J. Dissection of the qTGW1.1 region into two tightly-linked minor QTLs having stable effects for grain weight in rice. BMC Genet. 2016;17:98.
Ashikari M, Sakakibara H, Lin S, Yamamoto T, Takashi T, Nishimura A, et al. Cytokinin oxidase regulates rice grain production. Science. 2005;309:741–5.
Fan C, Xing Y, Mao H, Lu T, Han B, Xu C, 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:1164–71.
Mao H, Sun S, Yao J, Wang C, Yu S, Xu C, et al. Linking differential domain functions of the GS3 protein to natural variation of grain size in rice. Proc Natl Acad Sci U S A. 2010;107(45):19579–84.
Huang X, Qian Q, Liu Z, Sun H, He S, Luo D, et al. Natural variation at the DEP1 locus enhances grain yield in rice. Nat Genet. 2009;41:494–7.
Song X, Huang W, Shi M, Zhu M, Lin H. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat Genet. 2007;39:623–30.
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:1023–8.
Li Y, Fan C, Xing Y, Jiang Y, Luo L, 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–9.
Zhang X, Wang J, Huang J, Lan H, Wang C, Yin C, et al. Rare allele of OsPPKL1 associated with grain length causes extra-large grain and a significant yield increase in rice. Proc Natl Acad Sci U S A. 2012;109(52):21534–9.
Wang S, Wu K, Yuan Q, Liu X, Liu Z, Lin X, et al. Control of grain size, shape and quality by OsSPL16 in rice. Nat Genet. 2012;44:950–4.
Ishimaru K, Hirotsu N, Madoka Y, Murakami N, Hara N, Onodera H, et al. Loss of function of the IAA-glucose hydrolase gene TGW6 enhances rice grain weight and increases yield. Nat Genet. 2013;45:707–11.
Song X, Kuroha T, Ayano M, Furuta T, Nagai K, Komeda N, et al. Rare allele of a previously unidentified histone H4 acetyltransferase enhances grain weight, yield, and plant biomass in rice. Proc Natl Acad Sci U S A. 2015;112:76–81.
Li T, Jiang J, Zhang S, Shu H, Wang Y, Lai J, et al. OsAGSW1, an ABC1-like kinase gene, is involved in the regulation of grain size and weight in rice. J Exp Bot. 2015;66:5691–701.
Sun B, Zhan X, Lin Z, Wu W, Yu P, Zhang Y, et al. Fine mapping and candidate gene analysis of qHD5, a novel major QTL with pleiotropism for yield-related traits in rice (Oryza sativa L.). Theor Appl Genet. 2017;130:247–58.
Ma X, Feng F, Zhang Y, Elesawi IE, Xu K, Li T, et al. A novel rice grain size gene OsSNB was identified by genome-wide association study in natural population. PLoS Genet. 2019;15(5):e1008191.
Zhou Y, Tao Y, Zhu J, Miao J, Liu J, Liu Y, et al. GNS4, a novel allele of DWARF11, regulates grain number and grain size in a high-yield rice variety. Rice. 2017;10:34.
Huang K, Wang D, Duan P, Zhang B, Xu R, Li N, et al. WIDE AND THICK GRAIN 1, which encodes an otubain-like protease with deubiquitination activity, influences grain size and shape in rice. Plant J. 2017;91:849–60.
Zhou W, Wang X, Zhou D, Ouyang Y, Yao J. Overexpression of the 16-kDa a-amylase/trypsin inhibitor RAG2 improves grain yield and quality of rice. Plant Biotechnol J. 2017;15:568–80.
Hu Z, Lu S, Wang M, He H, Sun L, Wang H, et al. A novel QTL qTGW3 encodes the GSK3/SHAGGY-like kinase OsGSK5/OsSK41 that interacts with OsARF4 to negatively regulate grain size and weight in rice. Mol Plant. 2018;11:736–49.
Guo S, Qin F, Zhang D, Lin X. Characterization of interspecific hybrids and backcross progenies from a cross between Oryza minuta and Oryza sativa. Sci China Ser C-Life Sci. 2009;52:1148–55.
Rahman M, Chu S, Choi M, Qiao Y, Jiang W, Piao R, et al. Identification of QTLs for some agronomic traits in rice using an introgression line from Oryza minuta. Mol Cells. 2007;24:16–26.
Guo S. Development and characterization of substation lines from a cross of Oryza sativa and Oryza minuta. Wuhan: Dissertation thesis, Huazhong Agricultural University; 2009.
International Rice Genome Sequencing Project. The map-based sequence of the rice genome. Nature. 2005;436:793–800.
Zhu Z, Li X, WeiY GS, Sha A. Identification of a novel QTL for panicle length from wild rice (Oryza minuta) by Specific Locus Amplified Fragment Sequencing and high density genetic mapping. Front Plant Sci. 2018;9:1492.
Romero G, Amante-Bordeos A, Dalmacio R, Elloran R, Sitch L. Comparative studies of isozymes in Oryza sativa, O. minuta, and their interspecific derivatives: evidence for homoeology and recombination. Theor Appl Genet. 1993;87:609–15.
Sasaki K. Utilization of transcription factors for controlling floral morphogenesis in horticultural plants. Breed Sci. 2018;68:88–98.
Li Z, Pinson S, Park W, Paterson A, Stansel J. Epistasis for three grain yield components in rice (Oryza sativa L.). Genetics. 1997;145:453–65.
Parenicova L, de Folter S, Kieffer M, Horner D, Favalli C, Busscher J, et al. Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis: new openings to the MADS world. Plant Cell. 2003;15:1538–51.
de Folter S, Busscher J, Colombo L, Losa A, Angenent G. Transcript profiling of transcription factor genes during silique development in Arabidopsis. Plant Mol Biol. 2004;56:351–66.
Arora R, Agarwal P, Ray S, Kumar Singh A, Singh V, Tyagi A, et al. MADS-box gene family in rice: genome-wide identification, organization and expression profiling during reproductive development and stress. BMC Genomics. 2007;8:242.
Yu J, Miao J, Zhang Z, Xiong H, Zhu X, Sun X, et al. Alternative splicing of OsLG3b controls grain length and yield in japonica rice. Plant Biotechnol J. 2018;16(9):1667–78.
Murray M, Thompson W. Rapid isolation of high molecular-weight plant DNA. Nucleic Acids Res. 1980;8:4321–5.
Cho Y, Blair M, Panaud O, McCouch S. Cloning and mapping of variety-specific rice genomic DNA sequences: amplified fragment length polymorphisms (AFLP) from silver-stained polyacrylamide gels. Genome. 1996;39:373–8.
Lincoln S, Daly M, Lander E. Constructing a genetic linkage map with MAPMAKER/EXP v. 3.0: a tutorial and reference manual. Whitehead Institute Technical Report, Cambridge, Mass; 1993.
Yang J, Zhu J. Methods for predicting superior genotypes under multiple environments based on QTL effects. Theor Appl Genet. 2005;110(7):1268–74.
Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, et al. SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS One. 2013;8:e58700.
Sha A, Qi Y, Shan Z, Chen H, Yang Z, Qiu D, et al. Identifying patellin-like genes in Glycine max and elucidating their response to phosphorus starvation. Acta Physiol Plant. 2016;38:138.
This work was supported by grants from the National Science and Technology Major Project of China (2016ZX08001001), National Natural Science Fund of China (No. 31860368, 31771735, 31960402), Guangxi Natural Science Foundation of China (No. 2017GXNSFDA198039), and Guangxi Key Laboratory of Rice Genetics and Breeding Open Foundation of China (No. 160–380–16-1). The funders had no role in the study design, data collection, and analysis, decision to publish or preparation of the manuscript.
Ethics approval and consent to participate
The research was conducted in accordance with the IUCN Policy Statement on Research Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, X., Wei, Y., Li, J. et al. Identification of QTL TGW12 responsible for grain weight in rice based on recombinant inbred line population crossed by wild rice (Oryza minuta) introgression line K1561 and indica rice G1025. BMC Genet 21, 10 (2020). https://doi.org/10.1186/s12863-020-0817-x
- Wild rice
- Introgression lines
- Transcription factor