Next Article in Journal
Development of a Novel Tissue Blot Hybridization Chain Reaction for the Identification of Plant Viruses
Next Article in Special Issue
Wheat Omics: Advancements and Opportunities
Previous Article in Journal
Extraction Techniques and Analytical Methods for Isolation and Characterization of Lignans
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 17
10.3390/plants11172324
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of QTL under Brassinosteroid-Combined Cold Treatment at Seedling Stage in Rice Using Genotyping-by-Sequencing (GBS)

by
Zhifu Guo
1,†,
Jialu Yao
1,†,
Yishan Cheng
1,
Wenzhong Zhang
2,
Zhengjin Xu
2,
Maomao Li
3,
Jing Huang
4,
Dianrong Ma
2,* and
Minghui Zhao
2,*
1
Key Laboratory of Agricultural Biotechnology of Liaoning Province, College of Biosciences and Biotechnology, Shenyang Agricultural University, Shenyang 110866, China
2
Rice Research Institute, College of Agronomy, Shenyang Agricultural University, Shenyang 110866, China
3
Rice Research Institute, Jiangxi Academy of Agricultural Sciences, Nanchang 330200, China
4
Department of Agronomy, College of Agriculture, Purdue University, West Lafayette, IN 47907, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Plants 2022, 11(17), 2324; https://doi.org/10.3390/plants11172324
Submission received: 6 August 2022 / Revised: 21 August 2022 / Accepted: 1 September 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Applications and Challenges of ‘Omics’ Technologies for Generation of Climate-Resilient Crop Varieties II)
Browse Figures
Review Reports Versions Notes

Abstract

:
Cold stress is a major threat to the sustainability of rice yield. Brassinosteroids (BR) application can enhance cold tolerance in rice. However, the regulatory mechanism related to cold tolerance and the BR signaling pathway in rice has not been clarified. In the current study, the seedling shoot length (SSL), seedling root length (SRL), seedling dry weight (SDW), and seedling wet weight (SWW) were used as the indices for identifying cold tolerance under cold stress and BR-combined cold treatment in a backcross recombinant inbred lines (BRIL) population. According to the phenotypic characterization for cold tolerance and a high-resolution SNP genetic map obtained from the GBS technique, a total of 114 QTLs were identified, of which 27 QTLs were detected under cold stress and 87 QTLs under BR-combined cold treatment. Among them, the intervals of many QTLs were coincident under different treatments, as well as different traits. A total of 13 candidate genes associated with cold tolerance or BR pathway, such as BRASSINAZOLE RESISTANT1 (OsBZR1), OsWRKY77, AP2 domain-containing protein, zinc finger proteins, basic helix-loop-helix (bHLH) protein, and auxin-induced protein, were predicted. Among these, the expression levels of 10 candidate genes were identified under different treatments in the parents and representative BRIL individuals. These results were helpful in understanding the regulation relationship between cold tolerance and BR pathway in rice.

1. Introduction

Rice (Oryza sativa L.) is one of the most important food crops that grows in tropical to temperate regions worldwide [1]. Low temperature is a severe environmental restriction that most strongly impacts rice growth and development, especially at the early seedling and reproductive stages [2]. Cold stress at the early vegetative stage can result in stunted growth and increased seedling mortality, which in turn, leads to uneven seedling stand establishment, delayed panicle development, spikelet sterility, and subsequently, decreased yield of rice [3,4]. Therefore, it is a crucial object in rice cultivation and breeding to improve cold tolerance at seedling stage.
Cold tolerance of rice is a complex quantitative trait controlled by multiple quantitative trait loci (QTL), which is directly related to a large amount of physiological and biochemical processes and environmental factors [3,4]. To date, numerous backcross inbred lines, recombinant inbred lines, and near-isogenic lines have been developed using many cold tolerance varieties as donors, and a massive number of cold tolerance-related QTLs have been identified on 12 chromosomes at different development stages in rice [5,6]. However, these are not sufficient to mine the useful genes and reveal the molecular mechanisms of cold tolerance. In recent years, the methods to obtain a single nucleotide polymorphism (SNP) molecular marker based on the high-throughput sequencing technologies, such as genotyping-by-sequencing (GBS), genome-wide association studies (GWAS), and bulked segregant RNA-seq (BSR), are increasingly used to identify QTLs across the whole genome [7,8]. In particular, as GBS technology is rapid, simple, and easy to implement for QTL mapping, genomic selection, and genetic diversity, it has been successfully applied in QTL mapping for cold tolerance in rice [9,10].
Dongxiang wild rice (Oryza rufpogon Grif, DXWR) is a relative ancestor of the cultivated rice, which can be used as a donor of novel and favorable alleles for rice breeding [11,12,13]. DXWR has a high tolerance to cold stress at all growth stages, and its underground stem can tolerate temperatures up to −12.8 °C [14]. Many QTLs associated with cold tolerance have been identified in DXWR [15,16,17,18,19]. Therefore, DXWR is an elite material for enhancing the cold tolerance of rice by hybridization, backcrossing or genetic transformation.
Brassinosteroids (BR) are important steroid hormones, which play vital roles in the growth, development, and tolerance to various stresses in plants. The exogenous BR application can enhance low-temperature tolerance in plants [20,21,22]. It was shown that BRASSINOSTEROID INSENSITIVE 1 (BRI1) and BRASSINAZOLE RESISTANT1 (BZR1) positively regulate cold tolerance, whereas BRASSINOSTEROID INSENSITIVE 2 (BIN2) negatively regulates cold tolerance in Arabidopsis thaliana. In addition, several genes related to BR biosynthesis, BR6ox2, DWF4, and CPD, are rapidly downregulated under low-temperature stress [23,24,25,26]. In rice, the treatments with BR can also effectively improve cold tolerance at germination and seedling stages [27,28,29,30]. To date, however, few genes associated with both cold tolerance and BR treatments in rice have been studied. The regulatory relationship between cold tolerance and BR remain largely unclear.
In the current study, a backcross recombinant inbred lines (BRIL) population of 140 individuals, derived from a cross between the DXWR and a super rice SN265, was constructed. Based on this BRIL, a high-resolution genome-wide SNP genetic map was developed using the GBS technique. Furthermore, we identified some QTLs and candidate genes under cold stress and BR-combined cold treatment at the seedling stage. The results have an important significance for the discovery of the genes involved in cold tolerance and for understanding the molecular mechanism that regulates the cold tolerance related to the BR signaling pathway in rice.

2. Results

2.1. Phenotypic Characterization under Cold Stress and BR-Combined Cold Treatment

Cold tolerance of each BRIL individual was evaluated under cold stress and BR-combined cold treatment at the seedling stage. The BRIL population showed varying levels of cold tolerance and significant contrasting response in SSL, SRL, SDW, and SWW (Figure 1). SSL was in the range of 20–40 cm and the majority of BRIL individuals were distributed in the range of 30–35 cm under normal temperature. Meanwhile, SRL was in the range of 3–9 cm and the majority of BRIL individuals were concentrated in the range of 4–7 cm. After cold stress and BR-combined cold treatment, the SSL of the majority of BRIL individuals was significantly reduced. SSL was in the range of 10–35 cm and the majority of BRIL individuals were distributed in the range of 20–30 cm. Compared to cold stress, the number of BRIL individuals with SSL in the range of 25–30 cm significantly increased under the BR-combined cold treatment (Figure 1b). After cold stress and BR-combined cold treatment, the SRL of the majority of BRIL individuals also reduced; however, it was less pronounced than the SSL. The number of BRIL individuals with SRL in the range of 4–5 cm, 6–7 cm, and 8–9 cm decreased, respectively. In addition, the SRL difference between cold stress and BR-combined cold treatment was not significant (Figure 1c).
For SDW and SWW, both cold stress and BR-combined cold treatment resulted in a significant decrease in SDW and SWW of the majority of BRIL individuals. The change in the number of BRILs with SDW in the range of 0.01–0.03 g and SWW in the range of 0.1–0.2 g were remarkable, respectively. Compared to cold stress, the number of BRILs with SDW in the range of 0.02–0.03 g and SWW in the range of 0.15–0.2 g increased under the BR-combined cold treatment, respectively (Figure 1d,e).
The correlation among all four traits with different treatments in the BRIL population was analyzed. The results revealed that all traits showed continuous and approximately normal distributions. Most of the traits showed a significant positive correlation regardless of cold stress or BR-combined cold treatment (Figure 2). For example, the significant positive correlations were observed between SSL and SRL (r = 1.00, p ≤ 0.01), SSL and SDW (r = 0.29, p ≤ 0.01), SSL and SWW (r = 0.31, p ≤ 0.01), SRL and SDW (r = 0.29, p ≤ 0.01), SRL and SWW (r = 0.31, p ≤ 0.01), and SDW and SWW (r = 0.85, p ≤ 0.01) under cold stress, respectively. Under BR-combined cold treatment, the positive correlations between SSL and SRL (r = 1.00, p ≤ 0.01) as well as between SDW and SWW (r = 0.57, p ≤ 0.01) were shown similar to under cold stress. Additionally, under both cold stress and BR-combined cold treatment, all of the SSL, SRL, SDW, and SWW showed highly significant positive correlations with r = 0.19 (p ≤ 0.1), r = 0.19 (p ≤ 0.1), r = 0.43 (p ≤ 0.01), and r = 0.76 (p ≤ 0.01), respectively. These results suggested that all four traits showed continuous single-peak pattern distributions, which met the requirements for QTL mapping. Moreover, the treatment of BR improved the SSL, SDW, and SWW of some BRILs under cold stress.

