Nucleosome positioning and histone modifications define relationships between regulatory elements and nearby gene expression in breast epithelial cells

G-quadruplex is critical to epigenetic activation of the lncRNA HOTAIR in cancer cells

The cancer-promoting lncRNA HOTAIR has multiple isoforms. Which isoform of HOTAIR accounts for its expression and functions in cancer is unknown. Unlike HOTAIR's canonical intergenic isoform NR_003716 (HOTAIR-C), the novel isoform NR_047517 (HOTAIR-N) forms an overlapping antisense transcription locus with HOXC11. We identified HOTAIR-N as the dominant isoform that regulates the gene expression programs and networks for cell proliferation, survival, and death in cancer cells. The CpG island in the HOTAIR-N promoter was marked with epigenetic markers for active transcription. We identified a G-quadruplex (G4) motif rich region in the HOTAIR-N CpG island. Our findings indicate that G4s in HOTAIR-N CpG island is critical for expression of HOTAIR-N in cancer cells. Disruption of G4 may represent a novel therapeutic approach for cancer. The transcriptomes regulated by HOTAIR-N and Bloom in cancer cells as provided herein are important resources for the exploration of lncRNA, DNA helicases, and G4 in cancer.

Serine starvation silences estrogen receptor signaling through histone hypoacetylation

Loss of estrogen receptor (ER) pathway activity promotes breast cancer progression, yet how this occurs remains poorly understood. Here, we show that serine starvation, a metabolic stress often found in breast cancer, represses estrogen receptor alpha (ERα) signaling by reprogramming glucose metabolism and epigenetics. Using isotope tracing and time-resolved metabolomic analyses, we demonstrate that serine is required to maintain glucose flux through glycolysis and the TCA cycle to support acetyl-CoA generation for histone acetylation. Consequently, limiting serine depletes histone H3 lysine 27 acetylation (H3K27ac), particularly at the promoter region of ER pathway genes including the gene encoding ERα, ESR1. Mechanistically, serine starvation impairs acetyl-CoA-dependent gene expression by inhibiting the entry of glycolytic carbon into the TCA cycle and down-regulating the mitochondrial citrate exporter SLC25A1, a critical enzyme in the production of nucleocytosolic acetyl-CoA from glucose. Consistent with this model, total H3K27ac and ERα expression are suppressed by SLC25A1 inhibition and restored by acetate, an alternate source of acetyl-CoA, in serine-free conditions. We thus uncover an unexpected role for serine in sustaining ER signaling through the regulation of acetyl-CoA metabolism.

Characterizing chromatin interactions of regulatory elements and nucleosome positions, using Hi-C, Micro-C, and promoter capture Micro-C

BACKGROUND

Regulatory elements such as promoters, enhancers, and insulators interact each other to mediate molecular processes. To capture chromatin interactions of regulatory elements, 3C-derived methods such as Hi-C and Micro-C are developed. Here, we generated and analyzed Hi-C, Micro-C, and promoter capture Micro-C datasets with different sequencing depths to study chromatin interactions of regulatory elements and nucleosome positions in human prostate cancer cells.

RESULTS

Compared to Hi-C, Micro-C identifies more high-resolution loops, including ones around structural variants. By evaluating the effect of sequencing depth, we revealed that more than 2 billion reads of Micro-C are needed to detect chromatin interactions at 1 kb resolution. Moreover, we found that deep-sequencing identifies additional long-range loops that are longer than 1 Mb in distance. Furthermore, we found that more than 50% of the loops are involved in insulators while less than 10% of the loops are promoter-enhancer loops. To comprehensively capture chromatin interactions that promoters are involved in, we performed promoter capture Micro-C. Promoter capture Micro-C identifies loops near promoters with a lower amount of sequencing reads. Sequencing of 160 million reads of promoter capture Micro-C resulted in reaching a plateau of identifying loops. However, there was still a subset of promoters that are not involved in loops even after deep-sequencing. By integrating Micro-C with NOMe-seq and ChIP-seq, we found that active promoters involved in loops have a more accessible region with lower levels of DNA methylation and more highly phased nucleosomes, compared to active promoters that are not involved in loops.

CONCLUSION

We determined the required sequencing depth for Micro-C and promoter capture Micro-C to generate high-resolution chromatin interaction maps and loops. We also investigated the effect of sequencing coverage of Hi-C, Micro-C, and promoter capture Micro-C on detecting chromatin loops. Our analyses suggest the presence of distinct regulatory element groups, which are differently involved in nucleosome positions and chromatin interactions. This study does not only provide valuable insights on understanding chromatin interactions of regulatory elements, but also present guidelines for designing research projects on chromatin interactions among regulatory elements.

The Parkinson's disease variant rs356182 regulates neuronal differentiation independently from alpha-synuclein

One of the most significant risk variants for Parkinson's disease (PD), rs356182, is located at the PD-associated locus near the alpha-synuclein (α-syn) encoding gene, SNCA. SNCA-proximal variants, including rs356182, are thought to function in PD risk through enhancers via allele-specific regulatory effects on SNCA expression. However, this interpretation discounts the complex activity of genetic enhancers and possible non-conical functions of α-syn. Here we investigated a novel risk mechanism for rs356182. We use CRISPR-Cas9 in LUHMES cells, a model for dopaminergic midbrain neurons, to generate precise hemizygous lesions at rs356182. The PD-protective (A/-), PD-risk (G/-) and wild-type (A/G) clones were neuronally differentiated and then compared transcriptionally and morphologically. Among the affected genes was SNCA, whose expression was promoted by the PD-protective allele (A) and repressed in its absence. In addition to SNCA, hundreds of genes were differentially expressed and associated with neurogenesis and axonogenesis-an effect not typically ascribed to α-syn. We also found that the transcription factor FOXO3 specifically binds to the rs356182 A-allele in differentiated LUHMES cells. Finally, we compared the results from the rs356182-edited cells to our previously published knockouts of SNCA and found only minimal overlap between the sets of significant differentially expressed genes. Together, the data implicate a risk mechanism for rs356182 in which the risk-allele (G) is associated with abnormal neuron development, independent of SNCA expression. We speculate that these pathological effects manifest as a diminished population of dopaminergic neurons during development leading to the predisposition for PD later in life.

