Molecular Markers and Candidate Genes for Thermo-Sensitive Genic Male Sterile in Rice

Sudthana Khlaimongkhon, Sriprapai Chakhonkaen, Keasinee Pitngam, Khanittha Ditthab, Numphet Sangarwut,Natjaree Panyawut, Thiwawan Wasinanon, Chareerat Mongkolsiriwatana, Julapark Chunwongse, Amorntip Muangprom,

?

Molecular Markers and Candidate Genes for Thermo-Sensitive Genic Male Sterile in Rice

Sudthana Khlaimongkhon1, Sriprapai Chakhonkaen2, Keasinee Pitngam2, Khanittha Ditthab2, Numphet Sangarwut2,Natjaree Panyawut2, Thiwawan Wasinanon2, Chareerat Mongkolsiriwatana3, Julapark Chunwongse1, Amorntip Muangprom1, 2

(Center for Agricultural Biotechnology, Kasetsart University, Nakhon Pathom 73140, Thailand; National Center for Genetic Engineering and Biotechnology, National Science and Technology Development Agency, Pathum Thani 12120, Thailand; Division of Genetics, )

The discovery of thermo-sensitive genic male sterility (TGMS) has led to development of a simple and highly efficient two-line breeding system. In this study, genetic analysis was conducted using three F2populations derived from crosses between IR68301S, anTGMS rice line, and IR14632 (), Supanburi 91062 () and IR67966-188-2-2-1 (), respectively. Approximately 1:3 ratio between sterile and normal pollen of F2plants from the three populations revealed that TGMS is controlled by a single recessive gene. Bulked segregant analysis using simple sequence repeat (SSR) and insertion-deletion (InDel) markers were used to identify markers linked to thegene. The linkage analysis based on the three populations indicated that thelocus was located on chromosome 2 covering the same area. Using IR68301S × IR14632 F2population, the results showed that thelocus was located between SSR marker RM12676 and InDel marker 2gAP0050058. The genetic distance from thegene to these two flanking markers were 1.10 and 0.82 cM, respectively. InDel marker 2gAP004045 located between these two markers showed complete co-segregation with the TGMS phenotype. In addition, InDel marker vf0206114052 showed 2.94 cM linked to thegene using F2populations of IR68301S × Supanburi 91062. These markers are useful tool for developing new TGMS lines by marker-assisted selection. There were ten genes located between the two flanking markers RM12676 and 2gAP0050058. Using quantitative real-time PCR for expression analysis, 7 of the 10 genes showed expression in panicles, and response to temperatures. These genes could be the candidate gene controlling TGMS in IR68301S.

hybrid rice; thermo-sensitive genic male sterility; insertion/deletion; simple sequence repeat; marker-assisted selection

Rice (L.) is one of the most important global food crops.The world’s population has doubled since the early 1960s, therefore, existing yields of the major cereal crops will be insufficient to meet the food needs of the future (FAO, 2018). In addition, arable lands have declined worldwide, thus an increase in yield will be necessary to meet this demand (Bruinsma, 2003). Rice production in Thailand signifies a portion of the Thai economy and is a major export product.However,yield of Thai rice is less than that of the neighboring countries such as Indonesia, Vietnam, Malaysia, Laos and the Philippines (USDA, 2016). Yields have increased with the development of hybrid rice currently planted in many countries. The trend of hybrid rice development in both commercialization and research of private companies has been increasing. Hybrid rice technology has increased yields by 20%–30% under unchanged irrigation conditions (Yuan, 1998; Virmani, 2003). Hybrid rice seeds of large scale commercial production have been released to the market in China since 1976. The Ministry of Agriculture of China develops super rice breeding programs to increase rice yield per planting area, resulting thataverage rice yieldsin China rose from 1.89 t/hm2in 1949 to 6.71 t/hm2in 2012 (Cao and Zhan, 2014). In Thailand, average rice yield is only 2.89 t/hm2in 2018 (Thai Rice ExporterAssociation, 2018). Many countries which adopt hybrid rice technology are able to achieve higher yields.

