State-transition modeling of blood transcriptome predicts disease evolution and treatment response in chronic myeloid leukemia.

Chronic myeloid leukemia (CML) is initiated and maintained by BCR::ABL which is clinically targeted using tyrosine kinase inhibitors (TKIs). TKIs can induce long-term remission but are also not curative. Thus, CML is an ideal system to test our hypothesis that transcriptome-based state-transition models accurately predict cancer evolution and treatment response. We collected time-sequential blood samples from tetracycline-off (Tet-Off) BCR::ABL-inducible transgenic mice and wild-type controls. From the transcriptome, we constructed a CML state-space and a three-well leukemogenic potential landscape. The potential's stable critical points defined observable disease states. Early states were characterized by anti-CML genes opposing leukemia; late states were characterized by pro-CML genes. Genes with expression patterns shaped similarly to the potential landscape were identified as drivers of disease transition. Re-introduction of tetracycline to silence the BCR::ABL gene returned diseased mice transcriptomes to a near healthy state, without reaching it, suggesting parts of the transition are irreversible. TKI only reverted the transcriptome to an intermediate disease state, without approaching a state of health; disease relapse occurred soon after treatment. Using only the earliest time-point as initial conditions, our state-transition models accurately predicted both disease progression and treatment response, supporting this as a potentially valuable approach to time clinical intervention, before phenotypic changes become detectable.

State-transition Modeling of Blood Transcriptome Predicts Disease Evolution and Treatment Response in Chronic Myeloid Leukemia.

Chronic myeloid leukemia (CML) is initiated and maintained by BCR::ABL which is clinically targeted using tyrosine kinase inhibitors (TKIs). TKIs can induce long-term remission but are also not curative. Thus, CML is an ideal system to test our hypothesis that transcriptome-based state-transition models accurately predict cancer evolution and treatment response. We collected time-sequential blood samples from tetracycline-off (Tet-Off) BCR::ABL-inducible transgenic mice and wild-type controls. From the transcriptome, we constructed a CML state-space and a three-well leukemogenic potential landscape. The potential's stable critical points defined observable disease states. Early states were characterized by anti-CML genes opposing leukemia; late states were characterized by pro-CML genes. Genes with expression patterns shaped similarly to the potential landscape were identified as drivers of disease transition. Re-introduction of tetracycline to silence the BCR::ABL gene returned diseased mice transcriptomes to a near healthy state, without reaching it, suggesting parts of the transition are irreversible. TKI only reverted the transcriptome to an intermediate disease state, without approaching a state of health; disease relapse occurred soon after treatment. Using only the earliest time-point as initial conditions, our state-transition models accurately predicted both disease progression and treatment response, supporting this as a potentially valuable approach to time clinical intervention even before phenotypic changes become detectable.

Dynamic patterns of microRNA expression during acute myeloid leukemia state-transition.

MicroRNAs (miRNAs) have been shown to hold prognostic value in acute myeloid leukemia (AML); however, the temporal dynamics of miRNA expression in AML are poorly understood. Using serial samples from a mouse model of AML to generate time-series miRNA sequencing data, we are the first to show that the miRNA transcriptome undergoes state-transition during AML initiation and progression. We modeled AML state-transition as a particle undergoing Brownian motion in a quasi-potential and validated the AML state-space and state-transition model to accurately predict time to AML in an independent cohort of mice. The critical points of the model provided a framework to align samples from mice that developed AML at different rates. Our mathematical approach allowed discovery of dynamic processes involved during AML development and, if translated to humans, has the potential to predict an individual's disease trajectory.

Dietary omega-3 fatty acid intake impacts peripheral blood DNA methylation -anti-inflammatory effects and individual variability in a pilot study.

Omega-3 or n-3 polyunsaturated fatty acids (PUFAs) are widely studied for health benefits that may relate to anti-inflammatory activity. However, mechanisms mediating an anti-inflammatory response to n-3 PUFA intake are not fully understood. Of interest is the emerging role of fatty acids to impact DNA methylation (DNAm) and thereby modulate mediating inflammatory processes. In this pilot study, we investigated the impact of n-3 PUFA intake on DNAm in inflammation-related signaling pathways in peripheral blood mononuclear cells (PBMCs) of women at high risk of breast cancer. PBMCs of women at high risk of breast cancer (n=10) were obtained at baseline and after 6 months of n-3 PUFA (5 g/d EPA+DHA dose arm) intake in a previously reported dose finding trial. DNA methylation of PBMCs was assayed by reduced representation bisulfite sequencing (RRBS) to obtain genome-wide methylation profiles at the single nucleotide level. We examined the impact of n-3 PUFA on genome-wide DNAm and focused upon a set of candidate genes associated with inflammation signaling pathways and breast cancer. We identified 24,842 differentially methylated CpGs (DMCs) in gene promoters of 5507 genes showing significant enrichment for hypermethylation in both the candidate gene and genome-wide analyses. Pathway analysis identified significantly hypermethylated signaling networks after n-3 PUFA treatment, such as the Toll-like Receptor inflammatory pathway. The DNAm pattern in individuals and the response to n-3 PUFA intake are heterogeneous. PBMC DNAm profiling suggests a mechanism whereby n-3 PUFAs may impact inflammatory cascades associated with disease processes including carcinogenesis.

