Convergence of Goals: Phylogenetical,Morphological, and Physiological Characterization of Tolerance to Drought Stress in Tall Fescue (Festuca arundinacea Schreb.)

The aim of this study is to find Iranian tall fescue accessions that tolerate drought stress and investigation on phylogenetical, morphological, and physiological characterization of them. For this propose, inter-simple sequence repeats (ISSR) markers were used to examine the genetic variability of accessions from different provinces of Iran. Of 21 primers, 20 primers generated highly reproducible fragments. Using these primers, 390 discernible DNA fragments were produced with 367 (93.95 %) being polymorphic. The polymorphic information content (PIC) values ranged from 0.948 to 0.976, with a mean PIC value of 0.969. Probability identity (PI) and discriminating power (D = 1 ? PI) among the primers ranged from 0.001 to 0.004 and 0.998 to 0.995, respectively. A binary qualitative data matrix was constructed. Data analyses were performed using the NTSYS software and the similarity values were used to generate a dendrogram via UPGMA. To study the drought stress, plants were irrigated at 25 % FC condition for three times. Fresh leaves were collected to measure physiological characters including: superoxide dismutase, catalase, and peroxidase activities and proline and total chlorophyll content at two times, before and after stress application. Relative water content, fresh and dry weight ratio, survival percentage, and visual quality were evaluated after stress. Morphological and physiological characters were assessed in order to classify accessions as either tolerant or sensitive using Ward’s method of Hierarchical cluster analysis in SPSS software. The results of present study demonstrated that the ISSR markers are useful for studying tall fescue genetic diversity. Convergence of morphological and physiological characterizations during drought stress and phylogenetic relationship results showed that accessions can be grouped into four clusters; drought-tolerant accessions that collected from west of Iran, drought-tolerant accessions collected from northwest of Iran, drought semi-tolerant accessions collected from center of Iran, and drought-sensitive accessions collected from north of Iran. Data presented could be used to classify the tall fescue accessions based on suitability of cultivation in the regions studied or the regions with the similar environmental condition.     

Evaluation of genetic diversity amongst <Emphasis Type="Italic">Descurainia sophia</Emphasis> L. genotypes by inter-simple sequence repeat (ISSR) marker

Descurainia sophia is a valuable medicinal plant in family of Brassicaceae. To determine the range of diversity amongst D. sophia in Iran, 32 naturally distributed plants belonging to six natural populations of the Iranian plateau were investigated by inter-simple sequence repeat (ISSR) markers. The average percentage of polymorphism produced by 12 ISSR primers was 86 %. The PIC values for primers ranged from 0.22 to 0.40 and Rp values ranged between 6.5 and 19.9. The relative genetic diversity of the populations was not high (Gst =0.32). However, the value of gene flow revealed by the ISSR marker was high (Nm = 1.03). UPGMA clustering method based on Jaccard similarity coefficient grouped the genotypes into two major clusters. Graph results from Neighbor-Net Network generated after a 1000 bootstrap test using Jaccard coefficient, and STRUCTURE analysis confirmed the UPGMA clustering. The first three PCAs represented 57.31 % of the total variation. The high levels of genetic diversity were observed within populations, which is useful in breeding and conservation programs. ISSR is found to be an eligible marker to study genetic diversity of D. sophia.     

Genetic diversity analysis of Cynodon dactylon (bermudagrass) accessions and cultivars from different countries based on ISSR and SSR markers

ISSR and SSR markers were used to evaluate genetic diversity among 33 Cynodon dactylon accessions and 22 cultivars from four different countries in order to provide information on how to improve the utilization of bermudagrass germplasms. Eighty eight bands were amplified by nine SSR primer combinations and 236 bands were observed from 23 ISSR primers. The results showed that 97.7% of the SSR primers and 86.9% of the ISSR primers were polymorphic. The genetic similarity coefficients (GSC), gene diversity (He) and Shannon index (I) were 0.58–0.97, 0.27 and 0.41, respectively, for ISSR and 0.52–0.97, 0.29, and 0.43 for SSR. The UPGMA analysis clustered the 55 accessions (cultivars) into three groups. The cluster results produced by the ISSR data were close to the SSR data results. Analysis based on the combined ISSR and SSR data was more closely related to the geographical distribution of the tested germplasm.     

