APOBEC2 safeguards skeletal muscle cell fate through binding chromatin and regulating transcription of non-muscle genes during myoblast differentiation

The apolipoprotein B messenger RNA editing enzyme, catalytic polypeptide (APOBEC) family is composed of nucleic acid editors with roles ranging from antibody diversification to RNA editing. APOBEC2, a member of this family with an evolutionarily conserved nucleic acid-binding cytidine deaminase domain, has neither an established substrate nor function. Using a cellular model of muscle differentiation where APOBEC2 is inducibly expressed, we confirmed that APOBEC2 does not have the attributed molecular functions of the APOBEC family, such as RNA editing, DNA demethylation, and DNA mutation. Instead, we found that during muscle differentiation APOBEC2 occupied a specific motif within promoter regions; its removal from those regions resulted in transcriptional changes. Mechanistically, these changes reflect the direct interaction of APOBEC2 with histone deacetylase (HDAC) transcriptional corepressor complexes. We also found that APOBEC2 could bind DNA directly, in a sequence-specific fashion, suggesting that it functions as a recruiter of HDAC to specific genes whose promoters it occupies. These genes are normally suppressed during muscle cell differentiation, and their suppression may contribute to the safeguarding of muscle cell fate. Altogether, our results reveal a unique role for APOBEC2 within the APOBEC family.

The remission status of AML patients after allo-HCT is associated with a distinct single-cell bone marrow T-cell signature

ABSTRACT

Acute myeloid leukemia (AML) is a hematologic malignancy for which allogeneic hematopoietic cell transplantation (allo-HCT) often remains the only curative therapeutic approach. However, incapability of T cells to recognize and eliminate residual leukemia stem cells might lead to an insufficient graft-versus-leukemia (GVL) effect and relapse. Here, we performed single-cell RNA-sequencing (scRNA-seq) on bone marrow (BM) T lymphocytes and CD34+ cells of 6 patients with AML 100 days after allo-HCT to identify T-cell signatures associated with either imminent relapse (REL) or durable complete remission (CR). We observed a higher frequency of cytotoxic CD8+ effector and gamma delta (γδ) T cells in CR vs REL samples. Pseudotime and gene regulatory network analyses revealed that CR CD8+ T cells were more advanced in maturation and had a stronger cytotoxicity signature, whereas REL samples were characterized by inflammatory tumor necrosis factor/NF-κB signaling and an immunosuppressive milieu. We identified ADGRG1/GPR56 as a surface marker enriched in CR CD8+ T cells and confirmed in a CD33-directed chimeric antigen receptor T cell/AML coculture model that GPR56 becomes upregulated on T cells upon antigen encounter and elimination of AML cells. We show that GPR56 continuously increases at the protein level on CD8+ T cells after allo-HCT and confirm faster interferon gamma (IFN-γ) secretion upon re-exposure to matched, but not unmatched, recipient AML cells in the GPR56+ vs GPR56- CD8+ T-cell fraction. Together, our data provide a single-cell reference map of BM-derived T cells after allo-HCT and propose GPR56 expression dynamics as a surrogate for antigen encounter after allo-HCT.

Escape from X inactivation is directly modulated by levels of Xist non-coding RNA

In placental females, one copy of the two X chromosomes is largely silenced during a narrow developmental time window, in a process mediated by the non-coding RNA Xist1. Here, we demonstrate that Xist can initiate X-chromosome inactivation (XCI) well beyond early embryogenesis. By modifying its endogenous level, we show that Xist has the capacity to actively silence genes that escape XCI both in neuronal progenitor cells (NPCs) and in vivo, in mouse embryos. We also show that Xist plays a direct role in eliminating TAD-like structures associated with clusters of escapee genes on the inactive X chromosome, and that this is dependent on Xist's XCI initiation partner, SPEN2. We further demonstrate that Xist's function in suppressing gene expression of escapees and topological domain formation is reversible for up to seven days post-induction, but that sustained Xist up-regulation leads to progressively irreversible silencing and CpG island DNA methylation of facultative escapees. Thus, the distinctive transcriptional and regulatory topologies of the silenced X chromosome is actively, directly - and reversibly - controlled by Xist RNA throughout life.

GRaNIE and GRaNPA: inference and evaluation of enhancer-mediated gene regulatory networks

