Integrative analysis of drug response and clinical outcome in acute myeloid leukemia.

Acute myeloid leukemia (AML) is a cancer of myeloid-lineage cells with limited therapeutic options. We previously combined ex vivo drug sensitivity with genomic, transcriptomic, and clinical annotations for a large cohort of AML patients, which facilitated discovery of functional genomic correlates. Here, we present a dataset that has been harmonized with our initial report to yield a cumulative cohort of 805 patients (942 specimens). We show strong cross-cohort concordance and identify features of drug response. Further, deconvoluting transcriptomic data shows that drug sensitivity is governed broadly by AML cell differentiation state, sometimes conditionally affecting other correlates of response. Finally, modeling of clinical outcome reveals a single gene, PEAR1, to be among the strongest predictors of patient survival, especially for young patients. Collectively, this report expands a large functional genomic resource, offers avenues for mechanistic exploration and drug development, and reveals tools for predicting outcome in AML.

Genomic landscape of neutrophilic leukemias of ambiguous diagnosis.

Chronic neutrophilic leukemia (CNL), atypical chronic myeloid leukemia (aCML), and myelodysplastic/myeloproliferative neoplasms, unclassifiable (MDS/MPN-U) are a group of rare and heterogeneous myeloid disorders. There is strong morphologic resemblance among these distinct diagnostic entities as well as a lack of specific molecular markers and limited understanding of disease pathogenesis, which has made diagnosis challenging in certain cases. The treatment has remained empirical, resulting in dismal outcomes. We, therefore, performed whole-exome and RNA sequencing of these rare hematologic malignancies and present the most complete survey of the genomic landscape of these diseases to date. We observed a diversity of combinatorial mutational patterns that generally do not cluster within any one diagnosis. Gene expression analysis reveals enrichment, but not cosegregation, of clinical and genetic disease features with transcriptional clusters. In conclusion, these groups of diseases represent a continuum of related diseases rather than discrete diagnostic entities.

Clinical resistance to crenolanib in acute myeloid leukemia due to diverse molecular mechanisms.

FLT3 mutations are prevalent in AML patients and confer poor prognosis. Crenolanib, a potent type I pan-FLT3 inhibitor, is effective against both internal tandem duplications and resistance-conferring tyrosine kinase domain mutations. While crenolanib monotherapy has demonstrated clinical benefit in heavily pretreated relapsed/refractory AML patients, responses are transient and relapse eventually occurs. Here, to investigate the mechanisms of crenolanib resistance, we perform whole exome sequencing of AML patient samples before and after crenolanib treatment. Unlike other FLT3 inhibitors, crenolanib does not induce FLT3 secondary mutations, and mutations of the FLT3 gatekeeper residue are infrequent. Instead, mutations of NRAS and IDH2 arise, mostly as FLT3-independent subclones, while TET2 and IDH1 predominantly co-occur with FLT3-mutant clones and are enriched in crenolanib poor-responders. The remaining patients exhibit post-crenolanib expansion of mutations associated with epigenetic regulators, transcription factors, and cohesion factors, suggesting diverse genetic/epigenetic mechanisms of crenolanib resistance. Drug combinations in experimental models restore crenolanib sensitivity.

Functional genomic landscape of acute myeloid leukaemia.

The implementation of targeted therapies for acute myeloid leukaemia (AML) has been challenging because of the complex mutational patterns within and across patients as well as a dearth of pharmacologic agents for most mutational events. Here we report initial findings from the Beat AML programme on a cohort of 672 tumour specimens collected from 562 patients. We assessed these specimens using whole-exome sequencing, RNA sequencing and analyses of ex vivo drug sensitivity. Our data reveal mutational events that have not previously been detected in AML. We show that the response to drugs is associated with mutational status, including instances of drug sensitivity that are specific to combinatorial mutational events. Integration with RNA sequencing also revealed gene expression signatures, which predict a role for specific gene networks in the drug response. Collectively, we have generated a dataset-accessible through the Beat AML data viewer (Vizome)-that can be leveraged to address clinical, genomic, transcriptomic and functional analyses of the biology of AML.

Congenic mice provide in vivo evidence for a genetic locus that modulates intrinsic transforming growth factor ß1-mediated signaling and bone acquisition.

