Polymorphic short tandem repeats make widespread contributions to blood and serum traits.

Short tandem repeats (STRs) are genomic regions consisting of repeated sequences of 1-6 bp in succession. Single-nucleotide polymorphism (SNP)-based genome-wide association studies (GWASs) do not fully capture STR effects. To study these effects, we imputed 445,720 STRs into genotype arrays from 408,153 White British UK Biobank participants and tested for association with 44 blood phenotypes. Using two fine-mapping methods, we identify 119 candidate causal STR-trait associations and estimate that STRs account for 5.2%-7.6% of causal variants identifiable from GWASs for these traits. These are among the strongest associations for multiple phenotypes, including a coding CTG repeat associated with apolipoprotein B levels, a promoter CGG repeat with platelet traits, and an intronic poly(A) repeat with mean platelet volume. Our study suggests that STRs make widespread contributions to complex traits, provides stringently selected candidate causal STRs, and demonstrates the need to consider a more complete view of genetic variation in GWASs.

A novel quantitative trait locus implicates Msh3 in the propensity for genome-wide short tandem repeat expansions in mice.

Short tandem repeats (STRs) are a class of rapidly mutating genetic elements typically characterized by repeated units of 1-6 bp. We leveraged whole-genome sequencing data for 152 recombinant inbred (RI) strains from the BXD family of mice to map loci that modulate genome-wide patterns of new mutations arising during parent-to-offspring transmission at STRs. We defined quantitative phenotypes describing the numbers and types of germline STR mutations in each strain and performed quantitative trait locus (QTL) analyses for each of these phenotypes. We identified a locus on Chromosome 13 at which strains inheriting the C57BL/6J (B) haplotype have a higher rate of STR expansions than those inheriting the DBA/2J (D) haplotype. The strongest candidate gene in this locus is Msh3, a known modifier of STR stability in cancer and at pathogenic repeat expansions in mice and humans, as well as a current drug target against Huntington's disease. The D haplotype at this locus harbors a cluster of variants near the 5' end of Msh3, including multiple missense variants near the DNA mismatch recognition domain. In contrast, the B haplotype contains a unique retrotransposon insertion. The rate of expansion covaries positively with Msh3 expression-with higher expression from the B haplotype. Finally, detailed analysis of mutation patterns showed that strains carrying the B allele have higher expansion rates, but slightly lower overall total mutation rates, compared with those with the D allele, particularly at tetranucleotide repeats. Our results suggest an important role for inherited variants in Msh3 in modulating genome-wide patterns of germline mutations at STRs.

Mitotic H3K9ac is controlled by phase-specific activity of HDAC2, HDAC3, and SIRT1.

Histone acetylation levels are reduced during mitosis. To study the mitotic regulation of H3K9ac, we used an array of inhibitors targeting specific histone deacetylases. We evaluated the involvement of the targeted enzymes in regulating H3K9ac during all mitotic stages by immunofluorescence and immunoblots. We identified HDAC2, HDAC3, and SIRT1 as modulators of H3K9ac mitotic levels. HDAC2 inhibition increased H3K9ac levels in prophase, whereas HDAC3 or SIRT1 inhibition increased H3K9ac levels in metaphase. Next, we performed ChIP-seq on mitotic-arrested cells following targeted inhibition of these histone deacetylases. We found that both HDAC2 and HDAC3 have a similar impact on H3K9ac, and inhibiting either of these two HDACs substantially increases the levels of this histone acetylation in promoters, enhancers, and insulators. Altogether, our results support a model in which H3K9 deacetylation is a stepwise process-at prophase, HDAC2 modulates most transcription-associated H3K9ac-marked loci, and at metaphase, HDAC3 maintains the reduced acetylation, whereas SIRT1 potentially regulates H3K9ac by impacting HAT activity.

The use of base editing technology to characterize single nucleotide variants.