Enhancers play a vital role in gene regulation and are critical in mediating the impact of noncoding genetic variants associated with complex traits. Enhancer activity is a cell-type-specific process regulated by transcription factors (TFs), epigenetic mechanisms and genetic variants. Despite the strong mechanistic link between TFs and enhancers, we currently lack a framework for jointly analysing them in cell-type-specific gene regulatory networks (GRN). Equally important, we lack an unbiased way of assessing the biological significance of inferred GRNs since no complete ground truth exists. To address these gaps, we present GRaNIE (Gene Regulatory Network Inference including Enhancers) and GRaNPA (Gene Regulatory Network Performance Analysis). GRaNIE (https://git.embl.de/grp-zaugg/GRaNIE) builds enhancer-mediated GRNs based on covariation of chromatin accessibility and RNA-seq across samples (e.g. individuals), while GRaNPA (https://git.embl.de/grp-zaugg/GRaNPA) assesses the performance of GRNs for predicting cell-type-specific differential expression. We demonstrate their power by investigating gene regulatory mechanisms underlying the response of macrophages to infection, cancer and common genetic traits including autoimmune diseases. Finally, our methods identify the TF PURA as a putative regulator of pro-inflammatory macrophage polarisation.

MYT1L haploinsufficiency in human neurons and mice causes autism-associated phenotypes that can be reversed by genetic and pharmacologic intervention

MYT1L is an autism spectrum disorder (ASD)-associated transcription factor that is expressed in virtually all neurons throughout life. How MYT1L mutations cause neurological phenotypes and whether they can be targeted remains enigmatic. Here, we examine the effects of MYT1L deficiency in human neurons and mice. Mutant mice exhibit neurodevelopmental delays with thinner cortices, behavioural phenotypes, and gene expression changes that resemble those of ASD patients. MYT1L target genes, including WNT and NOTCH, are activated upon MYT1L depletion and their chemical inhibition can rescue delayed neurogenesis in vitro. MYT1L deficiency also causes upregulation of the main cardiac sodium channel, SCN5A, and neuronal hyperactivity, which could be restored by shRNA-mediated knockdown of SCN5A or MYT1L overexpression in postmitotic neurons. Acute application of the sodium channel blocker, lamotrigine, also rescued electrophysiological defects in vitro and behaviour phenotypes in vivo. Hence, MYT1L mutation causes both developmental and postmitotic neurological defects. However, acute intervention can normalise resulting electrophysiological and behavioural phenotypes in adulthood.

Convolutional networks for supervised mining of molecular patterns within cellular context

Cryo-electron tomograms capture a wealth of structural information on the molecular constituents of cells and tissues. We present DeePiCt (deep picker in context), an open-source deep-learning framework for supervised segmentation and macromolecular complex localization in cryo-electron tomography. To train and benchmark DeePiCt on experimental data, we comprehensively annotated 20 tomograms of Schizosaccharomyces pombe for ribosomes, fatty acid synthases, membranes, nuclear pore complexes, organelles, and cytosol. By comparing DeePiCt to state-of-the-art approaches on this dataset, we show its unique ability to identify low-abundance and low-density complexes. We use DeePiCt to study compositionally distinct subpopulations of cellular ribosomes, with emphasis on their contextual association with mitochondria and the endoplasmic reticulum. Finally, applying pre-trained networks to a HeLa cell tomogram demonstrates that DeePiCt achieves high-quality predictions in unseen datasets from different biological species in a matter of minutes. The comprehensively annotated experimental data and pre-trained networks are provided for immediate use by the community.

Identifying a novel role for the master regulator Tal1 in the Endothelial to Hematopoietic Transition

Progress in the generation of Hematopoietic Stem and Progenitor Cells (HSPCs) in vitro and ex vivo has been built on the knowledge of developmental hematopoiesis, underscoring the importance of understanding this process. HSPCs emerge within the embryonic vasculature through an Endothelial-to-Hematopoietic Transition (EHT). The transcriptional regulator Tal1 exerts essential functions in the earliest stages of blood development, but is considered dispensable for the EHT. Nevertheless, Tal1 is expressed with its binding partner Lmo2 and it homologous Lyl1 in endothelial and transitioning cells at the time of EHT. Here, we investigated the function of these genes using a mouse embryonic-stem cell (mESC)-based differentiation system to model hematopoietic development. We showed for the first time that the expression of TAL1 in endothelial cells is crucial to ensure the efficiency of the EHT process and a sustained hematopoietic output. Our findings uncover an important function of Tal1 during the EHT, thus filling the current gap in the knowledge of the role of this master gene throughout the whole process of hematopoietic development.

Operation of a TCA cycle subnetwork in the mammalian nucleus

Nucleic acid and histone modifications critically depend on the tricarboxylic acid (TCA) cycle for substrates and cofactors. Although a few TCA cycle enzymes have been reported in the nucleus, the corresponding pathways are considered to operate in mitochondria. Here, we show that a part of the TCA cycle is operational also in the nucleus. Using 13C-tracer analysis, we identified activity of glutamine-to-fumarate, citrate-to-succinate, and glutamine-to-aspartate routes in the nuclei of HeLa cells. Proximity labeling mass spectrometry revealed a spatial vicinity of the involved enzymes with core nuclear proteins. We further show nuclear localization of aconitase 2 and 2-oxoglutarate dehydrogenase in mouse embryonic stem cells. Nuclear localization of the latter enzyme, which produces succinyl-CoA, changed from pluripotency to a differentiated state with accompanying changes in the nuclear protein succinylation. Together, our results demonstrate operation of an extended metabolic pathway in the nucleus, warranting a revision of the canonical view on metabolic compartmentalization.

Drug-microenvironment perturbations reveal resistance mechanisms and prognostic subgroups in CLL

The tumour microenvironment and genetic alterations collectively influence drug efficacy in cancer, but current evidence is limited and systematic analyses are lacking. Using chronic lymphocytic leukaemia (CLL) as a model disease, we investigated the influence of 17 microenvironmental stimuli on 12 drugs in 192 genetically characterised patient samples. Based on microenvironmental response, we identified four subgroups with distinct clinical outcomes beyond known prognostic markers. Response to multiple microenvironmental stimuli was amplified in trisomy 12 samples. Trisomy 12 was associated with a distinct epigenetic signature. Bromodomain inhibition reversed this epigenetic profile and could be used to target microenvironmental signalling in trisomy 12 CLL. We quantified the impact of microenvironmental stimuli on drug response and their dependence on genetic alterations, identifying interleukin 4 (IL4) and Toll‐like receptor (TLR) stimulation as the strongest actuators of drug resistance. IL4 and TLR signalling activity was increased in CLL‐infiltrated lymph nodes compared with healthy samples. High IL4 activity correlated with faster disease progression. The publicly available dataset can facilitate the investigation of cell‐extrinsic mechanisms of drug resistance and disease progression.

Comparative chromatin accessibility upon BDNF stimulation delineates neuronal regulatory elements

Neuronal stimulation induced by the brain-derived neurotrophic factor (BDNF) triggers gene expression, which is crucial for neuronal survival, differentiation, synaptic plasticity, memory formation, and neurocognitive health. However, its role in chromatin regulation is unclear. Here, using temporal profiling of chromatin accessibility and transcription in mouse primary cortical neurons upon either BDNF stimulation or depolarization (KCl), we identify features that define BDNF-specific chromatin-to-gene expression programs. Enhancer activation is an early event in the regulatory control of BDNF-treated neurons, where the bZIP motif-binding Fos protein pioneered chromatin opening and cooperated with co-regulatory transcription factors (Homeobox, EGRs, and CTCF) to induce transcription. Deleting cis-regulatory sequences affect BDNF-mediated Arc expression, a regulator of synaptic plasticity. BDNF-induced accessible regions are linked to preferential exon usage by neurodevelopmental disorder-related genes and the heritability of neuronal complex traits, which were validated in human iPSC-derived neurons. Thus, we provide a comprehensive view of BDNF-mediated genome regulatory features using comparative genomic approaches to dissect mammalian neuronal stimulation.

CDK7/12/13 inhibition targets an oscillating leukemia stem cell network and synergizes with venetoclax in acute myeloid leukemia

The heterogeneous response of acute myeloid leukemia (AML) to current anti‐leukemic therapies is only partially explained by mutational heterogeneity. We previously identified GPR56 as a surface marker associated with poor outcome across genetic groups, which characterizes two leukemia stem cell (LSC)‐enriched compartments with different self‐renewal capacities. How these compartments self‐renew remained unclear. Here, we show that GPR56⁺ LSC compartments are promoted in a complex network involving epithelial‐to‐mesenchymal transition (EMT) regulators besides Rho, Wnt, and Hedgehog (Hh) signaling. Unexpectedly, Wnt pathway inhibition increased the more immature, slowly cycling GPR56⁺CD34⁺ fraction and Hh/EMT gene expression, while Wnt activation caused opposite effects. Our data suggest that the crucial role of GPR56 lies in its ability to co‐activate these opposing signals, thus ensuring the constant supply of both LSC subsets. We show that CDK7 inhibitors suppress both LSC‐enriched subsets in vivo and synergize with the Bcl‐2 inhibitor venetoclax. Our data establish reciprocal transition between LSC compartments as a novel concept underlying the poor outcome in GPR56^high AML and propose combined CDK7 and Bcl‐2 inhibition as LSC‐directed therapy in this disease.

Enhancers in disease: molecular basis and emerging treatment strategies

Enhancers are genomic sequences that play a key role in regulating tissue-specific gene expression levels. An increasing number of diseases are linked to impaired enhancer function through chromosomal rearrangement, genetic variation within enhancers, or epigenetic modulation. Here, we review how these enhancer disruptions have recently been implicated in congenital disorders, cancers, and common complex diseases and address the implications for diagnosis and treatment. Although further fundamental research into enhancer function, target genes, and context is required, enhancer-targeting drugs and gene editing approaches show great therapeutic promise for a range of diseases.

Transcription factors: Bridge between cell signaling and gene regulation

Transcription factors (TFs) are key regulators of intrinsic cellular processes, such as differentiation and development, and of the cellular response to external perturbation through signaling pathways. In this review we focus on the role of TFs as a link between signaling pathways and gene regulation. Cell signaling tends to result in the modulation of a set of TFs that then lead to changes in the cell's transcriptional program. We highlight the molecular layers at which TF activity can be measured and the associated technical and conceptual challenges. These layers include post-translational modifications (PTMs) of the TF, regulation of TF binding to DNA through chromatin accessibility and epigenetics, and expression of target genes. We highlight that a large number of TFs are understudied in both signaling and gene regulation studies, and that our knowledge about known TF targets has a strong literature bias. We argue that TFs serve as a perfect bridge between the fields of gene regulation and signaling, and that separating these fields hinders our understanding of cell functions. Multi-omics approaches that measure multiple dimensions of TF activity are ideally suited to study the interplay of cell signaling and gene regulation using TFs as the anchor to link the two fields.

Donor cell memory confers a metastable state of directly converted cells

Generation of induced oligodendrocyte progenitor cells (iOPCs) from somatic fibroblasts is a strategy for cell-based therapy of myelin diseases. However, iOPC generation is inefficient, and the resulting iOPCs exhibit limited expansion and differentiation competence. Here we overcome these limitations by transducing an optimized transcription factor combination into a permissive donor phenotype, the pericyte. Pericyte-derived iOPCs (PC-iOPCs) are stably expandable and functionally myelinogenic with high differentiation competence. Unexpectedly, however, we found that PC-iOPCs are metastable so that they can produce myelination-competent oligodendrocytes or revert to their original identity in a context-dependent fashion. Phenotypic reversion of PC-iOPCs is tightly linked to memory of their original transcriptome and epigenome. Phenotypic reversion can be disconnected from this donor cell memory effect, and in vivo myelination can eventually be achieved by transplantation of O4⁺ pre-oligodendrocytes. Our data show that donor cell source and memory can contribute to the fate and stability of directly converted cells.

Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Human Developmental Hematopoiesis

Regulation of hematopoiesis during human development remains poorly defined. Here we applied single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) to over 8,000 human immunophenotypic blood cells from fetal liver and bone marrow. We inferred their differentiation trajectory and identified three highly proliferative oligopotent progenitor populations downstream of hematopoietic stem cells (HSCs)/multipotent progenitors (MPPs). Along this trajectory, we observed opposing patterns of chromatin accessibility and differentiation that coincided with dynamic changes in the activity of distinct lineage-specific transcription factors. Integrative analysis of chromatin accessibility and gene expression revealed extensive epigenetic but not transcriptional priming of HSCs/MPPs prior to their lineage commitment. Finally, we refined and functionally validated the sorting strategy for the HSCs/MPPs and achieved around 90% enrichment. Our study provides a useful framework for future investigation of human developmental hematopoiesis in the context of blood pathologies and regenerative medicine.

