TXNIP loss expands Myc-dependent transcriptional programs by increasing Myc genomic binding.

The c-Myc protooncogene places a demand on glucose uptake to drive glucose-dependent biosynthetic pathways. To meet this demand, c-Myc protein (Myc henceforth) drives the expression of glucose transporters, glycolytic enzymes, and represses the expression of thioredoxin interacting protein (TXNIP), which is a potent negative regulator of glucose uptake. A Mychigh/TXNIPlow gene signature is clinically significant as it correlates with poor clinical prognosis in triple-negative breast cancer (TNBC) but not in other subtypes of breast cancer, suggesting a functional relationship between Myc and TXNIP. To better understand how TXNIP contributes to the aggressive behavior of TNBC, we generated TXNIP null MDA-MB-231 (231:TKO) cells for our study. We show that TXNIP loss drives a transcriptional program that resembles those driven by Myc and increases global Myc genome occupancy. TXNIP loss allows Myc to invade the promoters and enhancers of target genes that are potentially relevant to cell transformation. Together, these findings suggest that TXNIP is a broad repressor of Myc genomic binding. The increase in Myc genomic binding in the 231:TKO cells expands the Myc-dependent transcriptome we identified in parental MDA-MB-231 cells. This expansion of Myc-dependent transcription following TXNIP loss occurs without an apparent increase in Myc's intrinsic capacity to activate transcription and without increasing Myc levels. Together, our findings suggest that TXNIP loss mimics Myc overexpression, connecting Myc genomic binding and transcriptional programs to the nutrient and progrowth signals that control TXNIP expression.

Comprehensive evaluation and efficient classification of BRCA1 RING domain missense substitutions.

BRCA1 is a high-risk susceptibility gene for breast and ovarian cancer. Pathogenic protein-truncating variants are scattered across the open reading frame, but all known missense substitutions that are pathogenic because of missense dysfunction are located in either the amino-terminal RING domain or the carboxy-terminal BRCT domain. Heterodimerization of the BRCA1 and BARD1 RING domains is a molecularly defined obligate activity. Hence, we tested every BRCA1 RING domain missense substitution that can be created by a single nucleotide change for heterodimerization with BARD1 in a mammalian two-hybrid assay. Downstream of the laboratory assay, we addressed three additional challenges: assay calibration, validation thereof, and integration of the calibrated results with other available data, such as computational evidence and patient/population observational data to achieve clinically applicable classification. Overall, we found that 15%-20% of BRCA1 RING domain missense substitutions are pathogenic. Using a Bayesian point system for data integration and variant classification, we achieved clinical classification of 89% of observed missense substitutions. Moreover, among missense substitutions not present in the human observational data used here, we find an additional 45 with concordant computational and functional assay evidence in favor of pathogenicity plus 223 with concordant evidence in favor of benignity; these are particularly likely to be classified as likely pathogenic and likely benign, respectively, once human observational data become available.

MS4A3 promotes differentiation in chronic myeloid leukemia by enhancing common ß-chain cytokine receptor endocytosis.

The chronic phase of chronic myeloid leukemia (CP-CML) is characterized by the excessive production of maturating myeloid cells. As CML stem/progenitor cells (LSPCs) are poised to cycle and differentiate, LSPCs must balance conservation and differentiation to avoid exhaustion, similar to normal hematopoiesis under stress. Since BCR-ABL1 tyrosine kinase inhibitors (TKIs) eliminate differentiating cells but spare BCR-ABL1-independent LSPCs, understanding the mechanisms that regulate LSPC differentiation may inform strategies to eliminate LSPCs. Upon performing a meta-analysis of published CML transcriptomes, we discovered that low expression of the MS4A3 transmembrane protein is a universal characteristic of LSPC quiescence, BCR-ABL1 independence, and transformation to blast phase (BP). Several mechanisms are involved in suppressing MS4A3, including aberrant methylation and a MECOM-C/EBPε axis. Contrary to previous reports, we find that MS4A3 does not function as a G1/S phase inhibitor but promotes endocytosis of common ß-chain (ßc) cytokine receptors upon GM-CSF/IL-3 stimulation, enhancing downstream signaling and cellular differentiation. This suggests that LSPCs downregulate MS4A3 to evade ßc cytokine-induced differentiation and maintain a more primitive, TKI-insensitive state. Accordingly, knockdown (KD) or deletion of MS4A3/Ms4a3 promotes TKI resistance and survival of CML cells ex vivo and enhances leukemogenesis in vivo, while targeted delivery of exogenous MS4A3 protein promotes differentiation. These data support a model in which MS4A3 governs response to differentiating myeloid cytokines, providing a unifying mechanism for the differentiation block characteristic of CML quiescence and BP-CML. Promoting MS4A3 reexpression or delivery of ectopic MS4A3 may help eliminate LSPCs in vivo.