Treatment-induced arteriolar revascularization and miR-126 enhancement in bone marrow niche protect leukemic stem cells in AML.

BACKGROUND: During acute myeloid leukemia (AML) growth, the bone marrow (BM) niche acquires significant vascular changes that can be offset by therapeutic blast cytoreduction. The molecular mechanisms of this vascular plasticity remain to be fully elucidated. Herein, we report on the changes that occur in the vascular compartment of the FLT3-ITD+ AML BM niche pre and post treatment and their impact on leukemic stem cells (LSCs). METHODS: BM vasculature was evaluated in FLT3-ITD+ AML models (MllPTD/WT/Flt3ITD/ITD mouse and patient-derived xenograft) by 3D confocal imaging of long bones, calvarium vascular permeability assays, and flow cytometry analysis. Cytokine levels were measured by Luminex assay and miR-126 levels evaluated by Q-RT-PCR and miRNA staining. Wild-type (wt) and MllPTD/WT/Flt3ITD/ITD mice with endothelial cell (EC) miR-126 knockout or overexpression served as controls. The impact of treatment-induced BM vascular changes on LSC activity was evaluated by secondary transplantation of BM cells after administration of tyrosine kinase inhibitors (TKIs) to MllPTD/WT/Flt3ITD/ITD mice with/without either EC miR-126 KO or co-treatment with tumor necrosis factor alpha (TNFα) or anti-miR-126 miRisten. RESULTS: In the normal BM niche, CD31+Sca-1high ECs lining arterioles have miR-126 levels higher than CD31+Sca-1low ECs lining sinusoids. We noted that during FLT3-ITD+ AML growth, the BM niche lost arterioles and gained sinusoids. These changes were mediated by TNFα, a cytokine produced by AML blasts, which induced EC miR-126 downregulation and caused depletion of CD31+Sca-1high ECs and gain in CD31+Sca-1low ECs. Loss of miR-126high ECs led to a decreased EC miR-126 supply to LSCs, which then entered the cell cycle and promoted leukemia growth. Accordingly, antileukemic treatment with TKI decreased the BM blast-produced TNFα and increased miR-126high ECs and the EC miR-126 supply to LSCs. High miR-126 levels safeguarded LSCs, as shown by more severe disease in secondary transplanted mice. Conversely, EC miR-126 deprivation via genetic or pharmacological EC miR-126 knock-down prevented treatment-induced BM miR-126high EC expansion and in turn LSC protection. CONCLUSIONS: Treatment-induced CD31+Sca-1high EC re-vascularization of the leukemic BM niche may represent a LSC extrinsic mechanism of treatment resistance that can be overcome with therapeutic EC miR-126 deprivation.

Vascularity and Dynamic Contrast-Enhanced Breast Magnetic Resonance Imaging.

Angiogenesis is a key step in the initiation and progression of an invasive breast cancer. High microvessel density by morphological characterization predicts metastasis and poor survival in women with invasive breast cancers. However, morphologic characterization is subject to variability and only can evaluate a limited portion of an invasive breast cancer. Consequently, breast Magnetic Resonance Imaging (MRI) is currently being evaluated to assess vascularity. Recently, through the new field of radiomics, dynamic contrast enhanced (DCE)-MRI is being used to evaluate vascular density, vascular morphology, and detection of aggressive breast cancer biology. While DCE-MRI is a highly sensitive tool, there are specific features that limit computational evaluation of blood vessels. These include (1) DCE-MRI evaluates gadolinium contrast and does not directly evaluate biology, (2) the resolution of DCE-MRI is insufficient for imaging small blood vessels, and (3) DCE-MRI images are very difficult to co-register. Here we review computational approaches for detection and analysis of blood vessels in DCE-MRI images and present some of the strategies we have developed for co-registry of DCE-MRI images and early detection of vascularization.

State-Transition Analysis of Time-Sequential Gene Expression Identifies Critical Points That Predict Development of Acute Myeloid Leukemia.