Genetic Diversity in Various Accessions of Pineapple [<Emphasis Type="Italic">Ananas comosus</Emphasis> (L.) Merr.] Using ISSR and SSR Markers

Inter simple sequence repeat (ISSR) and simple sequence repeat (SSR) markers were used to assess the genetic diversity of 36 pineapple accessions that were introduced from 10 countries/regions. Thirteen ISSR primers amplified 96 bands, of which 91 (93.65%) were polymorphic, whereas 20 SSR primers amplified 73 bands, of which 70 (96.50%) were polymorphic. Nei’s gene diversity (h = 0.28), Shannon’s information index (I = 0.43), and polymorphism information content (PIC = 0.29) generated using the SSR primers were higher than that with ISSR primers (h =  0.23, I = 0.37, PIC = 0.24), thereby suggesting that the SSR system is more efficient than the ISSR system in assessing genetic diversity in various pineapple accessions. Mean genetic similarities were 0.74, 0.61, and 0.69, as determined using ISSR, SSR, and combined ISSR/SSR, respectively. These results suggest that the genetic diversity among pineapple accessions is very high. We clustered the 36 pineapple accessions into three or five groups on the basis of the phylogenetic trees constructed based on the results of ISSR, SSR, and combined ISSR/SSR analyses using the unweighted pair-group with arithmetic averaging (UPGMA) method. The results of principal components analysis (PCA) also supported the UPGMA clustering. These results will be useful not only for the scientific conservation and management of pineapple germplasm but also for the improvement of the current pineapple breeding strategies.     

Genetic diversity and structure of <Emphasis Type="Italic">Capparis spinosa</Emphasis> L. in Iran as revealed by ISSR markers

Capparis spinosa L. (caper bush) is an economically and ecologically important perennial shrub that grows across different regions of Iran. In this study, the genetic diversity and population structure of Iranian genepool of C. spinosa is evaluated using Inter Simple Sequence Repeat (ISSR) markers. Using 10 ISSR primers, 387 DNA fragments (bands) were amplified from the genomic DNA of 92 individuals belonging to twenty-one populations of C. spinosa, of which 378 (97.7%) were polymorphic. High level of genetic diversity (percentage of polymorphic loci = 98.2%, h = 0.1382, I = 0.243), high genetic differentiation (G_st = 0.5234) and low gene flow (Nm = 0.4553) among populations were observed. Caper bush populations were divided into 4 groups in the dendrogram, PCoA plot and Bayesian clustering results, mostly corresponded to their geographic regions. The results showed that there are value in sampling Iranian caper bush populations to look for valuable alleles for use in plant breeding programs.     

Evaluation of genetic diversity of <Emphasis Type="Italic">Clinacanthus nutans</Emphasis> (Acanthaceaea) using RAPD,ISSR and RAMP markers

Three polymerase chain reaction (PCR) techniques were compared to analyse the genetic diversity of Clinacanthus nutans eight populations in the northern region of Peninsular Malaysia. The PCR techniques were random amplified polymorphic deoxyribonucleic acids (RAPD), inter-simple sequence repeats (ISSR) and random amplified microsatellite polymorphisms (RAMP). Leaf genomic DNA was PCR amplified using 17 RAPD, 8 ISSR and 136 RAMP primers . However, only 10 RAPD primers, 5 ISSR primers and 37 RAMP primers produced reproducible bands. The results were evaluated for polymorphic information content (PIC), marker index (MI) and resolving power (RP). The RAMP marker was the most useful marker compared to RAPD and ISSR markers because it showed the highest average value of PIC (0.25), MI (11.36) and RP (2.86). The genetic diversity showed a high percentage of polymorphism at the species level compared to the population level. Furthermore, analysis of molecular variance revealed that the genetic diversity was higher within populations, as compared to among populations of C. nutans. From the results, the RAMP technique was recommended for the analysis of genetic diversity of C. nutans.     

Genetic diversity of some Iranian sweet cherry (<Emphasis Type="Italic">Prunus avium</Emphasis>) cultivars using microsatellite markers and morphological traits

