BACKGROUND: Polychromatic flow cytometry is a popular technique that has wide usage in the medical sciences, especially for studying phenotypic properties of cells. The high-dimensionality of data generated by flow cytometry usually makes it difficult to visualize. The naive solution of simply plotting two-dimensional graphs for every combination of observables becomes impractical as the number of dimensions increases. A natural solution is to project the data from the original high dimensional space to a lower dimensional space while approximately preserving the overall relationship between the data points. The expert can then easily visualize and analyze this low-dimensional embedding of the original dataset. RESULTS: This paper describes a new method, SANJAY, for visualizing high-dimensional flow cytometry datasets. This technique uses a decision procedure to automatically synthesize two-dimensional and three-dimensional projections of the original high-dimensional data while trying to minimize distortion. We compare SANJAY to the popular multidimensional scaling (MDS) approach for visualization of small data sets drawn from a representative set of benchmarks, and our experiments show that SANJAY produces distortions that are 1.44 to 4.15 times smaller than those caused due to MDS. Our experimental results show that SANJAY also outperforms the Random Projections technique in terms of the distortions in the projections. CONCLUSIONS: We describe a new algorithmic technique that uses a symbolic decision procedure to automatically synthesize low-dimensional projections of flow cytometry data that typically have a high number of dimensions. Our algorithm is the first application, to our knowledge, of using automated theorem proving for automatically generating highly-accurate, low-dimensional visualizations of high-dimensional data.

MOTIVATION: p38 mitogen-activated protein kinase activation plays an important role in resistance to chemotherapeutic cytotoxic drugs in treating multiple myeloma (MM). However, how the p38 mitogen-activated protein kinase signaling pathway is involved in drug resistance, in particular the roles that the various p38 isoforms play, remains largely unknown. METHOD: To explore the underlying mechanisms, we developed a novel systems biology approach by integrating liquid chromatography-mass spectrometry and reverse phase protein array data from human MM cell lines with computational pathway models in which the unknown parameters were inferred using a proposed novel algorithm called modularized factor graph. RESULTS: New mechanisms predicted by our models suggest that combined activation of various p38 isoforms may result in drug resistance in MM via regulating the related pathways including extracellular signal-regulated kinase (ERK) pathway and NFкB pathway. ERK pathway regulating cell growth is synergistically regulated by p38δ isoform, whereas nuclear factor kappa B (NFкB) pathway regulating cell apoptosis is synergistically regulated by p38α isoform. This finding that p38δ isoform promotes the phosphorylation of ERK1/2 in MM cells treated with bortezomib was validated by western blotting. Based on the predicted mechanisms, we further screened drug combinations in silico and found that a promising drug combination targeting ERK1/2 and NFκB might reduce the effects of drug resistance in MM cells. This study provides a framework of a systems biology approach to studying drug resistance and drug combination selection. AVAILABILITY AND IMPLEMENTATION: RPPA experimental Data and Matlab source codes of modularized factor graph for parameter estimation are freely available online at http://ctsb.is.wfubmc.edu/publications/modularized-factor-graph.php.

Cancer initiating cells have been documented in multiple myeloma and believed to be a key factor that initiates and drives tumor growth, differentiation,metastasis, and recurrence of the diseases. Although myeloma initiating cells (MICs) are likely to share many properties of normal stem cells, the underlying mechanisms regulating the fate of MICs are largely unknown. Studies designed to explore such communication are urgently needed to enhance our ability to predict the fate decisions of ICs (self-renewal, differentiation, and proliferation). In this study, we developed a novel system to understand the intercellular communication between MICs and their niche by seamlessly integrating experimental data and mathematical model. We first designed dynamic cell culture experiments and collected three types of cells (side population cells, progenitor cells, and mature myeloma cells) under various cultural conditions with flow cytometry. Then we developed a lineage model with ordinary differential equations by considering secreted factors, self-renewal, differentiation, and other biological functions of those cells, to model the cell–cell interactions among the three cell types. Particle swarm optimization was employed to estimate the model parameters by fitting the experimental data to the lineage model. The theoretical results show that the correlation coefficient analysis can reflect the feedback loops among the three cell types, the intercellular feedback signaling can regulate cell population dynamics, and the culture strategies can decide cell growth. This study provides a basic framework of studying cell–cell interactions in regulating MICs fate.