ZHX2 promotes HIF1α oncogenic signaling in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is an aggressive and highly lethal disease, which warrants the critical need to identify new therapeutic targets. We show that Zinc Fingers and Homeoboxes 2 (ZHX2) is amplified or overexpressed in TNBC cell lines and patients. Functionally, depletion of ZHX2 inhibited TNBC cell growth and invasion in vitro, orthotopic tumor growth, and spontaneous lung metastasis in vivo. Mechanistically, ZHX2 bound with hypoxia-inducible factor (HIF) family members and positively regulated HIF1α activity in TNBC. Integrated ChIP-seq and gene expression profiling demonstrated that ZHX2 co-occupied with HIF1α on transcriptionally active promoters marked by H3K4me3 and H3K27ac, thereby promoting gene expression. Among the identified ZHX2 and HIF1α coregulated genes, overexpression of AP2B1, COX20, KDM3A, or PTGES3L could partially rescue TNBC cell growth defect by ZHX2 depletion, suggested that these downstream targets contribute to the oncogenic role of ZHX2 in an accumulative fashion. Furthermore, multiple residues (R491, R581, and R674) on ZHX2 are important in regulating its phenotype, which correspond with their roles on controlling ZHX2 transcriptional activity in TNBC cells. These studies establish that ZHX2 activates oncogenic HIF1α signaling, therefore serving as a potential therapeutic target for TNBC.

BRD4 methylation by the methyltransferase SETD6 regulates selective transcription to control mRNA translation

The transcriptional coactivator BRD4 has a fundamental role in transcription regulation and thus became a promising epigenetic therapeutic candidate to target diverse pathologies. However, the regulation of BRD4 by posttranslational modifications has been largely unexplored. Here, we show that BRD4 is methylated on chromatin at lysine-99 by the protein lysine methyltransferase SETD6. BRD4 methylation negatively regulates the expression of genes that are involved in translation and inhibits total mRNA translation in cells. Mechanistically, we provide evidence that supports a model where BRD4 methylation by SETD6 does not have a direct role in the association with acetylated histone H4 at chromatin. However, this methylation specifically determines the recruitment of the transcription factor E2F1 to selected target genes that are involved in mRNA translation. Together, our findings reveal a previously unknown molecular mechanism for BRD4 methylation-dependent gene-specific targeting, which may serve as a new direction for the development of therapeutic applications.

TP63 links chromatin remodeling and enhancer reprogramming to epidermal differentiation and squamous cell carcinoma development

Squamous cell carcinoma (SCC) is an aggressive malignancy that can originate from various organs. TP63 is a master regulator that plays an essential role in epidermal differentiation. It is also a lineage-dependent oncogene in SCC. ΔNp63α is the prominent isoform of TP63 expressed in epidermal cells and SCC, and overexpression promotes SCC development through a variety of mechanisms. Recently, ΔNp63α was highlighted to act as an epidermal-specific pioneer factor that binds closed chromatin and enhances chromatin accessibility at epidermal enhancers. ΔNp63α coordinates chromatin-remodeling enzymes to orchestrate the tissue-specific enhancer landscape and three-dimensional high-order architecture of chromatin. Moreover, ΔNp63α establishes squamous-like enhancer landscapes to drive oncogenic target expression during SCC development. Importantly, ΔNp63α acts as an upstream regulator of super enhancers to activate a number of oncogenic transcripts linked to poor prognosis in SCC. Mechanistically, ΔNp63α activates genes transcription through physically interacting with a number of epigenetic modulators to establish enhancers and enhance chromatin accessibility. In contrast, ΔNp63α also represses gene transcription via interacting with repressive epigenetic regulators. ΔNp63α expression is regulated at multiple levels, including transcriptional, post-transcriptional, and post-translational levels. In this review, we summarize recent advances of p63 in epigenomic and transcriptional control, as well as the mechanistic regulation of p63.

TENET 2.0: Identification of key transcriptional regulators and enhancers in lung adenocarcinoma

Lung cancer is the leading cause of cancer-related death and lung adenocarcinoma is its most common subtype. Although genetic alterations have been identified as drivers in subsets of lung adenocarcinoma, they do not fully explain tumor development. Epigenetic alterations have been implicated in the pathogenesis of tumors. To identify epigenetic alterations driving lung adenocarcinoma, we used an improved version of the Tracing Enhancer Networks using Epigenetic Traits method (TENET 2.0) in primary normal lung and lung adenocarcinoma cells. We found over 32,000 enhancers that appear differentially activated between normal lung and lung adenocarcinoma. Among the identified transcriptional regulators inactivated in lung adenocarcinoma vs. normal lung, NKX2-1 was linked to a large number of silenced enhancers. Among the activated transcriptional regulators identified, CENPA, FOXM1, and MYBL2 were linked to numerous cancer-specific enhancers. High expression of CENPA, FOXM1, and MYBL2 is particularly observed in a subgroup of lung adenocarcinomas and is associated with poor patient survival. Notably, CENPA, FOXM1, and MYBL2 are also key regulators of cancer-specific enhancers in breast adenocarcinoma of the basal subtype, but they are associated with distinct sets of activated enhancers. We identified individual lung adenocarcinoma enhancers linked to CENPA, FOXM1, or MYBL2 that were associated with poor patient survival. Knockdown experiments of FOXM1 and MYBL2 suggest that these factors regulate genes involved in controlling cell cycle progression and cell division. For example, we found that expression of TK1, a potential target gene of a MYBL2-linked enhancer, is associated with poor patient survival. Identification and characterization of key transcriptional regulators and associated enhancers in lung adenocarcinoma provides important insights into the deregulation of lung adenocarcinoma epigenomes, highlighting novel potential targets for clinical intervention.

BRD4 promotes metastatic potential in oral squamous cell carcinoma through the epigenetic regulation of the MMP2 gene

BACKGROUND

Oral squamous cell carcinoma (OSCC) has increased morbidity, and its high metastatic potential affects patient survival. Bromodomain containing 4 (BRD4) is a chromatin protein that associates with acetylated histone lysines and facilitates transcription. BRD4 has been implicated in cell proliferation, metastasis, and prognosis in several types of cancer. However, the role of BRD4 in OSCC remains to be elucidated.

METHODS

We investigated the role of BRD4 and its potential utility as a therapeutic target in OSCC.

RESULTS

JQ1, the BRD4 inhibitor, suppressed the cell proliferation, migration, and invasion in the OSCC cell lines and in vivo. JQ1 reduced the expression levels of 15 metastasis genes in OSCC, including matrix metallopeptidase 2 (MMP2). Our chromatin immunoprecipitation assay showed that JQ1 reduced the BRD4 binding to the histone H3 lysine 27 acetylation-enriched sites in the MMP2 locus. Analyses of biopsy specimens from OSCC patients revealed that the BRD4 and MMP2 expression levels were correlated in the cancerous regions, and both were highly expressed in lymph node metastasis cases, including delayed metastasis.