Predictive features of gene expression variation reveal mechanistic link with differential expression

For most biological processes, organisms must respond to extrinsic cues, while maintaining essential gene expression programmes. Although studied extensively in single cells, it is still unclear how variation is controlled in multicellular organisms. Here, we used a machine-learning approach to identify genomic features that are predictive of genes with high versus low variation in their expression across individuals, using bulk data to remove stochastic cell-to-cell variation. Using embryonic gene expression across 75 Drosophila isogenic lines, we identify features predictive of expression variation (controlling for expression level), many of which are promoter-related. Genes with low variation fall into two classes reflecting different mechanisms to maintain robust expression, while genes with high variation seem to lack both types of stabilizing mechanisms. Applying this framework to humans revealed similar predictive features, indicating that promoter architecture is an ancient mechanism to control expression variation. Remarkably, expression variation features could also partially predict differential expression after diverse perturbations in both Drosophila and humans. Differential gene expression signatures may therefore be partially explained by genetically encoded gene-specific features, unrelated to the studied treatment.

Landscape of cohesin-mediated chromatin loops in the human genome

Physical interactions between distal regulatory elements have a key role in regulating gene expression, but the extent to which these interactions vary between cell types and contribute to cell-type-specific gene expression remains unclear. Here, to address these questions as part of phase III of the Encyclopedia of DNA Elements (ENCODE), we mapped cohesin-mediated chromatin loops, using chromatin interaction analysis by paired-end tag sequencing (ChIA-PET), and analysed gene expression in 24 diverse human cell types, including core ENCODE cell lines. Twenty-eight per cent of all chromatin loops vary across cell types; these variations modestly correlate with changes in gene expression and are effective at grouping cell types according to their tissue of origin. The connectivity of genes corresponds to different functional classes, with housekeeping genes having few contacts, and dosage-sensitive genes being more connected to enhancer elements. This atlas of chromatin loops complements the diverse maps of regulatory architecture that comprise the ENCODE Encyclopedia, and will help to support emerging analyses of genome structure and function.

