  Complex phenotypic differences among different acute leukemias cannot be
fully captured by analyzing the expression levels of one single molecule, such
as a miR, at a time, but requires systematic analysis of large sets of miRs.
While a popular approach for analysis of such datasets is principal component
analysis (PCA), this method is not designed to optimally discriminate different
phenotypes. Moreover, PCA and other low-dimensional representation methods
yield linear or non-linear combinations of all measured miRs. Global human miR
expression was measured in AML, B-ALL, and T-ALL cell lines and patient RNA
samples. By systematically applying support vector machines to all measured
miRs taken in dyad and triad groups, we built miR networks using cell line data
and validated our findings with primary patient samples. All the coordinately
transcribed members of the miR-23a cluster (which includes also miR-24 and
miR-27a), known to function as tumor suppressors of acute leukemias, appeared
in the AML, B-ALL and T-ALL centric networks. Subsequent qRT-PCR analysis
showed that the most connected miR in the B-ALL-centric network, miR-708, is
highly and specifically expressed in B-ALLs, suggesting that miR-708 might
serve as a biomarker for B-ALL. This approach is systematic, quantitative,
scalable, and unbiased. Rather than a single signature, our approach yields a
network of signatures reflecting the redundant nature of biological signaling
pathways. The network representation allows for visual analysis of all
signatures by an expert and for future integration of additional information.
Furthermore, each signature involves only small sets of miRs, such as dyads and
triads, which are well suited for in depth validation through laboratory
experiments such as loss- and gain-of-function assays designed to drive changes
in leukemia cell survival, proliferation and differentiation.
