Enhancer elements in the human genome control how genes are expressed in specific cell types and harbor thousands of genetic variants that influence risk for common diseases1-4. Yet, we still do not know how enhancers regulate specific genes, and we lack general rules to predict enhancer-gene connections across cell types5,6. We developed an experimental approach, CRISPRi-FlowFISH, to perturb enhancers in the genome, and we applied it to test >3,500 potential enhancer-gene connections for 30 genes. We found that a simple activity-by-contact model substantially outperformed previous methods at predicting the complex connections in our CRISPR dataset. This activity-by-contact model allows us to construct genome-wide maps of enhancer-gene connections in a given cell type, on the basis of chromatin state measurements. Together, CRISPRi-FlowFISH and the activity-by-contact model provide a systematic approach to map and predict which enhancers regulate which genes, and will help to interpret the functions of the thousands of disease risk variants in the noncoding genome.

/pdf/activity-by-contact-model-of-enhancer-promoter-regulation-1yqm2q76rp.pdf

Activity-by-contact model of enhancer-promoter regulation from thousands of CRISPR perturbations.

The recent development of experimental methods for measuring chromatin state at single-cell resolution has created a need for computational tools capable of analyzing these datasets. Here we developed Signac, a comprehensive toolkit for the analysis of single-cell chromatin data. Signac enables an end-to-end analysis of single-cell chromatin data, including peak calling, quantification, quality control, dimension reduction, clustering, integration with single-cell gene expression datasets, DNA motif analysis and interactive visualization. Through its seamless compatibility with the Seurat package, Signac facilitates the analysis of diverse multimodal single-cell chromatin data, including datasets that co-assay DNA accessibility with gene expression, protein abundance and mitochondrial genotype. We demonstrate scaling of the Signac framework to analyze datasets containing over 700,000 cells.

Single-cell chromatin state analysis with Signac.

Replication of the genome before mitotic cell division is a highly regulated process that ensures the fidelity of DNA duplication. DNA replication initiates at specific locations, termed origins of replication, and progresses in a defined temporal order during the S phase of the cell cycle. The relationship between replication timing and gene expression has been the subject of some speculation. A recent genome-wide analysis in Saccharomyces cerevisiae showed no association between replication timing and gene expression. In higher eukaryotes, the limited number of genomic loci analyzed has not permitted a firm conclusion regarding this association. To explore the relationship between DNA replication and gene expression in higher eukaryotes, we developed a strategy to measure the timing of DNA replication for thousands of genes in a single DNA array hybridization experiment. Using this approach, we generated a genome-wide map of replication timing for Drosophila melanogaster. Moreover, by surveying over 40% of all D. melanogaster genes, we found a strong correlation between DNA replication early in S phase and transcriptional activity. As this correlation does not exist in S. cerevisiae, this interplay between DNA replication and transcription may be a unique characteristic of higher eukaryotes.

Genome-wide DNA replication profile for Drosophila melanogaster: a link between transcription and replication timing: a link between transcription and replication timing

Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown1. The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer2-4 and coronary heart disease5-this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP)6. Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

/pdf/inherited-causes-of-clonal-haematopoiesis-in-97-691-whole-2ew4mk46pj.pdf

Inherited causes of clonal haematopoiesis in 97,691 whole genomes.

Genome-wide association studies (GWAS) have identified thousands of noncoding loci that are associated with human diseases and complex traits, each of which could reveal insights into the mechanisms of disease1. Many of the underlying causal variants may affect enhancers2,3, but we lack accurate maps of enhancers and their target genes to interpret such variants. We recently developed the activity-by-contact (ABC) model to predict which enhancers regulate which genes and validated the model using CRISPR perturbations in several cell types4. Here we apply this ABC model to create enhancer–gene maps in 131 human cell types and tissues, and use these maps to interpret the functions of GWAS variants. Across 72 diseases and complex traits, ABC links 5,036 GWAS signals to 2,249 unique genes, including a class of 577 genes that appear to influence multiple phenotypes through variants in enhancers that act in different cell types. In inflammatory bowel disease (IBD), causal variants are enriched in predicted enhancers by more than 20-fold in particular cell types such as dendritic cells, and ABC achieves higher precision than other regulatory methods at connecting noncoding variants to target genes. These variant-to-function maps reveal an enhancer that contains an IBD risk variant and that regulates the expression of PPIF to alter the membrane potential of mitochondria in macrophages. Our study reveals principles of genome regulation, identifies genes that affect IBD and provides a resource and generalizable strategy to connect risk variants of common diseases to their molecular and cellular functions. Mapping enhancer regulation across human cell types and tissues illuminates genome function and provides a resource to connect risk variants for common diseases to their molecular and cellular functions.

Genome-wide enhancer maps link risk variants to disease genes

Functional genomics assays based on high-throughput sequencing greatly expand our ability to understand the genome. Here, we define the ENCODE blacklist- a comprehensive set of regions in the human, mouse, worm, and fly genomes that have anomalous, unstructured, or high signal in next-generation sequencing experiments independent of cell line or experiment. The removal of the ENCODE blacklist is an essential quality measure when analyzing functional genomics data.

/pdf/the-encode-blacklist-identification-of-problematic-regions-4zhjjl4wv4.pdf

The ENCODE Blacklist: Identification of Problematic Regions of the Genome.

