Assessment of genetic variation of 15 Thai elite rice cultivars using InDel markers

: Many elite rice varieties have been developed in the last decades by Thai rice breeders. No information about the genetic gene pool of these elite rice varieties is available. This study evaluated the genetic relationship of 15 elite Thai rice cultivars using Indel markers. Of the 133 markers screened, 44 were polymorphic. The PIC scores ranged from 0.06 to-0.50. By the UPGMA clustering method, they were separated into three major groups. Several markers identified specific DNA bands to certain rice cultivars, e


INTRODUCTION
Rice (Oryza sativa), one of the most important cereal crops in the world, is a staple food and the main source of nutrition of Thai citizens.Thailand has proposed to expand the area of rice production, by adding 500,000 hectares to the currently 9.2 million hectares of rice-growing areas.From 2014 to 2017, Thailand was the second largest rice exporter after India on the world market (Thai Rice Exporters Association 2018).Although Thailand is one of the largest rice producers and exporters, the grain yield per area is almost the lowest of all rice-producing countries, due to many factors, e.g., diseases and environmental stress (Meyer and Prasertsri 2007).In the past three decades, Thailand has developed many elite rice varieties to increase the yield and quality of rice production.However, no information about the population gene pool of these elite rice varieties developed in Thailand is available.Many of them have the same parents, indicating that they share the same gene pool.It is very important to understand the genetic relationships of these developed elite rice varieties to ensure the sustainability of a rice breeding program.Genetic diversity plays a major role in plant breeding programs.Knowledge of the genetic relationships between the different accessions can be extremely useful in the development of efficient germplasm management and utilization strategies (Barrett et al. 1998.)The importance of this kind of information has been demonstrated in genetic fingerprinting and genotype selection programs (Fufa et al. 2005) studies.Insertion-Deletion (InDel) markers result from the insertion of transposable elements, slippage in simple sequence replication or unequal crossover events (Britten et al. 2003).InDel polymorphisms can be easily detected based on polymerase chain reaction (PCR) fragment length polymorphisms (Bhattramakki et al. 2002).InDels are also deemed a major source of gene defects and a significant source of evolutionary change (Britten et al. 2003).Because of their high abundance, reproducibility, co-dominant inheritance, and extensive genome coverage, InDels have been applied successfully for the analysis of genetic diversity and population structure in many plant species (García-Lor et al. 2012, Ollitrault et al. 2012, Wu et al. 2013).They have also been widely applied for genotyping, QTLs mapping, mapbased cloning, and even marker-assisted selection, including Arabidopsis (Hou et al. 2010), rice (Patel et al. 2014), wheat (Raman et al. 2006) and citrus (García-Lor et al. 2012).Our objective in this study was to use InDel markers distributed throughout the 12 rice chromosomes so as to analyze the genetic diversity and relationships of improved Thai rice cultivars.

Plant materials
Seeds of 15 improved Thai rice cultivars (Table 1) were collected and germinated in a tray filled with soil until the seedlings were 22 days old.Young and fresh leaf samples were collected for genomic DNA extraction.

DNA extraction
For each rice cultivar, young leaves of five 3-week-old plants were pooled and used for genomic DNA extraction according to the modified CTAB method (Doyle and Doyle 1987).The quality and concentration of the DNA samples were checked by 1% agarose gel electrophoresis and with a NanoDrop spectrophotometer.The DNA samples were diluted to a final concentration of 50 ng/μl and stored at -20 °C until use.

Data and statistical analysis
The presence (1) or absence (0) of amplified DNA fragments were scored.The discriminatory power of each InDel marker was determined by calculating the polymorphic information content (PIC).The PIC was calculated according to the formula PIC = 1 -Σ(fi 2 ) where fi represents the frequency of the ith allele of each InDel locus (Anderson et al. 1993).The genetic similarity coefficient, Nei-Li or Dice coefficient and Jaccard coefficient, with a bootstrap analysis of 1,000 for measuring pairwise band similarities between individuals was calculated, using the NTSYSpc version 2.01 (Rohlf 2000).Cluster analysis of the rice cultivars based on the genetic distance matrix was performed with the UPGMA (unweighted pair group method based on arithmetic averages) method (Sneath and Sokal 1973) using the SAHN function of the NTSYSpc version 2.01 to illustrate the genetic relationships.Tree Plot in Graphics was also used to draw a cluster analysis tree.Nei's gene diversity was calculated using POPGEN32 version 1.31 (Yeh and Yang 1999).Analysis of molecular variance (AMOVA) was performed to justify the suitability of group classification using program GenAlEx, version 6.5 (Peakal and Smouse 2012).