The aim of this study was to characterize 23 important Iranian sweet cherry (Prunus avium) cultivars collected from different provinces of Iran and 1 foreign cultivar, which was used as control, considered for breeding programs by using 21 microsatellite markers and 27 morphological traits. In sweet cherry (Prunus avium) accessions, leaf, fruit, and stone morphological characters were evaluated during two consecutive years. The study revealed a high variability in the set of evaluated sweet cherry accessions. The majority of important correlations were determined among variables representing fruit and leaf size and variables related to color. Cluster analysis distinguished sweet cherry accessions into two distinct groups. Principal component analysis (PCA) of qualitative and quantitative morphological parameters explained over 86.59% of total variability in the first seven axes. In PCA, leaf traits such as leaf length and width, and fruit traits such as length, width, and weight, and fruit flesh and juice color were predominant in the first two components, indicating that they were useful for the assessment of sweet cherry germplasm characterization. Out of 21 SSR markers, 16 were polymorphic, producing 177 alleles that varied from 4 to 16 alleles (9.35 on average) with a mean heterozygosity value of 0.82 that produced successful amplifications and revealed DNA polymorphisms. Allele size varied from 95 to 290 bp. Cluster analyses showed that the studied sweet cherry genotypes were classified into five main groups based mainly on their species characteristics and SSR data. In general, our results did not show a clear structuring of genetic variability within the Iranian diffusion area of sweet cherry, so it was not possible to draw any indications on regions of provenance delimitation. The results of this study contribute to a better understanding of sweet cherry genetic variations in Iran, thus making for more efficient programs aimed at preserving biodiversity and more rational planning of the management of reproductive material.     

Association of SSR Markers with Cold Tolerance Traits in Diverse Bermudagrass [Cynodon dactylon (L.) Pers.] Accessions

Bermudagrass [Cynodon dactylon (L.) Pers.] is a perennial and typical warm-season grass. It undergoes withering and severe damage under cold stress; thus, cold is considered as a key factor that restricts the widespread use in bermudagrass. Identification of association between molecular markers and cold tolerance-related traits would facilitate the efficient selection of cold tolerant bermudagrass cultivars. A total of 106 diverse bermudagrass accessions, including 4 commercial cultivars and 102 wild germplasms, were tested for cold tolerance and analyzed by 104 simple sequence repeat (SSR) markers. Cold significantly decreased transpiration rate, growth rate and turf quality. There were significant variations in these trait values among the accessions under cold conditions. Two subpopulations were detected in the panel of accessions based on the analysis of 1474 alleles with 104 SSR markers. Clustering analysis revealed that the genetic relationship was affected by the natural habitats. Thirty-four SSR markers were identified to be associated with two or three traits based on the corrected P values (P < 3.5 × 10⁻⁴). These markers can be used for genetic improvement of cold tolerance of bermudagrass after further validation.

     

Genetic diversity of <Emphasis Type="Italic">Opuntia</Emphasis> spp. varieties assessed by classical marker tools (RAPD and ISSR)

Opuntia, commonly named “nopal” in Mexico, is an important crop for its agronomical, economical, ecological and cultural value. Furthermore, it is known for its taxonomic complexity. In this paper, we report the genetic variability of 52 Opuntia cultivars with agronomic and economic importance, classified into 12 different species using random amplified polymorphic DNA (RAPD) and inter-simple sequence repeats (ISSR) markers. Ten primers, five for each marker type, were selected to assess their ability to detect polymorphisms in this plant accesions/varieties. Both marker systems generated a total of 307 bands, of which 50.8 % were polymorphic with an average of 15.6 polymorphic bands per primer. Thus, we assume that Mexican Opuntia varieties present broad genetic variation. Based on percentage of polymorphic bands; resolving power; polymorphic information content; and Marker Index, the K-12 (RAPD) and IS-06 (ISSR) primers were the most informative ones. Clusters obtained from RAPD, ISSR and a combination of both data sets did not match the actual taxonomic classification. On the other hand, the putative varieties currently classified in the same species were not located in the same cluster. Besides, the varieties included in O. ficus-indica, O. albicarpa and O. megacantha showed broad variation but were not well defined into separate clades; these cultivars possibly have common ancestry. The results presented here support our hypothesis about the existence of a smaller number of Opuntia species in accordance with those currently described, but with high intraspecific genetic variation.     

Morphological and molecular genetic variations of oat genotypes grown in Kermanshah,Iran

