Abstract
Reducing nitrogen (N) input is a key measure to achieve a sustainable rice production in China, especially in Jiangsu Province. Tiller is the basis for achieving panicle number that plays as a major factor in the yield determination. In actual production, excessive N is often applied in order to produce enough tillers in the early stages. Understanding how N regulates tillering in rice plants is critical to generate an integrative management to reduce N use and reaching tiller number target. Aiming at this objective, we utilized RNA sequencing and weighted gene co-expression network analysis (WGCNA) to compare the transcriptomes surrounding the shoot apical meristem of indica (Yangdao6, YD6) and japonica (Nipponbare, NPB) rice subspecies. Our results showed that N rate influenced tiller number in a different pattern between the two varieties, with NPB being more sensitive to N enrichment, and YD6 being more tolerant to high N rate. Tiller number was positively related to N content in leaf, culm and root tissue, but negatively related to the soluble carbohydrate content, regardless of variety. Transcriptomic comparisons revealed that for YD6 when N rate enrichment from low (LN) to medium (MN), it caused 115 DEGs (LN vs. MN), from MN to high level (HN) triggered 162 DEGs (MN vs. HN), but direct comparison of low with high N rate showed a 511 DEGs (LN vs. HN). These numbers of DEG in NPB were 87 (LN vs. MN), 40 (MN vs. HN), and 148 (LN vs. HN). These differences indicate that continual N enrichment led to a bumpy change at the transcription level. For the reported sixty-five genes which affect tillering, thirty-six showed decent expression in SAM at tiller starting phase, among them only nineteen being significantly influenced by N level, and two genes showed significant interaction between N rate and variety. Gene ontology analysis revealed that the majority of the common DEGs are involved in general stress responses, stimulus responses, and hormonal signaling process. WGCNA network identified twenty-two co-expressing gene modules and ten candidate hubgenes for each module. Several genes associated with tillering and N rate fall on the related modules. These indicate that there are more genes participating in tillering regulation in response to N enrichment.