InDel marker polymorphism
A total of 133 InDel markers distributed throughout 12 rice chromosomes were screened with 15 improved Thai rice DNA samples.Of these, 14 InDel markers were located on rice chromosome 1.Ten InDel markers were located on rice chromosome 2. Five, nine, nine, nine, eleven, eight, six, seven, seven and five markers were located on rice chromosome 3, 4, 5, 6, 7, 8, 9, 10, 11 and 12, respectively (Figure 1).Ninty-nine markers (74.44 %) successfully amplified rice genomic DNA and resulted in 151 clearly visible bands.Of 99 InDel markers, 44 (44.44%) were polymorphic and identified a total of 88 bands and on average 2 bands per marker (Table 2).Our result was similar to that reported by Patel et al (2014), with an average number of 4.42 DNA bands for colored and white rice.The InDel PIC scores ranged from 0.06 to 0.50, with an average of 0.30 and the genetic similarity values ranged from 0.64 to 1.00 (Table 2).The frequency of each amplified band ranged from 0.06 (present in 1 of 15 cultivars) to 0.93 (present in 14 of 15 cultivars) (Figure 2).Since rice is a self-pollinating species, most genes are normally in the homozygous form.Twelve of 99 InDel loci were found to contain heterozygous genotypes, which may be due to pollen contamination in the rice seeds.Among these 12 InDel loci, rice cultivar Chinat1 is heterozygous in 10 loci.This result suggested that the purity of Chinat1, used in this study, was unsatisfactory.Almost half (44/99) of the InDel markers used in this study detected polymorphic DNA bands that were suitable for diversity studies.The discriminating power of the PIC scores of InDel markers was considered high and the markers were informative to discriminate the samples from each other.There were 13 InDel markers that produced specific DNA bands to one rice cultivar.The InDel markers RD0102 and RD0410 identified specific DNA bands to Jao Hom Nin; RD0112 to RD8, RD0211; RD0305 to Chainat1; RD0504 to Supanburi3; RD0607 to RD41, RD0601, RD0602; RD0610 to Supanburi60; and RD0707 and RD1101 to Pathumtani1 (Figure 2).These markers can be used to identify a variety or for hybrid screening in a rice breeding program.

Genetic relationship analysis
Forty-four InDel markers, used to analyze the genetic relationships, allowed the discrimination and identification of 14 of 15 rice cultivars.Two cultivars (KDML105 and RD6) could not be distinguished from each other.One possible explanation for this is that RD6 was developed from KDML105 by gamma radiation mutation.Genetically they were more than 99% identical.RD6 has almost the same genetic background as KDML105, except that RD6 is glutinous and KDML105 non-glutinous rice.The waxy (Wx) gene regulates the quantitative level of the Wx protein and the amylose content.RD6 has a low (0%) but KDML a high amylose content (16%) (Wanchana et al. 2003).The genetic relationships among all rice cultivars were analyzed and observed by cluster analysis (UPGMA) based on Euclidean distances, by which dissimilarities between the studied rice cultivars were highlighted.Three major groups consisted of 2, 9 and 4 rice cultivars, separated at a similarity index value of 0.758 (Figure 3).The first major group consisted of two cultivars, namely Chainat1 and Suphanburi3, which were clearly divergent from the rest.The second major group consisted of nine cultivars with two sub-groups.Sub-group A consisted of three cultivars, i.e., Suphanburi1, RD31 and RD41.Sub-group B consisted of six cultivars, including Pissanulok2, RD49, RD47, Pathumthani1, Suphanburi60 and Jao Hom Nin.The third major group was formed by four cultivars, which were KDML105, RD6, RD8 and Niaw Ubon1.The rice cultivars RD8 and NeawUbon1 were placed close to each other, with a genetic similarity of 0.83.The reason is that these two rice cultivars were generated from the same parental cross, Niawsaphatong2 and IR262.Nei's gene diversity index among cultivars was 0.3542, which indicated the moderate level of gene diversity.The gene diversity within cultivars could not be calculated because each cultivar was represented by pooled DNA.

Figure 1 .
Figure 1.Linkage map showing the position of 99 InDel markers that can amplify the genomic DNA of rice cultivars used in this study (modified from Wu et.al. 2013).* polymorphic markers.

Figure 3 .
Figure 3. Dendrogram of the relationships of InDel markers of 15 Thai rice cultivars based on UPGMA cluster analysis.

Table 2 .
Polymorphism in 15 rice cultivars revealed by InDel primers