The genomic basis of childhood T-lineage acute lymphoblastic leukaemia.

T-lineage acute lymphoblastic leukaemia (T-ALL) is a high-risk tumour1 that has eluded comprehensive genomic characterization, which is partly due to the high frequency of noncoding genomic alterations that result in oncogene deregulation2,3. Here we report an integrated analysis of genome and transcriptome sequencing of tumour and remission samples from more than 1,300 uniformly treated children with T-ALL, coupled with epigenomic and single-cell analyses of malignant and normal T cell precursors. This approach identified 15 subtypes with distinct genomic drivers, gene expression patterns, developmental states and outcomes. Analyses of chromatin topology revealed multiple mechanisms of enhancer deregulation that involve enhancers and genes in a subtype-specific manner, thereby demonstrating widespread involvement of the noncoding genome. We show that the immunophenotypically described, high-risk entity of early T cell precursor ALL is superseded by a broader category of 'early T cell precursor-like' leukaemia. This category has a variable immunophenotype and diverse genomic alterations of a core set of genes that encode regulators of hematopoietic stem cell development. Using multivariable outcome models, we show that genetic subtypes, driver and concomitant genetic alterations independently predict treatment failure and survival. These findings provide a roadmap for the classification, risk stratification and mechanistic understanding of this disease.

Biologic and Clinical Analysis of Childhood Gamma Delta T-ALL Identifies LMO2/STAG2 Rearrangements as Extremely High Risk.

Acute lymphoblastic leukemia expressing the gamma delta T-cell receptor (γδ T-ALL) is a poorly understood disease. We studied 200 children with γδ T-ALL from 13 clinical study groups to understand the clinical and genetic features of this disease. We found age and genetic drivers were significantly associated with outcome. γδ T-ALL diagnosed in children under 3 years of age was extremely high-risk and enriched for genetic alterations that result in both LMO2 activation and STAG2 inactivation. Mechanistically, using patient samples and isogenic cell lines, we show that inactivation of STAG2 profoundly perturbs chromatin organization by altering enhancer-promoter looping, resulting in deregulation of gene expression associated with T-cell differentiation. High-throughput drug screening identified a vulnerability in DNA repair pathways arising from STAG2 inactivation, which can be targeted by poly(ADP-ribose) polymerase inhibition. These data provide a diagnostic framework for classification and risk stratification of pediatric γδ T-ALL. Significance: Patients with acute lymphoblastic leukemia expressing the gamma delta T-cell receptor under 3 years old or measurable residual disease ≥1% at end of induction showed dismal outcomes and should be classified as having high-risk disease. The STAG2/LMO2 subtype was enriched in this very young age group. STAG2 inactivation may perturb chromatin conformation and cell differentiation and confer vulnerability to poly(ADP-ribose) polymerase inhibition.

Redefining the biological basis of lineage-ambiguous leukemia through genomics: BCL11B deregulation in acute leukemias of ambiguous lineage.

Acute leukemias of ambiguous lineage (ALAL), including mixed phenotype acute leukemia (MPAL) and related entities such as early T-cell precursor acute leukemia (ETP-ALL), remain diagnostic and clinical challenges due to limited understanding of pathogenesis, reliance of immunophenotyping to classify disease, and the lack of a rational approach to guide selection of appropriate therapy. Recent studies utilizing genomic sequencing and complementary approaches have provided key insights that are changing the way in which such leukemias are classified, and potentially, treated. Several recurrent genomic alterations define leukemias that straddle immunophenotypic entities, such as ZNF384-rearranged childhood B-ALL and B/myeloid MPAL, and BCL11B-rearranged T/myeloid MPAL, ETP-ALL and AML. In contrast, some cases of MPAL represent canonical ALL/AML entities exhibiting lineage aberrancy. For many cases of ALAL, experimental approaches indicate lineage aberrancy arises from acquisition of a founding genetic alteration into a hematopoietic stem or progenitor cell. Determination of optimal therapeutic approach requires genomic characterization of uniformly treated ALAL patients in prospective studies, but several approaches, including kinase inhibitors and BH3 mimetics may be efficacious in subsets of ALAL.

A functional genomics pipeline identifies pleiotropy and cross-tissue effects within obesity-associated GWAS loci.

Genome-wide association studies (GWAS) have identified many disease-associated variants, yet mechanisms underlying these associations remain unclear. To understand obesity-associated variants, we generate gene regulatory annotations in adipocytes and hypothalamic neurons across cellular differentiation stages. We then test variants in 97 obesity-associated loci using a massively parallel reporter assay and identify putatively causal variants that display cell type specific or cross-tissue enhancer-modulating properties. Integrating these variants with gene regulatory information suggests genes that underlie obesity GWAS associations. We also investigate a complex genomic interval on 16p11.2 where two independent loci exhibit megabase-range, cross-locus chromatin interactions. We demonstrate that variants within these two loci regulate a shared gene set. Together, our data support a model where GWAS loci contain variants that alter enhancer activity across tissues, potentially with temporally restricted effects, to impact the expression of multiple genes. This complex model has broad implications for ongoing efforts to understand GWAS.