CONCLUSIONS

BRD4 contributes to metastasis in OSCC, through the epigenetic regulation of the MMP2 gene, and thus BRD4 may represent a therapeutic target and a novel prediction indicator for metastasis.

Gene-Specific Actions of Transcriptional Coregulators Facilitate Physiological Plasticity: Evidence for a Physiological Coregulator Code

The actions of transcriptional coregulators are highly gene-specific, that is, each coregulator is required only for a subset of the genes regulated by a specific transcription factor. These coregulator-specific gene subsets often represent selected physiological responses among multiple pathways targeted by a transcription factor. Regulating the activity of a coregulator via post-translational modifications would thus affect only a subset of the transcription factor's physiological actions. Using the context of transcriptional regulation by steroid hormone receptors, this review focuses on gene-specific actions of coregulators and evidence linking individual coregulators with specific physiological pathways. Such evidence suggests that there is a 'physiological coregulator code', which represents a fertile area for future research with important clinical implications.

Chromatin interactome mapping at 139 independent breast cancer risk signals

BACKGROUND

Genome-wide association studies have identified 196 high confidence independent signals associated with breast cancer susceptibility. Variants within these signals frequently fall in distal regulatory DNA elements that control gene expression.

RESULTS

We designed a Capture Hi-C array to enrich for chromatin interactions between the credible causal variants and target genes in six human mammary epithelial and breast cancer cell lines. We show that interacting regions are enriched for open chromatin, histone marks for active enhancers, and transcription factors relevant to breast biology. We exploit this comprehensive resource to identify candidate target genes at 139 independent breast cancer risk signals and explore the functional mechanism underlying altered risk at the 12q24 risk region.

CONCLUSIONS

Our results demonstrate the power of combining genetics, computational genomics, and molecular studies to rationalize the identification of key variants and candidate target genes at breast cancer GWAS signals.

CSA: a web service for the complete process of ChIP-Seq analysis

BACKGROUND

Chromatin immunoprecipitation sequencing (ChIP-seq) is a technology that combines chromatin immunoprecipitation (ChIP) with next generation of sequencing technology (NGS) to analyze protein interactions with DNA. At present, most ChIP-seq analysis tools adopt the command line, which lacks user-friendly interfaces. Although some web services with graphical interfaces have been developed for ChIP-seq analysis, these sites cannot provide a comprehensive analysis of ChIP-seq from raw data to downstream analysis.

RESULTS

In this study, we develop a web service for the whole process of ChIP-Seq Analysis (CSA), which covers mapping, quality control, peak calling, and downstream analysis. In addition, CSA provides a customization function for users to define their own workflows. And the visualization of mapping, peak calling, motif finding, and pathway analysis results are also provided in CSA. For the different types of ChIP-seq datasets, CSA can provide the corresponding tool to perform the analysis. Moreover, CSA can detect differences in ChIP signals between ChIP samples and controls to identify absolute binding sites.

CONCLUSIONS

The two case studies demonstrate the effectiveness of CSA, which can complete the whole procedure of ChIP-seq analysis. CSA provides a web interface for users, and implements the visualization of every analysis step. The website of CSA is available at http://CompuBio.csu.edu.cn.

CD44 Promotes PD-L1 Expression and Its Tumor-Intrinsic Function in Breast and Lung Cancers

The PD-L1 (CD274) immune-checkpoint ligand is often upregulated in cancers to inhibit T cells and elicit immunosuppression. Independent of this activity, PD-L1 has recently been shown to also exert a cancer cell-intrinsic function promoting tumorigenesis. Here, we establish this tumor-intrinsic role of PD-L1 in triple-negative breast cancer (TNBC) and non-small cell lung cancer (NSCLC). Using FACS-assisted shRNA screens, we identified the cell-surface adhesion receptor CD44 as a key positive regulator of PD-L1 expression in these cancers. Mechanistically, CD44 activated PD-L1 transcription in part through its cleaved intracytoplasmic domain (ICD), which bound to a regulatory region of the PD-L1 locus containing a consensus CD44-ICD binding site. Supporting this genetic interaction, CD44 positively correlated with PD-L1 expression at the mRNA and protein levels in primary tumor samples of TNBC and NSCLC patients. These data provide a novel basis for CD44 as a critical therapeutic target to suppress PD-L1 tumor-intrinsic function. SIGNIFICANCE: CD44 is a potential target to suppress PD-L1 function in TNBC. This finding has the potential to open a new area of therapy for TNBC.

A high-resolution 3D epigenomic map reveals insights into the creation of the prostate cancer transcriptome

To better understand the impact of chromatin structure on regulation of the prostate cancer transcriptome, we develop high-resolution chromatin interaction maps in normal and prostate cancer cells using in situ Hi-C. By combining the in situ Hi-C data with active and repressive histone marks, CTCF binding sites, nucleosome-depleted regions, and transcriptome profiling, we identify topologically associating domains (TADs) that change in size and epigenetic states between normal and prostate cancer cells. Moreover, we identify normal and prostate cancer-specific enhancer-promoter loops and involved transcription factors. For example, we show that FOXA1 is enriched in prostate cancer-specific enhancer-promoter loop anchors. We also find that the chromatin structure surrounding the androgen receptor (AR) locus is altered in the prostate cancer cells with many cancer-specific enhancer-promoter loops. This creation of 3D epigenomic maps enables a better understanding of prostate cancer biology and mechanisms of gene regulation.

In prostate cancer, chromatin structure can impact the transcriptome. Here, the authors develop high resolution chromatin interaction maps in prostate cancer cells using in situ Hi-C, revealing prostate cancer-specific TADs and enhancer-promoter loops surrounding the androgen receptor (AR) locus.

MCF-7 as a Model for Functional Analysis of Breast Cancer Risk Variants

BACKGROUND

Breast cancer genetic predisposition is governed by more than 142 loci as revealed by genome-wide association studies (GWAS). The functional contribution of these risk loci to breast cancer remains unclear, and additional post-GWAS analyses are required.

METHODS

We identified active regulatory elements (enhancers, promoters, and chromatin organizing elements) by histone H3K27 acetylation and CTCF occupancy and determined the enrichment of risk variants at these sites. We compared these results with previously published data and for other cell lines, including human mammary epithelial cells, and related these data to gene expression.