Single nucleotide variants (SNVs) represent the most common type of polymorphism in the human genome. However, in many cases the phenotypic impacts of such variants are not well understood. Intriguingly, while some SNVs cause debilitating diseases, other variants in the same gene may have no, or limited, impact. The mechanisms underlying these complex patterns are difficult to study at scale. Additionally, current data and research is mainly focused on European populations, and the mechanisms underlying genetic traits in other populations are poorly studied. Novel technologies may be able to mitigate this disparity and improve the applicability of personalized healthcare to underserved populations. In this review we discuss base editing technologies and their potential to accelerate progress in this field, particularly in combination with single-cell RNA sequencing. We further explore how base editing screens can help link SNVs to distinct disease phenotypes. We then highlight several studies that take advantage of single-cell RNA sequencing and CRISPR screens to emphasize the current limitations and future potential of this technique. Lastly, we consider the use of such approaches to potentially accelerate the study of genetic mechanisms in non-European populations.

Missing heritability may be hiding in repeats.

Unexplored variable number tandem repeats make a large contribution to complex traits.

A flexible ChIP-sequencing simulation toolkit.

BACKGROUND: A major challenge in evaluating quantitative ChIP-seq analyses, such as peak calling and differential binding, is a lack of reliable ground truth data. Accurate simulation of ChIP-seq data can mitigate this challenge, but existing frameworks are either too cumbersome to apply genome-wide or unable to model a number of important experimental conditions in ChIP-seq. RESULTS: We present ChIPs, a toolkit for rapidly simulating ChIP-seq data using statistical models of key experimental steps. We demonstrate how ChIPs can be used for a range of applications, including benchmarking analysis tools and evaluating the impact of various experimental parameters. ChIPs is implemented as a standalone command-line program written in C++ and is available from https://github.com/gymreklab/chips . CONCLUSIONS: ChIPs is an efficient ChIP-seq simulation framework that generates realistic datasets over a flexible range of experimental conditions. It can serve as an important component in various ChIP-seq analyses where ground truth data are needed.

The impact of short tandem repeat variation on gene expression.

Short tandem repeats (STRs) have been implicated in a variety of complex traits in humans. However, genome-wide studies of the effects of STRs on gene expression thus far have had limited power to detect associations and provide insights into putative mechanisms. Here, we leverage whole-genome sequencing and expression data for 17 tissues from the Genotype-Tissue Expression Project to identify more than 28,000 STRs for which repeat number is associated with expression of nearby genes (eSTRs). We use fine-mapping to quantify the probability that each eSTR is causal and characterize the top 1,400 fine-mapped eSTRs. We identify hundreds of eSTRs linked with published genome-wide association study signals and implicate specific eSTRs in complex traits, including height, schizophrenia, inflammatory bowel disease and intelligence. Overall, our results support the hypothesis that eSTRs contribute to a range of human phenotypes, and our data should serve as a valuable resource for future studies of complex traits.

The Histone Deacetylase SIRT6 Restrains Transcription Elongation via Promoter-Proximal Pausing.

Transcriptional regulation in eukaryotes occurs at promoter-proximal regions wherein transcriptionally engaged RNA polymerase II (Pol II) pauses before proceeding toward productive elongation. The role of chromatin in pausing remains poorly understood. Here, we demonstrate that the histone deacetylase SIRT6 binds to Pol II and prevents the release of the negative elongation factor (NELF), thus stabilizing Pol II promoter-proximal pausing. Genetic depletion of SIRT6 or its chromatin deficiency upon glucose deprivation causes intragenic enrichment of acetylated histone H3 at lysines 9 (H3K9ac) and 56 (H3K56ac), activation of cyclin-dependent kinase 9 (CDK9)-that phosphorylates NELF and the carboxyl terminal domain of Pol II-and enrichment of the positive transcription elongation factors MYC, BRD4, PAF1, and the super elongation factors AFF4 and ELL2. These events lead to increased expression of genes involved in metabolism, protein synthesis, and embryonic development. Our results identified SIRT6 as a Pol II promoter-proximal pausing-dedicated histone deacetylase.

H3K27me3-mediated silencing of structural genes is required for zebrafish heart regeneration.

Deciphering the genetic and epigenetic regulation of cardiomyocyte proliferation in organisms that are capable of robust cardiac renewal, such as zebrafish, represents an attractive inroad towards regenerating the human heart. Using integrated high-throughput transcriptional and chromatin analyses, we have identified a strong association between H3K27me3 deposition and reduced sarcomere and cytoskeletal gene expression in proliferative cardiomyocytes following cardiac injury in zebrafish. To move beyond an association, we generated an inducible transgenic strain expressing a mutant version of histone 3, H3.3K27M, that inhibits H3K27me3 catalysis in cardiomyocytes during the regenerative window. Hearts comprising H3.3K27M-expressing cardiomyocytes fail to regenerate, with wound edge cells showing heightened expression of structural genes and prominent sarcomeres. Although cell cycle re-entry was unperturbed, cytokinesis and wound invasion were significantly compromised. Collectively, our study identifies H3K27me3-mediated silencing of structural genes as requisite for zebrafish heart regeneration and suggests that repression of similar structural components in the border zone of an infarcted human heart might improve its regenerative capacity.