The lingering mysteries of metastatic recurrence in breast cancer.

Despite being the hallmark of cancer that is responsible for the highest number of deaths, very little is known about the biology of metastasis. Metastatic disease typically manifests after a protracted period of undetectable disease following surgery or systemic therapy, owing to relapse or recurrence. In the case of breast cancer, metastatic relapse can occur months to decades after initial diagnosis and treatment. In this review, we provide an overview of the known key factors that influence metastatic recurrence, with the goal of highlighting the critical unanswered questions that still need to be addressed to make a difference in the mortality of breast cancer patients.

Distinct properties of cell-type-specific and shared transcription factor binding sites.

Most human transcription factors bind a small subset of potential genomic sites and often use different subsets in different cell types. To identify mechanisms that govern cell-type-specific transcription factor binding, we used an integrative approach to study estrogen receptor α (ER). We found that ER exhibits two distinct modes of binding. Shared sites, bound in multiple cell types, are characterized by high-affinity estrogen response elements (EREs), inaccessible chromatin, and a lack of DNA methylation, while cell-specific sites are characterized by a lack of EREs, co-occurrence with other transcription factors, and cell-type-specific chromatin accessibility and DNA methylation. These observations enabled accurate quantitative models of ER binding that suggest tethering of ER to one-third of cell-specific sites. The distinct properties of cell-specific binding were also observed with glucocorticoid receptor and for ER in primary mouse tissues, representing an elegant genomic encoding scheme for generating cell-type-specific gene regulation.

TBCRC 002: a phase II, randomized, open-label trial of preoperative letrozole with or without bevacizumab in postmenopausal women with newly diagnosed stage 2/3 hormone receptor-positive and HER2-negative breast cancer.

BACKGROUND: In preclinical studies, the expression of vascular endothelial growth factor (VEGF) in hormone receptor-positive breast cancer is associated with estrogen-independent tumor growth and resistance to endocrine therapies. This study investigated whether the addition of bevacizumab, a monoclonal antibody against VEGF, to letrozole enhanced the antitumor activity of the letrozole in the preoperative setting. METHODS: Postmenopausal women with newly diagnosed stage 2 or 3 estrogen and/or progesterone receptor-positive, HER2-negative breast cancer were randomly assigned (2:1) between letrozole 2.5 mg PO daily plus bevacizumab 15 mg/kg IV every 3 weeks (Let/Bev) and letrozole 2.5 mg PO daily (Let) for 24 weeks prior to definitive surgery. Primary objective was within-arm pathologic complete remission (pCR) rate. Secondary objectives were safety, objective response, and downstaging rate. RESULTS: Seventy-five patients were randomized (Let/Bev n = 50, Let n = 25). Of the 45 patients evaluable for pathological response in the Let/Bev arm, 5 (11%; 95% CI, 3.7-24.1%) achieved pCR and 4 (9%; 95% CI, 2.5-21.2%) had microscopic residual disease; no pCRs or microscopic residual disease was seen in the Let arm (0%; 95% CI, 0-14.2%). The rates of downstaging were 44.4% (95% CI, 29.6-60.0%) and 37.5% (95% CI, 18.8-59.4%) in the Let/Bev and Let arms, respectively. Adverse events typically associated with letrozole (hot flashes, arthralgias, fatigue, myalgias) occurred in similar frequencies in the two arms. Hypertension, headache, and proteinuria were seen exclusively in the Let/Bev arm. The rates of grade 3 and 4 adverse events and discontinuation due to adverse events were 18% vs 8% and 16% vs none in the Let/Bev and Let arms, respectively. A small RNA-based classifier predictive of response to preoperative Let/Bev was developed and confirmed on an independent cohort. CONCLUSION: In the preoperative setting, the addition of bevacizumab to letrozole was associated with a pCR rate of 11%; no pCR was seen with letrozole alone. There was additive toxicity with the incorporation of bevacizumab. Responses to Let/Bev can be predicted from the levels of 5 small RNAs in a pretreatment biopsy. TRIAL REGISTRATION: This trial is registered with ClinicalTrials.gov (Identifier: NCT00161291), first posted on September 12, 2005, and is completed.

