Franck Rapaport

Comprehensive evaluation of differential gene expression analysis methods for RNA-seq data

A large number of computational methods have been developed for analyzing differential gene expression in RNA-seq data. We describe a comprehensive evaluation of common methods using the SEQC benchmark dataset and ENCODE data. We consider a number of key features, including normalization, accuracy of differential expression detection and differential expression analysis when one condition has no detectable expression. We find significant differences among the methods, but note that array-based methods adapted to RNA-seq data perform comparably to methods designed for RNA-seq. Our results demonstrate that increasing the number of replicate samples significantly improves detection power over increased sequencing depth.

Increased GVHD-related mortality with broad-spectrum antibiotic use after allogeneic hematopoietic stem cell transplantation in human patients and mice

Treating neutropenic fever with broad-spectrum antibiotics after allogeneic hematopoietic stem cell transplant is associated with an increase in graft-versus-host disease in mice and humans. Antibiotics for allogeneic transplant—A double-edged sword Patients undergoing allogeneic hematopoietic stem cell transplantation often receive antibiotics for infections, which can also unfortunately kill intestinal bacteria. These symbiotic bacteria in the gut generally do not cause disease and are thought to suppress inflammation. In a new study, Shono et al. examined the records of 857 transplant patients and found that certain antibiotics were linked with development of graft-versus-host disease (GVHD), which can cause severe intestinal inflammation. Using a mouse model, the authors showed that these antibiotics may select for bacteria that consume intestinal mucus and lead to loss of this important layer of protection for the gut, thus exacerbating GVHD in the intestine. This study suggests that not all antibiotic regimens are appropriate for treating transplant patients. Intestinal bacteria may modulate the risk of infection and graft-versus-host disease (GVHD) after allogeneic hematopoietic stem cell transplantation (allo-HSCT). Allo-HSCT recipients often develop neutropenic fever, which is treated with antibiotics that may target anaerobic bacteria in the gut. We retrospectively examined 857 allo-HSCT recipients and found that treatment of neutropenic fever with imipenem-cilastatin and piperacillin-tazobactam antibiotics was associated with increased GVHD-related mortality at 5 years (21.5% for imipenem-cilastatin–treated patients versus 13.1% for untreated patients, P = 0.025; 19.8% for piperacillin-tazobactam–treated patients versus 11.9% for untreated patients, P = 0.007). However, two other antibiotics also used to treat neutropenic fever, aztreonam and cefepime, were not associated with GVHD-related mortality (P = 0.78 and P = 0.98, respectively). Analysis of stool specimens from allo-HSCT recipients showed that piperacillin-tazobactam administration was associated with perturbation of gut microbial composition. Studies in mice demonstrated aggravated GVHD mortality with imipenem-cilastatin or piperacillin-tazobactam compared to aztreonam (P < 0.01 and P < 0.05, respectively). We found pathological evidence for increased GVHD in the colon of imipenem-cilastatin–treated mice (P < 0.05), but no difference in the concentration of short-chain fatty acids or numbers of regulatory T cells. Notably, imipenem-cilastatin treatment of mice with GVHD led to loss of the protective mucus lining of the colon (P < 0.01) and the compromising of intestinal barrier function (P < 0.05). Sequencing of mouse stool specimens showed an increase in Akkermansia muciniphila (P < 0.001), a commensal bacterium with mucus-degrading capabilities, raising the possibility that mucus degradation may contribute to murine GVHD. We demonstrate an underappreciated risk for the treatment of allo-HSCT recipients with antibiotics that may exacerbate GVHD in the colon.

Mutational Cooperativity Linked to Combinatorial Epigenetic Gain of Function in Acute Myeloid Leukemia

Specific combinations of acute myeloid leukemia (AML) disease alleles, including FLT3 and TET2 mutations, confer distinct biologic features and adverse outcome. We generated mice with mutations in Tet2 and Flt3, which resulted in fully penetrant, lethal AML. Multipotent Tet2(-/-);Flt3(ITD) progenitors (LSK CD48(+)CD150(-)) propagate disease in secondary recipients and were refractory to standard AML chemotherapy and FLT3-targeted therapy. Flt3(ITD) mutations and Tet2 loss cooperatively remodeled DNA methylation and gene expression to an extent not seen with either mutant allele alone, including at the Gata2 locus. Re-expression of Gata2 induced differentiation in AML stem cells and attenuated leukemogenesis. TET2 and FLT3 mutations cooperatively induce AML, with a defined leukemia stem cell population characterized by site-specific changes in DNA methylation and gene expression.

