  MicroRNAs (miRNAs) are small endogenous regulatory molecules that modulate
gene expression post-transcriptionally. Although differential expression of
miRNAs have been implicated in many diseases (including cancers), the
underlying mechanisms of action remain unclear. Because each miRNA can target
multiple genes, miRNAs may potentially have functional implications for the
overall behavior of entire pathways. Here we investigate the functional
consequences of miRNA dysregulation through an integrative analysis of miRNA
and mRNA expression data using a novel approach that incorporates pathway
information a priori. By searching for miRNA-pathway associations that differ
between healthy and tumor tissue, we identify specific relationships at the
systems-level which are disrupted in cancer. Our approach is motivated by the
hypothesis that if a miRNA and pathway are associated, then the expression of
the miRNA and the collective behavior of the genes in a pathway will be
correlated. As such, we first obtain an expression-based summary of pathway
activity using Isomap, a dimension reduction method which can articulate
nonlinear structure in high-dimensional data. We then search for miRNAs that
exhibit differential correlations with the pathway summary between phenotypes
as a means of finding aberrant miRNA-pathway coregulation in tumors. We apply
our method to cancer data using gene and miRNA expression datasets from The
Cancer Genome Atlas (TCGA) and compare ${\sim}10^5$ miRNA-pathway relationships
between healthy and tumor samples from four tissues (breast, prostate, lung,
and liver). Many of the flagged pairs we identify have a biological basis for
disruption in cancer.