Dynamic DNA methylation across diverse human cell lines and tissues.

As studies of DNA methylation increase in scope, it has become evident that methylation has a complex relationship with gene expression, plays an important role in defining cell types, and is disrupted in many diseases. We describe large-scale single-base resolution DNA methylation profiling on a diverse collection of 82 human cell lines and tissues using reduced representation bisulfite sequencing (RRBS). Analysis integrating RNA-seq and ChIP-seq data illuminates the functional role of this dynamic mark. Loci that are hypermethylated across cancer types are enriched for sites bound by NANOG in embryonic stem cells, which supports and expands the model of a stem/progenitor cell signature in cancer. CpGs that are hypomethylated across cancer types are concentrated in megabase-scale domains that occur near the telomeres and centromeres of chromosomes, are depleted of genes, and are enriched for cancer-specific EZH2 binding and H3K27me3 (repressive chromatin). In noncancer samples, there are cell-type specific methylation signatures preserved in primary cell lines and tissues as well as methylation differences induced by cell culture. The relationship between methylation and expression is context-dependent, and we find that CpG-rich enhancers bound by EP300 in the bodies of expressed genes are unmethylated despite the dense gene-body methylation surrounding them. Non-CpG cytosine methylation occurs in human somatic tissue, is particularly prevalent in brain tissue, and is reproducible across many individuals. This study provides an atlas of DNA methylation across diverse and well-characterized samples and enables new discoveries about DNA methylation and its role in gene regulation and disease.

Genistein and bisphenol A exposure cause estrogen receptor 1 to bind thousands of sites in a cell type-specific manner.

Endogenous estrogens that are synthesized in the body impact gene regulation by activating estrogen receptors in diverse cell types. Exogenous compounds that have estrogenic properties can also be found circulating in the blood in both children and adults. The genome-wide impact of these environmental estrogens on gene regulation is unclear. To obtain an integrated view of gene regulation in response to environmental and endogenous estrogens on a genome-wide scale, we performed ChIP-seq to identify estrogen receptor 1 (ESR1; previously estrogen receptor α) binding sites, and RNA-seq in endometrial cancer cells exposed to bisphenol A (BPA; found in plastics), genistein (GEN; found in soybean), or 17ß-estradiol (E2; an endogenous estrogen). GEN and BPA treatment induces thousands of ESR1 binding sites and >50 gene expression changes, representing a subset of E2-induced gene regulation changes. Genes affected by E2 were highly enriched for ribosome-associated proteins; however, GEN and BPA failed to regulate most ribosome-associated proteins and instead enriched for transporters of carboxylic acids. Treatment-dependent changes in gene expression were associated with treatment-dependent ESR1 binding sites, with the exception that many genes up-regulated by E2 harbored a BPA-induced ESR1 binding site but failed to show any expression change after BPA treatment. GEN and BPA exhibited a similar relationship to E2 in the breast cancer line T-47D, where cell type specificity played a much larger role than treatment specificity. Overall, both environmental estrogens clearly regulate gene expression through ESR1 on a genome-wide scale, although with lower potency resulting in less ESR1 binding sites and less gene expression changes compared to the endogenous estrogen, E2.

Transposase mediated construction of RNA-seq libraries.

