Co-Expression Network-based Analysis associated with potato initial resistance

RNA sequencing (RNA-seq) providing genome-wide expression datasets has been successfully used to study gene expression patterns and regulation mechanism among multiple samples. Gene co-expression networks (GCNs) studies within or across species showed that coordinated genes in expression patterns are often functionally related. For potatoes, a large amount of publicly available transcriptome datasets have been generated but an optimal GCN detecting expression patterns in different genotypes, tissues and environmental conditions, is lacking. We constructed a potato GCN using 16 published RNA-Seq datasets covering 11 cultivars from native habitat worldwide. The correlations of gene expression were assessed pair-wisely and biologically meaningful gene modules which are highly connected in GCN were identified. One of the primitively native-farmer-selected cultivars in the Andes, ssp.Andigena, had relative far distance in gene expression patterns with other modern varieties. GCN in further enriched 134 highly and specifically co-expressed genes in ssp.Andigena associated with potato disease and stress resistance, which underlying the dramatic shift in evolutionary pressures during potato artificial domestication. In total, the network was consisted of into 14 gene models that involves in a variety of plant processes, which sheds light on how gene modules organized intra- and inter-varieties in the context of evolutionary divergence and provides a basis of information resource for potato gene functional studies.