Jerry Davison

Genome-wide MyoD Binding in Skeletal Muscle Cells: A Potential for Broad Cellular Reprogramming

Recent studies have demonstrated that MyoD initiates a feed-forward regulation of skeletal muscle gene expression, predicting that MyoD binds directly to many genes expressed during differentiation. We have used chromatin immunoprecipitation and high-throughput sequencing to identify genome-wide binding of MyoD in several skeletal muscle cell types. As anticipated, MyoD preferentially binds to a VCASCTG sequence that resembles the in vitro-selected site for a MyoD:E-protein heterodimer, and MyoD binding increases during differentiation at many of the regulatory regions of genes expressed in skeletal muscle. Unanticipated findings were that MyoD was constitutively bound to thousands of additional sites in both myoblasts and myotubes, and that the genome-wide binding of MyoD was associated with regional histone acetylation. Therefore, in addition to regulating muscle gene expression, MyoD binds genome wide and has the ability to broadly alter the epigenome in myoblasts and myotubes.

Differential genomic targeting of the transcription factor TAL1 in alternate haematopoietic lineages

TAL1/SCL is a master regulator of haematopoiesis whose expression promotes opposite outcomes depending on the cell type: differentiation in the erythroid lineage or oncogenesis in the T‐cell lineage. Here, we used a combination of ChIP sequencing and gene expression profiling to compare the function of TAL1 in normal erythroid and leukaemic T cells. Analysis of the genome‐wide binding properties of TAL1 in these two haematopoietic lineages revealed new insight into the mechanism by which transcription factors select their binding sites in alternate lineages. Our study shows limited overlap in the TAL1‐binding profile between the two cell types with an unexpected preference for ETS and RUNX motifs adjacent to E‐boxes in the T‐cell lineage. Furthermore, we show that TAL1 interacts with RUNX1 and ETS1, and that these transcription factors are critically required for TAL1 binding to genes that modulate T‐cell differentiation. Thus, our findings highlight a critical role of the cellular environment in modulating transcription factor binding, and provide insight into the mechanism by which TAL1 inhibits differentiation leading to oncogenesis in the T‐cell lineage.

G9a/GLP-dependent histone H3K9me2 patterning during human hematopoietic stem cell lineage commitment

G9a and GLP are conserved protein methyltransferases that play key roles during mammalian development through mono- and dimethylation of histone H3 Lys 9 (H3K9me1/2), modifications associated with transcriptional repression. During embryogenesis, large H3K9me2 chromatin territories arise that have been proposed to reinforce lineage choice by affecting high-order chromatin structure. Here we report that in adult human hematopoietic stem and progenitor cells (HSPCs), H3K9me2 chromatin territories are absent in primitive cells and are formed de novo during lineage commitment. In committed HSPCs, G9a/GLP activity nucleates H3K9me2 marks at CpG islands and other genomic sites within genic regions, which then spread across most genic regions during differentiation. Immunofluorescence assays revealed the emergence of H3K9me2 nuclear speckles in committed HSPCs, consistent with progressive marking. Moreover, gene expression analysis indicated that G9a/GLP activity suppresses promiscuous transcription of lineage-affiliated genes and certain gene clusters, suggestive of regulation of HSPC chromatin structure. Remarkably, HSPCs continuously treated with UNC0638, a G9a/GLP small molecular inhibitor, better retain stem cell-like phenotypes and function during in vitro expansion. These results suggest that G9a/GLP activity promotes progressive H3K9me2 patterning during HSPC lineage specification and that its inhibition delays HSPC lineage commitment. They also inform clinical manipulation of donor-derived HSPCs.

Genome-wide CRISPR-Cas9 Screens Reveal Loss of Redundancy between PKMYT1 and WEE1 in Glioblastoma Stem-like Cells.

To identify therapeutic targets for glioblastoma (GBM), we performed genome-wide CRISPR-Cas9 knockout (KO) screens in patient-derived GBM stem-like cells (GSCs) and human neural stem/progenitors (NSCs), non-neoplastic stem cell controls, for genes required for their in vitro growth. Surprisingly, the vast majority GSC-lethal hits were found outside of molecular networks commonly altered in GBM and GSCs (e.g., oncogenic drivers). In vitro and in vivo validation of GSC-specific targets revealed several strong hits, including the wee1-like kinase, PKMYT1/Myt1. Mechanistic studies demonstrated that PKMYT1 acts redundantly with WEE1 to inhibit cyclin B-CDK1 activity via CDK1-Y15 phosphorylation and to promote timely completion of mitosis in NSCs. However, in GSCs, this redundancy is lost, most likely as a result of oncogenic signaling, causing GBM-specific lethality.

