Genome-wide RNAi screens in human brain tumor isolates reveal a novel viability requirement for PHF5A.

To identify key regulators of human brain tumor maintenance and initiation, we performed multiple genome-wide RNAi screens in patient-derived glioblastoma multiforme (GBM) stem cells (GSCs). These screens identified the plant homeodomain (PHD)-finger domain protein PHF5A as differentially required for GSC expansion, as compared with untransformed neural stem cells (NSCs) and fibroblasts. Given PHF5A's known involvement in facilitating interactions between the U2 snRNP complex and ATP-dependent helicases, we examined cancer-specific roles in RNA splicing. We found that in GSCs, but not untransformed controls, PHF5A facilitates recognition of exons with unusual C-rich 3' splice sites in thousands of essential genes. PHF5A knockdown in GSCs, but not untransformed NSCs, astrocytes, or fibroblasts, inhibited splicing of these genes, leading to cell cycle arrest and loss of viability. Notably, pharmacologic inhibition of U2 snRNP activity phenocopied PHF5A knockdown in GSCs and also in NSCs or fibroblasts overexpressing MYC. Furthermore, PHF5A inhibition compromised GSC tumor formation in vivo and inhibited growth of established GBM patient-derived xenograft tumors. Our results demonstrate a novel viability requirement for PHF5A to maintain proper exon recognition in brain tumor-initiating cells and may provide new inroads for novel anti-GBM therapeutic strategies.

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.

Transposon mutagenesis identifies candidate genes that cooperate with loss of transforming growth factor-beta signaling in mouse intestinal neoplasms.

Colorectal cancer (CRC) results from the accumulation of gene mutations and epigenetic alterations in colon epithelial cells, which promotes CRC formation through deregulating signaling pathways. One of the most commonly deregulated signaling pathways in CRC is the transforming growth factor ß (TGF-ß) pathway. Importantly, the effects of TGF-ß signaling inactivation in CRC are modified by concurrent mutations in the tumor cell, and these concurrent mutations determine the ultimate biological effects of impaired TGF-ß signaling in the tumor. However, many of the mutations that cooperate with the deregulated TGF-ß signaling pathway in CRC remain unknown. Therefore, we sought to identify candidate driver genes that promote the formation of CRC in the setting of TGF-ß signaling inactivation. We performed a forward genetic screen in mice carrying conditionally inactivated alleles of the TGF-ß receptor, type II (Tgfbr2) using Sleeping Beauty (SB) transposon mediated mutagenesis. We used TAPDANCE and Gene-centric statistical methods to identify common insertion sites (CIS) and, thus, candidate tumor suppressor genes and oncogenes within the tumor genome. CIS analysis of multiple neoplasms from these mice identified many candidate Tgfbr2 cooperating genes and the Wnt/ß-catenin, Hippo and MAPK pathways as the most commonly affected pathways. Importantly, the majority of candidate genes were also found to be mutated in human CRC. The SB transposon system provides an unbiased method to identify Tgfbr2 cooperating genes in mouse CRC that are functionally relevant and that may provide further insight into the pathogenesis of human CRC.

Prognostic methylation markers for overall survival in cytogenetically normal patients with acute myeloid leukemia treated on SWOG trials.

BACKGROUND: Aberrant DNA methylation is known to occur in patients with acute myeloid leukemia (AML), whereas methylation signatures and prognostic markers have been proposed. The objective of the current study was to evaluate all CpG sites of the genome and identify prognostic methylation markers for overall survival in patients with AML with normal karyotype (AML-NK). METHODS: AML-NK samples from 7 SWOG trials were analyzed using a novel genome-wide approach called "CHARMcox" (comprehensive high-throughput array-based relative methylation analysis combined with the Cox proportional hazards model) controlling for known clinical covariates. CHARMcox was applied to a phase 1 discovery cohort (72 patients) to identify survival-associated methylation regions (SAMRs). Subsequently, using bisulfite pyrosequencing, SAMRs were studied in phase 2 model-building (65 patients) and phase 3 validation (65 patients) cohorts. An independent external cohort from The Cancer Genome Atlas (TCGA) AML study (LAML) was used for further validation (93 patients). RESULTS: Two SAMRs, located at the CpG island shores of leucine zipper tumor suppressor 2 (LZTS2) and nuclear receptor subfamily 6 group a member 1 (NR6A1), respectively, were identified. Multivariable analyses demonstrated that hypomethylation of either LZTS2 or NR6A1 was associated with worse overall survival in the SWOG cohort (P<.001). The prognosis was validated in patients with AML-NK from the TCGA-LAML cohort. Methylation values below the median at both markers predicted worse overall survival (SWOG: hazard ratio, 1.89 [P<.001]; and TCGA-LAML: hazard ratio, 2.08 [P=.006]). The C-statistic was 0.71 for both cohorts, and the impact was independent of the Fms-related tyrosine kinase 3 internal tandem duplication (FLT3-ITD) status. CONCLUSIONS: The 2 methylation markers, measurable by clinically applicable assays such as bisulfite pyrosequencing, are promising for risk stratification among patients with AML-NK. Cancer 2017;123:2472-81. © 2017 American Cancer Society.

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.

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.