Hybrid rice seed production involves the use of male sterility systems. Two well established male sterility systems in rice are cytoplasmic genetic male sterility (CMS), a three-line system, and environmentally sensitive genic male sterility (EGMS), a two-line system including photoperiod-sensitive genic male sterility (PGMS) and thermo-sensitive genic male sterility (TGMS). The discovery of PGMS and TGMS has led to development of a simple and highly efficient two-line breeding system in hybrid rice seed production. Several EGMS-related genes have been mapped such as(Wang et al, 1995),(Yamaguchi et al, 1997; Pitnjam et al, 2008),(Subudhi et al, 1997;Lang et al, 1999),(Dong et al, 2000),(Wang et al, 2003;Nas et al, 2005; Jiang et al, 2006; Yang et al, 2007),(Lee et al, 2005),(t) (Li et al, 2005),(Hussain et al, 2011),(Sheng et al, 2013),(Qi et al, 2014),(Peng et al, 2010),(Liu et al, 2001; Zhou et al, 2011),(Zhanget al,1994),(Dinget al,2012a, b),(Xu et al, 2011),(Zhou et al, 2012),(Peng et al, 2008) and(Jia et al, 2001). Candidate genes for EGMS have been reported, Myb-like DNA-binding domain containing protein (Zhou et al, 2011), a nuclear ribonuclease Z gene (Xu et al, 2011; Zhou et al, 2014), ORMDL (Pitnjam et al, 2008; Chueasiri et al, 2014) and the putative pollen specific protein encoded by(Nas et al, 2005). In addition, several genes controlling male sterility phenotype have been reported in rice such as undeveloped() (Jung et al, 2005), tapetum degeneration retardation () (Li et al, 2006), CER-like protein causing wax-deficient() (Jung et al, 2006), and LRR receptor kinase affecting(Nonomura et al, 2003).

IR68301S is a TGMS rice line obtained from International Rice Research Institute (IRRI), the Philippines. Male sterility in IR68301S is stable in a high temperature (≥ 28oC), and it has been used for the development of two-line hybrid in Thailand. In this study, we mapped a gene controlling TGMS in IR68301S, and identified candidate genes using expression analysis. The knowledge in this study will be useful for rice breeding programs.

Materials and methods

Rice materials and phenotypic screening

An F2population, generated from a cross between IR68301S (B2), anTGMS mutant rice line, and IR14632 (B30), arice line, was used for identification of linked markers to thegene. These plants were obtained from IRRI. In addition, the other F2populations generated from crosses between B2 and Supanburi 91062 (B11), a Thairice line having high yield, and between B2 and IR67966-188-2-2-1(KM9), arice line, were also used. Theandcrosses were used to increase chances for polymorphic marker identification, and then these markers were used for linkage analysis. The F2populations were planted in paddy fields for observation on pollen sterility/fertility in summer of 2011, 2012 and 2013, in which the maximum temperatures were 33 oC–37 oC. Pollen sterility and fertility of the F2plants were evaluated on flowering day by observing anthers with the naked eyes and under a microscope after staining with 1% I2-KI solution. The χ2test was applied to determine the F2segregation ratio of sterility and fertility.

Genotyping and linkage analysis