Genomic and functional analysis of leukemic transformation of myeloproliferative neoplasms

Significance Myeloproliferative neoplasms (MPN) are chronic hematopoietic disorders characterized by clonal proliferation of mature myeloid elements. A subset of MPNs transforms to acute myeloid leukemia (AML). The mechanisms and pathways that contribute to transformation from MPN to AML have not been well delineated. We have characterized the somatic mutational spectrum of post-MPN AML and demonstrate that somatic tumor protein 53 (TP53) mutations are common in JAK2V617F-mutant, post-MPN AML but not in chronic-phase MPN. We demonstrate that expression of JAK2V617F combined with Tp53 loss in a murine model leads to fully penetrant AML in vivo. We have characterized this model and used it to test therapeutic strategies. These data reveal novel insights into the pathogenesis of, and potential therapeutic strategies for, leukemic transformation. Patients with myeloproliferative neoplasms (MPNs) are at significant, cumulative risk of leukemic transformation to acute myeloid leukemia (AML), which is associated with adverse clinical outcome and resistance to standard AML therapies. We performed genomic profiling of post-MPN AML samples; these studies demonstrate somatic tumor protein 53 (TP53) mutations are common in JAK2V617F-mutant, post-MPN AML but not in chronic-phase MPN and lead to clonal dominance of JAK2V617F/TP53-mutant leukemic cells. Consistent with these data, expression of JAK2V617F combined with Tp53 loss led to fully penetrant AML in vivo. JAK2V617F-mutant, Tp53-deficient AML was characterized by an expanded megakaryocyte erythroid progenitor population that was able to propagate the disease in secondary recipients. In vitro studies revealed that post-MPN AML cells were sensitive to decitabine, the JAK1/2 inhibitor ruxolitinib, or the heat shock protein 90 inhibitor 8-(6-iodobenzo[d][1.3]dioxol-5-ylthio)-9-(3-(isopropylamino)propyl)-9H-purine-6-amine (PU-H71). Treatment with ruxolitinib or PU-H71 improved survival of mice engrafted with JAK2V617F-mutant, Tp53-deficient AML, demonstrating therapeutic efficacy for these targeted therapies and providing a rationale for testing these therapies in post-MPN AML.

Distinct evolution and dynamics of epigenetic and genetic heterogeneity in acute myeloid leukemia

Genetic heterogeneity contributes to clinical outcome and progression of most tumors, but little is known about allelic diversity for epigenetic compartments, and almost no data exist for acute myeloid leukemia (AML). We examined epigenetic heterogeneity as assessed by cytosine methylation within defined genomic loci with four CpGs (epialleles), somatic mutations, and transcriptomes of AML patient samples at serial time points. We observed that epigenetic allele burden is linked to inferior outcome and varies considerably during disease progression. Epigenetic and genetic allelic burden and patterning followed different patterns and kinetics during disease progression. We observed a subset of AMLs with high epiallele and low somatic mutation burden at diagnosis, a subset with high somatic mutation and lower epiallele burdens at diagnosis, and a subset with a mixed profile, suggesting distinct modes of tumor heterogeneity. Genes linked to promoter-associated epiallele shifts during tumor progression showed increased single-cell transcriptional variance and differential expression, suggesting functional impact on gene regulation. Thus, genetic and epigenetic heterogeneity can occur with distinct kinetics likely to affect the biological and clinical features of tumors.

Dnmt3a regulates myeloproliferation and liver-specific expansion of hematopoietic stem and progenitor cells.

DNA methyltransferase 3A (DNMT3A) mutations are observed in myeloid malignancies, including myeloproliferative neoplasms (MPN), myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML). Transplantation studies have elucidated an important role for Dnmt3a in stem cell self-renewal and in myeloid differentiation. Here, we investigated the impact of conditional hematopoietic Dnmt3a loss on disease phenotype in primary mice. Mx1-Cre-mediated Dnmt3a ablation led to the development of a lethal, fully penetrant MPN with myelodysplasia (MDS/MPN) characterized by peripheral cytopenias and by marked, progressive hepatomegaly. We detected expanded stem/progenitor populations in the liver of Dnmt3a-ablated mice. The MDS/MPN induced by Dnmt3a ablation was transplantable, including the marked hepatomegaly. Homing studies showed that Dnmt3a-deleted bone marrow cells preferentially migrated to the liver. Gene expression and DNA methylation analyses of progenitor cell populations identified differential regulation of hematopoietic regulatory pathways, including fetal liver hematopoiesis transcriptional programs. These data demonstrate that Dnmt3a ablation in the hematopoietic system leads to myeloid transformation in vivo, with cell-autonomous aberrant tissue tropism and marked extramedullary hematopoiesis (EMH) with liver involvement. Hence, in addition to the established role of Dnmt3a in regulating self-renewal, Dnmt3a regulates tissue tropism and limits myeloid progenitor expansion in vivo.