Role of transcription complexes in the formation of the basal methylation pattern in early development.

Following erasure in the blastocyst, the entire genome undergoes de novo methylation at the time of implantation, with CpG islands being protected from this process. This bimodal pattern is then preserved throughout development and the lifetime of the organism. Using mouse embryonic stem cells as a model system, we demonstrate that the binding of an RNA polymerase complex on DNA before de novo methylation is predictive of it being protected from this modification, and tethering experiments demonstrate that the presence of this complex is, in fact, sufficient to prevent methylation at these sites. This protection is most likely mediated by the recruitment of enzyme complexes that methylate histone H3K4 over a local region and, in this way, prevent access to the de novo methylation complex. The topological pattern of H3K4me3 that is formed while the DNA is as yet unmethylated provides a strikingly accurate template for modeling the genome-wide basal methylation pattern of the organism. These results have far-reaching consequences for understanding the relationship between RNA transcription and DNA methylation.

Study of mitotic chromatin supports a model of bookmarking by histone modifications and reveals nucleosome deposition patterns.

Mitosis encompasses key molecular changes including chromatin condensation, nuclear envelope breakdown, and reduced transcription levels. Immediately after mitosis, the interphase chromatin structure is reestablished and transcription resumes. The reestablishment of the interphase chromatin is probably achieved by "bookmarking," i.e., the retention of at least partial information during mitosis. To gain a deeper understanding of the contribution of histone modifications to the mitotic bookmarking process, we merged proteomics, immunofluorescence, and ChIP-seq approaches. We focused on key histone modifications and employed HeLa-S3 cells as a model system. Generally, in spite of the general hypoacetylation observed during mitosis, we observed a global concordance between the genomic organization of histone modifications in interphase and mitosis, suggesting that the epigenomic landscape may serve as a component of the mitotic bookmarking process. Next, we investigated the nucleosome that enters nucleosome depleted regions (NDRs) during mitosis. We observed that in ∼60% of the NDRs, the entering nucleosome is distinct from the surrounding highly acetylated nucleosomes and appears to have either low levels of acetylation or high levels of phosphorylation in adjacent residues (since adjacent phosphorylation may interfere with the ability to detect acetylation). Inhibition of histone deacetylases (HDACs) by the small molecule TSA reverts this pattern, suggesting that these nucleosomes are specifically deacetylated during mitosis. Altogether, by merging multiple approaches, our study provides evidence to support a model where histone modifications may play a role in mitotic bookmarking and uncovers new insights into the deposition of nucleosomes during mitosis.

Postnatal DNA demethylation and its role in tissue maturation.

Development in mammals is accompanied by specific de novo and demethylation events that are thought to stabilize differentiated cell phenotypes. We demonstrate that a large percentage of the tissue-specific methylation pattern is generated postnatally. Demethylation in the liver is observed in thousands of enhancer-like sequences associated with genes that undergo activation during the first few weeks of life. Using. conditional gene ablation strategy we show that the removal of these methyl groups is stable and necessary for assuring proper hepatocyte gene expression and function through its effect on chromatin accessibility. These postnatal changes in methylation come about through exposure to hormone signaling. These results define the molecular rules of 5-methyl-cytosine regulation as an epigenetic mechanism underlying cellular responses to. changing environment.

Biases in the SMART-DNA library preparation method associated with genomic poly dA/dT sequences.