Genomic DNA was extracted from fresh leaves of individual F2sterile plants generated from all the three crosses using the cethy trimethylammonium bromide (CTAB) method (Murray and Thomson, 1980). Polymorphic markers between parents and bulked segregant analysis (BSA) using pooled DNA from 10 individuals each for fertile and sterile plants, were applied to identify linked markers togene. The resulting linked markers were used to determine the genotype of the individual F2plants. Primers for single nucleotide polymorphism (SNP), and insertion-deletion (InDel) markers were designed from sequences obtained from GRAMENE database (http://www.gramene.org) andRiceVarMap (http://ricevarmap.ncpgr.cn/).Amplification of DNA fragments using simple sequence repeat (SSR) and InDels was performed as previously described by McCouch et al (2002) and Pitnjam et al (2008). The resulting PCR products were electrophoretically separated on 3% agarose gel or 6% acrylamide gel, and DNA patterns were observed by ethidium bromide or silver staining.The recombination frequency () was calculated with the formula:

= (1+2)/

whereis the total number of sterile plants surveyed,1is the number of sterile individuals with the homozygous band of the fertile parent, and2is the number of individuals with heterozygous bands (Peng et al, 2010). Recombination frequency was converted into genetic distance (cM) for linkage analysis.

Expression analysis

Segregating wild type male fertile and TGMS male sterile BC3F5plants generated from Pathum Thani 1, a Thai elite line, crossed with B2 were used for expression analysis, These plants were grown under natural condition at National Center for Genetic Engineering and Biotechnology, National Science and Technology Development Agency (14o04′50.6′′N and 100o36′09.0′′E), Pathum Thani Province, Thailand. Pre-germinated grains wereseeded in a tray and then the seedlings were transplanted into plastic pots (without hole at the bottom, diameter of 20.5 cm and height of 16.0 cm) with one seedling per pot. About one month before flowering, each of at least five plants/line were moved to controlled growth rooms at temperatures(24 ± 2) oC and (32 ± 2) oC under 80% relative humidity, 12 h light/12 h darkness until flowering. For expression analysis, young panicles about 13?14 cm in length (dyad stage) were harvested from BC3F5wild type and TGMS mutant plants at the same day. Total RNA was extracted using TRIZOLTMreagent (Invitrogen, USA). Quantitative real-time PCR (qRT-PCR) experiments were performed with cDNA synthesis kit (Fermentus, Lithuania) using cDNAs transcribed from total RNA. Based on information in GRAMENE database, genes located between two flanking markers tightly linked togene were used for expression analysis by qRT-PCR.The qRT-PCR was performed using Bio-RAD iCycler iQ5 Machine (BioRAD, USA) and all reactions were conducted in 96-well plates (BioRad, USA). The qPCRBIO SyGreen Mix Lo-ROX (PCR Biosystems) assay was used in a total volume of 10 μL per reaction. Each reaction mixture contained 1 μL cDNA, 2μL distilled water, 1 μL (10 μmol/L) forward primer, 1 μL (10 μmol/L) reverse primer and 5 μL of 2× qPCRBIO SyGreen Mix. All qRT-PCR plates were carried out with following cycling condition: 95 oC for 3 min, following 35 cycles of 95 oC for 30 s, 57 oC to 63 oC for 30 s and 72 oC for 30 s, then 95 oC for 1 min and a melting curve from 60 oC to 95 oC in 0.5 oC increments.was used as an internal control. Three replications were performed for all reactions. The 2?ΔΔCtmethod was used to calculate relative gene expression(Livak and Schmittgen, 2001). The significantly of qRT-PCR were calculated byprogram (http://cran.r-project.org/)with ANOVA and Duncan’s new multiple range test.

Results

Pollen staining

Mature anthers of the TGMS mutant and wild type grown in controlled growth rooms under fertile and sterile conditions were observed. There was no pollen produced in anthers of TGMS mutant plants under sterile condition (32 ± 2) oC (Fig. 1-A), while underthe fertile condition of (24 ± 2) oC, pollen production was similar to that of wild type plants (Fig. 1-B). At these two conditions, wild type plants produced normal pollen (Fig. 1-C and -D).

Fig. 1. Pollen staining of thermo-sensitive genic male sterility (TGMS) mutant and wild type.

A, Abortive anther of the TGMS mutant grown under high temperature (32 ± 2)oC. B, Normal anther of the TGMS mutant grown under low temperature (24 ± 2) oC. C, Normal anther of wild type grown under high temperature (32 ± 2) oC. D, Normal anther of wild type grown under low temperature (24 ± 2) oC.

Genetic analysis of tms gene controlling TGMS in IR68301S

By growing the F2populations under sterile condition of(32 ± 2) oC, the segregation patterns of fertile to sterile plants in the F2populations were near to 3:1 ratio in all the three populations (Table 1). The results indicated that the TGMS of IR68301S was controlled by a single recessive gene.

Mapping of tms gene controlling TGMS in IR68301S

A total of 119 SSR, 10 SNP and 283 InDel markers were used to identify polymorphism between parents from the three crosses. The results showed that 53, 3 and 7 markers were polymorphic between the parents from B2 × B30, B2 × B11 and B2 × KM9, respectively.

The resulting polymorphic markers were applied to identify linked markers to thegene using BSA. And 24 out of the 53 markers showed linkage to thegene in F2male sterile plants from B2 × B30. In addition, three and four markers showed linkage to thegene in F2male sterile plants from B2 × B11 and B2 × KM9, respectively (Supplemental Table 1). According to GRAMENE and RiceVarMap databases, these markers are on chromosome 2 at the position of 4.19 to 7.44 Mb.

Table 1. Segregation patterns of F2 populations for fertility and sterility.

B2, IR68301S; B30, IR14632; B11, Supanburi 91062; KM9, IR67966-188-2-2-1.

The value of χ20.05is 3.84.

The linked markers identified from each cross were used to genotype the F2male sterile. These F2male sterile plants selected from flowering plants in summers (33 oC–37 oC) were completely sterile. A total of 217, 259 and 450 F2male sterile plants from B2 × B30, B2 × B11 and B2 × KM9 were used for genotyping, respectively. Segregation patterns were shown in Fig. 2. A linkage analysis based on the three F2populations showed that thegene in IR68301S was located on chromosome 2 (Fig. 3). Using F2sterile individuals of B2 × B30, the results showed that thegene was located between SSR marker RM12676 at 5.74 Mb and InDel marker 2gAP0050058 at 5.81 Mb. The genetic distance of thegene from the two markers was 1.10 and 0.82 cM, respectively. In addition, InDel marker 2gAP004045 showed complete co-segregation with the TGMS phenotype (Fig. 3-A). Using F2sterile individuals of B2 × B11, the results showed thatgene was located between RM126016 at 4.69 Mb and InDel marker vf0206114052 at 6.11 Mb. The genetic distance of thegene from the two markers was 14.97 and 2.94 cM, respectively (Fig. 3-B). Using F2population of B2 × KM9, the results showed that four SSR markers, RM126016, RM12649, RM6378 and RM12674 were linked to thegene. The genetic distance of thegene from these markers were 12.50, 8.75, 8.25 and 6.25 cM, respectively (Fig. 3-C).

Expression analysis of candidate genes located inside the flanking markers

Based on information in GRAMENE database, there are ten genes (Table 2) located between the two flanking markers, RM12676 and 2gAP0050058. Nine out of the ten genes were annotated as expressed proteins, and one was annotated as a conserved hypothetical protein. These genes were tested for expression in panicles of wild type and TGMS mutantplants under fertile and sterile conditions. The results showed that the expressions of seven genes were detected in the tested tissues, but those of,andcould not be detected (data not shown). All the seven expressed genes showed different expression levels under fertile and sterile conditions or between the TGMS mutant and the wild type. Genes,andshowed significantly different expression levelsbetween the two conditions in the TGMS rice plants but no significantly different expression levels were detected between the two conditions in the wild-type rice plants (Fig. 4).andshowed higher expression levelsunder fertile conditionthan sterile condition, whileshowed higher expression levelsunder sterile conditionthanthe fertile conditionsThese three genes showed higher expression levels in the TGMS mutants than in the wild type under both conditions. These two genesandshowed a similar pattern of expression both in the mutant and wild-type rice plants, by showing higher expression under lower temperature condition. These genesandshowed a similar pattern of expression by showing higher expression under lower temperature condition in the TGMS mutants, but higher expression under higher temperature condition in wild-type plants (Fig. 4).

Fig. 2. Samples of genotyping F2male sterile plants.

P1, Female parent (IR68301S, B2); P2, Male parent (IR14632, B30); 1–10, F2male sterile individuals of B2 × B30 (×) using Os02g12370 marker; P3, Male parent (Supanburi 91062, B11); 11–20, F2male sterile individuals of B2 × B11 (×) using vf0206114052 marker.

Fig. 3. Genetic linkage map of thermo-sensitive genic male sterility gene on chromosome 2.

The linkage maps were analyzed based on male sterile F2populations of IR68301S × IR14632 (A), IR68301S × Supanburi 91062 (B) and IR68301S × IR67966-188-2-2-1 (C). Distances of each marker in centiMorgans (cM) fromgene were given on the left side of the genetic map.

Table 2. Genes located between the two flanking markers RM12676 and 2gAP0050058.

Fig. 4. Relative expression levels of genes in panicles by quantitative real-time PCR.

HT, High temperature condition (32 oC); LT, Low temperature condition (24 oC); TGMS, Thermo-sensitive genic male sterility.

Data represent Mean ± SE, and the same lowercase letter(s) indicate no significant difference at< 0.05 by Duncan’s new multiple range test.

Discussion

The segregation patterns of fertile to sterile plants in the three F2populations followed the 3:1 ratio, which is typical of Mendelian low, indicating that the TGMS of anIR68301S is controlled by a single recessive gene, which is in agreement with other studies (Borkakati and Virmani, 1996; Lopez et al, 2003; Wang et al, 2003; Hussain et al, 2011; Qi et al, 2014).

Using the three mapping populations, thegene in IR68301S was mapped on chromosome 2. Using B2 × B30 F2population, thegene was located between SSR marker RM12676 and InDel marker 2gAP0050058. In addition, InDel marker 2gAP004045, located between these two markers, showed complete co-segregation with the TGMS phenotype. These markers are useful tool for developing new TGMS lines by marker-assisted selection (MAS) and identifying the TGMS individuals at earlier stages of line development in rice breeding program. Using F2population fromparents (B2 and B11),only two markers were linked to thegene. The genetic distance of thegene from the two markers was 14.97 and 7.50 cM. Since major subspecies of rice in Thailand is, the polymorphic and tightly linked markers generated fromparents will be more practical than the linked markers generated from×crosses. Therefore, several markers nearby the flanking region were designed and test for polymorphism between the twoparents (B2 and B11). The 30 InDel markers covering flanking region from 5.50 to 6.28 Mb were designed to test polymorphism between these two parents (Supplemental Table 2). The results showed that only vf0206114052 at 6.11 Mb showed polymorphism and linked (2.94 cM) to thegene. Therefore, this marker will be useful inbreeding programs. The two tightly linked InDel markers, vf0206114052 and 2gAP004045, will be useful for×and×rice breeding programs, respectively.

Candidate genes for EGMS have been reported. The locus of(photo period-sensitive male sterility) and(photo- or thermo-sensitive genic male sterility) represent the same locus on chromosome 12 conferring PGMS and TGMS traits (Zhou et al, 2012). This report revealed that a SNP C-to-G between the fertile lines and sterile lines led to increasing of methylation in putative promoter region of this non-coding gene, resulting in premature programmed cell death (PCD) in developing anthers, thus causing PGMS inand TGMS insubspecies.

Theis a candidate gene ofwhich locates on chromosome 2 at the position of 6.39 Mb (htpp://www.gramene.org) and this gene encodes RNase ZS1. Experimental results revealed that at high temperature RNase ZS1loses its function and causes defective pollen production (Zhou et al, 2014). This enzyme maintained mRNAs ofUbat normal level and led to male fertility.LOC_Os02g12290, a nuclear ribonuclease Z gene was identified as the candidate for thegene, located on chromosome 2 in Guangzhan 63S, aPTGMS line (Xu et al, 2011).

The gene controlling TGMS in IR68301S was located between the two flanking markers at the position of 5.74 to 5.81 Mb. There were no previously reported EGMS genes at this position. There were 10 genes located between the flanking markers, however, we detected the expressions level of only 7 genes in the TGMS mutant and the wild type under both conditions. All these genes responded to temperatures by showing different expression levels under fertile and sterile conditions in the TGMS rice plants.Threegenes,andshowed similar expression levelsunder fertile and sterile conditions in wild type rice plants but they showed significantly different levels of expression in the mutant plants under sterile condition. Interestingly,showed significantly higher expression level under sterile condition than under fertile condition in the TGMS rice plants, while it showed similar expression levelsunder both conditions in the wildtype rice plants. Theencodesa tetratricopeptide-like helical domain protein (TPR). Shin et al (2014) reported that the N-terminal tetratricopeptide repeat1 (TPR1) domain is essential for its interaction with pectate lyase-like proteins (PLLs) in petunia, maize and. Lacking the TPR1 domain, no interaction resulting in pollen that failed to germinate.pollen calmodulin-binding protein (OsPCBP) is a Ca2+-dependent calmodulin-binding (CaMBP) protein and it contains six TPR motifs that can interact with protein in pollen to regulate germination and also involve in starch accumulation. The absence or reduction of starch can cause abortive pollen (Zhang et al, 2012). This information indicated that the TPR domain is associated with pollen development.encodes a zinc ion binding protein andencodes a zinc finger protein. The zinc finger family proteins are transcription factors (TF) involved in abiotic stress. Bai et al (2015) reported that three zinc finger TF genes (LOC_Os06g14180, LOC_Os12g37800, LOC_Os11g14000) are up-regulated and two zinc finger TF genes (LOC_Os01g50750 and LOC_Os05g44550) are down-regulated following cold stress treatment in rice male sterile lines. In, the result of RNA sequencing showed that, a zinc finger TF protein gene associated with flowering time and tolerance to abiotic stress (Yang et al, 2014). Yang et al (2015) reported that zinc finger CCCH domain-containing proteins are two (Unigene43080_Zhong-531and Unigene43085_Zhong-531) of differential expression genes (DEGs) which respond to temperature interacting with nitrogen at meiosis stage of rice spikelet. In our study,andshowed higher expression levels under lower temperature in the TGMS plants, while in wild type rice plants,showed lower expression levelsunder lower temperature andshowed similar expression levels under both conditions.encodes an F-box domain and kelch repeat containing protein.F-box proteins are critical for degradation of cellular proteins. The F-box protein-encoding genes have specific and/or overlapping expression during floral transition as well as panicle and seed development. The F-box protein encoding genes are involved in different abiotic stress conditions (Jain et al, 2007; Chunthong et al, 2017). The TPR domains were predicted in rice F-box proteins responsible for processing and/or translation of mRNAs. DNA-binding domains such as zinc finger, ring finger and helix loop helix are also found in rice F-box proteins, which may be directly or indirectly involved in transcriptional regulation (Small and Peeters, 2000).

encodes a conserved hypothetical protein. It showed higher level of expression under the higher temperature condition than under the lower temperature condition in the wildtype rice plants, but in the TGMS rice plants the expression of this gene was similar under both conditions.encodes a putative uncharacterized protein. It showed similar expression patterns in the TGMS and wild type rice plants but the expression levels were higher in the TGMS than in wild type rice plants under both conditions. Although the functions of,andare not clear but their expressions in our study suggested that they could be involved in TGMS. Therefore, further study is needed to identify the candidate gene controlling TGMS in IR68301S rice line.

Acknowledgements

This study was supported by Center for Agricultural Biotechnology, Kasetsart University, Center of Excellence on Agricultural Biotechnology (AG-BIO/PERDO- CHE), Agricultural Research Development Agency (ARDA), and National Science and Technology Development Agency, Thailand. The authors thank Dr. Peera Jaruampornpan for her critical reading and comments on the manuscript.

supplemental datA

The following materials are available in the online version of this article at http://www.sciencedirect.com/science/ journal/16726308; http://www.ricescience.org.

Supplemental Table 1. Markers on chromosome 2 linked to thegene using F2male sterile plants from IR68301S × IR14632, IR68301S × Supanburi 91062 and IR68301S × IR67966-188-2-2-1.

Supplemental Table 2. InDel markers on chromosome 2 covering flanking region from 5.50 to 6.28 Mb.

Bai B, Wu J, Sheng WT, Zhou B, Zhou LJ, Zhuang W, Yao DP,Deng QY. 2015. Comparative analysis of anther transcriptome profiles of two different rice male sterile lines genotypes under cold stress., 16(5):11398–11416.

Borkakati R R, Virmani SS. 1996. Genetics of thermosensitive genic male sterility in rice.,88(1): 1–7.

Bruinsma J. 2003. World Agriculture: Towards 2015/2030. London, UK: Earthscan.

Cao L Y, Zhan X D. 2014. Chinese experiences in breeding three-line, two-line and super hybrid rice.: Yan W G, Bao J S. Rice: Germplasm, Genetics and Improvement. Intech: 279–308.

Chueasiri C, Chunthong K, Pitnjam K, Chakhonkaen S,Sangarwut N,Sangsawang K,Suksangpanomrung M,Michaelson LV,Napier JA,Muangprom A. 2014. Rice ORMDL controls sphingolipid homeostasis affecting fertility resulting from abnormal pollen development.,9(9): e106386.

Chunthong K, Pitnjam K, Chakhonkaen S, Sangarwut N, Panyawut N, Wasinanon T, Ukoskit K, Muangprom A. 2017. Differential drought responses in F-box gene expression and grain yield between two rice groups with contrasting drought tolerance.,36:970–982.

Ding J H, Lu Q, Ouyang Y D, Mao H L,Zhang P B,Yao J L,Xu C G,Li X H,Xiao J H,Zhang Q F. 2012a. A long noncoding RNA regulates photoperiod-sensitive male sterility, an essential component of hybrid rice.,109(7):2654–2659.

Ding J H, Shen J Q, Mao H L, Xie W B, Li X H,Zhang Q F. 2012b.RNA-directed DNA methylation is involved in regulating photoperiod-sensitive male sterility in rice.,5(6):1210–1216.

Dong NV, Subudhi PK, Luong PN, Quang VD, Quy TD, Zheng HG, Wang B, Nguyen HT. 2000. Molecular mapping of a rice gene conditioning thermo-sensitive genic male sterility using AFLP, RFLP and SSR techniques.,100(5):727–734.

FAO. 2018.FAO Statistical Pocketbook 2018: World Food and Agriculture. Rome: FAO.

Hussain A J, Ali J, Siddiq EA, Gupta VS, Reddy UK, Ranjekar PK. 2011. Mapping ofgene for temperature-sensitive genic male sterility (TGMS) in rice (L.).,131(1):42–47.

Jain M, Nijhawan A, Arora R, Agarwal P, Ray S, Sharma P, Kapoor S, Tyagi AK, Khurana JP. 2007. F-Box proteins in rice: Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress.,143(4):1467–1483.

Jia JH, Zhang DS, Li CY, Qu XP, Wang SW, Chamarerk V, Nguyen HT, Wang B Y. 2001. Molecular mapping of the reverse thermosensitive genic male-sterile gene () in rice.,103(4):607–612.

Jiang D G, Lu S, Zhou H, Wu X J, Zhuang C X, Liu Y G, Mei M T. 2006. Mapping of the rice (L.) thermo-sensitive genic male sterile genewith EST and SSR markers.,51(4):417–420.

Jung KH, Han MJ, Lee YS, Kim YW,Hwang I,Kim MJ,Kim YK,Nahm BH,An G. 2005. Rice undevelopedis a major regulator of early tapetum development.,17(10):2705–2722.

Jung KH, Han MJ, Lee DY, Lee YS, Schreiber L, Franke R, Faust A, Yephremov A, Saedler H, Kim YW, Hwang I, An G.2006. Wax-deficientis involved in cuticle and wax production in rice anther walls and is required for pollen development.,18(11):3015–3032.

Lang NT, Subudhi PK, Virmani SS, Brar DS, Khush GS, Li Z, Huang N. 1999. Development of PCR-based markers for thermosensitive genetic male sterility gene(t) in rice (L.).,131(2):121–127.

Lee DS, Chen LJ, Suh HS. 2005. Genetic characterization and fine mapping of a novel thermo-sensitive genic male-sterile genein rice (L.).,111(7):1271–1277.

Li R B, Pandey MP, Sharma P, Ballabh G. 2005. Inheritance of thermosensitive genic male sterility in rice (L.).,88(11):1809–1815.

Li N, Zhang DS, Liu HS, Yin CS,Li XX,Liang WQ,Yuan Z,Xu B,Chu HW,Wang J,Wen TQ,Huang H,Luo D,Ma H,Zhang DB. 2006. The rice tapetum degeneration retardation gene is required for tapetum degradation and anther development.,18(11):2999–3014.

Liu N, Shan Y, Wang F P, Xu CG, Peng KM, Li XH, Zhang QF. 2001. Identification of an 85-kb DNA fragment containing, a locus for photoperiod-sensitive genic male sterility in rice.,266(2):271–275.

Livak KJ, Schmittgen TD. 2001. Analysis of relative gene expression data using real-time quantitative PCR andthe 2 (-Delta Delta C(T)) method.,25(4):402–408.

Lopez MT, Toojinda T, Vanavichit A, Tragoonrung S. 2003. Microsatellite markers flanking thegene facilitated tropical TGMS rice line development.,43(6):2267–2271.

Matthayatthaworn W, Sripichitt P, Phumichai C, Rungmekarat S, Uckarach S, Sreewongchai T. 2011. Development of specific simple sequence repeat (SSR) markers for non-pollen type thermo-sensitive genic male sterile gene in rice (L.).,10:16437–16442.

McCouch SR,Teytelman L, Xu Y, Lobos KB,Clare K,Walton M,Fu B,Maghirang R,Li Z,Xing Y,Zhang Q, Kono I,Yano M,Fjellstrom R,DeClerck G,Schneider D,Cartinhour S,Ware D,Stein L. 2002. Development and mapping of 2240 new SSR markers for rice (L.).,9(6):199–207.

Murray MG, Thompson WF.1980. Rapid isolation of high molecular weight plant DNA.,8(19):4321–4325.

Nas TMS, Sanchez DL, Diaz GQ, Mendioro MS, Vermani SS. 2005. Pyramiding of thermosensitive genetic male sterility(TGMS) genes and identification of a candidategene in rice.,145:67–75.

Nonomura K I, Miyoshi K, Eiguchi M, Suzuki T, Miyao A, Hirochika H, Kurata N. 2003. Thegene is necessary to restrict the number of cells entering into male and female sporogenesis and to initiate anther wall formation in rice.,15(8):1728–1739.

Peng HF, Zhang ZF, Wu B, Chen XH, Zhang GQ, Zhang ZM, Wan BH, Lu YP. 2008. Molecular mapping of two reverse photoperiod-sensitive genic male sterility genes (and) in rice (L.).,118(1):77–83.

Peng HF, Chen XH, Lu YP, Peng YF,Wan BH,Chen ND,Wu B,Xin SP,Zhang GQ. 2010. Fine mapping of a gene for non-pollen type thermosensitive genic male sterility in rice (L.).,120(5):1013–1020.

Pitnjam K, Chakhonkaen S, Toojinda T, Muangprom A. 2008. Identification of a deletion inand development of gene-based markers for selection.,228(5):813–822.

Qi Y B, Liu Q L, Zhang L, Mao B Z, Yan D W, Jin Q S, He Z H. 2014. Fine mapping and candidate gene analysis of the novel thermo sensitive genic male sterility-gene in rice., 127(5): 1173–1182.

Sharma M, Pandey GK. 2016. Expansion and function of repeat domain proteins during stress and development in plants.,6:1218.

Sheng Z H, Wei X J, Shao G N, Chen M L, Song J, Tang S Q, Luo J, Hu Y C, Hu P S, Chen L Y. 2013. Genetic analysis and fine mapping of, a novel thermosensitive genic male-sterile gene in rice (L.).,132(2):159–164.

Shin SB, Golovkin M, Reddy ASN. 2014. A pollen-specific calmodulin-binding protein, NPG1, interacts with putative pectate lyases.,4:5263.

Small ID, Peeters N. 2000. The PPR motif:A TPR-related motif prevalent in plant organellar proteins.,25(2):45–47.

Subudhi PK, Borkakati R P, Virmani SS, Huang N. 1997. Molecular mapping of a thermo-sensitive genetic male-sterility gene in rice using bulked segregant analysis.,40(2):188–194.

Thai Rice Exporters Association. 2018. World Rice Production and Ending Stocks. http://www.thairiceexporters.or.th/default_eng.htm.

USDA. 2016.World Agricultural Production. http://usda.mannlib.cornell.edu/usda/fas/worldag-production//2010s/2016/worldag-production-11-09-2016.pdf.

Virmani SS. 2003. Advances in hybrid rice research and development in the tropics.: Virmani SS, Mao CX,Hardy B. Hybrid Rice for Food Security, Poverty Alleviation, and Environmental Protection. Proceedings of the 4th International Symposium on Hybrid Rice, Hanoi, Vietnam. Los Ba?os,the Philippines: International Rice Research Institute: 7–20.

Wang B Y, Xu WW, Wang JZ, Wu W, Zheng HG, Yang ZY, Ray JD, Nguyen HT. 1995. Tagging and mapping the thermo-sensitive genic male-sterile gene in rice () with molecular markers.,91:1111–1114.

Wang YG, Xing QH, Deng QY, Liang FS, Yuan LP, Weng ML, Wang B. 2003. Fine mapping of the rice thermo-sensitive genic male-sterile gene.,107(5):917–921.

Xu J J, Wang B H, Wu Y H, Du P N, Wang J, Wang M, Yi C D, Gu M H, Liang G H. 2011. Fine mapping and candidate gene analysis of, the photoperiod-thermo-sensitive genic male sterile gene in rice (L.).,122(2):365–372.

Yamaguchi Y, Ikeda R, Hirasawa H, Minami M, Ujihara A. 1997. Linkage analysis of thermo-sensitive genic male sterility genein rice (L.).,47:371–373.

Yang QK, Liang CY, Zhuang W, Li J, Deng HB, Deng QY, Wang B. 2007. Characterization and identification of the candidate gene of rice thermo-sensitive genic male sterile geneby mapping.,225(2):321–330.

Yang Y J, Ma C, Xu Y J, Wei Q, Imtiaz M, Lan H B, Gao S, Cheng L N, Wang M Y, Fei Z J, Hong B, Gao J P. 2014. A zinc finger protein regulates flowering time and abiotic stress tolerance in Chrysanthemum by modulating gibberellin biosynthesis.,26(5):2038–2054.

Yang J, Chen X R, Zhu C L, Peng X S, He X P, Fu J R, Ouyang L J, Bian J M, Hu L F, Sun X T, Xu J, He H H. 2015. RNA-Seq reveals differentially expressed genes of rice () spikelet in response to interacting with nitrogen at meiosis stage.,16:959.

Yuan LP. 1998. Hybrid rice breeding in China.: Virmani SS, Siddiq EA, Muralidharan K. Advances in Hybrid Rice Technology. Proceedings of the 3rd International Symposium on Hybrid Rice. Hyderabad, India. Manila, the Philippines: International Rice Research Institute: 27–33.

Zhang Q F, Shen BZ, Dai XK, Mei MH, Saghai Maroof MA, Li ZB. 1994. Using bulked extremes and recessive classes to map genes for photoperiod-sensitive genic male sterility in rice.,91(18):8675–8679.

Zhang Q S, Li Z, Yang J, Li S Q, Yang D C, Zhu Y G. 2012. A calmodulin-binding protein from rice is essential to pollen development.,55(1):8–14.

Zhou YF, Zhang XY, Xue QZ. 2011. Fine mapping and candidate gene prediction of the photoperiod and thermo-sensitive genic male sterile gene(t)in rice.,12(6):436–447.

Zhou H, Liu Q J, Li J, Jiang D G, Zhou L Y, Wu P, Lu S, Li F, Zhu L Y, Liu Z L, Chen L T, Liu Y G, Zhuang C X. 2012. Photoperiod- and thermo-sensitive genic male sterility in rice are caused by a point mutation in a novel noncoding RNA that produces a small RNA., 22(4):649–660.

Zhou H, Zhou M, Yang Y Z,Li J, Zhu L Y, Jiang D G, Dong J F, Liu Q J, Gu L F, Zhou L Y, Feng M J, Qin P, Hu X C, Song C L, Shi J F, Song X W, Ni E D, Wu X J, Deng Q Y, Liu Z L, Chen M S, Liu Y G, Cao X F, Zhuang C X. 2014. RNase ZS1processes UbL40 mRNAs and controls thermosensitive genic male sterility in rice.,5:4884.

Copyright ? 2019, China National Rice Research Institute. Hosting by Elsevier B V

This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/)

Peer review under responsibility of China National Rice Research Institute

http://dx.doi.org/10.1016/j.rsci.2018.08.006

9 July 2018;

24 August 2018

Amorntip Muangprom (Amorntip.mua@biotec.or.th)

(Managing Editor: Wang Caihong)