Integrative genetic analysis of mouse and human AML identifies cooperating disease alleles

Hatlen et al. provide an integrative analysis of the mutational landscape of mouse and human AML and identify functionally relevant cooperation between AML1-ETO and PTPN11 D61Y. Based on these findings, they generate a novel mouse model of t(8;21)+ AML.

Determining frequent patterns of copy number alterations in cancer.

Cancer progression is often driven by an accumulation of genetic changes but also accompanied by increasing genomic instability. These processes lead to a complicated landscape of copy number alterations (CNAs) within individual tumors and great diversity across tumor samples. High resolution array-based comparative genomic hybridization (aCGH) is being used to profile CNAs of ever larger tumor collections, and better computational methods for processing these data sets and identifying potential driver CNAs are needed. Typical studies of aCGH data sets take a pipeline approach, starting with segmentation of profiles, calls of gains and losses, and finally determination of frequent CNAs across samples. A drawback of pipelines is that choices at each step may produce different results, and biases are propagated forward. We present a mathematically robust new method that exploits probe-level correlations in aCGH data to discover subsets of samples that display common CNAs. Our algorithm is related to recent work on maximum-margin clustering. It does not require pre-segmentation of the data and also provides grouping of recurrent CNAs into clusters. We tested our approach on a large cohort of glioblastoma aCGH samples from The Cancer Genome Atlas and recovered almost all CNAs reported in the initial study. We also found additional significant CNAs missed by the original analysis but supported by earlier studies, and we identified significant correlations between CNAs.

Somatic mutations in leukocytes infiltrating primary breast cancers

Background:Malignant transformation requires the interaction of cancer cells with their microenvironment, including infiltrating leukocytes. However, somatic mutational studies have focused on alterations in cancer cells, assuming that the microenvironment is genetically normal. Because we hypothesized that this might not be a valid assumption, we performed exome sequencing and targeted sequencing to investigate for the presence of pathogenic mutations in tumor-associated leukocytes in breast cancers.Methods:We used targeted sequencing and exome sequencing to evaluate the presence of mutations in sorted tumor-infiltrating CD45-positive cells from primary untreated breast cancers. We used high-depth sequencing to determine the presence/absence of the mutations we identified in breast cancer-infiltrating leukocytes in purified tumor cells and in circulating blood cells.Results:Capture-based sequencing of 15 paired tumor-infiltrating leukocytes and matched germline DNA identified variants in known cancer genes in all 15 primary breast cancer patients in our cohort. We validated the presence of mutations identified by targeted sequencing in infiltrating leukocytes through orthogonal exome sequencing. Ten patients harbored alterations previously reported as somatically acquired variants, including in known leukemia genes (DNTM3A, TET2, and BCOR). One of the mutations observed in the tumor-infiltrating leukocytes was also detected in the circulating leukocytes of the same patients at a lower allele frequency than observed in the tumor-infiltrating cells.Conclusions:Here we show that somatic mutations, including mutations in known cancer genes, are present in the leukocytes infiltrating a subset of primary breast cancers. This observation allows for the possibility that the cancer cells interact with mutant infiltrating leukocytes, which has many potential clinical implications.

Cooperative Epigenetic Remodeling by TET2 Loss and NRAS Mutation Drives Myeloid Transformation and MEK Inhibitor Sensitivity

Mutations in epigenetic modifiers and signaling factors often co-occur in myeloid malignancies, including TET2 and NRAS mutations. Concurrent Tet2 loss and NrasG12D expression in hematopoietic cells induced myeloid transformation, with a fully penetrant, lethal chronic myelomonocytic leukemia (CMML), which was serially transplantable. Tet2 loss and Nras mutation cooperatively led to decrease in negative regulators of mitogen-activated protein kinase (MAPK) activation, including Spry2, thereby causing synergistic activation of MAPK signaling by epigenetic silencing. Tet2/Nras double-mutant leukemia showed preferential sensitivity to MAPK kinase (MEK) inhibition in both mouse model and patient samples. These data provide insights into how epigenetic and signaling mutations cooperate in myeloid transformation and provide a rationale for mechanism-based therapy in CMML patients with these high-risk genetic lesions.