Genomic Rewiring of SOX2 Chromatin Interaction Network during Differentiation of ESCs to Postmitotic Neurons

Cellular differentiation requires dramatic changes in chromatin organization, transcriptional regulation, and protein production. To understand the regulatory connections between these processes, we generated proteomic, transcriptomic, and chromatin accessibility data during differentiation of mouse embryonic stem cells (ESCs) into postmitotic neurons and found extensive associations between different molecular layers within and across differentiation time points. We observed that SOX2, as a regulator of pluripotency and neuronal genes, redistributes from pluripotency enhancers to neuronal promoters during differentiation, likely driven by changes in its protein interaction network. We identified ATRX as a major SOX2 partner in neurons, whose co-localization correlated with an increase in active enhancer marks and increased expression of nearby genes, which we experimentally confirmed for three loci. Collectively, our data provide key insights into the regulatory transformation of SOX2 during neuronal differentiation, and we highlight the significance of multi-omic approaches in understanding gene regulation in complex systems.

•

Complex interplay of RNA, protein, and chromatin during neuronal differentiation

•

Multi-omic profiling reveals divergent roles of SOX2 in stem cells and neurons

•

SOX2 on-chromatin interaction network changes from pluripotent to neuronal factors

•

ATRX interacts with SOX2 in neurons and co-binds highly expressed neuronal genes

