Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.

Response and resistance to BET bromodomain inhibitors in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is a heterogeneous and clinically aggressive disease for which there is no targeted therapy. BET bromodomain inhibitors, which have shown efficacy in several models of cancer, have not been evaluated in TNBC. These inhibitors displace BET bromodomain proteins such as BRD4 from chromatin by competing with their acetyl-lysine recognition modules, leading to inhibition of oncogenic transcriptional programs. Here we report the preferential sensitivity of TNBCs to BET bromodomain inhibition in vitro and in vivo, establishing a rationale for clinical investigation and further motivation to understand mechanisms of resistance. In paired cell lines selected for acquired resistance to BET inhibition from previously sensitive TNBCs, we failed to identify gatekeeper mutations, new driver events or drug pump activation. BET-resistant TNBC cells remain dependent on wild-type BRD4, which supports transcription and cell proliferation in a bromodomain-independent manner. Proteomic studies of resistant TNBC identify strong association with MED1 and hyper-phosphorylation of BRD4 attributable to decreased activity of PP2A, identified here as a principal BRD4 serine phosphatase. Together, these studies provide a rationale for BET inhibition in TNBC and present mechanism-based combination strategies to anticipate clinical drug resistance.

Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data.

MOTIVATION: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts. RESULTS: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays. AVAILABILITY: The R package abs filter for library clonality simulations and detection of amplification-biased sites is available from http://updepla1srv1.epfl.ch/waszaks/absfilter

Heterogeneity and transcriptional drivers of triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is a heterogeneous disease with limited treatment options. To characterize TNBC heterogeneity, we defined transcriptional, epigenetic, and metabolic subtypes and subtype-driving super-enhancers and transcription factors by combining functional and molecular profiling with computational analyses. Single-cell RNA sequencing revealed relative homogeneity of the major transcriptional subtypes (luminal, basal, and mesenchymal) within samples. We found that mesenchymal TNBCs share features with mesenchymal neuroblastoma and rhabdoid tumors and that the PRRX1 transcription factor is a key driver of these tumors. PRRX1 is sufficient for inducing mesenchymal features in basal but not in luminal TNBC cells via reprogramming super-enhancer landscapes, but it is not required for mesenchymal state maintenance or for cellular viability. Our comprehensive, large-scale, multiplatform, multiomics study of both experimental and clinical TNBC is an important resource for the scientific and clinical research communities and opens venues for future investigation.

YAP-Mediated Recruitment of YY1 and EZH2 Represses Transcription of Key Cell-Cycle Regulators.

The Hippo pathway regulates cell proliferation and organ size through control of the transcriptional regulators YAP (yes-associated protein) and TAZ. Upon extracellular stimuli such as cell-cell contact, the pathway negatively regulates YAP through cytoplasmic sequestration. Under conditions of low cell density, YAP is nuclear and associates with enhancer regions and gene promoters. YAP is mainly described as a transcriptional activator of genes involved in cell proliferation and survival. Using a genome-wide approach, we show here that, in addition to its known function as a transcriptional activator, YAP functions as a transcriptional repressor by interacting with the multifunctional transcription factor Yin Yang 1 (YY1) and Polycomb repressive complex member enhancer of zeste homologue 2 (EZH2). YAP colocalized with YY1 and EZH2 on the genome to transcriptionally repress a broad network of genes mediating a host of cellular functions, including repression of the cell-cycle kinase inhibitor p27, whose role is to functionally promote contact inhibition. This work unveils a broad and underappreciated aspect of YAP nuclear function as a transcriptional repressor and highlights how loss of contact inhibition in cancer is mediated in part through YAP repressive function. SIGNIFICANCE: This study provides new insights into YAP as a broad transcriptional repressor of key regulators of the cell cycle, in turn influencing contact inhibition and tumorigenesis.

Opposite effects of Mn2+ and Zn2+ on PsaR-mediated expression of the virulence genes pcpA, prtA, and psaBCA of Streptococcus pneumoniae.