RNA-seq has been widely adopted as a gene-expression measurement tool due to the detail, resolution, and sensitivity of transcript characterization that the technique provides. Here we present two transposon-based methods that efficiently construct high-quality RNA-seq libraries. We first describe a method that creates RNA-seq libraries for Illumina sequencing from double-stranded cDNA with only two enzymatic reactions. We generated high-quality RNA-seq libraries from as little as 10 pg of mRNA (∼1 ng of total RNA) with this approach. We also present a strand-specific RNA-seq library construction protocol that combines transposon-based library construction with uracil DNA glycosylase and endonuclease VIII to specifically degrade the second strand constructed during cDNA synthesis. The directional RNA-seq libraries maintain the same quality as the nondirectional libraries, while showing a high degree of strand specificity, such that 99.5% of reads map to the expected genomic strand. Each transposon-based library construction method performed well when compared with standard RNA-seq library construction methods with regard to complexity of the libraries, correlation between biological replicates, and the percentage of reads that align to the genome as well as exons. Our results show that high-quality RNA-seq libraries can be constructed efficiently and in an automatable fashion using transposition technology.

Effects of sequence variation on differential allelic transcription factor occupancy and gene expression.

A complex interplay between transcription factors (TFs) and the genome regulates transcription. However, connecting variation in genome sequence with variation in TF binding and gene expression is challenging due to environmental differences between individuals and cell types. To address this problem, we measured genome-wide differential allelic occupancy of 24 TFs and EP300 in a human lymphoblastoid cell line GM12878. Overall, 5% of human TF binding sites have an allelic imbalance in occupancy. At many sites, TFs clustered in TF-binding hubs on the same homolog in especially open chromatin. While genetic variation in core TF binding motifs generally resulted in large allelic differences in TF occupancy, most allelic differences in occupancy were subtle and associated with disruption of weak or noncanonical motifs. We also measured genome-wide differential allelic expression of genes with and without heterozygous exonic variants in the same cells. We found that genes with differential allelic expression were overall less expressed both in GM12878 cells and in unrelated human cell lines. Comparing TF occupancy with expression, we found strong association between allelic occupancy and expression within 100 bp of transcription start sites (TSSs), and weak association up to 100 kb from TSSs. Sites of differential allelic occupancy were significantly enriched for variants associated with disease, particularly autoimmune disease, suggesting that allelic differences in TF occupancy give functional insights into intergenic variants associated with disease. Our results have the potential to increase the power and interpretability of association studies by targeting functional intergenic variants in addition to protein coding sequences.

Widespread plasticity in CTCF occupancy linked to DNA methylation.

CTCF is a ubiquitously expressed regulator of fundamental genomic processes including transcription, intra- and interchromosomal interactions, and chromatin structure. Because of its critical role in genome function, CTCF binding patterns have long been assumed to be largely invariant across different cellular environments. Here we analyze genome-wide occupancy patterns of CTCF by ChIP-seq in 19 diverse human cell types, including normal primary cells and immortal lines. We observed highly reproducible yet surprisingly plastic genomic binding landscapes, indicative of strong cell-selective regulation of CTCF occupancy. Comparison with massively parallel bisulfite sequencing data indicates that 41% of variable CTCF binding is linked to differential DNA methylation, concentrated at two critical positions within the CTCF recognition sequence. Unexpectedly, CTCF binding patterns were markedly different in normal versus immortal cells, with the latter showing widespread disruption of CTCF binding associated with increased methylation. Strikingly, this disruption is accompanied by up-regulation of CTCF expression, with the result that both normal and immortal cells maintain the same average number of CTCF occupancy sites genome-wide. These results reveal a tight linkage between DNA methylation and the global occupancy patterns of a major sequence-specific regulatory factor.

CTCF/cohesin-mediated DNA looping is required for protocadherin α promoter choice.