Osteoporosis, the most common skeletal disorder, is characterized by low bone mineral density (BMD) and an increased risk of fragility fractures. BMD is the best clinical predictor of future osteoporotic fracture risk, but is a complex trait controlled by multiple environmental and genetic determinants with individually modest effects. Quantitative trait locus (QTL) mapping is a powerful method for identifying chromosomal regions encompassing genes involved in shaping complex phenotypes, such as BMD. Here we have applied QTL analysis to male and female genetically-heterogeneous F(2) mice derived from a cross between C57BL/6 and DBA/2 strains, and have identified 11 loci contributing to femoral BMD. Further analysis of a QTL on mouse chromosome 7 following the generation of reciprocal congenic strains has allowed us to determine that the high BMD trait, which tracks with the DBA/2 chromosome and exerts equivalent effects on male and female mice, is manifested by enhanced osteogenic differentiation of mesenchymal stem cells (MSCs) in vitro and by increased growth of metatarsal bones in short-term primary culture. An insertion/deletion DNA polymorphism in Ltbp4 exon 12 that causes the in-frame removal of 12 codons in the DBA/2-derived gene maps within 0.6 Mb of the marker most tightly linked to the QTL. LTBP4, one of four paralogous mouse proteins that modify the bioavailability of the transforming growth factor ß (TGF-ß) family of growth factors, is expressed in differentiating MSC-derived osteoblasts and in long bones, and reduced responsiveness to TGF-ß1 is observed in MSCs of mice homozygous for the DBA/2 chromosome 7. Taken together, our results identify a potential genetic and biochemical relationship between decreased TGF-ß1-mediated signaling and enhanced femoral BMD that may be regulated by a variant LTBP4 molecule.

Rare coding variants in ALPL are associated with low serum alkaline phosphatase and low bone mineral density.

Alkaline phosphatase (ALP) plays an essential role in the regulation of tissue mineralization, and its activity is highly heritable. Guided by genetic associations discovered in a murine model, we hypothesized a role for rare coding variants in determining serum ALP level and bone mineral density (BMD) in humans. We sequenced the coding regions of the ALP gene (ALPL) in men with low and normal serum ALP activity levels. Single-nucleotide ALPL variants, including 19 rare nonsynonymous variants (minor allele frequency <1%), were much more frequent among the low ALP group (33.8%) than the normal group (1.4%, p = 1 × 10(-11)). Within the low ALP group, men with a rare, nonsynonymous variant had 11.2% lower mean serum ALP (p = 3.9 × 10(-4)), 6.7% lower BMD (p = 0.03), and 11.1% higher serum phosphate (p = 0.002) than those without. In contrast, common nonsynonymous variants had no association with serum ALP, phosphate, or BMD. Multiple rare ALPL coding variants are present in the general population, and nonsynonymous coding variants may be responsible for heritable differences in mineralization and thus BMD.

Regulation of bone mass in mice by the lipoxygenase gene Alox15.

The development of osteoporosis involves the interaction of multiple environmental and genetic factors. Through combined genetic and genomic approaches, we identified the lipoxygenase gene Alox15 as a negative regulator of peak bone mineral density in mice. Crossbreeding experiments with Alox15 knockout mice confirmed that 12/15-lipoxygenase plays a role in skeletal development. Pharmacologic inhibitors of this enzyme improved bone density and strength in two rodent models of osteoporosis. These results suggest that drugs targeting the 12/15-lipoxygenase pathway merit investigation as a therapy for osteoporosis.

Mapping quantitative trait loci that influence femoral cross-sectional area in mice.

Size and shape are critical determinants of the mechanical properties of skeletal elements and can be anticipated to be highly heritable. Moreover, the genes responsible may be independent of those that regulate bone mineral density (BMD). To begin to identify the heritable determinants of skeletal geometry, we have examined femoral cross-sectional area (FCSA) in male and female mice from two inbred strains of mice with divergent FCSA (C57BL/6 [B6] and DBA/2 [D2]), a large genetically heterogeneous population (n = 964) of B6D2F2 mice and 18 BXD recombinant inbred (RI) strains derived from their F2 cross. Femora were harvested from 16-week-old mice and FCSA (bone and marrow space enclosed within the periosteum) was measured at the midshaft by digital image analysis. In all mouse populations examined, FCSA was positively correlated with body weight and weight-corrected FCSA (WC-FCSA) values were normally distributed in the BXD-RI and F2 populations, suggesting polygenic control of this trait. Genome-wide quantitative trait locus (QTL) analysis of the B6D2F2 population revealed regions on four different chromosomes that were very strongly linked to WC-FCSA (chromosomes 6, 8, 10, and X) in both genders. Evidence of gender-specific genetic influences on femoral geometry was also identified at three other chromosomal sites (chromosomes 2, 7, and 12). Supporting evidence for the WC-FCSA QTLs on chromosomes 2, 7, 8, 10, and 12 also was present in the RI strains. Interestingly, none of these WC-FCSA QTLs were identified in our previous QTL analysis of whole body BMD in the same B6D2F2 population. Thus, the genetic determinants of bone size appear to be largely, if not entirely, distinct from those that regulate BMD attainment. The identification of the genes responsible for geometric differences in bone development should reveal fundamentally important processes in the control of skeletal integrity.