RESULTS

In terms of mapping accuracy and resolution, our data augment previous annotations of the MCF-7 epigenome. After intersection with GWAS risk variants, we found 39 enhancers and 15 CTCF occupancy sites that, between them, overlapped 96 breast cancer credible risk variants at 42 loci. These risk enhancers likely regulate the expression of dozens of genes, which are enriched for GO categories, including estrogen and prolactin signaling.

CONCLUSIONS

Ten (of 142) breast cancer risk loci likely function via enhancers that are active in MCF-7 and are well suited to targeted manipulation in this system. In contrast, risk loci cannot be mapped to specific CTCF-binding sites, and the genes linked to risk CTCF sites did not show functional enrichment. The identity of risk enhancers and their associated genes suggests that some risk may function during later processes in cancer progression.

IMPACT

Here, we report how the ER⁺ cell line MCF-7 can be used to dissect risk mechanisms for breast cancer.

Using 3D epigenomic maps of primary olfactory neuronal cells from living individuals to understand gene regulation

As part of PsychENCODE, we developed a three-dimensional (3D) epigenomic map of primary cultured neuronal cells derived from olfactory neuroepithelium (CNON). We mapped topologically associating domains and high-resolution chromatin interactions using Hi-C and identified regulatory elements using chromatin immunoprecipitation and nucleosome positioning assays. Using epigenomic datasets from biopsies of 63 living individuals, we found that epigenetic marks at distal regulatory elements are more variable than marks at proximal regulatory elements. By integrating genotype and metadata, we identified enhancers that have different levels corresponding to differences in genetic variation, gender, smoking, and schizophrenia. Motif searches revealed that many CNON enhancers are bound by neuronal-related transcription factors. Last, we combined 3D epigenomic maps and gene expression profiles to predict enhancer-target gene interactions on a genome-wide scale. This study not only provides a framework for understanding individual epigenetic variation using a primary cell model system but also contributes valuable data resources for epigenomic studies of neuronal epithelium.

Beyond DNA: the Role of Epigenetics in the Premalignant Progression of Breast Cancer

Ductal Carcinoma in Situ (DCIS) is an early breast cancer lesion that is considered a nonobligate precursor to development of invasive ductal carcinoma (IDC). Although only a small subset of DCIS lesions are predicted to progress into a breast cancer, distinguishing innocuous from minacious DCIS lesions remains a clinical challenge. Thus, patients diagnosed with DCIS will undergo surgery with the potential for radiation and hormone therapy. This has led to a current state of overdiagnosis and overtreatment. Interrogating the transcriptome alone has yet to define clear functional determinants of progression from DCIS to IDC. Epigenetic changes, critical for imprinting and tissue specific development, in the incorrect context can lead to global signaling rewiring driving pathological phenotypes. Epigenetic signaling pathways, and the molecular players that interpret and sustain their signals, are critical to understanding the underlying pathology of breast cancer progression. The types of epigenetic changes, as well as the molecular players, are expanding. In addition to DNA methylation, histone modifications, and chromatin remodeling, we must also consider enhancers as well as the growing field of noncoding RNAs. Herein we will review the epigenetic interactions that have been uncovered in early stage lesions that impact breast cancer progression, and how these players may be utilized as biomarkers to mitigate overdiagnosis and overtreatment.

Genetics and Expression Profile of the Tubulin Gene Superfamily in Breast Cancer Subtypes and Its Relation to Taxane Resistance

Taxanes are a class of chemotherapeutic agents that inhibit cell division by disrupting the mitotic spindle through the stabilization of microtubules. Most breast cancer (BC) tumors show resistance against taxanes partially due to alterations in tubulin genes. In this project we investigated tubulin isoforms in BC to explore any correlation between tubulin alterations and taxane resistance. Genetic alteration and expression profiling of 28 tubulin isoforms in 6714 BC tumor samples from 4205 BC cases were analyzed. Protein-protein, drug-protein and alterations neighbor genes in tubulin pathways were examined in the tumor samples. To study correlation between promoter activity and expression of the tubulin isoforms in BC, we analyzed the ChIP-seq enrichment of active promoter histone mark H3K4me3 and mRNA expression profile of MCF-7, ZR-75-30, SKBR-3 and MDA-MB-231 cell lines. Potential correlation between tubulin alterations and taxane resistance, were investigated by studying the expression profile of taxane-sensitive and resistant BC tumors also the MDA-MB-231 cells acquired resistance to paclitaxel. All genomic data were obtained from public databases. Results showed that TUBD1 and TUBB3 were the most frequently amplified and deleted tubulin genes in the BC tumors respectively. The interaction analysis showed physical interactions of α-, β- and γ-tubulin isoforms with each other. The most of FDA-approved tubulin inhibitor drugs including taxanes target only β-tubulins. The analysis also revealed sex tubulin-interacting neighbor proteins including ENCCT3, NEK2, PFDN2, PTP4A3, SDCCAG8 and TBCE which were altered in at least 20% of the tumors. Three of them are tubulin-specific chaperons responsible for tubulin protein folding. Expression of tubulin genes in BC cell lines were correlated with H3K4me3 enrichment on their promoter chromatin. Analyzing expression profile of BC tumors and tumor-adjacent normal breast tissues showed upregulation of TUBA1A, TUBA1C, TUBB and TUBB3 and downregulation of TUBB2A, TUBB2B, TUBB6, TUBB7P pseudogene, and TUBGCP2 in the tumor tissues compared to the normal breast tissues. Analyzing taxane-sensitive versus taxane-resistant tumors revealed that expression of TUBB3 and TUBB6 was significantly downregulated in the taxane-resistant tumors. Our results suggest that downregulation of tumor βIII- and βV-tubulins is correlated with taxane resistance in BC. Based on our results, we conclude that aberrant protein folding of tubulins due to mutation and/or dysfunction of tubulin-specific chaperons may be potential mechanisms of taxane resistance. Thus, we propose studying the molecular pathology of tubulin mutations and folding in BC and their impacts on taxane resistance.

Genetics and Expression Profile of the Tubulin Gene Superfamily in Breast Cancer Subtypes and Its Relation to Taxane Resistance

