Transcriptomics profiling in response to cold stress in cultivated rice and weedy rice

Weedy rice is an important germplasm resource for rice improvement because it has useful genes for many abiotic stresses including cold tolerance. We identified the cold tolerance and cold sensitivity of two weedy rice lines (WR 03-35 and WR 03-26) and two cultivated rice lines (Kongyu 131 and 9311). During the seedling stage of these lines, we used RNA-seq to measure changes in weedy rice and cultivated rice whole-genome transcriptome before and after cold treatment. We identified 14,213 and 14,730 differentially expressed genes (DEGs) in cold-tolerant genotypes (WR 03-35, Kongyu 131), and 9219 and 720 DEGs were observed in two cold sensitive genotypes (WR 03-26, 9311). Many common and special DEGs were analyzed in cold-tolerant and cold sensitive genotypes, respectively. Some typical genes related to cold stress such as the basic helix-loop-helix (bHLH) gene and leucine-rich repeat (LRR) domain gene etc. The number of these DEGs in cold-tolerant genotypes is more than those found in cold-sensitive genotypes. The gene ontology (GO) enrichment analyses showed significantly enriched terms for biological processes, cellular components and molecular functions. In addition, some genes related to several plant hormones such as abscisic acid (ABA), gibberellic acid (GA), auxin and ethylene were identified. To confirm the RNA-seq data, semi-quantitative RT-PCR and qRT-PCR were performed on 12 randomly selected DEGs. The expression patterns of RNA-seq on these genes corresponded with the semi-quantitative RT-PCR and qRT-PCR method. This study suggests the gene resources related to cold stress from weedy rice could be valuable for understanding the mechanisms involved in cold stress and rice breeding for improving cold tolerance.