2.2. Construction of the Linkage Maps

Through the GBS approach, a total of 64.48 Gb high-quality filtered reads were filtered and 96.19% of these reads were mapped to the rice reference genome. Furthermore, a total of 10,836 SNPs were validated for the determination of recombinant events. The ratio of Q20 for BRILs was above 90% and the guanine-cytosine (GC) content was 43.55%, thus the quality of the data met the requirements for further analysis. A total of 10,836 unfiltered SNPs were validated for the determination of recombinant events. Finally, after filtration, followed by the ABH plugin, a total of 1145 bin markers were retained and used in the linkage map construction using the polymorphic markers to map all 12 rice chromosomes. All 12 linkage groups varied in the number and density of markers, with a total of 3188.33 cM in genetic distance and the individual linkage groups ranged from 186.74 to 346.46 cM in length. The interval of genetic distance of each polymorphic marker ranged from 1.9 to 146.3 cM and the average genetic distance between markers was 0.3 cM. The highest number of markers was 117 on chromosome 4, while the lowest number was 77 on chromosome 7 (Table 1).

2.3. Identification of QTLs for Cold Tolerance

The genetic linkage map was used to map the QTLs related to cold tolerance based on SSL, SRL, SDW, and SWW, respectively. A total of 114 QTLs were identified, of which 27 QTLs were detected under cold stress and 87 QTLs under BR-combined cold treatment. The maximum LOD score was 11.4737, while the maximum phenotypic variation explained (PVE) value was 18.2439% (Figure 3).
For SSL, a total of 9 QTLs with LOD scores ranging from 3.1103 to 5.0857 and PVE values ranging from 4.1682% to 10.7271% were identified on chromosomes 1, 2, 3, 4, 7, and 11, respectively. Among them, 3 QTLs (qSSL1-2, qSSL3-1, and qSSL11-1) were localized in cold stress. Under BR-combined cold treatment, 4 QTLs (qSSL1-1, qSSL2-1, qSSL4-1, and qSSL11-2) were identified using SSL under normal temperature condition as control, while 2 QTLs (qSSL2-2 and qSSL7-1) were localized using SSL under cold stress as control (Table 2 and Figure S1).
For SRL, only 2 QTLs (qSRL10-2 and qSRL11-2) were identified on chromosome 10 and 11 under cold stress. Under BR-combined cold treatment, only 2 QTLs (qSRL1-3 and qSRL8-1) were localized on chromosomes 1 and 8 using SRL under normal temperature condition as control, while 29 QTLs were localized on all chromosomes except for chromosome 8 with the control of SRL under cold stress. The LOD scores of all QTLs ranged from 3.0047 to 4.3086 and the PVE values ranged from 0.7550% to 9.4108% (Table 3 and Figure S2).
For SDW, a total of 49 QTLs with LOD scores ranging from 3.0094 to 11.4737 and PVE values ranging from 0.4007% to 6.5896% were localized on all chromosomes except for chromosome 7. Among these, 19 QTLs were localized under cold stress. Under BR-combined cold treatment, 21 QTLs were localized using SDW under normal temperature condition as control, while 9 QTLs were localized using SDW under normal temperature condition as control. Interestingly, the intervals of many QTLs under different treatments were coincident. For example, there were 12 pairs of QTLs between cold stress and BR-combined cold treatment with the control of SDW under normal temperature condition, such as qSDW1-3 and qSDW1-4, qSDW6-1 and qSDW6-2, qSDW6-11 and qSDW6-12 etc. Similarly, there were two pairs of QTLs between BR-combined cold treatment with the control of SDW under normal temperature condition and those using SDW under cold stress as control, including qSDW2-1 and qSDW2-2, as well as qSDW6-9 and qSDW6-10, respectively. In addition, even intervals of 3 QTLs under three different treatments were coincident, including qSDW6-1, qSDW6-2, and qSDW6-3, as well as qSDW6-4, qSDW6-5, and qSDW6-6, respectively (Table 4 and Figure S3).
For SWW, a total of 23 QTLs with LOD scores ranging from 3.1392 to 8.4626 and PVE values ranging from 0.7134% to 18.2439% were identified. Among them, only 3 QTLs (qSWW1-3, qSWW7-1, and qSWW12-2) were localized under cold stress. Under BR-combined cold treatment, 5 QTLs were localized with the control of SWW under normal temperature condition, while 15 QTLs were localized with the control of SWW under cold stress. Similar to SDW, the intervals between several pairs of QTLs under cold stress and BR-combined cold treatment were also coincident, including qSWW6-1 and qSWW6-2, qSWW6-3 and qSWW6-4, and qSWW6-6 and qSWW6-7, respectively (Table 5 and Figure S4).
Interestingly, intervals between some QTLs related to different traits were also coincident, including qSRL6-2 and qSDW6-4/qSDW6-6, qSRL6-1 and qSDW6-2, qSWW1-1/qSWW1-2 and qSDW1-1/qSDW1-3, qSWW6-2 and qSDW6-5/qSDW6-6, qSWW6-3 and qSDW6-7, qSWW6-6 and qSDW6-13, qSWW5-1 and qSDW5-1, qSWW12-4 and SDW12-5, qSWW9-1 and qSDW9-5, qSWW10-1 and SDW10-1, and qSDW11-1 and qSWW11-1, respectively.

2.4. Prediction of Candidate Genes

The QTL regions with coincident intervals and higher LOD (more than 4.0) were used to predict candidate genes based on the Rice Genome Annotation Project website. A total of 156 candidate genes were predicted, of which 65 were annotated with known function, while the other genes were identified as retrotransposon and transposon proteins (40), expressed proteins (44), and hypothetical proteins (7), respectively (Table S2).
Among the 65 annotated genes, 13 were possibly involved in cold resistance or BR signal transduction pathway (Table 6). They included several transcription factors, including OsWRKY77 (LOC_Os01g40260), AP2 domain-containing protein (OC_Os06g44750), MYB (LOC_Os11g45740 and LOC_Os05g51160), zinc finger proteins (C3HC4 type, LOC_Os05g39380; C2H2 type, LOC_Os09g27650; LSD1 type, LOC_Os12g41700), basic helix-loop-helix (bHLH) transcription factor (LOC_Os07g08440), and BRASSINAZOLE RESISTANT1 (OsBZR1, LOC_Os07g05805), and some other proteins, such as auxin-induced protein (LOC_Os06g01966), oxygen evolving enhancer protein (LOC_Os07g01480), and heat shock protein (LOC_Os12g41820). It should be mentioned that BZR1 is a critically important transcription factor in BR signal transduction pathway as well as cold response pathway in plants.

2.5. Validation of Candidate Genes

To further validate functions of the candidate genes, qRT-PCR was used to determine the expression levels of the 10 candidate genes under cold stress and BR-combined cold treatment in the parents and representative BRIL individuals (three excellent cold tolerance and three non-excellent cold tolerance).
The expression levels of all transcription factor genes, including LOC_Os01g40260 (OsWRKY77), LOC_Os11g45740 (MYB), LOC_Os07g05805 (OsBZR1), LOC_Os07g08440 (bHLH), LOC_Os06g44750 (AP2), and LOC_Os05g39380 (zinc finger), were higher under both cold stress and BR-combined cold treatment than those under normal temperature in all of the materials mentioned above. For different treatments, the expression levels of LOC_Os07g05805 (OsBZR1) and LOC_Os07g08440 (bHLH) under BR-combined cold treatment were higher than those under cold stress. For the maternal and paternal parents, the expression levels of LOC_Os11g45740 (MYB) and LOC_Os07g08440 (bHLH) were higher in DXWR than in SN265 under both cold stress and BR-combined cold treatment, whereas the opposite expression pattern was monitored for LOC_Os07g05805 (OsBZR1). For the representative BRIL individuals, the expression levels of LOC_Os11g45740 (MYB), LOC_Os07g05805 (OsBZR1), and LOC_Os07g08440 (bHLH) in three excellent BRIL individuals were higher than those in three non-excellent BRIL individuals under cold stress as well as under BR-combined cold treatment, whereas the opposite expression pattern was detected for LOC_Os01g40260 (OsWRKY77) (Figure 4). In addition to the candidate genes mentioned above, the expression levels of other genes LOC_Os06g01966 (auxin-induced protein), LOC_Os12g41820 (heat shock protein), LOC_Os07g01480 (oxygen evolving enhancer protein), and LOC_Os09g27660 (F-box protein) were not affected by cold stress and BR-combined cold treatment (Figure S5).