Lysine 4 of histone H3.3 is required for embryonic stem cell differentiation, histone enrichment at regulatory regions and transcription accuracy

Mutations in enzymes that modify histone H3 at lysine 4 (H3K4) or lysine 36 (H3K36) have been linked to human disease, yet the role of these residues in mammals is unclear. We mutated K4 or K36 to alanine in the histone variant H3.3 and showed that the K4A mutation in mouse embryonic stem cells (ESCs) impaired differentiation and induced widespread gene expression changes. K4A resulted in substantial H3.3 depletion, especially at ESC promoters; it was accompanied by reduced remodeler binding and increased RNA polymerase II (Pol II) activity. Regulatory regions depleted of H3.3K4A showed histone modification alterations and changes in enhancer activity that correlated with gene expression. In contrast, the K36A mutation did not alter H3.3 deposition and affected gene expression at the later stages of differentiation. Thus, H3K4 is required for nucleosome deposition, histone turnover and chromatin remodeler binding at regulatory regions, where tight regulation of Pol II activity is necessary for proper ESC differentiation.

Mechanistic insights into transcription factor cooperativity and its impact on protein-phenotype interactions

Recent high-throughput transcription factor (TF) binding assays revealed that TF cooperativity is a widespread phenomenon. However, a global mechanistic and functional understanding of TF cooperativity is still lacking. To address this, here we introduce a statistical learning framework that provides structural insight into TF cooperativity and its functional consequences based on next generation sequencing data. We identify DNA shape as driver for cooperativity, with a particularly strong effect for Forkhead-Ets pairs. Follow-up experiments reveal a local shape preference at the Ets-DNA-Forkhead interface and decreased cooperativity upon loss of the interaction. Additionally, we discover many functional associations for cooperatively bound TFs. Examination of the link between FOXO1:ETV6 and lymphomas reveals that their joint expression levels improve patient clinical outcome stratification. Altogether, our results demonstrate that inter-family cooperative TF binding is driven by position-specific DNA readout mechanisms, which provides an additional regulatory layer for downstream biological functions.

Although transcription factor (TF) cooperativity is widespread, a global mechanistic understanding of the role of TF cooperativity is still lacking. Here the authors introduce a statistical learning framework that provides structural insight into TF cooperativity and its functional consequences based on next generation sequencing data and provide mechanistic insights into TF cooperativity and its impact on protein-phenotype interactions.

TET2 binding to enhancers facilitates transcription factor recruitment in hematopoietic cells