Asthma-associated genetic variants induce IL33 differential expression through an enhancer-blocking regulatory region.

Genome-wide association studies (GWAS) have implicated the IL33 locus in asthma, but the underlying mechanisms remain unclear. Here, we identify a 5 kb region within the GWAS-defined segment that acts as an enhancer-blocking element in vivo and in vitro. Chromatin conformation capture showed that this 5 kb region loops to the IL33 promoter, potentially regulating its expression. We show that the asthma-associated single nucleotide polymorphism (SNP) rs1888909, located within the 5 kb region, is associated with IL33 gene expression in human airway epithelial cells and IL-33 protein expression in human plasma, potentially through differential binding of OCT-1 (POU2F1) to the asthma-risk allele. Our data demonstrate that asthma-associated variants at the IL33 locus mediate allele-specific regulatory activity and IL33 expression, providing a mechanism through which a regulatory SNP contributes to genetic risk of asthma.

Enhancer Hijacking Drives Oncogenic BCL11B Expression in Lineage-Ambiguous Stem Cell Leukemia.

Lineage-ambiguous leukemias are high-risk malignancies of poorly understood genetic basis. Here, we describe a distinct subgroup of acute leukemia with expression of myeloid, T lymphoid, and stem cell markers driven by aberrant allele-specific deregulation of BCL11B, a master transcription factor responsible for thymic T-lineage commitment and specification. Mechanistically, this deregulation was driven by chromosomal rearrangements that juxtapose BCL11B to superenhancers active in hematopoietic progenitors, or focal amplifications that generate a superenhancer from a noncoding element distal to BCL11B. Chromatin conformation analyses demonstrated long-range interactions of rearranged enhancers with the expressed BCL11B allele and association of BCL11B with activated hematopoietic progenitor cell cis-regulatory elements, suggesting BCL11B is aberrantly co-opted into a gene regulatory network that drives transformation by maintaining a progenitor state. These data support a role for ectopic BCL11B expression in primitive hematopoietic cells mediated by enhancer hijacking as an oncogenic driver of human lineage-ambiguous leukemia. SIGNIFICANCE: Lineage-ambiguous leukemias pose significant diagnostic and therapeutic challenges due to a poorly understood molecular and cellular basis. We identify oncogenic deregulation of BCL11B driven by diverse structural alterations, including de novo superenhancer generation, as the driving feature of a subset of lineage-ambiguous leukemias that transcend current diagnostic boundaries.This article is highlighted in the In This Issue feature, p. 2659.

Identification of atrial fibrillation associated genes and functional non-coding variants.

Disease-associated genetic variants that lie in non-coding regions found by genome-wide association studies are thought to alter the functionality of transcription regulatory elements and target gene expression. To uncover causal genetic variants, variant regulatory elements and their target genes, here we cross-reference human transcriptomic, epigenomic and chromatin conformation datasets. Of 104 genetic variant regions associated with atrial fibrillation candidate target genes are prioritized. We optimize EMERGE enhancer prediction and use accessible chromatin profiles of human atrial cardiomyocytes to more accurately predict cardiac regulatory elements and identify hundreds of sub-threshold variants that co-localize with regulatory elements. Removal of mouse homologues of atrial fibrillation-associated regions in vivo uncovers a distal regulatory region involved in Gja1 (Cx43) expression. Our analyses provide a shortlist of genes likely affected by atrial fibrillation-associated variants and provide variant regulatory elements in each region that link genetic variation and target gene regulation, helping to focus future investigations.

A promoter interaction map for cardiovascular disease genetics.

Over 500 genetic loci have been associated with risk of cardiovascular diseases (CVDs); however, most loci are located in gene-distal non-coding regions and their target genes are not known. Here, we generated high-resolution promoter capture Hi-C (PCHi-C) maps in human induced pluripotent stem cells (iPSCs) and iPSC-derived cardiomyocytes (CMs) to provide a resource for identifying and prioritizing the functional targets of CVD associations. We validate these maps by demonstrating that promoters preferentially contact distal sequences enriched for tissue-specific transcription factor motifs and are enriched for chromatin marks that correlate with dynamic changes in gene expression. Using the CM PCHi-C map, we linked 1999 CVD-associated SNPs to 347 target genes. Remarkably, more than 90% of SNP-target gene interactions did not involve the nearest gene, while 40% of SNPs interacted with at least two genes, demonstrating the importance of considering long-range chromatin interactions when interpreting functional targets of disease loci.