Morphological traits and molecular markers are two common methods for genetic variation studies. Molecular markers, morphological traits methods and relationship between the two were used to study genetic variation among 43 oat genotypes and varieties. For this purpose, an augmented design was conducted in three replicates at 2008–2009 cropping season in the experimental field of Campus of Agriculture and Natural Resources of Razi University, Kermanshah, Iran. Four wild oat accessions (Avena sterilis) were added to evaluated genotypes in molecular experiment. Results showed a significant variation among genotypes for all morphological traits and they were classified based on this variation in four groups by WARD cluster analysis. In molecular experiment, 28 inter simple sequence repeat (ISSR) primers amplified 206 polymorph bands. Based on Jaccard similarity matrix, similarity among genotypes was varied from 0.23 to 0.66 and cluster analysis classified genotypes in seven groups by complete linkage method. The correlation between ISSR marker and morphological traits classifications was not significant. ISSR showed to be a helpful marker for genotype identity and separation as it put wild accessions in a group.     

Estimation of genetic diversity in some Iranian cornelian cherries (Cornus mas L.) accessions using ISSR markers

Iran is an important producer of cornelian cherries (Cornus mas L.). Seed propagation and long term human selection have given rise to a great diversity of trees. In this study, ISSR marker was used for genetic diversity evaluation of 40 accessions of cornelian cherries (C. mas L.). With 20 ISSR primers, 171 polymorphic bands were detected with polymorphism ratio range of 80–100%. The average polymorphism information content (PIC) was 0.46, which shows that the majority of primers are informative. The average values of Effective Multiplex Ratio (EMR), Marker Index (MI), effective number of alleles (Ne), Nei's gene diversity (He) and Shannon's information index (I) for all primers were 8.0, 3.75, 1.756, 0.416 and 0.595 respectively. Based on these results, ISSR analysis can be used for the characterization and grouping of cornelian cherry accessions. This is the first study demonstrating that ISSR analyses can be used to differentiate and classify cornelian cherry accessions.     

Estimation of genetic variation in Cynodon dactylon accessions using the ISSR technique

Genetic variation in 55 accessions of Cynodon dactylon was estimated using inter-simple sequence repeat (ISSR) markers. The plant materials used in this study originated from 17 countries. A total of 236 ISSR fragments were generated with 14 primers. Fragment sizes ranged from 200 to 3000 bp. All scorable bands were polymorphic in nature and none of the primers used produced monomorphic bands, indicating a high level of genetic variation in this grass. The accessions were found to be clustered into eight major groups through the unweighted pair-group method with arithmetic averages. Genetic similarity coefficients (GSC) among the 55 accessions ranged from 0.52 to 0.95. The results clearly indicate that a high level of variation exists in Cynodon accessions. This study shows that the ISSR technique is a reliable tool for differentiating Cynodon accessions and for determining the genetic relationships among them.     

Assessment of genetic variability through ISSR and RAPD markers in <Emphasis Type="Italic">Commiphora wightii</Emphasis> (Arn.) Bhandari

Commiphora wightii (Arn.) Bhandari is a commercially, medicinally and traditionally important tropical shrub widely used to treat various ailments and disorders. Demand of this plant is increasing in the pharmaceutical and perfumery industries due to the presence of guggulsterone E and Z, two important isomers conferring lipid- and cholesterol-lowering, and anti-cancerous properties. Ruthless and unscientific harvesting of oleo-gum resin by local populations from the wild, with negligible conservation efforts has made this species endangered and led to its inclusion in the Red Data Book of IUCN. It is imperative to have broad information regarding the extent of genetic variability available in the species to accelerate the breeding and conservation programs. Therefore, the present study was undertaken to analyze the extent of genetic variability existing among the C. wightii germplasm collected from Rajasthan and Haryana, the diversity rich Indian states, using ISSR and RAPD markers. A total of 100 (50 each) RAPD and ISSR markers were screened of which 37 RAPD and 43 ISSR primers were able to amplify DNA fragments. RAPD markers were more efficient, detecting 74.16 % polymorphism, compared to ISSR which detected 62.52 % polymorphism. Also, the values of average number of polymorphic bands per assay, polymorphism information content (PIC), diversity index (DI) and marker index (MI) were more for RAPD (7.76, 0.19, 0.38 and 2.53, respectively) than for ISSR (7.02, 0.13, 0.32 and 1.88) markers. The UPGMA dendrogram constructed using individual as well as combined data of the two marker systems separated the collected accessions into two major clusters containing 47 and 4 accessions, respectively, while one accession from Bikaner was not included in any cluster. Genetic similarity values obtained from Jaccard’s coefficient using combined data of both the marker systems were between 0.50 and 0.97. These results indicated the existence of wide genetic variability within this species and can be used for further research in the area of germplasm conservation, population genetics and plant breeding.     

