Identification of candidate genes controlling cold tolerance at the early seedling stage from Dongxiang wild rice by QTL mapping, BSA-Seq and RNA-Seq
Authors: Zhou Shiqi , Wu Ting , Li Xia , Wang Shilin , Hu Biaolin 
Published: July 2024

Abstract
Background
The cold tolerance of rice is closely related to its production and geographic distribution. The identification of cold tolerance-related genes is of important significance for developing cold-tolerant rice. Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) is well-adapted to the cold climate of northernmost-latitude habitats ever found in the world, and is one of the most valuable rice germplasms for cold tolerance improvement.

Results
Transcriptome analysis revealed genes differentially expressed between Xieqingzao B (XB; a cold sensitive variety) and 19H19 (derived from an interspecific cross between DXWR and XB) in the room temperature (RT), low temperature (LT), and recovery treatments. The results demonstrated that chloroplast genes might be involved in the regulation of cold tolerance in rice. A high-resolution SNP genetic map was constructed using 120 BC5F2 lines derived from a cross between 19H19 and XB based on the genotyping-by-sequencing (GBS) technique. Two quantitative trait loci (QTLs) for cold tolerance at the early seedling stage (CTS), qCTS12 and qCTS8, were detected. Moreover, a total of 112 candidate genes associated with cold tolerance were identified based on bulked segregant analysis sequencing (BSA-seq). These candidate genes were divided into eight functional categories, and the expression trend of candidate genes related to ‘oxidation-reduction process’ and ‘response to stress’ differed between XB and 19H19 in the RT, LT and recovery treatments. Among these candidate genes, the expression level of LOC_Os12g18729 in 19H19 (related to ‘response to stress’) decreased in the LT treatment but restored and enhanced during the recovery treatment whereas the expression level of LOC_Os12g18729 in XB declined during recovery treatment. Additionally, XB contained a 42-bp deletion in the third exon of LOC_Os12g18729, and the genotype of BC5F2 individuals with a survival percentage (SP) lower than 15% was consistent with that of XB. Weighted gene coexpression network analysis (WGCNA) and modular regulatory network learning with per gene information (MERLIN) algorithm revealed a gene interaction/coexpression network regulating cold tolerance in rice. In the network, differentially expressed genes (DEGs) related to ‘oxidation-reduction process’, ‘response to stress’ and ‘protein phosphorylation’ interacted with LOC_Os12g18729. Moreover, the knockout mutant of LOC_Os12g18729 decreased cold tolerance in early rice seedling stage signifcantly compared with that of wild type.

Conclusions
In general, study of the genetic basis of cold tolerance of rice is important for the development of cold-tolerant rice varieties. In the present study, QTL mapping, BSA-seq and RNA-seq were integrated to identify two CTS QTLs qCTS8 and qCTS12. Furthermore, qRT-PCR, genotype sequencing and knockout analysis indicated that LOC_Os12g18729 could be the candidate gene of qCTS12. These results are expected to further exploration of the genetic mechanism of CTS in rice and improve cold tolerance of cultivated rice by introducing the cold tolerant genes from DXWR through marker-assisted selection.