Taxanes are a class of chemotherapeutic agents that inhibit cell division by disrupting the mitotic spindle through the stabilization of microtubules. Most breast cancer (BC) tumors show resistance against taxanes partially due to alterations in tubulin genes. In this project we investigated tubulin isoforms in BC to explore any correlation between tubulin alterations and taxane resistance. Genetic alteration and expression profiling of 28 tubulin isoforms in 6714 BC tumor samples from 4205 BC cases were analyzed. Protein-protein, drug-protein and alterations neighbor genes in tubulin pathways were examined in the tumor samples. To study correlation between promoter activity and expression of the tubulin isoforms in BC, we analyzed the ChIP-seq enrichment of active promoter histone mark H3K4me3 and mRNA expression profile of MCF-7, ZR-75-30, SKBR-3 and MDA-MB-231 cell lines. Potential correlation between tubulin alterations and taxane resistance, were investigated by studying the expression profile of taxane-sensitive and resistant BC tumors also the MDA-MB-231 cells acquired resistance to paclitaxel. All genomic data were obtained from public databases. Results showed that TUBD1 and TUBB3 were the most frequently amplified and deleted tubulin genes in the BC tumors respectively. The interaction analysis showed physical interactions of α-, β- and γ-tubulin isoforms with each other. The most of FDA-approved tubulin inhibitor drugs including taxanes target only β-tubulins. The analysis also revealed sex tubulin-interacting neighbor proteins including ENCCT3, NEK2, PFDN2, PTP4A3, SDCCAG8 and TBCE which were altered in at least 20% of the tumors. Three of them are tubulin-specific chaperons responsible for tubulin protein folding. Expression of tubulin genes in BC cell lines were correlated with H3K4me3 enrichment on their promoter chromatin. Analyzing expression profile of BC tumors and tumor-adjacent normal breast tissues showed upregulation of TUBA1A, TUBA1C, TUBB and TUBB3 and downregulation of TUBB2A, TUBB2B, TUBB6, TUBB7P pseudogene, and TUBGCP2 in the tumor tissues compared to the normal breast tissues. Analyzing taxane-sensitive versus taxane-resistant tumors revealed that expression of TUBB3 and TUBB6 was significantly downregulated in the taxane-resistant tumors. Our results suggest that downregulation of tumor βIII- and βV-tubulins is correlated with taxane resistance in BC. Based on our results, we conclude that aberrant protein folding of tubulins due to mutation and/or dysfunction of tubulin-specific chaperons may be potential mechanisms of taxane resistance. Thus, we propose studying the molecular pathology of tubulin mutations and folding in BC and their impacts on taxane resistance.

ZFX acts as a transcriptional activator in multiple types of human tumors by binding downstream from transcription start sites at the majority of CpG island promoters

High expression of the transcription factor ZFX is correlated with proliferation, tumorigenesis, and patient survival in multiple types of human cancers. However, the mechanism by which ZFX influences transcriptional regulation has not been determined. We performed ChIP-seq in four cancer cell lines (representing kidney, colon, prostate, and breast cancers) to identify ZFX binding sites throughout the human genome. We identified roughly 9000 ZFX binding sites and found that most of the sites are in CpG island promoters. Moreover, genes with promoters bound by ZFX are expressed at higher levels than genes with promoters not bound by ZFX. To determine if ZFX contributes to regulation of the promoters to which it is bound, we performed RNA-seq analysis after knockdown of ZFX by siRNA in prostate and breast cancer cells. Many genes with promoters bound by ZFX were down-regulated upon ZFX knockdown, supporting the hypothesis that ZFX acts as a transcriptional activator. Surprisingly, ZFX binds at +240 bp downstream from the TSS of the responsive promoters. Using Nucleosome Occupancy and Methylome Sequencing (NOMe-seq), we show that ZFX binds between the open chromatin region at the TSS and the first downstream nucleosome, suggesting that ZFX may play a critical role in promoter architecture. We have also shown that a closely related zinc finger protein ZNF711 has a similar binding pattern at CpG island promoters, but ZNF711 may play a subordinate role to ZFX. This functional characterization of ZFX provides important new insights into transcription, chromatin structure, and the regulation of the cancer transcriptome.

Identification of cis-regulatory mutations generating de novo edges in personalized cancer gene regulatory networks

The identification of functional non-coding mutations is a key challenge in the field of genomics. Here we introduce μ-cisTarget to filter, annotate and prioritize cis-regulatory mutations based on their putative effect on the underlying "personal" gene regulatory network. We validated μ-cisTarget by re-analyzing the TAL1 and LMO1 enhancer mutations in T-ALL, and the TERT promoter mutation in melanoma. Next, we re-sequenced the full genomes of ten cancer cell lines and used matched transcriptome data and motif discovery to identify master regulators with de novo binding sites that result in the up-regulation of nearby oncogenic drivers. μ-cisTarget is available from http://mucistarget.aertslab.org .

Systematic functional perturbations uncover a prognostic genetic network driving human breast cancer

Prognostic classifiers conceivably comprise biomarker genes that functionally contribute to the oncogenic and metastatic properties of cancer, but this has not been investigated systematically. The transcription factor Fra-1 not only has an essential role in breast cancer, but also drives the expression of a highly prognostic gene set. Here, we systematically perturbed the function of 31 individual Fra-1-dependent poor-prognosis genes and examined their impact on breast cancer growth in vivo. We find that stable shRNA depletion of each of nine individual signature genes strongly inhibits breast cancer growth and aggressiveness. Several factors within this nine-gene set regulate each others expression, suggesting that together they form a network. The nine-gene set is regulated by estrogen, ERBB2 and EGF signaling, all established breast cancer factors. We also uncover three transcription factors, MYC, E2F1 and TP53, which act alongside Fra-1 at the core of this network. ChIP-Seq analysis reveals that a substantial number of genes are bound, and regulated, by all four transcription factors. The nine-gene set retains significant prognostic power and includes several potential therapeutic targets, including the bifunctional enzyme PAICS, which catalyzes purine biosynthesis. Depletion of PAICS largely cancelled breast cancer expansion, exemplifying a prognostic gene with breast cancer activity. Our data uncover a core genetic and prognostic network driving human breast cancer. We propose that pharmacological inhibition of components within this network, such as PAICS, may be used in conjunction with the Fra-1 prognostic classifier towards personalized management of poor prognosis breast cancer.

Characterizing Genetic Susceptibility to Breast Cancer in Women of African Ancestry