Temporal dynamics of gene expression inform cellular and molecular perturbations associated with disease development and evolution. Given the complexity of high-dimensional temporal genomic data, an analytic framework guided by a robust theory is needed to interpret time-sequential changes and to predict system dynamics. Here we model temporal dynamics of the transcriptome of peripheral blood mononuclear cells in a two-dimensional state-space representing states of health and leukemia using time-sequential bulk RNA-seq data from a murine model of acute myeloid leukemia (AML). The state-transition model identified critical points that accurately predict AML development and identifies stepwise transcriptomic perturbations that drive leukemia progression. The geometry of the transcriptome state-space provided a biological interpretation of gene dynamics, aligned gene signals that are not synchronized in time across mice, and allowed quantification of gene and pathway contributions to leukemia development. Our state-transition model synthesizes information from multiple cell types in the peripheral blood and identifies critical points in the transition from health to leukemia to guide interpretation of changes in the transcriptome as a whole to predict disease progression. SIGNIFICANCE: These findings apply the theory of state transitions to model the initiation and development of acute myeloid leukemia, identifying transcriptomic perturbations that accurately predict time to disease development.See related commentary by Kuijjer, p. 3072 GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/15/3157/F1.large.jpg.

Methylomic correlates of autophagy activity in cystic fibrosis.

Autophagy is a highly regulated, biological process that provides energy during periods of stress and starvation. This conserved process also acts as a defense mechanism and clears microbes from the host cell. Autophagy is impaired in Cystic Fibrosis (CF) patients and CF mice, as their cells exhibit low expression levels of essential autophagy molecules. The genetic disorder in CF is due to mutations in the cystic fibrosis transmembrane conductance regulator (cftr) gene that encodes for a chloride channel. CF patients are particularly prone to infection by pathogens that are otherwise cleared by autophagy in healthy immune cells including Burkholderia cenocepacia (B. cenocepacia). The objective of this study is to determine the mechanism underlying weak autophagic activity in CF macrophages and find therapeutic targets to correct it. Using reduced representation bisulfite sequencing (RRBS) to determine DNA methylation profile, we found that the promoter regions of Atg12 in CF macrophages are significantly more methylated than in the wild-type (WT) immune cells, accompanied by low protein expression. The natural product epigallocatechin-3-gallate (EGCG) significantly reduced the methylation of Atg12 promoter improving its expression. Accordingly, EGCG restricted B. cenocepacia replication within CF mice and their derived macrophages by improving autophagy and preventing dissemination. In addition, EGCG improved the function of CFTR protein. Altogether, utilizing RRBS for the first time in the CF field revealed a previously unrecognized mechanism for reduced autophagic activity in CF. Our data also offers a mechanism by which EGCG exerts its positive effects in CF.

FuSpot: a web-based tool for visual evaluation of fusion candidates.

BACKGROUND: Gene fusions often occur in cancer cells and in some cases are the main driver of oncogenesis. Correct identification of oncogenic gene fusions thus has implications for targeted cancer therapy. Recognition of this potential has led to the development of a myriad of sequencing-based fusion detection tools. However, given the same input, many of these detectors will find different fusion points or claim different sets of supporting data. Furthermore, the rate at which these tools falsely detect fusion events in data varies greatly. This discrepancy between tools underscores the fact that computation algorithms still cannot perfectly evaluate evidence; especially when provided with small amounts of supporting data as is typical in fusion detection. We assert that when evidence is provided in an easily digestible form, humans are more proficient in identifying true positives from false positives. RESULTS: We have developed a web tool that, given the genomic coordinates of a candidate fusion breakpoint, will extract fusion and non-fusion reads adjacent to the fusion point from partner transcripts, and color code reads by transcript origin and read orientation for ease of intuitive inspection by the user. Fusion partner transcript read alignments are performed using a novel variant of the Smith-Waterman algorithm. CONCLUSIONS: Combined with dynamic filtering parameters, the visualization provided by our tool introduces a powerful new investigative step that allows researchers to comprehensively evaluate fusion evidence. Additionally, this allows quick identification of false positives that may deceive most fusion detectors, thus eliminating unnecessary gene fusion validation. We apply our visualization tool to publicly available datasets and provide examples of true as well as false positives reported by open source fusion detection tools.

Statistical methods for detecting differentially methylated regions based on MethylCap-seq data.