3. Discussion

The rich genetic diversity is an important basis for improving the agronomic trait and stresses tolerance in plants. Previous studies regarding evolution and genetic divergence showed that cultivated rice was domesticated from wild rice species. However, during the domestication process with thousands of years of evolution and natural selection, the diversity of many morphological traits was reduced, and many gene resources were lost in cultivated rice [14,31,32,33]. As the relative ancestor of the cultivated rice, common wild rice has high diversity of desirable genes, which are very important genetic resources for improving the agronomic trait and stresses tolerance in cultivated rice [34,35,36]. DXWR (O. rufipogon) was discovered in Dongxiang County, Jiangxi Province, China, which has the northernmost distribution in all species of wild rice (28°14′ N). Since DWXR plays important roles in basic rice research and industrial development, it is known as “the panda of wild plants” [17,18,37]. DXWR possesses abundant genetic resources associated with wide cross-compatibility, fertility restoration, cytoplasmic male sterility, high grain yield, and resistance to a range of biotic and abiotic stresses, especially low-temperature tolerance.
Although wild rice and cultivated rice have the same genome type, they exhibit many differences in their genome sequences. If the temporary genetic populations, such as F2 or BC1 populations and the permanent primary genetic populations, such as recombinant inbred line were constructed and used for QTL mapping, they could show not only a low QTL detection efficiency, but also poor stability. Therefore, a stable and reliable mapping population is essential for utilizing wild rice [32,33,38]. To date, multiple types of mapping populations have been constructed using DXWR as paternal parent, while a large number of QTLs associated with cold tolerance at the germination, seedling, booting, and flowering stages have been mapped to all chromosomes [16,17,18,19,30]. In this study, the female parent SN265 was used to backcross for an additional four times after a cross between the female parent SN265 and female parent DXWR (BC4F1). Furthermore, the BC4F1 population was self-crossed for eight times to obtain advanced BRIL populations (BC4F8). The back- and self-crossing for multiple generations reduced the genetic differences caused by different genetic bases between cultivated rice and wild rice, which made the result of QTL mapping more stable and reliable.
It is well known that the low-temperature tolerance of plants is a complicated quantitative trait controlled by polygenes. To cope with low-temperature stress, plants modify their physiology, metabolism, and growth by means of signal transduction and expression regulation of many genes associated with cold tolerance when plants are exposed to low-temperature stress [5]. A substantial number of genes that facilitate cold signaling and control the expression of cold regulons have been identified in many plants. A cascade signaling pathway, ICE–CBF–COR, is the most intensively studied, and is thought to be pretty important for cold tolerance in plants. This pathway contains the core components CBF/DREB (C-repeat binding factors/dehydration-responsive element-binding proteins) transcriptional factor, ICE (inducer of CBF expression) activated factor, and diverse downstream functional proteins called cold-regulated (COR) proteins. The CBFs belong to the AP2/ERF (apetala 2/ethylene response factor) family of transcription factors, which are the most important in the ICE–CBF–COR pathway. The CBFs control the expression of COR genes in response to cold stress, while the ICE acts as a positive upstream regulator of CBF genes [39,40,41].
In rice, 10 CBF/DREB homologous genes (OsDREB1A-OsDREB1J) have been identified. Among them, OsDREB1A, OsDREB1B, OsDREB1D, OsDREB1F, and OsDREB1G enhanced cold tolerance in Arabidopsis or rice [42,43]. In the case of QTL mapping studies, OsDREB1G and OsDREB1J were identified as the candidate genes related to cold tolerance for three QTLs, qLOP2, qPSR2-1 and qSR8-3, using the DXWR populations as the donor parent, respectively [16,30]. In addition, an ICE1-like gene, OrbHLH001, was isolated from DXWR, which enhanced cold tolerance when expressed in transgenic Arabidopsis [44]. Apart from these core components, many other activators or repressors are related to this pathway that directly or indirectly affect the low-temperature tolerance of plants. A MYB transcriptional factor, MYB15, interacted with ICE1 and was bound to MYB recognition sequences in the promoters of CBF genes, which repressed the expression of CBF genes. ICE1 activated the transcription of CBF3 by binding to MYC recognition elements in the promoter [39,45]. OsbHLH002, a homolog of ICE1 in rice, positively regulated cold tolerance by promoting the expression of OsTPP1, which encoded a key enzyme for trehalose biosynthesis [46]. Two F-box proteins, EIN3-binding F-box 1/2 (EBF1/2), positively regulated the expression of CBF genes by regulating the degradation of EIN3 and PIF3 [47,48]. In the current study, we identified some genes similar to those described above, such as bHLH transcription factor (LOC_Os07g08440), F-box and another domain-containing protein (LOC_Os09g27660 and LOC_Os10g04590), MYB family transcription factor (LOC_Os11g45740), and AP2 domain-containing protein (LOC_Os06g44750). The transcript levels of these genes except for LOC_Os09g27660 and LOC_Os10g04590 (F-box) were upregulated both under cold stress and BR-combined cold treatment. In the representative BRILs, LOC_Os07g08440 (bHLH), LOC_Os11g45740 (MYB), and LOC_Os06g44750 (AP2) showed higher transcript levels in the excellent BRILs than the non-excellent BRILs under cold stress and BR-combined cold treatment.
Previous studies have shown that the exogenous application of BR can improve low-temperature tolerance at different growth stages in various plants, including rice [27,28], maize [49], wheat [50], winter rye [51], and Arabidopsis [24,26], etc. To date, the molecular mechanisms regarding both cold stress and BR signaling pathway were only intensively studied in Arabidopsis. Moreover, the cold tolerance in BR signaling pathway was involved in ICE–CBF–COR signaling pathway. Under cold stress, BR directed a bHLH transcription factor CESTA (CES) to regulate the expression of the CBF and downstream COR genes in Arabidopsis [45]. BZR1 was a very important transcription factor in BR signaling pathway, acting upstream of CBF1 and CBF2 to directly enhance low-temperature tolerance in Arabidopsis. Interestingly, OsBZR1 (LOC_Os07g05805) was predicted for qSSL7-1 under BR-combined cold treatment in the current study. The expression level of OsBZR1 was upregulated both under cold stress and BR-combined cold treatment, and it showed higher transcript levels in the excellent BRILs than the non-excellent BRILs under cold stress and BR-combined cold treatment. Moreover, BZR1 regulated other genes uncoupled with CBFs, such as WKRY6, PYL6, and SOC1, etc., to improve cold tolerance in Arabidopsis [24,26]. As a negative modulator in BR signaling pathway, the protein kinase BIN2 interacted with and phosphorylated ICE1 under cold stress, which facilitated the interaction between ICE1 and the E3 ubiquitin ligase HIGH EXPRESSION OF OSMOTICALLY RESPONSIVE GENE1 (HOS1), and thereby promoted ICE1 degradation. It was suggested that BIN2 mainly downregulated ICE1 abundance when the expression levels of CBF genes were attenuated [26]. In this study, OsWRKY77 (LOC_Os01g40260), a WRKY transcription factor, was also identified for qSSL1-2 under BR-combined cold treatment. It was shown that the expression level of OsWRKY77 was upregulated in the sensitive genotype, suggesting its negative role in cold tolerance at the germination stage in rice [52]. The expression levels of OsWRKY77 were upregulated both under cold stress and BR-combined cold treatment, and its expression levels in the excellent BRILs were lower than the non-excellent BRILs under cold stress and BR-combined cold treatment. In addition to the transcription factors and genes related to cold tolerance and BR signaling pathway mentioned above, some other genes related to cold tolerance were identified, such as zinc finger domain containing protein, auxin-induced protein, oxygen evolving enhancer protein, heat shock protein, etc. The precise function and regulatory mechanism of these genes still needs further investigation.
In conclusion, many QTLs and candidate genes related to cold tolerance or BR pathway were identified by mapping analysis and qRT-PCR. This study provided a basis for further mining the genes involved in low-temperature tolerance or BR signaling pathway and investigating the mechanism regulating low-temperature tolerance in rice. Further studies were needed to thoroughly investigate whether these genes are associated with both cold tolerance and BR signaling pathway in rice.

4. Materials and Methods

4.1. Plant Materials and Population Development

A super rice variety SN265 was used as the recipient parent to cross with the donor parent DXWR. F1 plants were backcrossed with the recipient parent SN265 for 4 times to construct the BC4F1 generation. Then, the BC4F1 population was used to develop a BRIL population with 140 individuals in F8 generation (BC4F8).

4.2. Phenotypic Evaluation for Cold Tolerance