BACKGROUND: Gene fusions, which result from abnormal chromosome rearrangements, are a pathogenic factor in cancer development. The emerging RNA-Seq technology enables us to detect gene fusions and profile their features. RESULTS: In this paper, we proposed a novel fusion detection tool, FusionQ, based on paired-end RNA-Seq data. This tool can detect gene fusions, construct the structures of chimerical transcripts, and estimate their abundances. To confirm the read alignment on both sides of a fusion point, we employed a new approach, "residual sequence extension", which extended the short segments of the reads by aggregating their overlapping reads. We also proposed a list of filters to control the false-positive rate. In addition, we estimated fusion abundance using the Expectation-Maximization algorithm with sparse optimization, and further adopted it to improve the detection accuracy of the fusion transcripts. Simulation was performed by FusionQ and another two stated-of-art fusion detection tools. FusionQ exceeded the other two in both sensitivity and specificity, especially in low coverage fusion detection. Using paired-end RNA-Seq data from breast cancer cell lines, FusionQ detected both the previously reported and new fusions. FusionQ reported the structures of these fusions and provided their expressions. Some highly expressed fusion genes detected by FusionQ are important biomarkers in breast cancer. The performances of FusionQ on cancel line data still showed better specificity and sensitivity in the comparison with another two tools. CONCLUSIONS: FusionQ is a novel tool for fusion detection and quantification based on RNA-Seq data. It has both good specificity and sensitivity performance. FusionQ is free and available at http://www.wakehealth.edu/CTSB/Software/Software.htm.

BACKGROUND: The development of techniques for sequencing the messenger RNA (RNA-Seq) enables it to study the biological mechanisms such as alternative splicing and gene expression regulation more deeply and accurately. Most existing methods employ RNA-Seq to quantify the expression levels of already annotated isoforms from the reference genome. However, the current reference genome is very incomplete due to the complexity of the transcriptome which hiders the comprehensive investigation of transcriptome using RNA-Seq. Novel study on isoform inference and estimation purely from RNA-Seq without annotation information is desirable. RESULTS: A Nonnegativity and Sparsity constrained Maximum APosteriori (NSMAP) model has been proposed to estimate the expression levels of isoforms from RNA-Seq data without the annotation information. In contrast to previous methods, NSMAP performs identification of the structures of expressed isoforms and estimation of the expression levels of those expressed isoforms simultaneously, which enables better identification of isoforms. In the simulations parameterized by two real RNA-Seq data sets, more than 77% expressed isoforms are correctly identified and quantified. Then, we apply NSMAP on two RNA-Seq data sets of myelodysplastic syndromes (MDS) samples and one normal sample in order to identify differentially expressed known and novel isoforms in MDS disease. CONCLUSIONS: NSMAP provides a good strategy to identify and quantify novel isoforms without the knowledge of annotated reference genome which can further realize the potential of RNA-Seq technique in transcriptome analysis. NSMAP package is freely available at https://sites.google.com/site/nsmapforrnaseq.

BACKGROUND: DNA copy number aberration (CNA) is very important in the pathogenesis of tumors and other diseases. For example, CNAs may result in suppression of anti-oncogenes and activation of oncogenes, which would cause certain types of cancers. High density single nucleotide polymorphism (SNP) array data is widely used for the CNA detection. However, it is nontrivial to detect the CNA automatically because the signals obtained from high density SNP arrays often have low signal-to-noise ratio (SNR), which might be caused by whole genome amplification, mixtures of normal and tumor cells, experimental noise or other technical limitations. With the reduction in SNR, many false CNA regions are often detected and the true CNA regions are missed. Thus, more sophisticated statistical models are needed to make the CNAs detection, using the low SNR signals, more robust and reliable. RESULTS: This paper presents a conditional random pattern (CRP) model for CNA detection where much contextual cues are explored to suppress the noise and improve CNA detection accuracy. Both simulated and the real data are used to evaluate the proposed model, and the validation results show that the CRP model is more robust and reliable in the presence of noise for CNA detection using high density SNP array data, compared to a number of widely used software packages. CONCLUSIONS: The proposed conditional random pattern (CRP) model could effectively detect the CNA regions in the presence of noise.

Flow has been widely used for the diagnosis of various diseases. We proposed an automated gating approach to analyze and cluster flow data and provide researchers and physicians a 3-D way of viewing flow data. This automated approach was designed to give reproducible results, which avoid the subjective and time-consuming human manipulation. Doublets problem will be eliminated by our tool, and the result were compared with the manual gating. In addition, a 3D visualization will not only give a user-friendly representation but also avert misclassification would occur in lower dimension. We demonstrated how the 3D approach can be used to augment the K-means method in classifying the data. To validate the feasibility of our proposed automatic approach, our result is compared with the result done at the Methodist Hospital in Houston.

MOTIVATION: Loss of heterozygosity (LOH) is one of the most important mechanisms in the tumor evolution. LOH can be detected from the genotypes of the tumor samples with or without paired normal samples. In paired sample cases, LOH detection for informative single nucleotide polymorphisms (SNPs) is straightforward if there is no genotyping error. But genotyping errors are always unavoidable, and there are about 70% non-informative SNPs whose LOH status can only be inferred from the neighboring informative SNPs. RESULTS: This article presents a novel LOH inference and segmentation algorithm based on the conditional random pattern (CRP) model. The new model explicitly considers the distance between two neighboring SNPs, as well as the genotyping error rate and the heterozygous rate. This new method is tested on the simulated and real data of the Affymetrix Human Mapping 500K SNP arrays. The experimental results show that the CRP method outperforms the conventional methods based on the hidden Markov model (HMM). AVAILABILITY: Software is available upon request.