Large-scale polymorphism of heterochromatic repeats in the DNA of Arabidopsis thaliana.

BACKGROUND: The composition of the individual eukaryote's genome and its variation within a species remain poorly defined. Even for a sequenced genome such as that of the model plant Arabidopsis thaliana accession Col-0, the large arrays of heterochromatic repeats are incompletely sequenced, with gaps of uncertain size persisting in them. RESULTS: Using geographically separate populations of A. thaliana, we assayed variation in the heterochromatic repeat arrays using two independent methods and identified significant polymorphism among them, with variation by as much as a factor of two in the centromeric 180 bp repeat, in the 45S rDNA arrays and in the Athila retroelements. In the accession with highest genome size as measured by flow cytometry, Loh-0, we found more than a two-fold increase in 5S RNA gene copies relative to Col-0; results from fluorescence in situ hybridization with 5S probes were consistent with the existence of size polymorphism between Loh-0 and Col-0 at the 5S loci. Comparative genomic hybridization results of Loh-0 and Col-0 did not support contiguous variation in copy number of protein-coding genes on the scale needed to explain their observed genome size difference. We developed a computational data model to test whether the variation we measured in the repeat fractions could account for the different genome sizes determined with flow cytometry, and found that this proposed relationship could account for about 50% of the variance in genome size among the accessions. CONCLUSION: Our analyses are consistent with substantial repeat number polymorphism for 5S and 45S ribosomal genes among accession of A. thaliana. Differences are also suggested for centromeric and pericentromeric repeats. Our analysis also points to the difficulties in measuring the repeated fraction of the genome and suggests that independent validation of genome size should be sought in addition to flow cytometric measurements.

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

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.

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.

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.

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, p16(Ink4a), 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.

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

BACKGROUND: 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. METHODS: DNA was extracted from frozen tumor sections of 21 breast carcinoma specimens and analyzed with a proprietary 50K XbaI SNP array. Copy number and LOH probability values were derived for each sample. Data were analyzed using bioinformatics and computational software, and permutation tests were used to estimate the significance of these values. RESULTS: There was a global increase in CNAs and LOH in ER-negative relative to ER-positive cancers. Gain of the long arm of chromosome 1 (1q) and 8q were the most obvious changes common in both subtypes: An increase in the chromosome 1 short arm (1p)/1q ratio was observed in ER-negative samples, and an increased 16p/16q ratio was observed in ER-positive samples. Significant CNAs (adjusted P<.05) in ER-negative relative to ER-positive tumors included 5q deletion, loss of 15q, and gain of 2p and 21q. Copy-neutral LOH (cnLOH) common to both ER-positive and ER-negative samples included 9p21, the p16 tumor suppressor locus, and 4q13, the RCHY1 (ring finger and CHY zinc finger domain-containing 1) oncogene locus. Of particular interest was an enrichment of 17q LOH among the ER-negative tumors, potentially suggesting breast cancer 1 gene (BRCA1) mutations. CONCLUSIONS: SNP array detected both genetic imbalances and cnLOH and was capable of discriminating ER-negative breast cancer from ER-positive breast cancer.

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.

Remodeling of DNA methylation and phenotypic and transcriptional changes in synthetic Arabidopsis allotetraploids.

The joining of different genomes in allotetraploids played a major role in plant evolution, but the molecular implications of this event are poorly understood. In synthetic allotetraploids of Arabidopsis and Cardaminopsis arenosa, we previously demonstrated the occurrence of frequent gene silencing. To explore the involvement of epigenetic phenomena, we investigated the occurrence and effects of DNA methylation changes. Changes in DNA methylation patterns were more frequent in synthetic allotetraploids than in the parents. Treatment with 5-aza-2'-deoxycytidine, an inhibitor of DNA methyltransferase, resulted in the development of altered morphologies in the synthetic allotetraploids, but not in the parents. We profiled mRNAs in control and 5-aza-2'-deoxycytidine-treated parents and allotetraploids by amplified fragment length polymorphism-cDNA. We show that DNA demethylation induced and repressed two different transcriptomes. Our results are consistent with the hypothesis that synthetic allotetraploids have compromised mechanisms of epigenetic gene regulation.