MYC-driven tumorigenesis is inhibited by WRN syndrome gene deficiency

MYC-induced DNA damage is exacerbated in WRN-deficient cells, leading to replication stress and accelerated cellular senescence. To determine whether WRN deficiency impairs MYC-driven tumor development, we used both xenograft and autochthonous tumor models. Conditional silencing of WRN expression in c-MYC overexpressing non–small cell lung cancer xenografts impaired both tumor establishment and tumor growth. This inhibitory effect of WRN knockdown was accompanied by increased DNA damage, decreased proliferation, and tumor necrosis. In the Eμ-Myc mouse model of B-cell lymphoma, a germline mutation in the helicase domain of Wrn (WrnΔhel/Δhel) resulted in a significant delay in emergence of lethal lymphomas, extending tumor-free survival by more than 30%. Analysis of preneoplastic B cells from Eμ-Myc Wrn mutant mice revealed increased DNA damage, elevation of senescence markers, and decreased proliferation in comparison with cells from age-matched Eμ-Myc mice. Immunohistochemical and global gene expression analysis of overt Eμ-Myc WrnΔhel/Δhel lymphomas showed a marked increase in expression of the CDK inhibitor, p16Ink4a, as well as elevation of TAp63, a known mediator of senescence. Collectively, these studies show that in the context of Myc-associated tumorigenesis, loss of Wrn amplifies the DNA damage response, both in preneoplastic and neoplastic tissue, engaging activation of tumor suppressor pathways. This leads to inhibition of tumor growth and prolonged tumor-free survival. Targeting WRN or its enzymatic function could prove to be an effective strategy in the treatment of MYC-associated cancers. Mol Cancer Res; 10(4); 535–45. ©2012 AACR.

Gene expression changes in adipose tissue with diet- and/or exercise-induced weight loss

Adipose tissue plays a role in obesity-related cancers via increased production of inflammatory factors, steroid hormones, and altered adipokines. The impact of weight loss on adipose tissue gene expression may provide insights into pathways linking obesity with cancer risk. We conducted an ancillary study within a randomized trial of diet, exercise, or combined diet + exercise versus control among overweight/obese postmenopausal women. In 45 women, subcutaneous adipose tissue biopsies were conducted at baseline and after 6 months, and changes in adipose tissue gene expression were determined by microarray with an emphasis on prespecified candidate pathways as well as by unsupervised clustering of more than 37,000 transcripts (Illumina). Analyses were conducted first by randomization group and then by degree of weight change at 6-months in all women combined. At 6 months, diet, exercise, and diet + exercise participants lost a mean of 8.8, 2.5, and 7.9 kg (all P < 0.05 vs. no change in controls). There was no significant change in candidate gene expression by intervention group. In analysis by weight change category, greater weight loss was associated a decrease in 17β-hydroxysteroid dehydrogenase-1 (HSD17B1, Ptrend < 0.01) and leptin (LEP, Ptrend < 0.01) expression, and marginally significant increased expression of estrogen receptor-1 (ESR1, Ptrend = 0.08) and insulin-like growth factor–binding protein-3 (IGFBP3, Ptrend = 0.08). Unsupervised clustering revealed 83 transcripts with statistically significant changes. Multiple gene expression changes correlated with changes in associated serum biomarkers. Weight loss was associated with changes in adipose tissue gene expression after 6 months, particularly in two pathways postulated to link obesity and cancer, that is, steroid hormone metabolism and IGF signaling. Cancer Prev Res; 6(3); 217–31. ©2013 AACR.

Genomic differences between estrogen receptor (ER)‐positive and ER‐negative human breast carcinoma identified by single nucleotide polymorphism array comparative genome hybridization analysis

Estrogen receptor (ER) remains one of the most important biomarkers for breast cancer subtyping and prognosis, and comparative genome hybridization has greatly contributed to the understanding of global genetic imbalance. The authors used single‐nucleotide polymorphism (SNP) arrays to compare overall copy number aberrations (CNAs) as well as loss of heterozygosity (LOH) of the entire human genome in ER‐positive and ER‐negative breast carcinomas.

Xbp1 directs global repression of budding yeast transcription during the transition to quiescence and is important for the longevity and reversibility of the quiescent state.