DNA methylation is a well-established epigenetic mark, whose pattern throughout the genome, especially in the promoter or CpG islands, may be modified in a cell at a disease stage. Recently developed probabilistic approaches allow distributing methylation signals at nucleotide resolution from MethylCap-seq data. Standard statistical methods for detecting differential methylation suffer from 'curse of dimensionality' and sparsity in signals, resulting in high false-positive rates. Strong correlation of signals between CG sites also yields spurious results. In this article, we review applicability of high-dimensional mean vector tests for detection of differentially methylated regions (DMRs) and compare and contrast such tests with other methods for detecting DMRs. Comprehensive simulation studies are conducted to highlight the performance of these tests under different settings. Based on our observation, we make recommendations on the optimal test to use. We illustrate the superiority of mean vector tests in detecting cancer-related canonical gene pathways, which are significantly enriched for acute myeloid leukemia and ovarian cancer.

Preclinical Investigation of the Novel Histone Deacetylase Inhibitor AR-42 in the Treatment of Cancer-Induced Cachexia.

BACKGROUND: Cancer cachexia is a debilitating condition that impacts patient morbidity, mortality, and quality of life and for which effective therapies are lacking. The anticachectic activity of the novel HDAC inhibitor AR-42 was investigated in murine models of cancer cachexia. METHODS: The effects of AR-42 on classic features of cachexia were evaluated in the C-26 colon adenocarcinoma and Lewis lung carcinoma (LLC) models. Effects on survival in comparison with approved HDAC inhibitors (vorinostat, romidepsin) were determined. The muscle metabolome and transcriptome (by RNA-seq), as well as serum cytokine profile, were evaluated. Data were analyzed using mixed effects models, analysis of variance, or log-rank tests. All statistical tests were two-sided. RESULTS: In the C-26 model, orally administered AR-42 preserved body weight (23.9±2.6 grams, AR-42-treated; 20.8±1.3 grams, vehicle-treated; P = .005), prolonged survival (P < .001), prevented reductions in muscle and adipose tissue mass, muscle fiber size, and muscle strength and restored intramuscular mRNA expression of the E3 ligases MuRF1 and Atrogin-1 to basal levels (n = 8). This anticachectic effect, confirmed in the LLC model, was not observed after treatment with vorinostat and romidepsin. AR-42 suppressed tumor-induced changes in inflammatory cytokine production and multiple procachexia drivers (IL-6, IL-6Rα, leukemia inhibitory factor, Foxo1, Atrogin-1, MuRF1, adipose triglyceride lipase, uncoupling protein 3, and myocyte enhancer factor 2c). Metabolomic analysis revealed cachexia-associated changes in glycolysis, glycogen synthesis, and protein degradation in muscle, which were restored by AR-42 to a state characteristic of tumor-free mice. CONCLUSIONS: These findings support further investigation of AR-42 as part of a comprehensive therapeutic strategy for cancer cachexia.

Quality Control for RNA-Seq (QuaCRS): An Integrated Quality Control Pipeline.

QuaCRS (Quality Control for RNA-Seq) is an integrated, simplified quality control (QC) system for RNA-seq data that allows easy execution of several open-source QC tools, aggregation of their output, and the ability to quickly identify quality issues by performing meta-analyses on QC metrics across large numbers of samples in different studies. It comprises two main sections. First is the QC Pack wrapper, which executes three QC tools: FastQC, RNA-SeQC, and selected functions from RSeQC. Combining these three tools into one wrapper provides increased ease of use and provides a much more complete view of sample data quality than any individual tool. Second is the QC database, which displays the resulting metrics in a user-friendly web interface. It was designed to allow users with less computational experience to easily generate and view QC information for their data, to investigate individual samples and aggregate reports of sample groups, and to sort and search samples based on quality. The structure of the QuaCRS database is designed to enable expansion with additional tools and metrics in the future. The source code for not-for-profit use and a fully functional sample user interface with mock data are available at http://bioserv.mps.ohio-state.edu/QuaCRS/.

PrEMeR-CG: inferring nucleotide level DNA methylation values from MethylCap-seq data.

MOTIVATION: DNA methylation is an epigenetic change occurring in genomic CpG sequences that contribute to the regulation of gene transcription both in normal and malignant cells. Next-generation sequencing has been used to characterize DNA methylation status at the genome scale, but suffers from high sequencing cost in the case of whole-genome bisulfite sequencing, or from reduced resolution (inability to precisely define which of the CpGs are methylated) with capture-based techniques. RESULTS: Here we present a computational method that computes nucleotide-resolution methylation values from capture-based data by incorporating fragment length profiles into a model of methylation analysis. We demonstrate that it compares favorably with nucleotide-resolution bisulfite sequencing and has better predictive power with respect to a reference than window-based methods, often used for enrichment data. The described method was used to produce the methylation data used in tandem with gene expression to produce a novel and clinically significant gene signature in acute myeloid leukemia. In addition, we introduce a complementary statistical method that uses this nucleotide-resolution methylation data for detection of differentially methylated features.