Homeostasis of Zn(2+) and Mn(2+) is important for the physiology and virulence of the human pathogen Streptococcus pneumoniae. Here, transcriptome analysis was used to determine the response of S. pneumoniae D39 to a high concentration of Zn(2+). Interestingly, virulence genes encoding the choline binding protein PcpA, the extracellular serine protease PrtA, and the Mn(2+) uptake system PsaBC(A) were strongly upregulated in the presence of Zn(2+). Using random mutagenesis, a previously described Mn(2+)-responsive transcriptional repressor, PsaR, was found to mediate the observed Zn(2+)-dependent derepression. In addition, PsaR is also responsible for the Mn(2+)-dependent repression of these genes. Subsequently, we investigated how these opposite effects are mediated by the same regulator. In vitro binding of purified PsaR to the prtA, pcpA, and psaB promoters was stimulated by Mn(2+), whereas Zn(2+) destroyed the interaction of PsaR with its target promoters. Mutational analysis of the pcpA promoter demonstrated the presence of a PsaR operator that mediates the transcriptional effects. In conclusion, PsaR is responsible for the counteracting effects of Mn(2+) and Zn(2+) on the expression of several virulence genes in S. pneumoniae, suggesting that the ratio of these metal ions exerts an important influence on pneumococcal pathogenesis.

TRPS1 Is a Lineage-Specific Transcriptional Dependency in Breast Cancer.

Perturbed epigenomic programs play key roles in tumorigenesis, and chromatin modulators are candidate therapeutic targets in various human cancer types. To define singular and shared dependencies on DNA and histone modifiers and transcription factors in poorly differentiated adult and pediatric cancers, we conducted a targeted shRNA screen across 59 cell lines of 6 cancer types. Here, we describe the TRPS1 transcription factor as a strong breast cancer-specific hit, owing largely to lineage-restricted expression. Knockdown of TRPS1 resulted in perturbed mitosis, apoptosis, and reduced tumor growth. Integrated analysis of TRPS1 transcriptional targets, chromatin binding, and protein interactions revealed that TRPS1 is associated with the NuRD repressor complex. These findings uncover a transcriptional network that is essential for breast cancer cell survival and propagation.

Genes with monoallelic expression contribute disproportionately to genetic diversity in humans.

An unexpectedly large number of human autosomal genes are subject to monoallelic expression (MAE). Our analysis of 4,227 such genes uncovers surprisingly high genetic variation across human populations. This increased diversity is unlikely to reflect relaxed purifying selection. Remarkably, MAE genes exhibit an elevated recombination rate and an increased density of hypermutable sequence contexts. However, these factors do not fully account for the increased diversity. We find that the elevated nucleotide diversity of MAE genes is also associated with greater allelic age: variants in these genes tend to be older and are enriched in polymorphisms shared by Neanderthals and chimpanzees. Both synonymous and nonsynonymous alleles of MAE genes have elevated average population frequencies. We also observed strong enrichment of the MAE signature among genes reported to evolve under balancing selection. We propose that an important biological function of widespread MAE might be the generation of cell-to-cell heterogeneity; the increased genetic variation contributes to this heterogeneity.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications.

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this "neighboring effect", we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples. We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype. GEO SERIES ACCESSION NUMBER: GSE33784, GSE33867.

Coordinated effects of sequence variation on DNA binding, chromatin structure, and transcription.

DNA sequence variation has been associated with quantitative changes in molecular phenotypes such as gene expression, but its impact on chromatin states is poorly characterized. To understand the interplay between chromatin and genetic control of gene regulation, we quantified allelic variability in transcription factor binding, histone modifications, and gene expression within humans. We found abundant allelic specificity in chromatin and extensive local, short-range, and long-range allelic coordination among the studied molecular phenotypes. We observed genetic influence on most of these phenotypes, with histone modifications exhibiting strong context-dependent behavior. Our results implicate transcription factors as primary mediators of sequence-specific regulation of gene expression programs, with histone modifications frequently reflecting the primary regulatory event.