AFLP,RAPD, and ISSR analysis of intraspecific polymorphism and interspecific differences of allotetraploid species <Emphasis Type="Italic">Aegilops kotschyi</Emphasis> Boiss. and <Emphasis Type="Italic">Aegilops variabilis</Emphasis> Eig

To evaluate genetic variation, 27 accessions of allotetraploid species Aegilops kotschyi and Ae. variabilis with the US genome were analyzed using the AFLP, RAPD, and ISSR methods. A total of 316 polymorphic RAPD fragments, 750 polymorphic AFLP fragments, and 234 polymorphic ISSR fragments were obtained. It was demonstrated that the analyzed species were characterized by a considerable level of nuclear genome variation. According to the data of ISSR and RAPD analysis, the average value of the Jaccard similarity coefficient for the accessions of Ae. variabilis from different geographical regions was slightly lower than that for the accessions of Ae. kotschyi. At the same time, AFLP analysis showed no considerable differences in the levels of intraspecific variation of the studied species. Analysis of the summarized RAPD, ISSR, and AFLP marking data in the Structure software program showed that most of the analyzed accessions with high degree of probability could be assigned to one of two groups, the first of which corresponded to Ae. kotschyi and the second corresponded to Ae. variabilis, thereby confirming the species independence of Ae. kotschyi and Ae. variabilis. Accessions k900, k907, k908, and v90 could not with a sufficiently high degree of probability be assigned to one of the species, which possibly was the result of interspecific hybridization. Analysis of the species diversity using different molecular markers made it possible to identify the accessions that were notably different from other accessions of its species.     

Assessment of genetic diversity in Ethiopian cowpea [<Emphasis Type="Italic">Vigna unguiculata</Emphasis> (L.) Walp.] germplasm using simple sequence repeat markers

The genetic diversity of cowpea (Vigna unguiculata L. Walp.) in Ethiopia was analyzed using 19 uniform accessions, 62 variable accessions (yielding 185 sub-types), and two mungbean (Vigna radiata) accessions (four subtypes) as outgroup. A set of 23 polymorphic simple sequence repeat (SSR) markers was identified, and polymorphism in the various accessions was scored by determining amplicon variability. Allele frequency, genetic diversity, and polymorphism information content (PIC) were determined for each SSR marker, and a neighbor joining dendrogram was generated to show the genetic relationship among the individual accessions. A total of 75 allelic variants was defined, with the average number of alleles per locus calculated to be three. The average genetic diversity (D) was 0.47, and PIC was 0.4. Three main clusters were identified by phylogenetic analysis, and the clusters and sub-grouping were supported by STRUCTURE and principal component analysis. This grouping had a moderate fixation index value of 0.075 and gene flow (Nm) of 3.176, indicating that the accessions possess wide diversity within individuals and populations. The accessions showed no clustering by geographical origins. Three well-characterized molecular markers (SSR1, C42-2B, and 61RM2) for race specific resistance to Striga gesnerioides in the cowpea cultivar B301 were used to evaluate the accessions for their potential for use in genetic improvement against this pest. Based on this analysis, only two accessions, 222890–2 from Gambela and 286–2 from the Southern Nations, Nationalities, and Peoples (SNNP) region, were found to cluster with B301 and contain the SSR1 resistance allele. These findings will assist in germplasm conservation efforts by the Institute of Biodiversity and Conservation of Ethiopia, and contribute to future studies aimed at the genetic improvement of local germplasm for improved overall agronomic performance as well as Striga resistance in particular.     

Cytotype-specific ISSR profiles and karyotypes in the Neotropical genus <Emphasis Type="Italic">Eigenmannia</Emphasis> (Teleostei: Gymnotiformes)