The closely linked human protocadherin (Pcdh) α, ß, and γ gene clusters encode 53 distinct protein isoforms, which are expressed in a combinatorial manner to generate enormous diversity on the surface of individual neurons. This diversity is a consequence of stochastic promoter choice and alternative pre-mRNA processing. Here, we show that Pcdhα promoter choice is achieved by DNA looping between two downstream transcriptional enhancers and individual promoters driving the expression of alternate Pcdhα isoforms. In addition, we show that this DNA looping requires specific binding of the CTCF/cohesin complex to two symmetrically aligned binding sites in both the transcriptionally active promoters and in one of the enhancers. These findings have important implications regarding enhancer/promoter interactions in the generation of complex Pcdh cell surface codes for the establishment of neuronal identity and self-avoidance in individual neurons.

Recurrent read-through fusion transcripts in breast cancer.

Read-through fusion transcripts that result from the splicing of two adjacent genes in the same coding orientation are a recently discovered type of chimeric RNA. We sought to determine if read-through fusion transcripts exist in breast cancer. We performed paired-end RNA-seq of 168 breast samples, including 28 breast cancer cell lines, 42 triple negative breast cancer primary tumors, 42 estrogen receptor positive (ER+) breast cancer primary tumors, and 56 non-malignant breast tissue samples. We analyzed the sequencing data to identify breast cancer associated read-through fusion transcripts. We discovered two recurrent read-through fusion transcripts that were identified in breast cancer cell lines, confirmed across breast cancer primary tumors, and were not detected in normal tissues (SCNN1A-TNFRSF1A and CTSD-IFITM10). Both fusion transcripts use canonical splice sites to join the last splice donor of the 5' gene to the first splice acceptor of the 3' gene, creating an in-frame fusion transcript. Western blots indicated that the fusion transcripts are translated into fusion proteins in breast cancer cells. Custom small interfering RNAs targeting the CTSD-IFITM10 fusion junction reduced expression of the fusion transcript and reduced breast cancer cell proliferation. Read-through fusion transcripts between adjacent genes with different biochemical functions represent a new type of recurrent molecular defect in breast cancer that warrant further investigation as potential biomarkers and therapeutic targets. Both breast cancer associated fusion transcripts identified in this study involve membrane proteins (SCNN1A-TNFRSF1A and CTSD-IFITM10), which raises the possibility that they could be breast cancer-specific cell surface markers.

Analysis of DNA methylation in a three-generation family reveals widespread genetic influence on epigenetic regulation.

The methylation of cytosines in CpG dinucleotides is essential for cellular differentiation and the progression of many cancers, and it plays an important role in gametic imprinting. To assess variation and inheritance of genome-wide patterns of DNA methylation simultaneously in humans, we applied reduced representation bisulfite sequencing (RRBS) to somatic DNA from six members of a three-generation family. We observed that 8.1% of heterozygous SNPs are associated with differential methylation in cis, which provides a robust signature for Mendelian transmission and relatedness. The vast majority of differential methylation between homologous chromosomes (>92%) occurs on a particular haplotype as opposed to being associated with the gender of the parent of origin, indicating that genotype affects DNA methylation of far more loci than does gametic imprinting. We found that 75% of genotype-dependent differential methylation events in the family are also seen in unrelated individuals and that overall genotype can explain 80% of the variation in DNA methylation. These events are under-represented in CpG islands, enriched in intergenic regions, and located in regions of low evolutionary conservation. Even though they are generally not in functionally constrained regions, 22% (twice as many as expected by chance) of genes harboring genotype-dependent DNA methylation exhibited allele-specific gene expression as measured by RNA-seq of a lymphoblastoid cell line, indicating that some of these events are associated with gene expression differences. Overall, our results demonstrate that the influence of genotype on patterns of DNA methylation is widespread in the genome and greatly exceeds the influence of imprinting on genome-wide methylation patterns.

3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression.

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1 , ETV1 , CRKL , and ID4 . Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC , on the same complex amplicons such as ecDNA. We characterized a MYC - ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC 's enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

IS-seq: a novel high throughput survey of in vivo IS6110 transposition in multiple Mycobacterium tuberculosis genomes.

