A Parallelizable Model for Analyzing Cancer Tissue Heterogeneity.

In a cancer study, the heterogeneous nature of a cell population creates a lot of challenges. Efficient determination of the compositional breakup of a cell population, from gene expression measurements, is critical to the success in a cancer study. This paper presents a new model for analyzing heterogeneity in cancer tissue using Markov chain Monte Carlo (MCMC) algorithms; we aim to compute the proportion wise breakup of the cell population on a GPU. We also show that the model computation time does not depend on the input data size, because the computation required to estimate the compositional breakup are parallelized. This model uses qPCR (quantitative polymerase chain reaction) gene expression data to determine compositional breakup in the heterogeneous cell population. We test this model on synthetic data and real-world data collected from fibroblasts. We also show how well this model scales to hundreds of gene expression data.

novoCaller: a Bayesian network approach for de novo variant calling from pedigree and population sequence data.

MOTIVATION: De novo mutations (i.e. newly occurring mutations) are a pre-dominant cause of sporadic dominant monogenic diseases and play a significant role in the genetics of complex disorders. De novo mutation studies also inform population genetics models and shed light on the biology of DNA replication and repair. Despite the broad interest, there is room for improvement with regard to the accuracy of de novo mutation calling. RESULTS: We designed novoCaller, a Bayesian variant calling algorithm that uses information from read-level data both in the pedigree and in unrelated samples. The method was extensively tested using large trio-sequencing studies, and it consistently achieved over 97% sensitivity. We applied the algorithm to 48 trio cases of suspected rare Mendelian disorders as part of the Brigham Genomic Medicine gene discovery initiative. Its application resulted in a significant reduction in the resources required for manual inspection and experimental validation of the calls. Three de novo variants were found in known genes associated with rare disorders, leading to rapid genetic diagnosis of the probands. Another 14 variants were found in genes that are likely to explain the phenotype, and could lead to novel disease-gene discovery. AVAILABILITY AND IMPLEMENTATION: Source code implemented in C++ and Python can be downloaded from https://github.com/bgm-cwg/novoCaller. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A Conjugate Exponential Model for Cancer Tissue Heterogeneity.

The diagnosis and treatment of cancer is made difficult by the heterogeneous nature of the cell population. Determining its compositional breakup from measurements of various measurable traits (such as gene expression measurements) is an important problem in the field of cancer diagnosis and treatment. In addition, the computational aspect of the problem also needs attention. The processing of the collected data must be done as efficiently as possible in terms of computational speed and memory requirements. The use of Markov chain Monte Carlo methods is time consuming, and hence, other methods need to be used for the analysis. In this paper, we develop a model for heterogeneous cancer tissue, which uses quantitative polymerase chain reaction gene expression data to determine the compositional breakup of cell populations in the heterogeneous tissue. We develop a fast algorithm for the model using variational methods and demonstrate its use on synthetic and real-world gene expression data collected from fibroblasts and compare the performance of the algorithm with other methods such as Markov chain Monte Carlo and expectation maximization.

A model for cancer tissue heterogeneity.

An important problem in the study of cancer is the understanding of the heterogeneous nature of the cell population. The clonal evolution of the tumor cells results in the tumors being composed of multiple subpopulations. Each subpopulation reacts differently to any given therapy. This calls for the development of novel (regulatory network) models, which can accommodate heterogeneity in cancerous tissues. In this paper, we present a new approach to model heterogeneity in cancer. We model heterogeneity as an ensemble of deterministic Boolean networks based on prior pathway knowledge. We develop the model considering the use of qPCR data. By observing gene expressions when the tissue is subjected to various stimuli, the compositional breakup of the tissue under study can be determined. We demonstrate the viability of this approach by using our model on synthetic data, and real-world data collected from fibroblasts.