CHPC Summer Presentation Schedule Now Available
- Intro to Python
- part 1: Thur, June 20
- part 2: Tue, June 25
- part 3 (Numpy 1) Thur, June 27
- part 4 (Numpy 2) Tue, July 2
Registration Open for CHPC's Build a Cluster Boot Camp
Posted June 3rd, 2019
ATTENTION Owners of Virtual Machines in the General Environment
Posted May 23rd, 2019
Notchpeak Outage: notch046, notch065-089
May 21st, from 1pm until 4pm
New Scratch File System Available in the General Environment
Posted: May 15, 2019
New GPU nodes added to Notchpeak
Posted: May 8, 2019
New Enforcement of Interactive Node Usage Policy
Posted: April 4, 2019
CHPC Spring 2019 Newsletter
News History...
Uncertainty Quantification of RNA-Seq Co-expression Networks
By Lance Pflieger and Julio Facelli, Department of Biomedical Informatics
Systems biology utilizes the complex and copious data originating from the “omics” fields to increase understanding of biology by studying interactions among biological entities. Gene co-expression network analysis is a systems biology technique derived from graph theory that uses RNA expression data to infer functional similar genes or regulatory pathways. Gene co-expression network analysis is a computationally intensive process that requires matrix operations on tens-of-thousands of genes/transcripts. This technique has been useful in drug discovery, functional annotation of a gene and insight into disease pathology.
To assess the effect of uncertainty inherent with gene expression data, our group utilized CHPC resources to characterize variation in gene expression estimates and simulate a large quantity of co-expression networks based on this variation. The figure shown is a representation of network generated using WGCNA and expression data from the disease Spinocerebellar Type 2 (SCA2). The colors represent highly connected subnetworks of genes which are used to correlate similar gene clusters with a phenotypic trait. Our results show that uncertainty has a large effect on downstream results including subnetwork structure, hub genes identification and enrichment analysis. For instance, we find that the number of subnetworks correlating with the SCA2 phenotype varies from 1 to 6 subnetworks. While a small gene co-expression network analysis can be performed using only modest computation resources, the scale of resources required to perform uncertainty quantification (UQ) using Monte Carlo ensemble methods is several orders of magnitude larger, which are only available at CHPC.
System Status
General Environment
| General Nodes | ||
|---|---|---|
| system | cores | % util. |
| ember | 960/960 | 100% |
| kingspeak | 848/848 | 100% |
| notchpeak | 1212/1212 | 100% |
| lonepeak | 1100/1100 | 100% |
| Owner/Restricted Nodes | ||
| system | cores | % util. |
| ash | 5240/7300 | 71.78% |
| notchpeak | 1720/1720 | 100% |
| ember | 708/1220 | 58.03% |
| kingspeak | 5832/5916 | 98.58% |
| lonepeak | 400/400 | 100% |
Protected Environment
| General Nodes | ||
|---|---|---|
| system | cores | % util. |
| redwood | 60/500 | 12% |
| Owner/Restricted Nodes | ||
| system | cores | % util. |
| redwood | 2192/3232 | 67.82% |