Background: Genome-wide association studies have identified approximately 100 common genetic variants associated with breast cancer risk, the majority of which were discovered in women of European ancestry. Because of different patterns of linkage disequilibrium, many of these genetic markers may not represent signals in populations of African ancestry.Methods: We tested 74 breast cancer risk variants and conducted fine-mapping of these susceptibility regions in 6,522 breast cancer cases and 7,643 controls of African ancestry from three genetic consortia (AABC, AMBER, and ROOT).Results: Fifty-four of the 74 variants (73%) were found to have ORs that were directionally consistent with those previously reported, of which 12 were nominally statistically significant (P < 0.05). Through fine-mapping, in six regions (3p24, 12p11, 14q13, 16q12/FTO, 16q23, 19p13), we observed seven markers that better represent the underlying risk variant for overall breast cancer or breast cancer subtypes, whereas in another two regions (11q13, 16q12/TOX3), we identified suggestive evidence of signals that are independent of the reported index variant. Overlapping chromatin features and regulatory elements suggest that many of the risk alleles lie in regions with biological functionality.Conclusions: Through fine-mapping of known susceptibility regions, we have revealed alleles that better characterize breast cancer risk in women of African ancestry.Impact: The risk alleles identified represent genetic markers for modeling and stratifying breast cancer risk in women of African ancestry. Cancer Epidemiol Biomarkers Prev; 26(7); 1016-26. ©2017 AACR.

Priming of transcriptional memory responses via the chromatin accessibility landscape in T cells

Memory T cells exhibit transcriptional memory and “remember” their previous pathogenic encounter to increase transcription on re-infection. However, how this transcriptional priming response is regulated is unknown. Here we performed global FAIRE-seq profiling of chromatin accessibility in a human T cell transcriptional memory model. Primary activation induced persistent accessibility changes, and secondary activation induced secondary-specific opening of previously less accessible regions associated with enhanced expression of memory-responsive genes. Increased accessibility occurred largely in distal regulatory regions and was associated with increased histone acetylation and relative H3.3 deposition. The enhanced re-stimulation response was linked to the strength of initial PKC-induced signalling, and PKC-sensitive increases in accessibility upon initial stimulation showed higher accessibility on re-stimulation. While accessibility maintenance was associated with ETS-1, accessibility at re-stimulation-specific regions was linked to NFAT, especially in combination with ETS-1, EGR, GATA, NFκB, and NR4A. Furthermore, NFATC1 was directly regulated by ETS-1 at an enhancer region. In contrast to the factors that increased accessibility, signalling from bHLH and ZEB family members enhanced decreased accessibility upon re-stimulation. Interplay between distal regulatory elements, accessibility, and the combined action of sequence-specific transcription factors allows transcriptional memory-responsive genes to “remember” their initial environmental encounter.

Small genomic insertions form enhancers that misregulate oncogenes

The non-coding regions of tumour cell genomes harbour a considerable fraction of total DNA sequence variation, but the functional contribution of these variants to tumorigenesis is ill-defined. Among these non-coding variants, somatic insertions are among the least well characterized due to challenges with interpreting short-read DNA sequences. Here, using a combination of Chip-seq to enrich enhancer DNA and a computational approach with multiple DNA alignment procedures, we identify enhancer-associated small insertion variants. Among the 102 tumour cell genomes we analyse, small insertions are frequently observed in enhancer DNA sequences near known oncogenes. Further study of one insertion, somatically acquired in primary leukaemia tumour genomes, reveals that it nucleates formation of an active enhancer that drives expression of the LMO2 oncogene. The approach described here to identify enhancer-associated small insertion variants provides a foundation for further study of these abnormalities across human cancers.

Sequencing initiatives have detected multiple types of mutations in cancer. Here the authors, analysing enhancer-targeting sequence data, show that small insertions in transcriptional enhancers are frequently found near oncogenes, and demonstrate how one mutation deregulates expression of LMO2 in leukemia cells.

Identification of activated enhancers and linked transcription factors in breast, prostate, and kidney tumors by tracing enhancer networks using epigenetic traits

Background

Although technological advances now allow increased tumor profiling, a detailed understanding of the mechanisms leading to the development of different cancers remains elusive. Our approach toward understanding the molecular events that lead to cancer is to characterize changes in transcriptional regulatory networks between normal and tumor tissue. Because enhancer activity is thought to be critical in regulating cell fate decisions, we have focused our studies on distal regulatory elements and transcription factors that bind to these elements.

Results

Using DNA methylation data, we identified more than 25,000 enhancers that are differentially activated in breast, prostate, and kidney tumor tissues, as compared to normal tissues. We then developed an analytical approach called Tracing Enhancer Networks using Epigenetic Traits that correlates DNA methylation levels at enhancers with gene expression to identify more than 800,000 genome-wide links from enhancers to genes and from genes to enhancers. We found more than 1200 transcription factors to be involved in these tumor-specific enhancer networks. We further characterized several transcription factors linked to a large number of enhancers in each tumor type, including GATA3 in non-basal breast tumors, HOXC6 and DLX1 in prostate tumors, and ZNF395 in kidney tumors. We showed that HOXC6 and DLX1 are associated with different clusters of prostate tumor-specific enhancers and confer distinct transcriptomic changes upon knockdown in C42B prostate cancer cells. We also discovered de novo motifs enriched in enhancers linked to ZNF395 in kidney tumors.

Conclusions

Our studies characterized tumor-specific enhancers and revealed key transcription factors involved in enhancer networks for specific tumor types and subgroups. Our findings, which include a large set of identified enhancers and transcription factors linked to those enhancers in breast, prostate, and kidney cancers, will facilitate understanding of enhancer networks and mechanisms leading to the development of these cancers.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0102-4) contains supplementary material, which is available to authorized users.

Functional characterization of open chromatin in bidirectional promoters of rice

Bidirectional gene pairs tend to be highly coregulated and function in similar biological processes in eukaryotic genomes. Structural features and functional consequences of bidirectional promoters (BDPs) have received considerable attention among diverse species. However, the underlying mechanisms responsible for the bidirectional transcription and coexpression of BDPs remain poorly understood in plants. In this study, we integrated DNase-seq, RNA-seq, ChIP-seq and MNase-seq data and investigated the effect of physical DNase I hypersensitive site (DHS) positions on the transcription of rice BDPs. We found that the physical position of a DHS relative to the TSS of bidirectional gene pairs can affect the expression of the corresponding genes: the closer a DHS is to the TSS, the higher is the expression level of the genes. Most importantly, we observed that the distribution of DHSs plays a significant role in the regulation of transcription and the coexpression of gene pairs, which are possibly mediated by orchestrating the positioning of histone marks and canonical nucleosomes around BDPs. Our results demonstrate that the combined actions of chromatin structures with DHSs, which contain functional cis-elements for interaction with transcriptional machinery, may play an important role in the regulation of the bidirectional transcription or coexpression in rice BDPs. Our findings may help to enhance the understanding of DHSs in the regulation of bidirectional gene pairs.

