Network preservation analysis to identify transcriptional biomarkers related to flowering in Crocus sativus

Crocus sativus L. is known as an ornamental geophyte and a source of valuable spice and secondary metabolites. Network preservation module analysis is one of the best approaches to revealing special features of different conditions. It can determine patterns of divergence and conservation between transcriptome data. Herein, we explored the regulatory genes of the flowering process by RNA-Seq data containing flowering and non-flowering samples in gene expression profiles. Persevered module analysis revealed three significant non-persevered modules related to the flowering process, namely pink, green, and blue. Several hub genes associated with non-preserved modules such as PIA1, NAC90, ALY3, Sus3, MYB31, ARF5/MP, MYB31, HD-ZIP, SEP3d, OR_B, AGL6a, bZIP(TGA1) and GRAS were identified. These candidate genes can be considered key diagnostic biomarkers for the flowering process. Here, we also compared two approaches, WGCNA and NetRep for module preservation analysis. The results of these methods were consistent with non-preserved modules. NetRep was a faster (11 times) and more efficient (run more than 10000 permutations for each comparison) method than WGCNA module preservation. Differential expression genes (DEGs) screening showed that many hub genes were downregulated in non-flowering than flowering samples. Our finding revealed regulatory mechanisms of the flowering process in C. sativus as can be developed transcriptional biomarkers which could pave the way for promoting saffron yield via flowering induction.