Evaluation of cold tolerance at seedling stage was performed according to the previous studies with minor changes [6,29,30]. For breaking dormancy, seeds were placed in a drying oven at 50 °C for 72 h. Then, the surface of the seeds, sterilized in 1% NaClO solution, was germinated and grew in a growth chamber with 75% humidity and 12 h light/12 h dark conditions. When the seedlings grew to the third leaf stage, the seedlings were watered and sprayed with sterilized distilled water (for cold treatment) and 0.1 μmol/L BR solution (for BR-combined cold treatment) for 24 h before the cold treatment, respectively. Furthermore, the seedlings were stressed at 4 °C for 7 days, and then moved to a greenhouse at 25 °C for 7 days to allow the seedlings to resume normal growth. Cold tolerance was evaluated based on the phenotypic changes of seedling shoot length (SSL), seedling root length (SRL), seedling dry weight (SDW), and seedling wet weight (SWW), respectively. SSL, SRL, SDW, and SWW of the abovementioned seedlings after recovering for 7 days were measured. The average of the three replicates of 10 seedlings for each treatment was analyzed. For the cold treatment, the data obtained under normal temperature (25 °C) were used as control. Cold tolerance scores based on the reduction rate of SSL, SRL, SDW, and SWW were calculated as follows: Reduction rate = ((data under normal temperature−data under cold treatment)/data under normal temperature) ×100%. For BR-combined cold treatment, both data obtained under normal temperature and cold treatment were used as controls, respectively. Cold tolerance scores were based on the reduction rate of SSL, SRL, SDW, and SWW = ((data under normal temperature−data under BR-combined cold treatment)/data under normal temperature) × 100% or ((data under BR-combined cold treatment−data under cold treatment)/data under BR-combined cold treatment) × 100%, respectively.
All of the experiments for evaluating cold tolerance were repeated three times under the same conditions, and the average cold tolerance scores from three replicates were used for the QTL mapping analysis. R software 3.5.1 (R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing) was used to analyze the frequency distributions and correlations of all 4 phenotypic traits (www.R-project.org, accessed on 1 July 2022).

4.3. Genotyping-by-Sequencing and SNP Identification

Young leaves of 140 BRIL individuals and their parental lines were collected separately for genomic DNA extraction using a DNA extraction kit (Aidlab, Beijing, China). DNA libraries from BRILs were constructed using the GBS technology. The restriction enzymes PstI and MspI and T4 ligase were used for digestion and ligation, respectively. The DNA libraries were enriched and sequenced using Illumina HiSeq 4000 instrument (Huada Gene Technology Co., Ltd., Shenzhen, China). The Tassel software was used to analyze the raw data and high-quality SNPs between parents were termed by alignment with Nipponbare reference genome [53].

4.4. Construction of Linkage Map and QTL Analysis

The Illumina data of BRIL population was used to construct a linkage map using the QTL IciMapping sofware v4.163 (Meng et al. 2015, Beijing, China) [54]. The genotypic maps were aligned and split into recombination bins according to the recombination breakpoints, with the parameter of window size of 10 cm and walk speed of 1 cm.
QTLs were mapped using the inclusive composite interval mapping of the additive (ICIM-ADD) mapping method. QTLs were computed by a permutation test involving 1000 runs at a significance level of p = 0.05. The threshold for the logarithm of odds (LOD) scores was set to 3.0. QTL nomenclature was followed by the method of McCouch et al. (2008) [55].

4.5. Prediction and Validation of Candidate Genes

