Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

This article is part of the supplement: UT-ORNL-KBRIN Bioinformatics Summit 2010

Open Access Oral presentation

Inferring gene coexpression networks for low dose ionizing radiation using graph theoretical algorithms and systems genetics

Sudhir Naswa1, Gary L Rogers1, Rachel M Lynch2, Stephen A Kania3, Suchita Das2, Elissa J Chesler24, Arnold M Saxton5, Brynn H Voy25 and Michael A Langston1*

Author affiliations

1 Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, TN 37996, USA

2 Systems Genetics Group, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA

3 Department of Comparative Medicine, University of Tennessee, Knoxville, TN 37996, USA

4 The Jackson Laboratory, Bar Harbor, ME 04609, USA

5 Department of Animal Science, University of Tennessee, Knoxville, TN 37996, USA

For all author emails, please log on.

Citation and License

BMC Bioinformatics 2010, 11(Suppl 4):O5  doi:10.1186/1471-2105-11-S4-O5

The electronic version of this article is the complete one and can be found online at:

Published:23 July 2010

© 2010 Langston et al; licensee BioMed Central Ltd.


Biological data generated through large scale -omics technologies have resulted in a new paradigm in the study of biological systems. Instead of focusing on individual genes or proteins these technologies enable us to extract biological networks using powerful computing and statistical algorithms that are scalable to very large datasets.

Materials and methods

We have developed a tool chain using novel graph algorithms to extract gene coexpression networks from microarray data. We highlight implementation of our tool chain to investigate the effects of in vivo low dose ionizing radiation treatments on mice. We are using systems genetics approach to investigate the biological effects of low dose (10 cGy) ionizing radiation. We measured the base line gene expression profile from spleen tissue of BXD recombinant inbred mice using Illumina microarrays. The data was filtered using coefficient of variance after robust spline normalization and variance stabilizing transformation. A graph was then derived from this data, with probes as vertices and edges between them representing correlations. The graph was analyzed using our toolkit to find the size and number of maximal cliques. We deployed another tool called paraclique that relaxes clique’s requirement that every edge be present between all vertices. Paraclique enables us to account for inherent noise in the microarray data and stochastic nature of biological processes. Using immunophenotype data from the baseline BXD mice, we employed biclique analysis to determine interactions between genotypes and immunophenotypes (%CD4, %CD3, LN T:B, %CD8, and LN CD4:CD8). We also extracted eQTLs from BXD data using QTL-Reaper from base line gene expression profiles. 1881 transcripts were associated with 686 loci. The eQTLs were classified as cis or trans according to their genomic positions. Besides population level studies we also investigated the differential effect of low dose and high dose (1Gy) of ionizing radiations on spleen gene expression in inbred parental strains (C57BL/6J and DBA/2J) of BXD recombinant inbred mice as well as BALB/c mice, a known radiation-sensitive strain.