Introduction

Rice (Oryza sativa L.) is one of the major food crops and a staple food for over 50% of the world population. Asia is responsible for more than 80% of global rice consumption, and China is the largest (30%) rice producer and consumer in the world (http://ricepedia.org). The northeast region is one of the leading rice-producing regions in China. Additionally, it is also classified as a high-latitude cold-region rice cultivation area along with Hokkaido, Japan. “Dry breeding at lower planting density” is the dominant plantation strategy in cold areas that grow rice at high latitudes. Using this strategy, strong seedlings were cultivated that reduced the damage caused by weeds (Farooq et al. 2011; Olajumoke et al. 2016; Singh et al. 2001). The economic development of China has reduced the rural labor force and increased labor costs, which have dramatically increased the cost of growing rice through the “Dry breeding at lower planting density” planting strategy. Contrastingly, the direct seeding method reduces seedling field, water, labor, and effort (Feng et al. 2020; Ya et al. 2020; Zhang et al. 2015; Zhen 2015). This method not only significantly reduces the cost of rice cultivation but also significantly alleviates labor shortages. However, due to frequent occurrences of low-temperature damage in the spring in cold regions at high latitudes, rice that can germinate and grow at low temperatures is required for direct seeding. Unfortunately, the existing rice varieties cannot germinate well at low temperatures, which adversely affects the development of direct-seeded rice (Cao et al. 2018; Chen et al. 2006).

Many low-temperature germination-related QTLs are located on the 12 chromosomes; however, only a few can be used for molecular-assisted breeding. Currently, only qLTG-3 has been cloned (Os03 g0103300), which is the only gene related to low-temperature germination (Fujino et al. 2008). In this study, 94 recombinant inbred lines (RILs) from Heigu and Ha 9366 were used to perform the QTL analysis on low temperature-resistant germination traits, which might provide a theoretical basis for further fine positioning and the development of molecular markers for molecular marker-assisted breeding (MAS).

The cold tolerance of rice at the germination stage is mainly identified by its ability to germinate at low temperatures, which is evaluated by identifying the seedling vigor and viability at low temperatures. The low-temperature germination of crops is a quantitative trait controlled by multiple genes. The QTLs associated with low-temperature germination have been detected on all 12 chromosomes. There are 97 QTLs associated with cold tolerance at the germination stage. More QTLs were detected on chromosomes 3, 5, 7, and 11, and fewer on chromosomes 2, 9, and 12 (Sheng et al., 2001; Fujino et al. 2004; Chen and Li 2005; Qiao et al. 2005; Han et al. 2007; Li et al., 2008; Gong et al. 2009; Li et al. 2011; Yang et al. 2016; Jiang et al. 2017; Yang et al. 2018) (Table 1).

Table 1 QTLs related to cold tolerance by germination stage in rice

Although many QTLs associated with low-temperature germination have been localized on the 12 chromosomes, the number of molecular markers, as well as, the information that can facilitate the screening of suitable direct-seeded rice varieties with cold tolerance, is insufficient. In this study, 94 recombinant inbred lines (RILs) and 32 advanced generation lines of Heigu and Ha 9366 were used for QTL analysis and to evaluate the number of QTLs associated with cold tolerance during germination. We aimed to provide a theoretical basis for further fine localization and development of molecular markers for molecular marker-assisted breeding (MAB).

Materials and methods

Experimental materials

The LTH/H9366 recombinant inbred line population, consisting of 94 lines, was obtained by the single seed transmission method using the low-temperature germination variety Heigu (LTH) and the conventional high-latitude cold rice variety Ha 9366 (H9366) as parents. The advanced generation lines are rice varieties suitable for cultivation in the Heilongjiang Province, China, above generation F9.

In 2017, two parents and 94 recombinant inbred lines were planted in the paddy field base of China Heilongjiang Academy of Agricultural Sciences Democracy Park (E: 126.8, N: 45.4). The planting density is 4 inches *8 inches, sowing in mid-April, harvesting in late September. After harvest, seeds were germinated at low temperature.

Low-temperature germination test

After the parents and population materials were harvested at maturity, we interrupted dormancy at 50 °C for seven days. Fifty filled seeds were randomly selected from each strain and placed in a 9 cm Petri dish lined with two layers of filter paper after disinfection. Each Petri dish was filled with 10 mL of sterile water and incubated at 15 °C in an incubator for germination. Three replications were performed for each line and parent. The number of germinated grains was recorded, and the germination rates were calculated on the fifth, sixth, and seventh day after germination, and the germination rate was calculated. Germination rate (%) = (total number of germinations 5, 6, and 7 days later)/50 grains. Successful germination was determined by the radicle or germ length >1 mm. After seven days of low-temperature treatment, the seeds were transferred to an incubator at 30 °C, and the germination rate was tested after four days. A seed germination rate above 80% indicated that the low-temperature germination ability was not affected by secondary dormancy (refer to partial method modifications by Jang et al., 2006).

Marker analysis

The CTAB method introduced by Murray and Thompson (1980) was used as a reference to extract the DNA of each strain. Parent screening was conducted for 350 SSR markers, and 125 SSR markers were finally selected for marker analysis. The PCR reaction volume was 20 µL, which included 1 µL of 5 µmol/L primers, 1 µL of DNA, and 2× TaqMix. The PCR conditions were: pre-denaturation at 94 °C for 5 min, denaturation at 94 °C for 30 s, renaturation at 52–60 °C for 30 s, and extension at 72 °C for 1 min, with a total of 35 cycles. Then, a final extension was performed at 72 °C for 5 min, and the products were stored at 4 °C. The PCR products were electrophoretically separated on a 4% agarose gel and then detected using an imaging system.

Construction of linkage mapping and QTL analysis

Mapmaker 3.0 and WinQTLCart 2.5 software were used for the QTL analysis of low-temperature germination traits of rice. Mapmarker 3.0 was used to construct a genetic linkage map and obtain the genetic linkage centimorgan (cM) between markers (Fig. 3); Win-QTLCart 2.5 was used to generate LOD curves to determine QTL loci and their distribution on the chromosomes. During the analysis, the effects of the main-effect markers and the interaction markers were controlled simultaneously. The likelihood ratio test and the t-test were used in combination to test the significance of the effect of QTL (P ≤ 0.005) and interaction effect (P ≤ 0.001), as well as for calculating the contribution rate.

Analysis of the utilization value of the QTL locus

Thirty-five rice lines were genotyped using SSR primers near the QTL. The method was similar to that described in “Marker analysis” section. The cold tolerance at the germination stage was determined as described in “Low temperature germination test” section. The cold tolerance at the germination stage was investigated using plants that had grown to the two-leaf stage by treating them at 5 °C for five days, and after treatment, the phenotype was investigated after five days of recovery under ambient conditions. The phenotypes were classified as grades 0 (cold sensitive), 1, 2, 3, and 4 (cold tolerance).

Results

In this study, parental Heigu (LTH) and Ha9366 plants were tested for germination at 15 °C. Heigu started to germinate on the fourth day; the germination rate reached 30.4%, 77.7%, and 88.9% on the fifth, sixth, and seventh day, respectively. Ha068143 started germinating on the fifth day at 15 °C, and the germination rates reached 4.5% and 30.2% on the sixth and seventh day, respectively. The low-temperature germination rate of the parents reached a very significant level (Fig. 1). The low-temperature germination of the recombinant inbred line population at 15 °C had a distribution range of 0–90.8% on the fifth day, with an average value of 24.1%, and a distribution range between 2.0 and 95.9% on the sixth day, with an average value of 49.9%, while the distribution range was 13.7-99.0%, with an average value of 69.2%. Overall, it showed an approximately normal distribution and a strong separation of super-parents (Fig. 2).

Fig. 1
figure 1

Differences in germination rates between parental lines at 15 °C

Fig. 2
figure 2

Frequency distributions for low-temperature germination in RILs. A Experimental treatment for 5 days. B Experimental treatment for 6 days. C Experimental treatment for 7 days. Black and white arrows indicate the mean values of LTH and H9366, respectively

Fig. 3
figure 3

The linkage mapping of recombinant inbred line populations between Heigu and Ha 9366

QTL distribution and effects of cold-tolerance germination in populations

A total of four QTLs for cold tolerance of rice at the germination stage were detected by the MapQTL 6.0 software on the fifth, sixth, and seventh day which were distributed on chromosomes 3 and 12 of rice. The QTLs of qLTG-12.5d, qLTG-12.6d, and qLTG-12.7d were detected between the marker range of RM415 and RM453 on chromosome 12, which improved cold tolerance at the germination stage on the fifth, sixth, and seventh days. The interval size was 2.48 cM, and the LOD values were 4.22, 4.05, and 4.12, respectively. The PVE values reached 25.9%, 24.9%, and 21.3%, respectively. This QTL enhanced the low-temperature tolerance of rice during germination when it was present in the Ha 9366 allele. On the eighth day of low-temperature treatment, we detected a QTL on chromosome 3, with LOD of 4.09, PVE of 22.1%, and improved the cold tolerance of rice at the germination stage when present in the Heigu allele (Table 1; Figs. 4 and 5).

Fig. 4
figure 4

The germination stage tolerance gene QTL-qLTG-3 is located on chromosome 3

Fig. 5
figure 5

Germination stage tolerance gene QTL-qLTG-12 is located on chromosome 12

Table 2 Results of MOM mapping for Low-Temperature Germination (LTG) in the RILs

The solid line represents the LOD value of germination rate on the fifth day at 15 °C; the double line represents the LOD value of germination rate on the sixth day at 15 °C; the dashed line represents the LOD value of germination rate on the seventh day at 15 °C; the dotted line represents the standard LOD value calculated by permeation test (P ≤ 0.005).

The solid line represents the LOD value of germination rate on the fifth day at 15 °C; the double line represents the LOD value of germination rate on the sixth day at 15 °C; the dashed line represents the LOD value of germination rate on the seventh day at 15 °C; the dotted line represents the standard LOD value calculated by permeation test (P ≤ 0.005).

Validation of cold tolerance QTL and analysis of interaction effects during germination

We validated the detected cold tolerance QTLs during the germination period and evaluated the interaction effects of these QTLs. First, the genotypes of qLTG-3 and qLTG-12, which govern cold tolerance at the germination stage, were distinguished, and a combination of two QTL sites and four alleles were selected from the recombinant inbred lines (qLTG-3 is the allele of H9366, qLTG-12 is also the allele of LTH; b: qLTG-3 is the allele of LTH, qLTG-12 is the allele of LTH; c: qLTG-3 is the allele of H9366, qLTG-12 is also the allele of H9366; d: qLTG-3 is the allele of LTH allele, and qLTG-12 is also the allele of H9366). Among the 94 recombinant inbred lines, there were 36 strains in category a, 17 strains in category b, 29 strains in category c, and 12 strains in category d. The qLTG-12 locus exhibited excellent low-temperature tolerance on the fifth, sixth, and seventh days of germination. The germination rates at different time points (five, six, and seven days) are shown in Fig. 6. The germination rate of category c was 16–27% higher than that of category a, and that of category d was 20–36% higher than category b, both reached a significant difference. Moreover, the results of the two-way ANOVA showed that the F values of qLTG-12 reached 10.78 to 14.81 on the fifth, sixth, and seventh days and the significance reached over 0.1% (Fig. 6; Table 2). The qLTG-3 locus showed no significant effect on the germination rate on day 5 and day 6, and only showed a significant difference for the germination rate on day 7; the F value and significance were 4.52 and 5%, respectively. The HTL allelic type effectively increased the germination rate from 4 to 8% under low-temperature conditions in this population (Fig. 6; Table 2). The above results suggested that the interaction effects between the QTLs, which were detected for the cold-tolerant recombinant inbred lines during germination, were mainly additive.

Fig. 6
figure 6

Interactions between two QTLs related to low-temperature germination on chromosomes 3 (qLTG-3) and 12 (qLTG-12). (I) The 5-day germination rate at 15 °C, (II) The 6-day germination rate at 15 °C, (III) The 7-day germination rate at 15 °C. The alleles of the QTLs were estimated by their flanking markers (Table 2). L and H represent alleles of LTH and Ha9366, respectively. The bar plots show germination rate (%) ±SD. Different letters above the bars indicate a significant difference based on the Tukey–Kramer test

Table 3 Two-way ANOVA analysis for the effects of one QTL and one gene for Low-Temperature Germination (LTG) in the RILs

Verification of cold tolerance of qLTG-3 and qLTG-12 at the germination stage

Thirty-two rice lines suitable for cultivation in the Heilongjiang Province, China, were genotyped using the qLTG-3 and qLTG-12-related quantitative SSR primers RM232 and RM20A. Four allelic combinations were obtained (a: qLTG-3 was the H9366 allele and qLTG-12 was also the LTH allele; b: qLTG-3 was the LTH allele and qLTG-12 was the LTH allele; c: qLTG-3 was the H9366 allele and qLTG-12 was the H9366 allele; d: qLTG-3 was the LTH allele and qLTG-12 was the H9366 allele) (Fig. 7). The cold resistance test at the germination stage was also conducted for the 32 lines. The 32 lines were distributed from grades 0 to 4 for cold resistance. According to the genotype classification, grade d had the best cold tolerance, with a mean value of 3.7, which was significantly higher than grades a, b, and c. The mean values for cold tolerance in grades a and b were 0.7 and 1.4, respectively, and the difference between the two grades (0.7) was not significant (Fig. 8).

Fig. 7
figure 7

Genotypic analysis of 32 rice lines using qLTG-3 and qLTG-12

Fig. 8
figure 8

Interactions between two QTLs for the low-temperature seeding stage on chromosomes 3 (qLTG-3) and 12 (qLTG-12). L and H indicate the LTH and Ha9366 alleles, respectively. The bar plots show the germination rate (%) ±SD. Different letters above the bars indicate a significant difference based on the Tukey-Kramer test

Discussion

It is well-known that direct-sown rice can reduce the time and cost of growing rice, but it is difficult to apply this method in the northernmost areas of rice cultivation, such as northern China and Hokkaido, Japan, where the temperature is a limiting factor. Many QTLs related to low temperature at the germination stage have been reported. However, there is a lack of information on the actual molecular marker-assisted breeding among the QTLs.

Two QTLs associated with low-temperature germination, i.e., qLTG-3 and qLTG-12, were detected in this study. The qLTG-3 locus was located near the qLTG-3.1 locus, as reported by Funino (2008). Thus, qLTG-3, detected on chromosome 3 in this study, was probably influenced by the same gene as qLTG-3.1. The qLTG-3.1 locus is the main QTL controlling germination at low temperatures and can substantially increase the germination rate under low-temperature conditions. The main QTL locus detected in this study was qLTG-12. Chromosome 12 is the least reported of the many QTLs for low-temperature germination, with only three studies providing evidence of its role (Ji et al. 2007; Hou et al. 2004; Yao et al. 2020). The qLTG-12 locus detected in this study was different from the three previously reported loci and can be considered to be a novel QTL for the study on cold tolerance at the germination stage in rice. Reciprocal effects are important and can improve molecular breeding in crops, while research on the supernumerary reciprocity in low-temperature germination of rice is lacking (Zhang et al. 2005; Ji et al. 2008). In this study, we detected a QTL (qLTG-12) that possessed a strong epistatic interaction effect with the main low-temperature tolerance germination QTL (qLTG-3.1). A QTL locus with cold tolerance was reported near qLTG-12. The gene controlling qCTS-12 was OsGSTZ2, and Kim et al. (2011) found that the SNP locus in the OsGSTZ2 gene resulted in amino acid differences that were significantly associated with low-temperature sensitivity in rice seedlings. Additionally, Takesawa et al. (2020) reported that OsGSTZ2, when overexpressed, improved germination under low-temperature conditions and low-temperature tolerance at the seedling stage. In this study, qLTG-3 and qLTG-12 effectively improved cold tolerance at the seedling stage. Therefore, it is likely that the qLTG-12 locus detected in this study was due to the effect of OsGSTZ2.

In this study, a QTL (qLTG-12) was detected on chromosome 12 with a reciprocal effect with the main QTL (qTLG-3.1) for cold tolerance at germination. The double locus improved cold tolerance at the germination stage by 22–27%, overall. Also, qLTG-12 improved the cold tolerance of rice at the germination stage and the seedling stage. The results of this study were therefore well suited to molecular marker-assisted breeding and provided the materials and the molecular markers for the germination traits at low temperatures.

Conclusions

We used recombinant inbred lines LTH and H9366 to locate the quantitative trait loci related to low-temperature germination of rice. The two QTL sites (qLTG-3 and qLTG-12) were located on chromosomes 3 and 12. QTL accounted for 21.3-25.9% of the phenotypic variations at three germination time points within the two QTLs. The reciprocal effects between the two QTL sites were also detected. In this population, double QTLs increased the germination rate by 22 –27%. Moreover, qLTG-12 enhanced cold tolerance at the germination stage and the seedling stage. This study provided the materials and molecular markers for molecular marker-assisted breeding of cold tolerance.