To predict potential candidate genes within QTL intervals, the physical positions of SNP markers flanking the QTLs were searched in the Rice Genome Annotation Project website (http://rice.plantbiology.msu.edu, accessed on 1 July 2022). The seedlings at the third leaf stage from 2 parental lines, 3 excellent cold tolerance BRIL individuals, and 3 non-excellent cold tolerance BRIL individuals at normal temperature (25 °C), for 12 h under cold stress (4 °C), and for 12 h under BR-combined cold treatment (pre-watering and pre-spraying with 0.1 μmol/L BR solution for 24 h, 4 °C) were collected to extract total RNA using the Trizol reagent. Quantitative real-time PCR (qRT-PCR) was used to investigate the expression patterns of 10 representative candidate genes. Moreover, qRT-PCR was carried out using a q225 Real-Time PCR System (Monad, China) under the following conditions: 95 °C for 5 min, 30 cycles of 95 °C for 15 s, 60 °C for 30 s, and 72 °C for 1 min, and 1 cycle of 72 °C for 15 min. Each sample was analyzed in three biological and three technical replicates. The relative expression levels were calculated via the ∆Ct method [56]. The actin gene of rice (LOC_Os03g50885) was used as a reference gene. All of the primer sequences are listed in Table S1.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11172324/s1. Figure S1: QTL mapping based on seedling shoot length (SSL); Figure S2: QTL mapping based on seedling root length (SRL); Figure S3: QTL mapping based on seedling dry weight (SDW); Figure S4: QTL mapping based on seedling wet weight (SWW); Figure S5: Expression patterns of some candidate genes (LOC_Os12g41820, LOC_Os06g01966, LOC_Os07g01480, and LOC_Os09g27660); Table S1: Details of primers used for this study; Table S2: Summary of all candidate genes.

Author Contributions

Z.G., D.M. and M.Z. designed the experiments and wrote the manuscript. M.Z. and M.L. constructed the BRIL population. Z.G., Y.C. and J.Y. performed the experiments. D.M. and W.Z. performed the candidate gene association. J.H., D.M. and Z.X. contributed to the revision of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Liaoning Revitalization Talents Program (XLYC2008025, M.Z.), the Natural Science Foundation of Liaoning Province (2019JH3/10300327, Z.G.), the Liaoning BaiQianWan Talents Program (2020921092, Z.G.), and the Basic Research Project of Education Department of Liaoning Province (LSNJC202022, Z.G.).

Data Availability Statement

The raw data supporting the conclusions of this manuscript can be found in the National Center for Biotechnology Information (NCBI) BioProject, https://www.ncbi.nlm.nih.gov/bioproject (accessed on 10 May 2022), PRJNA836720.

Acknowledgments

We express our gratitude to the reviewers for their helpful comments to improve the manuscript. We thank Ao Zhang from Shenyang Agricultural University for help with handling the data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cheng, S.H.; Cao, L.Y.; Zhuang, J.Y.; Chen, S.G.; Zhan, X.D.; Fan, Y.Y.; Min, S.K. Super hybrid rice breeding in China: Achievements and prospects. J. Integr. Plant Biol. 2007, 49, 805–810. [Google Scholar] [CrossRef]
  2. Fujino, K.; Sekiguchi, H.; Sato, T.; Kiuchi, H.; Nonoue, Y.; Takeuchi, Y.; Ando, T.; Lin, S.Y.; Yano, M. Mapping of Quantitative Trait Loci Controlling Low-Temperature Germinability in Rice (Oryza sativa L.). Theor. Appl. Genet. 2004, 108, 794–799. [Google Scholar] [CrossRef]
  3. Nakagahra, M.; Okuno, K.; Vaughan, D. Rice Genetic Resources: History, Conservation, Investigative Characterization and Use in Japan. Plant Mol. Biol. 1997, 35, 69–77. [Google Scholar] [CrossRef] [PubMed]
  4. Najeeb, S.; Ali, J.; Mahender, A.; Pang, Y.L.; Zilhas, J.; Murugaiyan, V.; Vemireddy, L.R.; Li, Z. Identification of Main-Effect Quantitative Trait Loci (QTLs) for Low-Temperature Stress Tolerance Germination-and Early Seedling Vigor-Related Traits in Rice (Oryza sativa L.). Mol. Breed. 2020, 40, 10. [Google Scholar] [CrossRef] [PubMed]
  5. Miura, K.; Lin, S.Y.; Yano, M.; Nagamine, T. Mapping Quantitative Trait Loci Controlling Low Temperature Germinability in Rice (Oryza sativa L.). Breed. Sci. 2001, 51, 293–299. [Google Scholar] [CrossRef]
  6. Yang, L.; Wang, J.; Han, Z.; Lei, L.; Liu, H.L.; Zheng, H.; Xin, W.; Zou, D. Combining QTL-seq and Linkage Mapping to Fine Map a Candidate Gene in qCTS6 for Cold Tolerance at the Seedling Stage in Rice. BMC Plant Biol. 2021, 21, 278. [Google Scholar] [CrossRef]
  7. Mielczarek, M.; Szyda, J. Review of Alignment and Snp Calling Algorithms for Next-Generation Sequencing Data. J. Appl. Genet. 2016, 57, 71–79. [Google Scholar] [CrossRef] [PubMed]
  8. Jiang, Z.; Wang, H.; Michal, J.J.; Zhou, X.; Liu, B.; Woods, L.C.; Fuchs, R.A. Genome Wide Sampling Sequencing for Snp Genotyping: Methods, Challenges and Future Development. Int. J. Biol. Sci. 2016, 12, 100–108. [Google Scholar] [CrossRef]
  9. Spindel, J.; Wright, M.; Chen, C.; Cobb, J.; Gage, J.; Harrington, S.; Lorieux, M.; Ahmadi, N.; McCouch, S. Bridging the Genotyping Gap: Using Genotyping by Sequencing (GBS) to Add High-Density Snp Markers and New Value to Traditional Bi-Parental Mapping and Breeding Populations. Theor. Appl. Genet. 2013, 126, 2699–2716. [Google Scholar] [CrossRef]
  10. Arbelaez, J.D.; Moreno, L.T.; Singh, N.; Tung, C.W.; Maron, L.G.; Ospina, Y.; Martinez, C.P.; Grenier, C.; Lorieux, M.; McCouch, S. Development and GBS-Genotyping of Introgression Lines (ILS) Using Two Wild Species of Rice, O. meridionalis and O. rufipogon, in a Common Recurrent Parent, O. sativa cv. Curinga. Mol. Breed. 2015, 35, 81. [Google Scholar] [CrossRef] [Green Version]
  11. Atwell, B.J.; Wang, H.; Scafaro, A.P. Could Abiotic Stress Tolerance in Wild Relatives of Rice Be Used to Improve Oryza Sativa? Plant Sci. 2014, 215–216, 48–58. [Google Scholar] [CrossRef]
  12. Zhang, K.; Su, J.; Xu, M.; Zhou, Z.; Zhu, X.; Ma, X.; Hou, J.; Tan, L.; Zhu, Z.; Cai, H.; et al. A Common Wild Rice-Derived Boc1 Allele Reduces Callus Browning in Indica Rice Transformation. Nat. Commun. 2020, 11, 443. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, Q.; Chen, Q.; Wang, S.; Hong, Y.; Wang, Z. Rice and Cold Stress: Methods for Its Evaluation and Summary of Cold Tolerance-Related Quantitative Trait Loci. Rice 2014, 7, 24. [Google Scholar] [CrossRef] [PubMed]
  14. Liu, F.; Sun, C.; Tan, L.; Fu, Y.; Li, D.; Wang, X. Identification and Mapping of Quantitative Trait Loci Controlling Cold-Tolerance of Chinese Common Wild Rice (O. rufipogon Griff.) at Booting to Flowering Stages. Chin. Sci. Bull. 2003, 48, 2068–2071. [Google Scholar] [CrossRef]
  15. Tian, F.; Li, D.J.; Fu, Q.; Zhu, Z.F.; Fu, Y.C.; Wang, X.K.; Sun, C.Q. Construction of Introgression Lines Carrying Wild Rice (Oryza Rufipogon Griff.) Segments in Cultivated Rice (Oryza sativa L.) Background and Characterization of Introgressed Segments Associated with Yield-Related Traits. Theor. Appl. Genet. 2006, 112, 570–580. [Google Scholar] [CrossRef]
  16. Xiao, N.; Huang, W.N.; Li, A.H.; Gao, Y.; Li, Y.H.; Pan, C.H.; Ji, H.; Zhang, X.X.; Dai, Y.; Dai, Z.Y.; et al. Fine Mapping of the qLOP2 and qPSR2-1 Loci Associated with Chilling Stress Tolerance of Wild Rice Seedlings. Theor. Appl. Genet. 2015, 128, 173–185. [Google Scholar] [CrossRef]
  17. Liang, Y.; Yan, C.; Zheng, J.; Nan, W.; Qin, X.; Zhang, H. Locating Qtl Associated with Spike Traits of Dongxiang Wild Rice (Oryza rufipogon Griff.). Euphytica 2019, 215, 26. [Google Scholar] [CrossRef]
  18. Liang, Y.; Zheng, J.; Yan, C.; Li, X.; Liu, S.; Zhou, J.; Qin, X.; Nan, W.; Yang, Y.; Zhang, H. Locating Qtls Controlling Overwintering Trait in Chinese Perennial Dongxiang Wild Rice. Mol. Genet. Genom. 2018, 293, 81–93. [Google Scholar] [CrossRef]
  19. Bai, L.W.; Liu, J.; Dai, L.F.; Deng, Q.W.; Chen, Y.L.; Xie, J.K.; Luo, X.D. Identification and Characterisation of Cold Stress-Related Proteins in Oryza rufipogon at the Seedling Stage Using Label-Free Quantitative Proteomic Analysis. Funct. Plant Biol. 2021, 48, 542–555. [Google Scholar] [CrossRef]
  20. Sharma, I.; Kaur, N.; Pati, P.K. Brassinosteroids: A Promising Option in Deciphering Remedial Strategies for Abiotic Stress Tolerance in Rice. Front. Plant Sci. 2017, 8, 2151. [Google Scholar] [CrossRef]
  21. Manghwar, H.; Hussain, A.; Ali, Q.; Liu, F. Brassinosteroids (BRs) Role in Plant Development and Coping with Different Stresses. Int. J. Mol. Sci. 2022, 23, 1012. [Google Scholar] [CrossRef]
  22. Basit, F.; Liu, J.; An, J.; Chen, M.; He, C.; Zhu, X.; Li, Z.; Hu, J.; Guan, Y. Brassinosteroids as a Multidimensional Regulator of Plant Physiological and Molecular Responses under Various Environmental Stresses. Environ. Sci. Pollut. Res. Int. 2021, 28, 44768–44779. [Google Scholar] [CrossRef] [PubMed]
  23. Sadura, I.; Janeczko, A. Brassinosteroids and the Tolerance of Cereals to Low and High Temperature Stress: Photosynthesis and the Physicochemical Properties of Cell Membranes. Int. J. Mol. Sci. 2021, 23, 342. [Google Scholar] [CrossRef] [PubMed]
  24. Li, H.; Ye, K.; Shi, Y.; Cheng, J.; Zhang, X.; Yang, S. BZR1 Positively Regulates Freezing Tolerance Via CBF-Dependent and CBF-Independent Pathways in Arabidopsis. Mol. Plant 2017, 10, 545–559. [Google Scholar] [CrossRef] [PubMed]
  25. Peres, A.; Soares, J.S.; Tavares, R.G.; Righetto, G.; Zullo, M.A.T.; Mandava, N.B.; Menossi, M. Brassinosteroids, the Sixth Class of Phytohormones: A Molecular View from the Discovery to Hormonal Interactions in Plant Development and Stress Adaptation. Int. J. Mol. Sci. 2019, 20, 331. [Google Scholar] [CrossRef]
  26. Ye, K.; Li, H.; Ding, Y.; Shi, Y.; Song, C.; Gong, Z.; Yang, S. Brassinosteroid-Insensitive2 Negatively Regulates the Stability of Transcription Factor Ice1 in Response to Cold Stress in Arabidopsis. Plant Cell 2019, 31, 2682–2696. [Google Scholar] [CrossRef]
  27. Hwang, O.J.; Back, K. Melatonin Deficiency Confers Tolerance to Multiple Abiotic Stresses in Rice Via Decreased Brassinosteroid Levels. Int. J. Mol. Sci. 2019, 20, 5173. [Google Scholar] [CrossRef]
  28. Wu, W.; Zhou, W.; Tang, C.; Liu, K.; Zeng, H.; Wang, Y. Effects of Exogenous 2,4-Epibrassinolide on Germination and Physiological Characteristics of Rice Seeds under Chilling Stress. Mol. Plant Breed. 2020, 18, 4427–4434. [Google Scholar] [CrossRef]
  29. Guo, H.; Wu, T.; Li, S.; He, Q.; Yang, Z.; Zhang, W.; Gan, Y.; Sun, P.; Xiang, G.; Zhang, H.; et al. The Methylation Patterns and Transcriptional Responses to Chilling Stress at the Seedling Stage in Rice. Int. J. Mol. Sci. 2019, 20, 5089. [Google Scholar] [CrossRef]
  30. Guo, Z.F.; Wang, H.T.; Cheng, Y.S.; Guo, S.Q.; Zhang, W.Z.; Xu, Z.J.; Li, M.M.; Huang, J.; Zhao, M.H. QTL Mapping Analysis for Cold Tolerance under Cold Stress and Brassinosteroids combined Cold Treatment at Germination and Bud Burst Stages in Rice. Front. Plant Sci. 2022, 13, 938339. [Google Scholar] [CrossRef]
  31. Hamzelou, S.; Kamath, K.S.; Masoomi-Aladizgeh, F.; Johnsen, M.M.; Atwell, B.J.; Haynes, P.A. Wild and Cultivated Species of Rice Have Distinctive Proteomic Responses to Drought. Int. J. Mol. Sci. 2020, 21, 5980. [Google Scholar] [CrossRef] [PubMed]
  32. Shi, C.; Li, W.; Zhang, Q.J.; Zhang, Y.; Tong, Y.; Li, K.; Liu, Y.L.; Gao, L.Z. The Draft Genome Sequence of an Upland Wild Rice Species, Oryza granulata. Sci. Data 2020, 7, 131. [Google Scholar] [CrossRef] [PubMed]
  33. Qiao, W.; Wang, Y.; Xu, R.; Yang, Z.; Sun, Y.; Su, L.; Zhang, L.; Wang, J.; Huang, J.; Zheng, X.; et al. A Functional Chromogen Gene C from Wild Rice Is Involved in a Different Anthocyanin Biosynthesis Pathway in Indica and Japonica. Theor. Appl. Genet. 2021, 134, 1531–1543. [Google Scholar] [CrossRef] [PubMed]
  34. Xing, J.; Zhang, D.; Yin, F.; Zhong, Q.; Wang, B.; Xiao, S.; Ke, X.; Wang, L.; Zhang, Y.; Zhao, C.; et al. Identification and Fine-Mapping of a New Bacterial Blight Resistance Gene, Xa47(T), in G252, an Introgression Line of Yuanjiang Common Wild Rice (Oryza rufipogon). Plant Dis. 2021, 105, 4106–4112. [Google Scholar] [CrossRef]
  35. Liu, W.; Ghouri, F.; Yu, H.; Li, X.; Yu, S.; Shahid, M.Q.; Liu, X. Genome Wide Re-Sequencing of Newly Developed Rice Lines from Common Wild Rice (Oryza rufipogon Griff.) for the Identification of Nbs-Lrr Genes. PLoS ONE 2017, 12, e0180662. [Google Scholar] [CrossRef]
  36. Yuan, R.; Zhao, N.; Usman, B.; Luo, L.; Liao, S.; Qin, Y.; Nawaz, G.; Li, R. Development of Chromosome Segment Substitution Lines (Cssls) Derived from Guangxi Wild Rice (Oryza rufipogon Griff.) under Rice (Oryza sativa L.) Background and the Identification of Qtls for Plant Architecture, Agronomic Traits and Cold Tolerance. Genes 2020, 11, 980. [Google Scholar] [CrossRef]
  37. Yun, Y.T.; Chung, C.T.; Lee, Y.J.; Na, H.J.; Lee, J.C.; Lee, S.G.; Lee, K.W.; Yoon, Y.H.; Kang, J.W.; Lee, H.S.; et al. Qtl Mapping of Grain Quality Traits Using Introgression Lines Carrying Oryza rufipogon Chromosome Segments in Japonica Rice. Rice 2016, 9, 62. [Google Scholar] [CrossRef]
  38. Cen, W.; Liu, J.; Lu, S.; Jia, P.; Yu, K.; Han, Y.; Li, R.; Luo, J. Comparative Proteomic Analysis of Qtl Cts-12 Derived from Wild Rice (Oryza rufipogon Griff.), in the Regulation of Cold Acclimation and De-Acclimation of Rice (Oryza sativa L.) in Response to Severe Chilling Stress. BMC Plant Biol. 2018, 18, 163. [Google Scholar] [CrossRef]
  39. Wang, D.Z.; Jin, Y.N.; Ding, X.H.; Wang, W.J.; Zhai, S.S.; Bai, L.P.; Guo, Z.F. Gene Regulation and Signal Transduction in the ICE-CBF-COR Signaling Pathway During Cold Stress in Plants. Biochemistry 2017, 82, 1103–1117. [Google Scholar] [CrossRef]
  40. Huang, G.T.; Ma, S.L.; Bai, L.P.; Zhang, L.; Ma, H.; Jia, P.; Liu, J.; Zhong, M.; Guo, Z.F. Signal Transduction During Cold, Salt, and Drought Stresses in Plants. Mol. Biol. Rep. 2012, 39, 969–987. [Google Scholar] [CrossRef]
  41. Hwarari, D.; Guan, Y.; Ahmad, B.; Movahedi, A.; Min, T.; Hao, Z.; Lu, Y.; Chen, J.; Yang, L. ICE-CBF-COR Signaling Cascade and Its Regulation in Plants Responding to Cold Stress. Int. J. Mol. Sci. 2022, 23, 1549. [Google Scholar] [CrossRef]
  42. Dubouzet, J.G.; Sakuma, Y.; Ito, Y.; Kasuga, M.; Dubouzet, E.G.; Miura, S.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Osdreb Genes in Rice, Oryza sativa L., Encode Transcription Activators That Function in Drought-, High-Salt- and Cold-Responsive Gene Expression. Plant J. 2003, 33, 751–763. [Google Scholar] [CrossRef] [PubMed]
  43. Huang, J.; Sun, S.; Xu, D.; Lan, H.; Sun, H.; Wang, Z.; Bao, Y.; Wang, J.; Tang, H.; Zhang, H. A Tfiiia-Type Zinc Finger Protein Confers Multiple Abiotic Stress Tolerances in Transgenic Rice (Oryza sativa L.). Plant Mol. Biol. 2012, 80, 337–350. [Google Scholar] [CrossRef]
  44. Li, F.; Guo, S.; Zhao, Y.; Chen, D.; Chong, K.; Xu, Y. Overexpression of a Homopeptide Repeat-Containing Bhlh Protein Gene (Orbhlh001) from Dongxiang Wild Rice Confers Freezing and Salt Tolerance in Transgenic Arabidopsis. Plant Cell Rep. 2010, 29, 977–986. [Google Scholar] [CrossRef] [PubMed]
  45. Eremina, M.; Unterholzner, S.J.; Rathnayake, A.I.; Castellanos, M.; Khan, M.; Kugler, K.G.; May, S.T.; Mayer, K.F.; Rozhon, W.; Poppenberger, B. Brassinosteroids Participate in the Control of Basal and Acquired Freezing Tolerance of Plants. Proc. Natl. Acad. Sci. USA 2016, 113, e5982–e5991. [Google Scholar] [CrossRef] [PubMed]
  46. Zhang, Z.; Li, J.; Li, F.; Liu, H.; Yang, W.; Chong, K.; Xu, Y. OsMAPK3 Phosphorylates OsbHLH002/OsICE1 and Inhibits Its Ubiquitination to Activate OsTPP1 and Enhances Rice Chilling Tolerance. Dev. Cell 2017, 43, 731–743. [Google Scholar] [CrossRef]
  47. Jiang, B.; Shi, Y.; Zhang, X.; Xin, X.; Qi, L.; Guo, H.; Li, J.; Yang, S. PIF3 Is a Negative Regulator of the Cbf Pathway and Freezing Tolerance in Arabidopsis. Proc. Natl. Acad. Sci. USA 2017, 114, e6695–e6702. [Google Scholar] [CrossRef]
  48. Shi, Y.; Tian, S.; Hou, L.; Huang, X.; Zhang, X.; Guo, H.; Yang, S. Ethylene Signaling Negatively Regulates Freezing Tolerance by Repressing Expression of CBF and Type-A ARR Genes in Arabidopsis. Plant Cell 2012, 24, 2578–2595. [Google Scholar] [CrossRef]
  49. Singh, I.; Kumar, U.; Singh, S.K.; Gupta, C.; Singh, M.; Kushwaha, S.R. Physiological and Biochemical Effect of 24-Epibrassinoslide on Cold Tolerance in Maize Seedlings. Physiol. Mol. Biol. Plants 2012, 18, 229–236. [Google Scholar] [CrossRef]
  50. Janeczko, A.; Pociecha, E.; Dziurka, M.; Jurczyk, B.; Libik-Konieczny, M.; Oklestkova, J.; Novák, O.; Pilarska, M.; Filek, M.; Rudolphi-Skórska, E.; et al. Changes in Content of Steroid Regulators During Cold Hardening of Winter Wheat-Steroid Physiological/Biochemical Activity and Impact on Frost Tolerance. Plant Physiol. Biochem. 2019, 139, 215–228. [Google Scholar] [CrossRef]
  51. Pociecha, E.; Dziurka, M.; Oklestkova, J.; Janeczko, A. Brassinosteroids Increase Winter Survival of Winter Rye (Secale Cereale L.) by Affecting Photosynthetic Capacity and Carbohydrate Metabolism During the Cold Acclimation Process. Plant Growth Regul. 2016, 80, 127–135. [Google Scholar] [CrossRef]
  52. Viana, V.E.; Carlos da Maia, L.; Busanello, C.; Pegoraro, C.; Costa de Oliveira, A. When Rice Gets the Chills: Comparative Transcriptome Profiling at Germination Shows Wrky Transcription Factor Responses. Plant Biol. 2021, 23, 100–112. [Google Scholar] [CrossRef] [PubMed]
  53. Kawahara, Y.; de la Bastide, M.; Hamilton, J.P.; Kanamori, H.; McCombie, W.R.; Ouyang, S.; Schwartz, D.C.; Tanaka, T.; Wu, J.; Zhou, S.; et al. Improvement of the Oryza Sativa Nipponbare Reference Genome Using Next Generation Sequence and Optical Map Data. Rice 2013, 6, 4. [Google Scholar] [CrossRef] [PubMed]
  54. Meng, L.; Li, H.H.; Zhang, L.Y.; Wang, J.K. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015, 3, 269–283. [Google Scholar] [CrossRef]
  55. McCouch, S.R. Gene Nomenclature System for Rice. Rice 2008, 1, 72–84. [Google Scholar] [CrossRef]
  56. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative Pcr and the 2-DDCt Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Plants 11 02324 g001 550
Figure 1. Phenotypic characterization under cold stress and BR-combined cold treatment at seedling stage in BRILs. (a) Phenotypes of the representative BRILs with excellent, moderate, and non-excellent cold tolerance; (b) ranges of the seedling shoot length (SSL) in BRILs; (c) ranges of the seedling root length (SRL) in BRILs; (d) ranges of the seedling dry weight (SDW) in BRILs; (e) ranges of the seedling wet weight (SWW) in BRILs. NT: Normal temperature; C: Cold stress; C + B: BR-combined cold treatment.
Figure 1. Phenotypic characterization under cold stress and BR-combined cold treatment at seedling stage in BRILs. (a) Phenotypes of the representative BRILs with excellent, moderate, and non-excellent cold tolerance; (b) ranges of the seedling shoot length (SSL) in BRILs; (c) ranges of the seedling root length (SRL) in BRILs; (d) ranges of the seedling dry weight (SDW) in BRILs; (e) ranges of the seedling wet weight (SWW) in BRILs. NT: Normal temperature; C: Cold stress; C + B: BR-combined cold treatment.
Plants 11 02324 g001
Plants 11 02324 g002 550
Figure 2. Correlation analysis of SSL, SRL, SDW, and SWW under cold stress and BR-combined cold treatment. Plots on the diagonal line show the phenotypic distribution of each trait as indicated; values above the diagonal line are Pearson’s correlation coefficients between traits; plots below the diagonal line are scatter plots of compared traits. ** p ≤ 0.01; * p ≤ 0.05; · p ≤ 0.1. SSL: Seedling shoot length; SRL: Seedling root length; SDW: Seedling dry weight; and SWW: Seedling wet weight; C: Cold stress; C + B: BR-combined cold treatment.
Figure 2. Correlation analysis of SSL, SRL, SDW, and SWW under cold stress and BR-combined cold treatment. Plots on the diagonal line show the phenotypic distribution of each trait as indicated; values above the diagonal line are Pearson’s correlation coefficients between traits; plots below the diagonal line are scatter plots of compared traits. ** p ≤ 0.01; * p ≤ 0.05; · p ≤ 0.1. SSL: Seedling shoot length; SRL: Seedling root length; SDW: Seedling dry weight; and SWW: Seedling wet weight; C: Cold stress; C + B: BR-combined cold treatment.
Plants 11 02324 g002
Plants 11 02324 g003 550
Figure 3. Molecular genetic map showing the positions of QTLs for four traits investigated under cold stress and BR-combined cold treatment.
Figure 3. Molecular genetic map showing the positions of QTLs for four traits investigated under cold stress and BR-combined cold treatment.
Plants 11 02324 g003
Plants 11 02324 g004 550
Figure 4. Expression patterns of some candidate genes. Biological triplicates were averaged and statistically analyzed via Student’s t-test (* p < 0.05; ** p < 0.01; *** p < 0.001). NT: Normal temperature; C: Cold stress; C + B: BR-combined cold treatment.
Figure 4. Expression patterns of some candidate genes. Biological triplicates were averaged and statistically analyzed via Student’s t-test (* p < 0.05; ** p < 0.01; *** p < 0.001). NT: Normal temperature; C: Cold stress; C + B: BR-combined cold treatment.
Plants 11 02324 g004
Table
Table 1. Summary of genetic linkage map characteristics in BRILs.
Table 1. Summary of genetic linkage map characteristics in BRILs.
ChrNumber of Bin MarkersGenetic Distance (cM)Ave. Genetic Distance between Markers (cM)Ave. Interval (cM)Interval Range (cM)
1104268.200.425.85.1–120.5
2112288.600.325.84.0–74.5
3116314.040.327.14.9–73.6
4117346.460.329.67.5–88.7
581229.410.428.33.6–127.8
6100293.940.429.41.9–98.6
777186.740.424.34.0–75.6
891271.540.529.86.3–146.3
979227.080.328.76.2–121.3
1085230.990.327.25.4–112.6
11105303.460.328.92.0–126.8
1278227.870.329.22.0–137.7
All11453188.330.3 (Ave.)27.8 (Ave.)1.9–146.3
Chr: Chromosome number; Ave.: Average.
Table
Table 2. QTL summary under different conditions for SSL.
Table 2. QTL summary under different conditions for SSL.
ConditionQTLChrIntervalLODPVE (%)Add
ColdqSSL1-2122,514,210–23,217,5914.74699.89180.0402
qSSL3-1311,339,312–12,047,0813.38507.7212−0.0398
qSSL11-1116,640,029–6,662,6254.81628.6890−0.0356
BR + C-NqSSL1-1112,150,755–12,275,9564.39926.3653−0.0364
qSSL2-128,052,187–8,229,9913.618410.72710.0374
qSSL4-1435,027,717–35,296,1195.08579.47700.0401
qSSL11-21116,545,655–16,884,9153.91689.1290−0.0356
BR + C-CqSSL2-2234,425,783–34,742,6713.925815.0490−0.0537
qSSL7-1723,483,968–23,485,5313.11034.16820.0259
Chr: Chromosome number; LOD: Logarithm of odds value; PVE: Phenotypic variance explained; Add: Additive effect; Cold: Cold stress; Cold + BR-N: Data under normal temperature condition were used as the control; Cold + BR-C: Data under cold treatment condition were used as the control.
Table
Table 3. QTL summary under different conditions for SRL.
Table 3. QTL summary under different conditions for SRL.
ConditionQTLChrIntervalLODPVE (%)Add
ColdqSRL10-21020,367,970–20,569,7863.01027.72000.1398
qSRL11-21116,889,698–17,163,4583.14697.6157−0.1713
BR + C-NqSRL1-3128,250,887–28,381,9103.03849.10640.1096
qSRL8-1823,572,813–23,607,9993.68789.41080.1096
BR + C-CqSRL1-119,795,818–9,824,2824.01581.6769−0.2551
qSRL1-2127,714,994–27,766,9013.11171.1912−0.1597
qSRL1-4136,461,005–37,047,6784.00280.75500.1129
qSRL1-5140,759,127–40,984,4573.03190.95150.9515
qSRL2-1218,887,339–19,398,5803.30441.8057−0.2096
qSRL3-137,023,285–7,063,0903.05001.6192−0.2303
qSRL4-1435,487,338–35,496,2663.26741.7908−0.2234
qSRL5-151,622,816–2,206,6453.17171.7711−0.2284
qSRL5-2515,157,600–16,159,8263.25481.7136−0.2395
qSRL5-3523,101,672–23,110,0604.19141.7021−0.2340
qSRL5-4529,195,759–29,365,2393.00471.7353−0.2379
qSRL6-1634,503–49,7103.23851.4471−0.2590
qSRL6-2649,710–75,0333.73941.5946−0.2537
qSRL6-362,026,942–2,573,6483.31841.6972−0.2214
qSRL6-4626,319,487–27,021,0924.11961.6335−0.2495
qSRL6-5627,021,092–27,066,7253.61841.6633−0.2488
qSRL7-172,206,752–2,342,1903.10011.25490.1573
qSRL7-2729,374,392–29,637,6533.33821.5548 −0.2212
qSRL9-192,137,063–2,395,1823.34161.6710−0.2369
qSRL9-292,395,182–2,547,9663.89061.7777−0.2301
qSRL9–3911,240,451–11,542,5983.58381.6649−0.2410
qSRL10-1105,080,068–5,747,4383.03491.6997−0.2413
qSRL11-1118,584,535–9,024,1403.47371.31060.2003
qSRL12-112978,610–3,267,4943.23281.6911−0.2408
qSRL12-21218,216,503–18,482,2343.14961.8460−0.1963
qSRL12-31218,482,234–19,002,4364.30861.6868−0.2264
qSRL12-41220,179,243–20,499,1133.15601.6292−0.2544
qSRL12-51220,499,113–20,522,9903.05511.6346−0.2470
qSRL12-61225,801,108–25,921,2384.05881.5263−0.2428
Chr: Chromosome number; LOD: Logarithm of odds value; PVE: phenotypic variance explained; Add: Additive effect; Cold: Cold stress; Cold + BR-N: Data under normal temperature condition were used as the control; Cold + BR-C: Data under cold treatment condition were used as the control.
Table
Table 4. QTL summary under different conditions for SDW.
Table 4. QTL summary under different conditions for SDW.
ConditionQTLChrIntervalLODPVE (%)Add
ColdqSDW1-118,430,459–9,795,8183.00941.6880−0.4911
qSDW1-319,795,818–9,824,2823.61081.7076−0.4813
qSDW2-3223,029,608–23,036,1933.28481.8389−0.3107
qSDW4-142,035,471–2,255,9313.66310.84530.1214
qSDW5-1528,253,886–28,304,1363.21731.8218−0.3146
qSDW6-1634,503–49,7105.27511.8124−0.3869
qSDW6-5649,710–75,0334.68021.7782−0.3912
qSDW6-8675,033–86,2573.99711.7254 −0.3676
qSDW6-116543,717–961,4034.59501.7380−0.3096
qSDW6-136918,728–8,551,0703.12241.7444−0.3801
qSDW6-15626,319,487–27,021,0923.57231.7885 −0.3040
qSDW6-18627,021,092–27,066,7253.04431.7837−0.3076
qSDW8-4822,126,649–23,227,3403.43431.9619−0.2550
qSDW8-6827,136,983–27,371,9183.55041.8118−0.3283
qSDW9-295,908,433–5,989,5123.25791.5819−0.2097
qSDW9-4911,240,451–11,542,5983.16151.8326−0.2532
qSDW12-112978,610–3,267,4943.73061.8958−0.3331
qSDW12-41220,179,243–20,499,1133.06381.8157−0.3170
qSDW12-51220,499,113–20,522,9903.10361.8231−0.3093
BR + C-NqSDW1-218,430,459–9,795,8184.11631.0370−0.4150
qSDW1-419,795,818–9,824,2829.61581.6340−1.3208
qSDW2-2214,192,829–15,359,1689.48881.6340 1.3208
qSDW6-2634,503–49,71010.49201.6340−1.3208
qSDW6-6649,710–75,03311.47371.6340−1.3208
qSDW6-7675,033–86,2579.88761.6340−1.3208
qSDW6-96536,139–543,7179.11141.6340−1.3208
qSDW6-126543,717–961,4039.19931.6340−1.3208
qSDW6-1468,551,070–9,187,28710.43301.6340−1.3208
qSDW6-16626,319,487–27,021,0923.30650.9447−0.3687
qSDW6-17627,021,092–27,066,7258.01831.6340−1.3208
qSDW8-1820,095,951–20,541,2824.64131.63401.3208
qSDW8-5827,136,983–27,371,9188.64421.6340−1.3208
qSDW9-192,395,182–2,547,9667.97561.6340−1.3208
qSDW9-395,908,433–5,989,5123.05300.4007−0.1190
qSDW9-5916,774,761–16,829,1258.60801.6340−1.3208
qSDW10-1102,188,025–2,329,7348.80751.6340−1.3208
qSDW11-21126,184,650–26,729,0258.89931.6340−1.3208
qSDW12-212978,610–3,267,4948.31681.6340−1.3208
qSDW12-31220,179,243–20,499,1138.35371.6340−1.3208
qSDW12-61220,499,113–20,522,9908.08661.6340−1.3208
BR + C-CqSDW2-1214,192,829–15,359,1689.22106.58960.9641
qSDW3-1328,456,345–29,627,2343.95771.7956−0.1035
qSDW6-3634,503–49,7106.68534.5867−0.8043
qSDW6-4649,710–75,0336.52554.5867−0.8043
qSDW6-106536,139–543,7175.28864.5867−0.8043
qSDW8-2820,095,951–20,541,2823.98326.08080.9261
qSDW8-3820,095,951–20,541,2823.62781.64610.1350
qSDW11-11123,793,431–23,795,2474.94214.5867−0.8043
qSDW12-71225,921,238–26,742,3064.03961.85630.1470
Chr: Chromosome number; LOD: Logarithm of odds value; PVE: Phenotypic variance explained; Add: Additive effect; Cold: Cold stress; Cold + BR-N: Data under normal temperature condition were used as the control; Cold + BR-C: Data under cold treatment condition were used as the control.
Table
Table 5. QTL summary under different conditions for SWW.
Table 5. QTL summary under different conditions for SWW.
ConditionQTLChrIntervalLODPVE (%)Add
ColdqSWW1-3117,059,460–18,491,4568.462618.24390.1201
qSWW7-17187,721–321,1994.15118.7113−0.1427
qSWW12-21219,002,436–19,612,4273.92557.57600.0745
BR + C-NqSWW3-1317,877,339–18,608,2363.75392.37490.1061
qSWW6-2649,710–75,0335.23045.4240−0.7790
qSWW6-3675,033–86,2573.68945.4240−0.7790
qSWW6-668,551,070–9,187,2874.68275.4241−0.7790
qSWW11-21127,344,371–27,803,8344.92433.23510.1113
BR + C-CqSWW1-118,430,459–9,795,8183.94050.7179−0.6502
qSWW1-219,795,818–9,824,2824.29510.7179−0.6501
qSWW4-1430,444,673–30,569,8203.19410.7601−0.4591
qSWW5-1528,253,886–28,304,1363.13920.8057−0.3553
qSWW6-1649,710–75,0335.61210.7489−0.4876
qSWW6-4675,033–86,2574.05250.7134−0.2935
qSWW6-56536,139–543,7173.30690.7941−0.3672
qSWW6-768,551,070–9,187,2875.49220.7179−0.6503
qSWW6-869,187,287–9,300,7703.57350.7532−0.4835
qSWW9-1916,774,761–16,829,1254.14660.8111−0.3491
qSWW10-1102,188,025–2,329,7343.41120.7576−0.4731
qSWW11-11123,793,431–23,795,2473.47240.7573−0.4714
qSWW12-11218,482,234–19,002,4364.06640.7574−0.4714
qSWW12-31220,179,243–20,499,1133.53130.7576−0.4694
qSWW12-41220,499,113–20,522,9903.35910.7582−0.4672
Chr: Chromosome number; LOD: Logarithm of odds value; PVE: Phenotypic variance explained; Add: Additive effect; Cold: Cold stress; Cold + BR-N: Data under normal temperature condition were used as the control; Cold + BR-C: Data under cold treatment condition were used as the control.
Table
Table 6. Summary of the annotated candidate genes possibly associated with cold tolerance or BR signaling pathway.
Table 6. Summary of the annotated candidate genes possibly associated with cold tolerance or BR signaling pathway.
ConditionQTLChrLocus NameGene Coordinates Gene Product
CqSSL1-21LOC_Os01g4026022,731,943–22,733,237OsWRKY77- Superfamily of TFs with WRKY and zinc finger domains
C + BRqSWW11-211LOC_Os11g4574027,670,321–27,673,334MYB family transcription factor
C + BR-NqSSL7-17LOC_Os07g0580523,483,968–23,485,531OsBZR1, transcription factor, Brassinosteroid (BR)-regulated growth response
CqSSL7-37LOC_Os07g084404,338,514–4,342,219bHLH transcription factor
C/C + BR-NqSDW6-18/176LOC_Os06g4475027,025,437–27,029,339AP2 domain-containing protein, expressed
C + BR-N/C + BR-CqSDW6-9/106LOC_Os06g01966543,057–545,780auxin-induced protein 5NG4, putative, expressed
C + BR-CqSRL12-612LOC_Os12g4182025,901,456–25,907,573heat shock protein DnaJ, putative, expressed
CqSSL7-17LOC_Os07g01480306,009–307,555oxygen evolving enhancer protein 3 domain-containing protein, expressed
C + BR-CqSRL12-612LOC_Os12g4170025,815,277–25,818,874LSD1, zinc finger domain-containing protein, expressed
C + BR-CqSRL5-35LOC_Os05g3938023,101,671–23,104,153zinc finger, C3HC4 type domain-containing protein, expressed
C + BR-CqSWW9-19LOC_Os09g2765016,822,234–16,825,686ZOS9-14, C2H2 zinc finger protein, expressed
C + BR-NqSDW9-59LOC_Os09g2766016,829,409–16,837,307OsFBO21, F-box, and other domain-containing protein, expressed
C + BR-NqSDW10-110LOC_Os10g045902,185,226–2,188,547OsFBX358, F-box domain-containing protein, expressed
Chr: Chromosome number; C: Cold treatment; C + BR-N: Data under normal temperature condition were used as the control; C + BR-C: Data under cold treatment condition were used as the control.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, Z.; Yao, J.; Cheng, Y.; Zhang, W.; Xu, Z.; Li, M.; Huang, J.; Ma, D.; Zhao, M. Identification of QTL under Brassinosteroid-Combined Cold Treatment at Seedling Stage in Rice Using Genotyping-by-Sequencing (GBS). Plants 2022, 11, 2324. https://doi.org/10.3390/plants11172324

AMA Style

Guo Z, Yao J, Cheng Y, Zhang W, Xu Z, Li M, Huang J, Ma D, Zhao M. Identification of QTL under Brassinosteroid-Combined Cold Treatment at Seedling Stage in Rice Using Genotyping-by-Sequencing (GBS). Plants. 2022; 11(17):2324. https://doi.org/10.3390/plants11172324

Chicago/Turabian Style

Guo, Zhifu, Jialu Yao, Yishan Cheng, Wenzhong Zhang, Zhengjin Xu, Maomao Li, Jing Huang, Dianrong Ma, and Minghui Zhao. 2022. "Identification of QTL under Brassinosteroid-Combined Cold Treatment at Seedling Stage in Rice Using Genotyping-by-Sequencing (GBS)" Plants 11, no. 17: 2324. https://doi.org/10.3390/plants11172324

APA Style

Guo, Z., Yao, J., Cheng, Y., Zhang, W., Xu, Z., Li, M., Huang, J., Ma, D., & Zhao, M. (2022). Identification of QTL under Brassinosteroid-Combined Cold Treatment at Seedling Stage in Rice Using Genotyping-by-Sequencing (GBS). Plants, 11(17), 2324. https://doi.org/10.3390/plants11172324

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop