Inferring the transcriptional landscape of bovine skeletal muscle by integrating co-expression networks

Nicholas J. Hudson, Antonio Reverter, YongHong Wang, Paul L. Greenwood and Brian P. Dalrymple
CSIRO Livestock Industries, 306 Carmody Road, St. Lucia, Qld 4067, Australia


Background: Despite modern technologies and novel computational approaches, decoding causal transcriptional regulation remains challenging. This is particularly true for less well studied organisms and when only gene expression data is available. In muscle a small number of well characterised transcription factors are proposed to regulate development. Therefore, muscle appears to be a tractable system for proposing new computational approaches.

Methodology / Principal Findings: Here we report a simple algorithm that asks which transcriptional regulator has the highest average absolute co-expression correlation to the genes in a co-expression module? It correctly infers a number of known causal regulators of fundamental biological processes, including cell cycle activity (E2F1), glycolysis (HLF), mitochondrial transcription (TFB2M), adipogenesis (PIAS1), neuronal development (TLX3), immune function (IRF1) and vasculogenesis (SOX17), within a skeletal muscle context. However, none of the canonical pro-myogenic transcription factors (MYOD1, MYOG, MYF5, MYF6 and MEF2C) were linked to muscle structural gene expression modules. Co-expression values were computed using developing bovine muscle from 60 days post conception (early foetal) to 30 months post natal (adulthood) for two breeds of cattle, in addition to a nutritional comparison with a third breed. A number of transcriptional landscapes were constructed and integrated into an always correlated landscape. One notable feature was a metabolic axis formed from glycolysis genes at one end, nuclear-encoded mitochondrial protein genes at the other, and centrally tethered by mitochondrially-encoded mitochondrial protein genes.

Conclusions / Significance: The new module-to-regulator algorithm complements our recently described Regulatory Impact Factor analysis. Together with a simple examination of a co-expression modules contents, these three gene expression approaches are starting to illuminate the in vivo transcriptional regulation of skeletal muscle development.


An image of the bovine muscle always correlated co-expression landsape is displayed below, available here. The Cytoscape network file can be found here



Note: Internet Explorer users right click and "Save as" with "cys" format to avoid IE guessing this is a zip file.