A Comparative Analysis of Genetic and Epigenetic Events of Breast and Ovarian Cancer Related to Tumorigenesis

Breast cancer persists as the most common cause of cancer death in women worldwide. Ovarian cancer is also a significant source of morbidity and mortality, as the fifth leading cause of cancer death among women. This reflects the continued need for further understanding and innovation in cancer treatment. Though breast and ovarian cancer usually present as distinct clinical entities, the recent explosion of large-scale -omics research has uncovered many overlaps, particularly with respect to genetic and epigenetic alterations. We compared genetic, microenvironmental, stromal, and epigenetic changes common between breast and ovarian cancer cells, as well as the clinical relevance of these changes. Some of the most striking commonalities include genetic alterations of BRCA1 and 2, TP53, RB1, NF1, FAT3, MYC, PTEN, and PIK3CA; down regulation of miRNAs 9, 100, 125a, 125b, and 214; and epigenetic alterations such as H3K27me3, H3K9me2, H3K9me3, H4K20me3, and H3K4me. These parallels suggest shared features of pathogenesis. Furthermore, preliminary evidence suggests a shared epigenetic mechanism of oncogenesis. These similarities, warrant further investigation in order to ultimately inform development of more effective chemotherapeutics, as well as strategies to circumvent drug resistance.

Identification of four novel susceptibility loci for oestrogen receptor negative breast cancer

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

Identification of chromatin accessibility domains in human breast cancer stem cells

Epithelial-to-mesenchymal transition (EMT) is physiological in embryogenesis and wound healing but also associated with the formation of cancer stem cells (CSCs). Many EMT signaling pathways are implicated in CSC formation, but the precise underlying mechanisms of CSC formation remain elusive. We have previously demonstrated that PKC is critical for EMT induction and CSC formation in inducible breast EMT/CSC models. Here, we used formaldehyde-assisted isolation of regulatory elements-sequencing (FAIRE-seq) to investigate DNA accessibility changes after PKC activation and determine how they influence EMT and CSC formation. During EMT, DNA accessibility principally increased in regions distant from transcription start sites, low in CpG content, and enriched with chromatin enhancer marks. ChIP-sequencing revealed that a subset of these regions changed from poised to active enhancers upon stimulation, with some even more acteylated in CSCs. While regions with increased accessibility were enriched for FOX, AP-1, TEAD, and TFAP2 motifs, those containing FOX and AP-1 motif were associated with increased expression of CSC-associated genes, while those with TFAP2 were associated with genes with increased expression in non-CSCs. Silencing of 2 members of the FOX family, FOXN2 and FOXQ1, repressed CSCs and the mesenchymal phenotype and inhibited the CSC gene signature. These novel, PKC-induced DNA accessibility regions help explain how the epigenomic plasticity of cells undergoing EMT leads to CSC gene activation.

Integration of multiethnic fine-mapping and genomic annotation to prioritize candidate functional SNPs at prostate cancer susceptibility regions

Interpretation of biological mechanisms underlying genetic risk associations for prostate cancer is complicated by the relatively large number of risk variants (n = 100) and the thousands of surrogate SNPs in linkage disequilibrium. Here, we combined three distinct approaches: multiethnic fine-mapping, putative functional annotation (based upon epigenetic data and genome-encoded features), and expression quantitative trait loci (eQTL) analyses, in an attempt to reduce this complexity. We examined 67 risk regions using genotyping and imputation-based fine-mapping in populations of European (cases/controls: 8600/6946), African (cases/controls: 5327/5136), Japanese (cases/controls: 2563/4391) and Latino (cases/controls: 1034/1046) ancestry. Markers at 55 regions passed a region-specific significance threshold (P-value cutoff range: 3.9 × 10(-4)-5.6 × 10(-3)) and in 30 regions we identified markers that were more significantly associated with risk than the previously reported variants in the multiethnic sample. Novel secondary signals (P < 5.0 × 10(-6)) were also detected in two regions (rs13062436/3q21 and rs17181170/3p12). Among 666 variants in the 55 regions with P-values within one order of magnitude of the most-associated marker, 193 variants (29%) in 48 regions overlapped with epigenetic or other putative functional marks. In 11 of the 55 regions, cis-eQTLs were detected with nearby genes. For 12 of the 55 regions (22%), the most significant region-specific, prostate-cancer associated variant represented the strongest candidate functional variant based on our annotations; the number of regions increased to 20 (36%) and 27 (49%) when examining the 2 and 3 most significantly associated variants in each region, respectively. These results have prioritized subsets of candidate variants for downstream functional evaluation.

Genome-wide association between YAP/TAZ/TEAD and AP-1 at enhancers drives oncogenic growth

YAP/TAZ are nuclear effectors of the Hippo pathway regulating organ growth and tumorigenesis. Yet, their function as transcriptional regulators remains underinvestigated. By ChIP-seq analyses in breast cancer cells, we discovered that the YAP/TAZ transcriptional response is pervasively mediated by a dual element: TEAD factors, through which YAP/TAZ bind to DNA, co-occupying chromatin with activator protein-1 (AP-1, dimer of JUN and FOS proteins) at composite cis-regulatory elements harbouring both TEAD and AP-1 motifs. YAP/TAZ/TEAD and AP-1 form a complex that synergistically activates target genes directly involved in the control of S-phase entry and mitosis. This control occurs almost exclusively from distal enhancers that contact target promoters through chromatin looping. YAP/TAZ-induced oncogenic growth is strongly enhanced by gain of AP-1 and severely blunted by its loss. Conversely, AP-1-promoted skin tumorigenesis is prevented in YAP/TAZ conditional knockout mice. This work highlights a new layer of signalling integration, feeding on YAP/TAZ function at the chromatin level.

Comprehensive assembly of novel transcripts from unmapped human RNA-Seq data and their association with cancer

