Abstract: Systems and methods of selecting combinatorial coordinately dysregulated biomarker subnetworks are provided. In one embodiment, a method comprises comparing the normalized gene expression data to a predetermined threshold to provide binary gene expression data associated with phenotype samples and control samples, analyzing subnetwork states of the binary gene expression data associated with phenotype samples and control samples to identify gene expression patterns that occur in phenotype samples and do not occur in control samples and identifying a subnetwork that provides gene expression patterns indicative of a sample being a phenotype sample.