The genus Eigenmannia (Teleostei: Gymnotiformes), a widely distributed fish genus from the Neotropical region, presents very complex morphological patterns and many taxonomic problems. It is suggested that this genus harbors a species complex that is hard to differentiate using only morphological characteristics. As a result, many species of Eigenmannia may be currently gathered under a common name. With the objective of providing new tools for species characterization in this group, an analysis of the polymorphism of DNA inter-simple sequence repeats (ISSR), obtained by single primer amplification reaction (SPAR), combined with karyotype identification, was carried out in specimens sampled from populations of the Upper Paraná, São Francisco and Amazon river basins (Brazil). Specific ISSR patterns generated by primers (AAGC)₄ and (GGAC)₄ were found to characterize the ten cytotypes analyzed, even though the cytotypes 2n = 38 and 2n = 38 XX:XY, from the Upper Paraná basin, share some ISSR amplification patterns. The geographical distribution of all Eigenmannia specimens sampled was inferred, showing the cytotype 2n = 31/2n = 32 as the most frequent and largely distributed in the Upper Paraná basin. The cytotype 2n = 34 was reported for the first time in the genus Eigenmania, restricted to the São Francisco basin. Polymorphic ISSR patterns were also detected for each cytotype. Considering our results and the data reported previously in the literature, it is suggested that many of the forms of Eigenmannia herein analyzed might be regarded as different species. This work reinforces the importance of employing diverse approaches, such as molecular and cytogenetic characterization, to address taxonomic and evolutionary issues.     

Simple Sequence Repeat and <Emphasis Type="Italic">S</Emphasis>-locus Genotyping to Explore Genetic Variability in Polyploid <Emphasis Type="Italic">Prunus spinosa</Emphasis> and <Emphasis Type="Italic">P. insititia</Emphasis>

Polyploid Prunus spinosa (2n = 4×) and P. insititia (2n = 6×) represent enormous genetic potential in Central Europe, which can be exploited in breeding programmes. In Hungary, 17 cultivar candidates were selected from wild-growing populations including 10 P. spinosa, 4 P. insititia and three P. spinosa × P. domestica hybrids (2n = 5×). Their taxonomic classification was based on their phenotypic characteristics. Six simple sequence repeats (SSRs) and the multiallelic S-locus genotyping were used to characterize genetic variability and reliable identification of the tested accessions. A total of 98 SSR alleles were identified, which presents 19.5 average allele number per locus, and each of the 17 genotypes could be discriminated based on unique SSR fingerprints. A total of 23 S-RNase alleles were identified. The complete and partial S-genotype was determined for 8 and 9 accessions, respectively. The identification of a cross-incompatible pair of cultivar candidates and several semi-compatible combinations help maximize fruit set in commercial orchards. Our results indicate that the S-allele pools of wild-growing P. spinosa and P. insititia are overlapping in Hungary. A phylogenetic and principal component analysis confirmed the high level of diversity and genetic differentiation present within the analysed genotypes and helped clarify doubtful taxonomic identities. Our data confirm that S-locus genotyping is suitable for diversity studies in polyploid Prunus species. The analysed accessions represent huge genetic potential that can be exploited in commercial cultivation.     

Identification and pyramiding of QTLs for cold tolerance at the bud bursting and the seedling stages by use of single segment substitution lines in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

Cold stress is one of the main constraints in rice production, and damage from cold can occur at different developmental stages in rice. Understanding the genetic basis of cold tolerance is the key for breeding cold-tolerant variety. In this study, we used single segment substitution lines (SSSLs) derived from a cross between cold-tolerant japonica variety “Nan-yang-zhan” and a popular indica variety “Hua-jing-xian 74” to detect and pyramid QTLs for cold tolerance at the bud bursting and the seedling stages. Evaluation of cold tolerance of these SSSLs and their recurrent parent helped identify two cold-tolerant QTLs (qCTBB-5 and qCTBB-6) at the bud bursting stage and two cold-tolerant QTLs (qCTS-6 and qCTS-12) at the seedling stage. The SSSLs carrying these QTLs showed stronger cold tolerance than their recurrent parent HJX74 did in three independent experiments. The qCTBB-6 and qCTS-6 were mapped to the same chromosomal region. QTL pyramiding was performed by intercrossing of SSSLs carrying the respective QTLs for cold tolerance at the bud bursting stage and the seedling stage and marker-assisted selection (MAS). The selected pyramiding line SC1-1 with different cold-tolerant QTLs showed cumulative effects on cold tolerance. Our results suggest that different genes (QTLs) control cold tolerance at bud bursting and seedling stages, and pyramiding of stable expression QTLs for cold tolerance at different developmental stages through MAS is a good strategy to prevent cold damage in rice.