Crucial parts of the genome including genes encoding microRNAs and noncoding RNAs went unnoticed for years, and even now, despite extensive annotation and assembly of the human genome, RNA-sequencing continues to yield millions of unmappable and thus uncharacterized reads. Here, we examined > 300 billion reads from 536 normal donors and 1,873 patients encompassing 21 cancer types, identified ∼300 million such uncharacterized reads, and using a distinctive approach de novo assembled 2,550 novel human transcripts, which mainly represent long noncoding RNAs. Of these, 230 exhibited relatively specific expression or non-expression in certain cancer types, making them potential markers for those cancers, whereas 183 exhibited tissue specificity. Moreover, we used lentiviral-mediated expression of three selected transcripts that had higher expression in normal than in cancer patients and found that each inhibited the growth of HepG2 cells. Our analysis provides a comprehensive and unbiased resource of unmapped human transcripts and reveals their associations with specific cancers, providing potentially important new genes for therapeutic targeting.

Cell-type-specific enrichment of risk-associated regulatory elements at ovarian cancer susceptibility loci

Understanding the regulatory landscape of the human genome is a central question in complex trait genetics. Most single-nucleotide polymorphisms (SNPs) associated with cancer risk lie in non-protein-coding regions, implicating regulatory DNA elements as functional targets of susceptibility variants. Here, we describe genome-wide annotation of regions of open chromatin and histone modification in fallopian tube and ovarian surface epithelial cells (FTSECs, OSECs), the debated cellular origins of high-grade serous ovarian cancers (HGSOCs) and in endometriosis epithelial cells (EECs), the likely precursor of clear cell ovarian carcinomas (CCOCs). The regulatory architecture of these cell types was compared with normal human mammary epithelial cells and LNCaP prostate cancer cells. We observed similar positional patterns of global enhancer signatures across the three different ovarian cancer precursor cell types, and evidence of tissue-specific regulatory signatures compared to non-gynecological cell types. We found significant enrichment for risk-associated SNPs intersecting regulatory biofeatures at 17 known HGSOC susceptibility loci in FTSECs (P = 3.8 × 10(-30)), OSECs (P = 2.4 × 10(-23)) and HMECs (P = 6.7 × 10(-15)) but not for EECs (P = 0.45) or LNCaP cells (P = 0.88). Hierarchical clustering of risk SNPs conditioned on the six different cell types indicates FTSECs and OSECs are highly related (96% of samples using multi-scale bootstrapping) suggesting both cell types may be precursors of HGSOC. These data represent the first description of regulatory catalogues of normal precursor cells for different ovarian cancer subtypes, and provide unique insights into the tissue specific regulatory variation with respect to the likely functional targets of germline genetic susceptibility variants for ovarian cancer.

Lipid-induced epigenomic changes in human macrophages identify a coronary artery disease-associated variant that regulates PPAP2B Expression through Altered C/EBP-beta binding

Genome-wide association studies (GWAS) have identified over 40 loci that affect risk of coronary artery disease (CAD) and the causal mechanisms at the majority of loci are unknown. Recent studies have suggested that many causal GWAS variants influence disease through altered transcriptional regulation in disease-relevant cell types. We explored changes in transcriptional regulation during a key pathophysiological event in CAD, the environmental lipid-induced transformation of macrophages to lipid-laden foam cells. We used a combination of open chromatin mapping with formaldehyde-assisted isolation of regulatory elements (FAIRE-seq) and enhancer and transcription factor mapping using chromatin immuno-precipitation (ChIP-seq) in primary human macrophages before and after exposure to atherogenic oxidized low-density lipoprotein (oxLDL), with resultant foam cell formation. OxLDL-induced foam cell formation was associated with changes in a subset of open chromatin and active enhancer sites that strongly correlated with expression changes of nearby genes. OxLDL-regulated enhancers were enriched for several transcription factors including C/EBP-beta, which has no previously documented role in foam cell formation. OxLDL exposure up-regulated C/EBP-beta expression and increased genomic binding events, most prominently around genes involved in inflammatory response pathways. Variants at CAD-associated loci were significantly and specifically enriched in the subset of chromatin sites altered by oxLDL exposure, including rs72664324 in an oxLDL-induced enhancer at the PPAP2B locus. OxLDL increased C/EBP beta binding to this site and C/EBP beta binding and enhancer activity were stronger with the protective A allele of rs72664324. In addition, expression of the PPAP2B protein product LPP3 was present in foam cells in human atherosclerotic plaques and oxLDL exposure up-regulated LPP3 in macrophages resulting in increased degradation of pro-inflammatory mediators. Our results demonstrate a genetic mechanism contributing to CAD risk at the PPAP2B locus and highlight the value of studying epigenetic changes in disease processes involving pathogenic environmental stimuli.

Coronary artery disease is a complex disease where over 40 genomic loci contributing to genetic risk have been identified. However, identifying the precise variants, genomic elements and genes that mediate this risk at each locus has proved challenging. We hypothesized that some genetic risk variants may influence a key step in development of coronary artery disease, which occurs when macrophages encounter environmentally-derived lipid. These cells take up lipid and accumulate in atherosclerotic plaques in the walls of blood vessels where they contribute to the inflammatory atherosclerotic disease process. Therefore, we studied the effects of this lipid exposure on the genomic activity of these cells. Environmental lipid exposure triggered changes in transcriptional regulation and gene expression. Variants at coronary artery disease risk loci were enriched for genomic regions altered by lipid exposure. We studied one such risk variant rs72664324 in detail and found that it altered binding of the C/EBP-beta transcription factor and altered expression of the PPAP2B gene. PPAP2B encodes an enzyme that degrades pro-inflammatory substances. Our study demonstrates a hitherto unknown genetic mechanism underlying atherosclerotic heart disease and demonstrates the value of studying changes in transcriptional regulation in key disease processes involving environmental influences.

Disturbing the histone code in leukemia: translocations and mutations affecting histone methyl transferases

Leukemia is characterized by increased numbers of blasts originating from transformed early hematopoietic stem and progenitor cells. Genetic alterations are widely recognized as the main drivers of oncogenic transformation. Of considerable interest are mutations affecting the writers of epigenetic marks. In this review, we focus on histone methyltransferases--enzymes that catalyze the methylation of lysine residues in core histones. Histone methylation is a tightly controlled mechanism that is responsible for both activating as well as repressing gene expression in a site-specific manner, depending on which lysine residue is methylated. Histone methyltransferases, including MLL1, DOT1L, EZH2, and SETD2 are recurrently deregulated in human leukemia, either directly by gene mutations or balanced translocations, or indirectly as components of protein complexes that are disturbed in leukemia due to alterations of the other components in these complexes. Several small molecule inhibitors of histone methyltransferases are currently being clinically evaluated for their therapeutic potential in human leukemia. These drugs reverse some of the adverse effects of aberrant histone methylation, and can induce differentiation and cell death in leukemic blasts.