Avoiding biases in next generation sequencing (NGS) library preparation is crucial for obtaining reliable sequencing data. Recently, a new library preparation method has been introduced which has eliminated the need for the ligation step. This method, termed SMART (switching mechanism at the 5' end of the RNA transcript), is based on template switching reverse transcription. To date, there has been no systematic analysis of the additional biases introduced by this method. We analysed the genomic distribution of sequenced reads prepared from genomic DNA using the SMART methodology and found a strong bias toward long (≥12bp) poly dA/dT containing genomic loci. This bias is unique to the SMART-based library preparation and does not appear when libraries are prepared with conventional ligation based methods. Although this bias is obvious only when performing paired end sequencing, it affects single end sequenced samples as well. Our analysis demonstrates that sequenced reads originating from SMART-DNA libraries are heavily skewed toward genomic poly dA/dT tracts. This bias needs to be considered when deciding to use SMART based technology for library preparation.

Systematic comparison of monoclonal versus polyclonal antibodies for mapping histone modifications by ChIP-seq.

BACKGROUND: The robustness of ChIP-seq datasets is highly dependent upon the antibodies used. Currently, polyclonal antibodies are the standard despite several limitations: They are non-renewable, vary in performance between lots and need to be validated with each new lot. In contrast, monoclonal antibody lots are renewable and provide consistent performance. To increase ChIP-seq standardization, we investigated whether monoclonal antibodies could replace polyclonal antibodies. We compared monoclonal antibodies that target five key histone modifications (H3K4me1, H3K4me3, H3K9me3, H3K27ac and H3K27me3) to their polyclonal counterparts in both human and mouse cells. RESULTS: Overall performance was highly similar for four monoclonal/polyclonal pairs, including when we used two distinct lots of the same monoclonal antibody. In contrast, the binding patterns for H3K27ac differed substantially between polyclonal and monoclonal antibodies. However, this was most likely due to the distinct immunogen used rather than the clonality of the antibody. CONCLUSIONS: Altogether, we found that monoclonal antibodies as a class perform equivalently to polyclonal antibodies for the detection of histone post-translational modifications in both human and mouse. Accordingly, we recommend the use of monoclonal antibodies in ChIP-seq experiments.

SIRT6 Suppresses Pancreatic Cancer through Control of Lin28b.

Chromatin remodeling proteins are frequently dysregulated in human cancer, yet little is known about how they control tumorigenesis. Here, we uncover an epigenetic program mediated by the NAD(+)-dependent histone deacetylase Sirtuin 6 (SIRT6) that is critical for suppression of pancreatic ductal adenocarcinoma (PDAC), one of the most lethal malignancies. SIRT6 inactivation accelerates PDAC progression and metastasis via upregulation of Lin28b, a negative regulator of the let-7 microRNA. SIRT6 loss results in histone hyperacetylation at the Lin28b promoter, Myc recruitment, and pronounced induction of Lin28b and downstream let-7 target genes, HMGA2, IGF2BP1, and IGF2BP3. This epigenetic program defines a distinct subset with a poor prognosis, representing 30%-40% of human PDAC, characterized by reduced SIRT6 expression and an exquisite dependence on Lin28b for tumor growth. Thus, we identify SIRT6 as an important PDAC tumor suppressor and uncover the Lin28b pathway as a potential therapeutic target in a molecularly defined PDAC subset. PAPERCLIP.

Oct1 and OCA-B are selectively required for CD4 memory T cell function.

Epigenetic changes are crucial for the generation of immunological memory. Failure to generate or maintain these changes will result in poor memory responses. Similarly, augmenting or stabilizing the correct epigenetic states offers a potential method of enhancing memory. Yet the transcription factors that regulate these processes are poorly defined. We find that the transcription factor Oct1 and its cofactor OCA-B are selectively required for the in vivo generation of CD4(+) memory T cells. More importantly, the memory cells that are formed do not respond properly to antigen reencounter. In vitro, both proteins are required to maintain a poised state at the Il2 target locus in resting but previously stimulated CD4(+) T cells. OCA-B is also required for the robust reexpression of multiple other genes including Ifng. ChIPseq identifies ∼50 differentially expressed direct Oct1 and OCA-B targets. We identify an underlying mechanism involving OCA-B recruitment of the histone lysine demethylase Jmjd1a to targets such as Il2, Ifng, and Zbtb32. The findings pinpoint Oct1 and OCA-B as central mediators of CD4(+) T cell memory.

Single-cell ChIP-seq reveals cell subpopulations defined by chromatin state.

Chromatin profiling provides a versatile means to investigate functional genomic elements and their regulation. However, current methods yield ensemble profiles that are insensitive to cell-to-cell variation. Here we combine microfluidics, DNA barcoding and sequencing to collect chromatin data at single-cell resolution. We demonstrate the utility of the technology by assaying thousands of individual cells and using the data to deconvolute a mixture of ES cells, fibroblasts and hematopoietic progenitors into high-quality chromatin state maps for each cell type. The data from each single cell are sparse, comprising on the order of 1,000 unique reads. However, by assaying thousands of ES cells, we identify a spectrum of subpopulations defined by differences in chromatin signatures of pluripotency and differentiation priming. We corroborate these findings by comparison to orthogonal single-cell gene expression data. Our method for single-cell analysis reveals aspects of epigenetic heterogeneity not captured by transcriptional analysis alone.

The histone deacetylase SIRT6 controls embryonic stem cell fate via TET-mediated production of 5-hydroxymethylcytosine.

How embryonic stem cells (ESCs) commit to specific cell lineages and yield all cell types of a fully formed organism remains a major question. ESC differentiation is accompanied by large-scale histone and DNA modifications, but the relations between these epigenetic categories are not understood. Here we demonstrate the interplay between the histone deacetylase sirtuin 6 (SIRT6) and the ten-eleven translocation enzymes (TETs). SIRT6 targets acetylated histone H3 at Lys 9 and 56 (H3K9ac and H3K56ac), while TETs convert 5-methylcytosine into 5-hydroxymethylcytosine (5hmC). ESCs derived from Sirt6 knockout (S6KO) mice are skewed towards neuroectoderm development. This phenotype involves derepression of OCT4, SOX2 and NANOG, which causes an upregulation of TET-dependent production of 5hmC. Genome-wide analysis revealed neural genes marked with 5hmC in S6KO ESCs, thereby implicating TET enzymes in the neuroectoderm-skewed differentiation phenotype. We demonstrate that SIRT6 functions as a chromatin regulator safeguarding the balance between pluripotency and differentiation through Tet-mediated production of 5hmC.

Pluripotency transcription factor Oct4 mediates stepwise nucleosome demethylation and depletion.

The mechanisms whereby the crucial pluripotency transcription factor Oct4 regulates target gene expression are incompletely understood. Using an assay system based on partially differentiated embryonic stem cells, we show that Oct4 opposes the accumulation of local H3K9me2 and subsequent Dnmt3a-mediated DNA methylation. Upon binding DNA, Oct4 recruits the histone lysine demethylase Jmjd1c. Chromatin immunoprecipitation (ChIP) time course experiments identify a stepwise Oct4 mechanism involving Jmjd1c recruitment and H3K9me2 demethylation, transient FACT (facilitates chromatin transactions) complex recruitment, and nucleosome depletion. Genome-wide and targeted ChIP confirms binding of newly synthesized Oct4, together with Jmjd1c and FACT, to the Pou5f1 enhancer and a small number of other Oct4 targets, including the Nanog promoter. Histone demethylation is required for both FACT recruitment and H3 depletion. Jmjd1c is required to induce endogenous Oct4 expression and fully reprogram fibroblasts to pluripotency, indicating that the assay system identifies functional Oct4 cofactors. These findings indicate that Oct4 sequentially recruits activities that catalyze histone demethylation and depletion.

Dissecting neural differentiation regulatory networks through epigenetic footprinting.

Models derived from human pluripotent stem cells that accurately recapitulate neural development in vitro and allow for the generation of specific neuronal subtypes are of major interest to the stem cell and biomedical community. Notch signalling, particularly through the Notch effector HES5, is a major pathway critical for the onset and maintenance of neural progenitor cells in the embryonic and adult nervous system. Here we report the transcriptional and epigenomic analysis of six consecutive neural progenitor cell stages derived from a HES5::eGFP reporter human embryonic stem cell line. Using this system, we aimed to model cell-fate decisions including specification, expansion and patterning during the ontogeny of cortical neural stem and progenitor cells. In order to dissect regulatory mechanisms that orchestrate the stage-specific differentiation process, we developed a computational framework to infer key regulators of each cell-state transition based on the progressive remodelling of the epigenetic landscape and then validated these through a pooled short hairpin RNA screen. We were also able to refine our previous observations on epigenetic priming at transcription factor binding sites and suggest here that they are mediated by combinations of core and stage-specific factors. Taken together, we demonstrate the utility of our system and outline a general framework, not limited to the context of the neural lineage, to dissect regulatory circuits of differentiation.