The Saccharomyces cerevisiae ribosomal DNA (rDNA) locus is known to exhibit greater instability relative to the rest of the genome. However, wild-type cells preferentially maintain a stable number of rDNA copies, suggesting underlying genetic control of the size of this locus. We performed a screen of a subset of the Yeast Knock-Out (YKO) single gene deletion collection to identify genetic regulators of this locus and to determine if rDNA copy number correlates with yeast replicative lifespan. While we found no correlation between replicative lifespan and rDNA size, we identified 64 candidate strains with significant rDNA copy number differences. However, in the process of validating candidate rDNA variants, we observed that independent isolates of our de novo gene deletion strains had unsolicited but significant changes in rDNA copy number. Moreover, we were not able to recapitulate rDNA phenotypes from the YKO yeast deletion collection. Instead, we found that the standard lithium acetate transformation protocol is a significant source of rDNA copy number variation, with lithium acetate exposure being the treatment causing variable rDNA copy number events after transformation. As the effects of variable rDNA copy number are being increasingly reported, our finding that rDNA is affected by lithium acetate exposure suggested that rDNA copy number variants may be influential passenger mutations in standard strain construction in S. cerevisiae.

rDNA Copy Number Variants Are Frequent Passenger Mutations in Saccharomyces cerevisiae Deletion Collections and de Novo Transformants

A form of dwarfism known as Meier-Gorlin syndrome (MGS) is caused by recessive mutations in one of six different genes (ORC1, ORC4, ORC6, CDC6, CDT1, and MCM5). These genes encode components of the pre-replication complex, which assembles at origins of replication prior to S phase. Also, variants in two additional replication initiation genes have joined the list of causative mutations for MGS (Geminin and CDC45). The identity of the causative MGS genetic variants strongly suggests that some aspect of replication is amiss in MGS patients; however, little evidence has been obtained regarding what aspect of chromosome replication is faulty. Since the site of one of the missense mutations in the human ORC4 alleles is conserved between humans and yeast, we sought to determine in what way this single amino acid change affects the process of chromosome replication, by introducing the comparable mutation into yeast (orc4Y232C). We find that yeast cells with the orc4Y232C allele have a prolonged S-phase, due to compromised replication initiation at the ribosomal DNA (rDNA) locus located on chromosome XII. The inability to initiate replication at the rDNA locus results in chromosome breakage and a severely reduced rDNA copy number in the survivors, presumably helping to ensure complete replication of chromosome XII. Although reducing rDNA copy number may help ensure complete chromosome replication, orc4Y232C cells struggle to meet the high demand for ribosomal RNA synthesis. This finding provides additional evidence linking two essential cellular pathways-DNA replication and ribosome biogenesis.

Defective replication initiation results in locus specific chromosome breakage and a ribosomal RNA deficiency in yeast

Genome architecture has proven to be critical in determining gene regulation across almost all domains of life. While many of the key components and mechanisms of eukaryotic genome organization have been described, the interplay between bacterial DNA organization and gene regulation is only now being fully appreciated. An increasing pool of evidence has demonstrated that the bacterial chromosome can reasonably be thought of as chromatin, and that bacterial chromosomes contain transcriptionally silent and transcriptionally active regions analogous to heterochromatin and euchromatin, respectively. The roles played by histones in eukaryotic systems appear to be shared across a range of nucleoid-associated proteins (NAPs) in bacteria, which function to compact, structure, and regulate large portions of bacterial chromosomes. The broad range of extant NAPs, and the extent to which they differ from species to species, has raised additional challenges in identifying and characterizing their roles in all but a handful of model bacteria. Here we review the regulatory roles played by NAPs in several well-studied bacteria and use the resulting state of knowledge to provide a working definition for NAPs, based on their function, binding pattern, and expression levels. We present a screening procedure which can be applied to any species for which transcriptomic data are available. Finally, we note that NAPs tend to play two major regulatory roles - xenogeneic silencers and developmental regulators - and that many unrecognized potential NAPs exist in each bacterial species examined.

Nucleoid-associated proteins shape chromatin structure and transcriptional regulation across the bacterial kingdom.

Free living bacteria adapt to changes in the environment by reprogramming gene expression through precise interactions of hundreds of DNA-binding proteins. Technologies such as ChIP-seq enable targeted characterization of regulatory interactions for individual DNA-binding proteins. However, in order to understand the cell’s global regulatory logic, we need to simultaneously monitor all such interactions in response to diverse genetic and environmental perturbations. To address this challenge, we have developed high-resolution in vivo protein occupancy display (IPOD-HR), a technology that enables rapid, quantitative, and comprehensive monitoring of DNA-protein interactions across a bacterial chromosome. IPOD-HR enables simultaneous activity profiling of all known sequence specific transcription factors, discovery of novel condition-dependent DNA-binding proteins, and systematic inference of binding specificity models for all bound transcription factors. IPOD-HR also reveals many large domains of extended protein occupancy in Escherichia coli that define relatively stable, transcriptionally silent regions with unique sequence and gene functional features.

https://www.biorxiv.org/content/biorxiv/early/2020/08/10/2020.01.29.924811.full.pdf

Haley M. Amemiya

Papers

The ENCODE Blacklist: Identification of Problematic Regions of the Genome.

rDNA Copy Number Variants Are Frequent Passenger Mutations in Saccharomyces cerevisiae Deletion Collections and de Novo Transformants

Defective replication initiation results in locus specific chromosome breakage and a ribosomal RNA deficiency in yeast

Nucleoid-associated proteins shape chromatin structure and transcriptional regulation across the bacterial kingdom.

Dynamic landscape of protein occupancy across the Escherichia coli chromosome