The epigenetic regulator TET2 is frequently mutated in hematological diseases. Mutations have been shown to arise in hematopoietic stem cells early in disease development and lead to altered DNA methylation landscapes and an increased risk of hematopoietic malignancy. Here, we show by genome-wide mapping of TET2 binding sites in different cell types that TET2 localizes to regions of open chromatin and cell-type-specific enhancers. We find that deletion of Tet2 in native hematopoiesis as well as fully transformed acute myeloid leukemia (AML) results in changes in transcription factor (TF) activity within these regions, and we provide evidence that loss of TET2 leads to attenuation of chromatin binding of members of the basic helix-loop-helix (bHLH) TF family. Together, these findings demonstrate that TET2 activity shapes the local chromatin environment at enhancers to facilitate TF binding and provides an example of how epigenetic dysregulation can affect gene expression patterns and drive disease development.

Haplotype-specific MAPT exon 3 expression regulated by common intronic polymorphisms associated with Parkinsonian disorders

BACKGROUND

Genome wide association studies have identified microtubule associated protein tau (MAPT) H1 haplotype single nucleotide polymorphisms (SNPs) as leading common risk variants for Parkinson's disease, progressive supranuclear palsy and corticobasal degeneration. The MAPT risk variants fall within a large 1.8 Mb region of high linkage disequilibrium, making it difficult to discern the functionally important risk variants. Here, we leverage the strong haplotype-specific expression of MAPT exon 3 to investigate the functionality of SNPs that fall within this H1 haplotype region of linkage disequilibrium.

METHODS

In this study, we dissect the molecular mechanisms by which haplotype-specific SNPs confer allele-specific effects on the alternative splicing of MAPT exon 3. Firstly, we use haplotype-hybrid whole-locus genomic MAPT vectors studies to identify functional SNPs. Next, we characterise the RNA-protein interactions at two loci by mass spectrometry. Lastly, we knockdown candidate splice factors to determine their effect on MAPT exon 3 using a novel allele-specific qPCR assay.

RESULTS

Using whole-locus genomic DNA expression vectors to express MAPT haplotype variants, we demonstrate that rs17651213 regulates exon 3 inclusion in a haplotype-specific manner. We further investigated the functionality of this region using RNA-electrophoretic mobility shift assays to show differential RNA-protein complex formation at the H1 and H2 sequence variants of SNP rs17651213 and rs1800547 and subsequently identified candidate trans-acting splicing factors interacting with these functional SNPs sequences by RNA-protein pull-down experiment and mass spectrometry. Finally, gene knockdown of candidate splice factors identified by mass spectrometry demonstrate a role for hnRNP F and hnRNP Q in the haplotype-specific regulation of exon 3 inclusion.

CONCLUSIONS

We identified common splice factors hnRNP F and hnRNP Q regulating the haplotype-specific splicing of MAPT exon 3 through intronic variants rs1800547 and rs17651213. This work demonstrates an integrated approach to characterise the functionality of risk variants in large regions of linkage disequilibrium.

Backmasking in the yeast genome: encoding overlapping information for protein-coding and RNA degradation

Backmasking is a recording technique used to hide a sound or message in a music track in reverse, meaning that it is only audible when the record is played backwards. Analogously, the compact yeast genome encodes for diverse sources of information such as overlapping coding and non-coding transcripts, and protein-binding sites on the two complementary DNA strands. Examples are the consensus binding site sequences of the RNA-binding proteins Nrd1 and Nab3 that target non-coding transcripts for degradation. Here, by examining the overlap of stable (SUTs, stable unannotated transcripts) and unstable (CUTs, cryptic unstable transcripts) transcripts with protein-coding genes, we show that the predicted Nrd1 and Nab3-binding site sequences occur at differing frequencies. They are always depleted in the sense direction of protein-coding genes, thus avoiding degradation of the transcript. However in the antisense direction, predicted binding sites occur at high frequencies in genes with overlapping unstable ncRNAs (CUTs), so limiting the availability of non-functional transcripts. In contrast they are depleted in genes with overlapping stable ncRNAs (SUTs), presumably to avoid degrading the non-coding transcript. The protein-coding genes maintain similar amino-acid contents, but they display distinct codon usages so that Nrd1 and Nab3-binding sites can arise at differing frequencies in antisense depending on the overlapping transcript type. Our study demonstrates how yeast has evolved to encode multiple layers of information-protein-coding genes in one strand and the relative chance of degrading antisense RNA in the other strand-in the same regions of a compact genome.

SNPhood: investigate, quantify and visualise the epigenomic neighbourhood of SNPs using NGS data

MOTIVATION

The vast majority of the many thousands of disease-associated single nucleotide polymorphisms (SNPs) lie in the non-coding part of the genome. They are likely to affect regulatory elements, such as enhancers and promoters, rather than the function of a protein. To understand the molecular mechanisms underlying genetic diseases, it is therefore increasingly important to study the effect of a SNP on nearby molecular traits such as chromatin or transcription factor binding.

RESULTS

