Defining super-enhancers by highly ranked histone H4 multi-acetylation levels identifies transcription factors associated with glioblastoma stem-like properties.

BACKGROUND: Super-enhancers (SEs), which activate genes involved in cell-type specificity, have mainly been defined as genomic regions with top-ranked enrichment(s) of histone H3 with acetylated K27 (H3K27ac) and/or transcription coactivator(s) including a bromodomain and extra-terminal domain (BET) family protein, BRD4. However, BRD4 preferentially binds to multi-acetylated histone H4, typically with acetylated K5 and K8 (H4K5acK8ac), leading us to hypothesize that SEs should be defined by high H4K5acK8ac enrichment at least as well as by that of H3K27ac. RESULTS: Here, we conducted genome-wide profiling of H4K5acK8ac and H3K27ac, BRD4 binding, and the transcriptome by using a BET inhibitor, JQ1, in three human glial cell lines. When SEs were defined as having the top ranks for H4K5acK8ac or H3K27ac signal, 43% of H4K5acK8ac-ranked SEs were distinct from H3K27ac-ranked SEs in a glioblastoma stem-like cell (GSC) line. CRISPR-Cas9-mediated deletion of the H4K5acK8ac-preferred SEs associated with MYCN and NFIC decreased the stem-like properties in GSCs. CONCLUSIONS: Collectively, our data highlights H4K5acK8ac's utility for identifying genes regulating cell-type specificity.

A flexible cross-platform single-cell data processing pipeline.

Single-cell RNA-sequencing analysis to quantify the RNA molecules in individual cells has become popular, as it can obtain a large amount of information from each experiment. We introduce UniverSC ( https://github.com/minoda-lab/universc ), a universal single-cell RNA-seq data processing tool that supports any unique molecular identifier-based platform. Our command-line tool, docker image, and containerised graphical application enables consistent and comprehensive integration, comparison, and evaluation across data generated from a wide range of platforms. We also provide a cross-platform application to run UniverSC via a graphical user interface, available for macOS, Windows, and Linux Ubuntu, negating one of the bottlenecks with single-cell RNA-seq analysis that is data processing for researchers who are not bioinformatically proficient.

Integrative analysis of scRNA-seq and scATAC-seq revealed transit-amplifying thymic epithelial cells expressing autoimmune regulator.

Medullary thymic epithelial cells (mTECs) are critical for self-tolerance induction in T cells via promiscuous expression of tissue-specific antigens (TSAs), which are controlled by the transcriptional regulator, AIRE. Whereas AIRE-expressing (Aire+) mTECs undergo constant turnover in the adult thymus, mechanisms underlying differentiation of postnatal mTECs remain to be discovered. Integrative analysis of single-cell assays for transposase-accessible chromatin (scATAC-seq) and single-cell RNA sequencing (scRNA-seq) suggested the presence of proliferating mTECs with a specific chromatin structure, which express high levels of Aire and co-stimulatory molecules, CD80 (Aire+CD80hi). Proliferating Aire+CD80hi mTECs detected using Fucci technology express a minimal number of Aire-dependent TSAs and are converted into quiescent Aire+CD80hi mTECs expressing high levels of TSAs after a transit amplification. These data provide evidence for the existence of transit-amplifying Aire+mTEC precursors during the Aire+mTEC differentiation process of the postnatal thymus.

E-cadherin-deficient cells have synthetic lethal vulnerabilities in plasma membrane organisation, dynamics and function.

BACKGROUND: The E-cadherin gene (CDH1) is frequently mutated in diffuse gastric cancer and lobular breast cancer, and germline mutations predispose to the cancer syndrome Hereditary Diffuse Gastric Cancer. We are taking a synthetic lethal approach to identify druggable vulnerabilities in CDH1-mutant cancers. METHODS: Density distributions of cell viability data from a genome-wide RNAi screen of isogenic MCF10A and MCF10A-CDH1-/- cells were used to identify protein classes affected by CDH1 mutation. The synthetic lethal relationship between selected protein classes and E-cadherin was characterised by drug sensitivity assays in both the isogenic breast MCF10A cells and CDH1-isogenic gastric NCI-N87. Endocytosis efficiency was quantified using cholera toxin B uptake. Pathway metagene expression of 415 TCGA gastric tumours was statistically correlated with CDH1 expression. RESULTS: MCF10A-CDH1-/- cells showed significantly altered sensitivity to RNAi inhibition of groups of genes including the PI3K/AKT pathway, GPCRs, ion channels, proteosomal subunit proteins and ubiquitinylation enzymes. Both MCF10A-CDH1-/- and NCI-N87-CDH1-/- cells were more sensitive than wild-type cells to compounds that disrupted plasma membrane composition and trafficking, but showed contrasting sensitivities to inhibitors of actin polymerisation and the chloride channel inhibitor NS3728. The MCF10A-CDH1-/- cell lines showed reduced capacity to endocytose cholera toxin B. Pathway metagene analysis identified 20 Reactome pathways that were potentially synthetic lethal in tumours. Genes involved in GPCR signalling, vesicle transport and the metabolism of PI3K and membrane lipids were strongly represented amongst the candidate synthetic lethal genes. CONCLUSIONS: E-cadherin loss leads to disturbances in receptor signalling and plasma membrane trafficking and organisation, creating druggable vulnerabilities.

Population-genetic models of sex-limited genomic imprinting.

Genomic imprinting is a form of epigenetic modification involving parent-of-origin-dependent gene expression, usually the inactivation of one gene copy in some tissues, at least, for some part of the diploid life cycle. Occurring at a number of loci in mammals and flowering plants, this mode of non-Mendelian expression can be viewed more generally as parentally-specific differential gene expression. The effects of natural selection on genetic variation at imprinted loci have previously been examined in a several population-genetic models. Here we expand the existing one-locus, two-allele population-genetic models of viability selection with genomic imprinting to include sex-limited imprinting, i.e., imprinted expression occurring only in one sex, and differential viability between the sexes. We first consider models of complete inactivation of either parental allele and these models are subsequently generalized to incorporate differential expression. Stable polymorphic equilibrium was possible without heterozygote advantage as observed in some prior models of imprinting in both sexes. In contrast to these latter models, in the sex-limited case it was critical whether the paternally inherited or maternally inherited allele was inactivated. The parental origin of inactivated alleles had a different impact on how the population responded to the different selection pressures between the sexes. Under the same fitness parameters, imprinting in the other sex altered the number of possible equilibrium states and their stability. When the parental origin of imprinted alleles and the sex in which they are inactive differ, an allele cannot be inactivated in consecutive generations. The system dynamics became more complex with more equilibrium points emerging. Our results show that selection can interact with epigenetic factors to maintain genetic variation in previously unanticipated ways.