BACKGROUND: The insertion element IS6110 is one of the main sources of genomic variability in Mycobacterium tuberculosis, the etiological agent of human tuberculosis. Although IS 6110 has been used extensively as an epidemiological marker, the identification of the precise chromosomal insertion sites has been limited by technical challenges. Here, we present IS-seq, a novel method that combines high-throughput sequencing using Illumina technology with efficient combinatorial sample multiplexing to simultaneously probe 519 clinical isolates, identifying almost all the flanking regions of the element in a single experiment. RESULTS: We identified a total of 6,976 IS6110 flanking regions on the different isolates. When validated using reference strains, the method had 100% specificity and 98% positive predictive value. The insertions mapped to both coding and non-coding regions, and in some cases interrupted genes thought to be essential for virulence or in vitro growth. Strains were classified into families using insertion sites, and high agreement with previous studies was observed. CONCLUSIONS: This high-throughput IS-seq method, which can also be used to map insertions in other organisms, extends previous surveys of in vivo interrupted loci and provides a baseline for probing the consequences of disruptions in M. tuberculosis strains.

Quantification of rare allelic variants from pooled genomic DNA.

We report a targeted, cost-effective method to quantify rare single-nucleotide polymorphisms from pooled human genomic DNA using second-generation sequencing. We pooled DNA from 1,111 individuals and targeted four genes to identify rare germline variants. Our base-calling algorithm, SNPSeeker, derived from large deviation theory, detected single-nucleotide polymorphisms present at frequencies below the raw error rate of the sequencing platform.

Tet2 coordinates with Foxo1 and Runx1 to balance T follicular helper cell and T helper 1 cell differentiation.

In response to various types of infection, naïve CD4+ T cells differentiate into diverse helper T cell subsets; however, the epigenetic programs that regulate differentiation in response to viral infection remain poorly understood. Demethylation of CpG dinucleotides by Tet methylcytosine dioxygenases is a key component of epigenetic programing that promotes specific gene expression, cellular differentiation, and function. We report that following viral infection, Tet2-deficient CD4+ T cells preferentially differentiate into highly functional germinal center T follicular helper (TFH) cells that provide enhanced help for B cells. Using genome-wide DNA methylation and transcription factor binding analyses, we find that Tet2 coordinates with multiple transcription factors, including Foxo1 and Runx1, to mediate the demethylation and expression of target genes, including genes encoding repressors of TFH differentiation. Our findings establish Tet2 as an important regulator of TFH cell differentiation and reveal pathways that could be targeted to enhance immune responses against infectious disease.

Loss of G0/G1 switch gene 2 (G0S2) promotes disease progression and drug resistance in chronic myeloid leukaemia (CML) by disrupting glycerophospholipid metabolism.

Tyrosine kinase inhibitors (TKIs) targeting BCR::ABL1 have turned chronic myeloid leukaemia (CML) from a fatal disease into a manageable condition for most patients. Despite improved survival, targeting drug-resistant leukaemia stem cells (LSCs) remains a challenge for curative CML therapy. Aberrant lipid metabolism can have a large impact on membrane dynamics, cell survival and therapeutic responses in cancer. While ceramide and sphingolipid levels were previously correlated with TKI response in CML, the role of lipid metabolism in TKI resistance is not well understood. We have identified downregulation of a critical regulator of lipid metabolism, G0/G1 switch gene 2 (G0S2), in multiple scenarios of TKI resistance, including (1) BCR::ABL1 kinase-independent TKI resistance, (2) progression of CML from the chronic to the blast phase of the disease, and (3) in CML versus normal myeloid progenitors. Accordingly, CML patients with low G0S2 expression levels had a worse overall survival. G0S2 downregulation in CML was not a result of promoter hypermethylation or BCR::ABL1 kinase activity, but was rather due to transcriptional repression by MYC. Using CML cell lines, patient samples and G0s2 knockout (G0s2-/- ) mice, we demonstrate a tumour suppressor role for G0S2 in CML and TKI resistance. Our data suggest that reduced G0S2 protein expression in CML disrupts glycerophospholipid metabolism, correlating with a block of differentiation that renders CML cells resistant to therapy. Altogether, our data unravel a new role for G0S2 in regulating myeloid differentiation and TKI response in CML, and suggest that restoring G0S2 may have clinical utility.