We developed SNPhood, a user-friendly Bioconductor R package to investigate, quantify and visualise the local epigenetic neighbourhood of a set of SNPs in terms of chromatin marks or TF binding sites using data from NGS experiments.

AVAILABILITY AND IMPLEMENTATION

SNPhood is publicly available and maintained as an R Bioconductor package at http://bioconductor.org/packages/SNPhood/

CONTACT

judith.zaugg@embl.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genetic Control of Chromatin States in Humans Involves Local and Distal Chromosomal Interactions

Deciphering the impact of genetic variants on gene regulation is fundamental to understanding human disease. Although gene regulation often involves long-range interactions, it is unknown to what extent non-coding genetic variants influence distal molecular phenotypes. Here, we integrate chromatin profiling for three histone marks in lymphoblastoid cell lines (LCLs) from 75 sequenced individuals with LCL-specific Hi-C and ChIA-PET-based chromatin contact maps to uncover one of the largest collections of local and distal histone quantitative trait loci (hQTLs). Distal QTLs are enriched within topologically associated domains and exhibit largely concordant variation of chromatin state coordinated by proximal and distal non-coding genetic variants. Histone QTLs are enriched for common variants associated with autoimmune diseases and enable identification of putative target genes of disease-associated variants from genome-wide association studies. These analyses provide insights into how genetic variation can affect human disease phenotypes by coordinated changes in chromatin at interacting regulatory elements.

Influence of the bleaching interval on the luminosity of long-term discolored enamel-dentin discs

OBJECTIVES

The aim of this study is to investigate the influence of changing the sodium perborate-tetrahydrate (PBS-4) at a 4-day interval versus no change after 16 days of internal bleaching.

MATERIALS AND METHODS

Two hundred and ten bovine enamel-dentin discs were discolored for 3.5 years with 14 different endodontic materials. All groups with a discoloring index of ∆E (mean) ≥ 5.5 were included in the present investigation: ApexCal (APCA), MTA white + blood (WMTA+BL), Portland cement + blood (PC+BL), blood (BL), MTA gray (GMTA), MTA gray + blood (GMTA+BL), Ledermix (LED), and triple antibiotic paste containing minocycline (3Mix). Fourteen specimens of each group were randomly assigned into two treatment groups: (1) no change of the PBS-4 (n = 7); (2) change of the PBS-4 every 4 days (n = 7). Color measurements were taken at 10 different time intervals and the L*a*b* values were recorded with a spectrophotometer (VITA Easyshade® compact).

RESULTS

In the group 3Mix, significantly better results were achieved by changing the bleaching agent every 4 days (P = 0.0049; q = 0.04), while the group WMTA+BL indicated better results by no change of the bleaching agent (P = 0.0222, q = 0.09). All remaining groups showed no statistical difference between the two treatment procedures.

CONCLUSIONS

Moderate discolorations can be successfully treated without changing the bleaching agent over a period of 16 days. Changing the sodium perborate-tetrahydrate every 4 days is preferred in case of severe discolored enamel-dentin discs only.

CLINICAL RELEVANCE

This approach may offer a reduced number of clinical appointments and a secondary cost reduction to the patient.

Role of histone modifications and early termination in pervasive transcription and antisense-mediated gene silencing in yeast

Most genomes, including yeast Saccharomyces cerevisiae, are pervasively transcribed producing numerous non-coding RNAs, many of which are unstable and eliminated by nuclear or cytoplasmic surveillance pathways. We previously showed that accumulation of PHO84 antisense RNA (asRNA), in cells lacking the nuclear exosome component Rrp6, is paralleled by repression of sense transcription in a process dependent on the Hda1 histone deacetylase (HDAC) and the H3K4 histone methyl transferase Set1. Here we investigate this process genome-wide and measure the whole transcriptome of various histone modification mutants in a Δrrp6 strain using tiling arrays. We confirm widespread occurrence of potentially antisense-dependent gene regulation and identify three functionally distinct classes of genes that accumulate asRNAs in the absence of Rrp6. These classes differ in whether the genes are silenced by the asRNA and whether the silencing is HDACs and histone methyl transferase-dependent. Among the distinguishing features of asRNAs with regulatory potential, we identify weak early termination by Nrd1/Nab3/Sen1, extension of the asRNA into the open reading frame promoter and dependence of the silencing capacity on Set1 and the HDACs Hda1 and Rpd3 particularly at promoters undergoing extensive chromatin remodelling. Finally, depending on the efficiency of Nrd1/Nab3/Sen1 early termination, asRNA levels are modulated and their capability of silencing is changed.

Extensive variation in chromatin states across humans