Pure populations of quiescent yeast can be obtained from stationary phase cultures that have ceased proliferation after exhausting glucose and other carbon sources from their environment. They are uniformly arrested in the G1 phase of the cell cycle, and display very high thermo-tolerance and longevity. We find that G1 arrest is initiated before all the glucose has been scavenged from the media. Maintaining G1 arrest requires transcriptional repression of the G1 cyclin, CLN3, by Xbp1. Xbp1 is induced as glucose is depleted and it is among the most abundant transcripts in quiescent cells. Xbp1 binds and represses CLN3 transcription and in the absence of Xbp1, or with extra copies of CLN3, cells undergo ectopic divisions and produce very small cells. The Rad53-mediated replication stress checkpoint reinforces the arrest and becomes essential when Cln3 is overproduced. The XBP1 transcript also undergoes metabolic oscillations under glucose limitation and we identified many additional transcripts that oscillate out of phase with XBP1 and have Xbp1 binding sites in their promoters. Further global analysis revealed that Xbp1 represses 15% of all yeast genes as they enter the quiescent state and over 500 of these transcripts contain Xbp1 binding sites in their promoters. Xbp1-repressed transcripts are highly enriched for genes involved in the regulation of cell growth, cell division and metabolism. Failure to repress some or all of these targets leads xbp1 cells to enter a permanent arrest or senescence with a shortened lifespan.

Methylation profiling identified novel differentially methylated markers including OPCML and FLRT2 in prostate cancer

ABSTRACT To develop new methods to distinguish indolent from aggressive prostate cancers (PCa), we utilized comprehensive high-throughput array-based relative methylation (CHARM) assay to identify differentially methylated regions (DMRs) throughout the genome, including both CpG island (CGI) and non-CGI regions in PCa patients based on Gleason grade. Initially, 26 samples, including 8 each of low [Gleason score (GS) 6] and high (GS ≥7) grade PCa samples and 10 matched normal prostate tissues, were analyzed as a discovery cohort. We identified 3,567 DMRs between normal and cancer tissues, and 913 DMRs distinguishing low from high-grade cancers. Most of these DMRs were located at CGI shores. The top 5 candidate DMRs from the low vs. high Gleason comparison, including OPCML, ELAVL2, EXT1, IRX5, and FLRT2, were validated by pyrosequencing using the discovery cohort. OPCML and FLRT2 were further validated in an independent cohort consisting of 20 low-Gleason and 33 high-Gleason tissues. We then compared patients with biochemical recurrence (n=70) vs. those without (n=86) in a third cohort, and they showed no difference in methylation at these DMR loci. When GS 3+4 cases and GS 4+3 cases were compared, OPCML-DMR methylation showed a trend of lower methylation in the recurrence group (n=30) than in the no-recurrence (n=52) group. We conclude that whole-genome methylation profiling with CHARM revealed distinct patterns of differential DNA methylation between normal prostate and PCa tissues, as well as between different risk groups of PCa as defined by Gleason scores. A panel of selected DMRs may serve as novel surrogate biomarkers for Gleason score in PCa.

Identification of differentially methylated markers among cytogenetic risk groups of acute myeloid leukemia

Aberrant DNA methylation is known to occur in cancer, including hematological malignancies such as acute myeloid leukemia (AML). However, less is known about whether specific methylation profiles characterize specific subcategories of AML. We examined this issue by using comprehensive high-throughput array-based relative methylation analysis (CHARM) to compare methylation profiles among patients in different AML cytogenetic risk groups. We found distinct profiles in each group, with the high-risk group showing overall increased methylation compared with low- and mid-risk groups. The differentially methylated regions (DMRs) distinguishing cytogenetic risk groups of AML were enriched in the CpG island shores. Specific risk-group associated DMRs were located near genes previously known to play a role in AML or other malignancies, such as MN1, UHRF1, HOXB3, and HOXB4, as well as TRIM71, the function of which in cancer is not well characterized. These findings were verified by quantitative bisulfite pyrosequencing and by comparison with results available at the TCGA cancer genome browser. To explore the potential biological significance of the observed methylation changes, we correlated our findings with gene expression data available through the TCGA database. The results showed that decreased methylation at HOXB3 and HOXB4 was associated with increased gene expression of both HOXB genes specific to the mid-risk AML, while increased DNA methylation at DCC distinctive to the high-risk AML was associated with increased gene expression. Our results suggest that the differential impact of cytogenetic changes on AML prognosis may, in part, be mediated by changes in methylation.