An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia.

Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

Chromatin Interaction Analysis with Paired-End Tag Sequencing (ChIA-PET) for mapping chromatin interactions and understanding transcription regulation.

Genomes are organized into three-dimensional structures, adopting higher-order conformations inside the micron-sized nuclear spaces (7, 2, 12). Such architectures are not random and involve interactions between gene promoters and regulatory elements (13). The binding of transcription factors to specific regulatory sequences brings about a network of transcription regulation and coordination (1, 14). Chromatin Interaction Analysis by Paired-End Tag Sequencing (ChIA-PET) was developed to identify these higher-order chromatin structures (5,6). Cells are fixed and interacting loci are captured by covalent DNA-protein cross-links. To minimize non-specific noise and reduce complexity, as well as to increase the specificity of the chromatin interaction analysis, chromatin immunoprecipitation (ChIP) is used against specific protein factors to enrich chromatin fragments of interest before proximity ligation. Ligation involving half-linkers subsequently forms covalent links between pairs of DNA fragments tethered together within individual chromatin complexes. The flanking MmeI restriction enzyme sites in the half-linkers allow extraction of paired end tag-linker-tag constructs (PETs) upon MmeI digestion. As the half-linkers are biotinylated, these PET constructs are purified using streptavidin-magnetic beads. The purified PETs are ligated with next-generation sequencing adaptors and a catalog of interacting fragments is generated via next-generation sequencers such as the Illumina Genome Analyzer. Mapping and bioinformatics analysis is then performed to identify ChIP-enriched binding sites and ChIP-enriched chromatin interactions (8). We have produced a video to demonstrate critical aspects of the ChIA-PET protocol, especially the preparation of ChIP as the quality of ChIP plays a major role in the outcome of a ChIA-PET library. As the protocols are very long, only the critical steps are shown in the video.

Extensive promoter-centered chromatin interactions provide a topological basis for transcription regulation.

Higher-order chromosomal organization for transcription regulation is poorly understood in eukaryotes. Using genome-wide Chromatin Interaction Analysis with Paired-End-Tag sequencing (ChIA-PET), we mapped long-range chromatin interactions associated with RNA polymerase II in human cells and uncovered widespread promoter-centered intragenic, extragenic, and intergenic interactions. These interactions further aggregated into higher-order clusters, wherein proximal and distal genes were engaged through promoter-promoter interactions. Most genes with promoter-promoter interactions were active and transcribed cooperatively, and some interacting promoters could influence each other implying combinatorial complexity of transcriptional controls. Comparative analyses of different cell lines showed that cell-specific chromatin interactions could provide structural frameworks for cell-specific transcription, and suggested significant enrichment of enhancer-promoter interactions for cell-specific functions. Furthermore, genetically-identified disease-associated noncoding elements were found to be spatially engaged with corresponding genes through long-range interactions. Overall, our study provides insights into transcription regulation by three-dimensional chromatin interactions for both housekeeping and cell-specific genes in human cells.

Dynamic changes in the human methylome during differentiation.

DNA methylation is a critical epigenetic regulator in mammalian development. Here, we present a whole-genome comparative view of DNA methylation using bisulfite sequencing of three cultured cell types representing progressive stages of differentiation: human embryonic stem cells (hESCs), a fibroblastic differentiated derivative of the hESCs, and neonatal fibroblasts. As a reference, we compared our maps with a methylome map of a fully differentiated adult cell type, mature peripheral blood mononuclear cells (monocytes). We observed many notable common and cell-type-specific features among all cell types. Promoter hypomethylation (both CG and CA) and higher levels of gene body methylation were positively correlated with transcription in all cell types. Exons were more highly methylated than introns, and sharp transitions of methylation occurred at exon-intron boundaries, suggesting a role for differential methylation in transcript splicing. Developmental stage was reflected in both the level of global methylation and extent of non-CpG methylation, with hESC highest, fibroblasts intermediate, and monocytes lowest. Differentiation-associated differential methylation profiles were observed for developmentally regulated genes, including the HOX clusters, other homeobox transcription factors, and pluripotence-associated genes such as POU5F1, TCF3, and KLF4. Our results highlight the value of high-resolution methylation maps, in conjunction with other systems-level analyses, for investigation of previously undetectable developmental regulatory mechanisms.

GRIM-19, a death-regulatory gene product, suppresses Stat3 activity via functional interaction.

Signal transducer and activator of transcription 3 (Stat3) is a latent cytoplasmic transcription factor that can be activated by cytokines and growth factors. Stat3 plays important roles in cell growth, anti-apoptosis and cell transformation, and is constitutively active in various cancers. We examined its potential regulators by yeast two-hybrid screening. GRIM-19, a gene product related to interferon-beta- and retinoic acid-induced cancer cell death, was identified and demonstrated to interact with Stat3 in various cell types. The interaction is specific for Stat3, but not for Stat1 and Stat5a. The interaction regions in both proteins were mapped, and the cellular localization of the interaction was examined. GRIM-19 itself co-localizes with mitochondrial markers, and forms aggregates at the perinulear region with co-expressed Stat3, which inhibits Stat3 nuclear translocation stimulated by epidermal growth factor (EGF). GRIM-19 represses Stat3 transcriptional activity and its target gene expression, and also suppresses cell growth in Src-transformed cells and a Stat3-expressing cell line. Our data suggest that GRIM-19 is a novel negative regulator of Stat3.

Interdomain interaction of Stat3 regulates its Src homology 2 domain-mediated receptor binding activity.

Activation of Stat proteins by cytokines is initiated by their Src homology 2 (SH2) domain-mediated association with the cytokine receptors. Previously, we identified an essential role of the coiled-coil domain of Stat3 in binding of the receptor peptides derived from the interleukin-6 receptor subunit, gp130. In this study, we further investigated the molecular basis of this regulation. We found that the C-terminal domain of Stat3 negatively regulates its receptor binding activity only in the absence of the first alpha-helix of the coiled-coil domain, which leads to a hypothesis of intramolecular interaction. Physical interactions between the coiled-coil domain and the C-terminal domain, as well as the SH2 domain, were indeed detected. Furthermore, a sub-region of the C-terminal domain (amino acids 720-740), which is also involved in the interaction with the coiled-coil domain, was demonstrated to be critical for the regulation of the receptor binding. Correspondingly, phosphorylation on Ser-727 within this region inhibits this interaction. In agreement with the peptide binding results, both the coiled-coil domain and the C-terminal sub-region are necessary for the functional recruitment of Stat3 to the cellular gp130 in response to interleukin-6, suggesting that the interdomain interaction is a prerequisite for the SH2-mediated receptor binding in interleukin-6 signaling.