The majority of disease-associated variants lie outside protein-coding regions, suggesting a link between variation in regulatory regions and disease predisposition. We studied differences in chromatin states using five histone modifications, cohesin, and CTCF in lymphoblastoid lines from 19 individuals of diverse ancestry. We found extensive signal variation in regulatory regions, which often switch between active and repressed states across individuals. Enhancer activity is particularly diverse among individuals, whereas gene expression remains relatively stable. Chromatin variability shows genetic inheritance in trios, correlates with genetic variation and population divergence, and is associated with disruptions of transcription factor binding motifs. Overall, our results provide insights into chromatin variation among humans.

Gene loops enhance transcriptional directionality

Eukaryotic genomes are extensively transcribed, forming both messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). ncRNAs made by RNA polymerase II often initiate from bidirectional promoters (nucleosome-depleted chromatin) that synthesize mRNA and ncRNA in opposite directions. We demonstrate that, by adopting a gene-loop conformation, actively transcribed mRNA encoding genes restrict divergent transcription of ncRNAs. Because gene-loop formation depends on a protein factor (Ssu72) that coassociates with both the promoter and the terminator, the inactivation of Ssu72 leads to increased synthesis of promoter-associated divergent ncRNAs, referred to as Ssu72-restricted transcripts (SRTs). Similarly, inactivation of individual gene loops by gene mutation enhances SRT synthesis. We demonstrate that gene-loop conformation enforces transcriptional directionality on otherwise bidirectional promoters.

A genomic model of condition-specific nucleosome behavior explains transcriptional activity in yeast

Nucleosomes play an important role in gene regulation. Molecular studies observed that nucleosome binding in promoters tends to be repressive. In contrast, genomic studies have delivered conflicting results: An analysis of yeast grown on diverse carbon sources reported that nucleosome occupancies remain largely unchanged between conditions, whereas a study of the heat-shock response suggested that nucleosomes get evicted at promoters of genes with increased expression. Consequently, there are few general principles that capture the relationship between chromatin organization and transcriptional regulation. Here, we present a qualitative model for nucleosome positioning in Saccharomyces cerevisiae that helps explain important properties of gene expression. By integrating publicly available data sets, we observe that promoter-bound nucleosomes assume one of four discrete configurations that determine the active and silent transcriptional states of a gene, but not its expression level. In TATA-box-containing promoters, nucleosome architecture indicates the amount of transcriptional noise. We show that >20% of genes switch promoter states upon changes in cellular conditions. The data suggest that DNA-binding transcription factors together with chromatin-remodeling enzymes are primarily responsible for the nucleosome architecture. Our model for promoter nucleosome architecture reconciles genome-scale findings with molecular studies; in doing so, we establish principles for nucleosome positioning and gene expression that apply not only to individual genes, but across the entire genome. The study provides a stepping stone for future models of transcriptional regulation that encompass the intricate interplay between cis- and trans-acting factors, chromatin, and the core transcriptional machinery.

N5-CAIR mutase: role of a CO2 binding site and substrate movement in catalysis

N5-Carboxyaminoimidazole ribonucleotide mutase (N5-CAIR mutase or PurE) from Escherichia coli catalyzes the reversible interconversion of N5-CAIR to carboxyaminoimidazole ribonucleotide (CAIR) with direct CO2 transfer. Site-directed mutagenesis, a pH-rate profile, DFT calculations, and X-ray crystallography together provide new insight into the mechanism of this unusual transformation. These studies suggest that a conserved, protonated histidine (His45) plays an essential role in catalysis. The importance of proton transfers is supported by DFT calculations on CAIR and N5-CAIR analogues in which the ribose 5'-phosphate is replaced with a methyl group. The calculations suggest that the nonaromatic tautomer of CAIR (isoCAIR) is only 3.1 kcal/mol higher in energy than its aromatic counterpart, implicating this species as a potential intermediate in the PurE-catalyzed reaction. A structure of wild-type PurE cocrystallized with 4-nitroaminoimidazole ribonucleotide (NO2-AIR, a CAIR analogue) and structures of H45N and H45Q PurEs soaked with CAIR have been determined and provide the first insight into the binding of an intact PurE substrate. A comparison of 19 available structures of PurE and PurE mutants in apo and nucleotide-bound forms reveals a common, buried carboxylate or CO2 binding site for CAIR and N5-CAIR in a hydrophobic pocket in which the carboxylate or CO2 interacts with backbone amides. This work has led to a mechanistic proposal in which the carboxylate orients the substrate for proton transfer from His45 to N5-CAIR to form an enzyme-bound aminoimidazole ribonucleotide (AIR) and CO2 intermediate. Subsequent movement of the aminoimidazole moiety of AIR reorients it for addition of CO2 at C4 to generate isoCAIR. His45 is now in a position to remove a C4 proton to produce CAIR.