Borbala Mifsud

Mapping long-range promoter contacts in human cells with high-resolution capture Hi-C

Transcriptional control in large genomes often requires looping interactions between distal DNA elements, such as enhancers and target promoters. Current chromosome conformation capture techniques do not offer sufficiently high resolution to interrogate these regulatory interactions on a genomic scale. Here we use Capture Hi-C (CHi-C), an adapted genome conformation assay, to examine the long-range interactions of almost 22,000 promoters in 2 human blood cell types. We identify over 1.6 million shared and cell type–restricted interactions spanning hundreds of kilobases between promoters and distal loci. Transcriptionally active genes contact enhancer-like elements, whereas transcriptionally inactive genes interact with previously uncharacterized elements marked by repressive features that may act as long-range silencers. Finally, we show that interacting loci are enriched for disease-associated SNPs, suggesting how distal mutations may disrupt the regulation of relevant genes. This study provides new insights and accessible tools to dissect the regulatory interactions that underlie normal and aberrant gene regulation.

The pluripotent regulatory circuitry connecting promoters to their long-range interacting elements

The mammalian genome harbors up to one million regulatory elements often located at great distances from their target genes. Long-range elements control genes through physical contact with promoters and can be recognized by the presence of specific histone modifications and transcription factor binding. Linking regulatory elements to specific promoters genome-wide is currently impeded by the limited resolution of high-throughput chromatin interaction assays. Here we apply a sequence capture approach to enrich Hi-C libraries for >22,000 annotated mouse promoters to identify statistically significant, long-range interactions at restriction fragment resolution, assigning long-range interacting elements to their target genes genome-wide in embryonic stem cells and fetal liver cells. The distal sites contacting active genes are enriched in active histone modifications and transcription factor occupancy, whereas inactive genes contact distal sites with repressive histone marks, demonstrating the regulatory potential of the distal elements identified. Furthermore, we find that coregulated genes cluster nonrandomly in spatial interaction networks correlated with their biological function and expression level. Interestingly, we find the strongest gene clustering in ES cells between transcription factor genes that control key developmental processes in embryogenesis. The results provide the first genome-wide catalog linking gene promoters to their long-range interacting elements and highlight the complex spatial regulatory circuitry controlling mammalian gene expression.

Polycomb repressive complex PRC1 spatially constrains the mouse embryonic stem cell genome.

The Polycomb repressive complexes PRC1 and PRC2 maintain embryonic stem cell (ESC) pluripotency by silencing lineage-specifying developmental regulator genes. Emerging evidence suggests that Polycomb complexes act through controlling spatial genome organization. We show that PRC1 functions as a master regulator of mouse ESC genome architecture by organizing genes in three-dimensional interaction networks. The strongest spatial network is composed of the four Hox gene clusters and early developmental transcription factor genes, the majority of which contact poised enhancers. Removal of Polycomb repression leads to disruption of promoter-promoter contacts in the Hox gene network. In contrast, promoter-enhancer contacts are maintained in the absence of Polycomb repression, with accompanying widespread acquisition of active chromatin signatures at network enhancers and pronounced transcriptional upregulation of network genes. Thus, PRC1 physically constrains developmental transcription factor genes and their enhancers in a silenced but poised spatial network. We propose that the selective release of genes from this spatial network underlies cell fate specification during early embryonic development.

Genome-wide association analysis identifies novel blood pressure loci and offers biological insights into cardiovascular risk.

Elevated blood pressure is the leading heritable risk factor for cardiovascular disease worldwide. We report genetic association of blood pressure (systolic, diastolic, pulse pressure) among UK Biobank participants of European ancestry with independent replication in other cohorts, and robust validation of 107 independent loci. We also identify new independent variants at 11 previously reported blood pressure loci. In combination with results from a range of in silico functional analyses and wet bench experiments, our findings highlight new biological pathways for blood pressure regulation enriched for genes expressed in vascular tissues and identify potential therapeutic targets for hypertension. Results from genetic risk score models raise the possibility of a precision medicine approach through early lifestyle intervention to offset the impact of blood pressure–raising genetic variants on future cardiovascular disease risk.

Flipping between Polycomb repressed and active transcriptional states introduces noise in gene expression

Polycomb repressive complexes (PRCs) are important histone modifiers, which silence gene expression; yet, there exists a subset of PRC-bound genes actively transcribed by RNA polymerase II (RNAPII). It is likely that the role of Polycomb repressive complex is to dampen expression of these PRC-active genes. However, it is unclear how this flipping between chromatin states alters the kinetics of transcription. Here, we integrate histone modifications and RNAPII states derived from bulk ChIP-seq data with single-cell RNA-sequencing data. We find that Polycomb repressive complex-active genes have greater cell-to-cell variation in expression than active genes, and these results are validated by knockout experiments. We also show that PRC-active genes are clustered on chromosomes in both two and three dimensions, and interactions with active enhancers promote a stabilization of gene expression noise. These findings provide new insights into how chromatin regulation modulates stochastic gene expression and transcriptional bursting, with implications for regulation of pluripotency and development.Polycomb repressive complexes modify histones but it is unclear how changes in chromatin states alter kinetics of transcription. Here, the authors use single-cell RNAseq and ChIPseq to find that actively transcribed genes with Polycomb marks have greater cell-to-cell variation in expression.

Thirty loci identified for heart rate response to exercise and recovery implicate autonomic nervous system

Impaired capacity to increase heart rate (HR) during exercise (ΔHRex), and a reduced rate of recovery post-exercise (ΔHRrec) are associated with higher cardiovascular mortality rates. Currently, the genetic basis of both phenotypes remains to be elucidated. We conduct genome-wide association studies (GWASs) for ΔHRex and ΔHRrec in ~40,000 individuals, followed by replication in ~27,000 independent samples, all from UK Biobank. Six and seven single-nucleotide polymorphisms for ΔHRex and ΔHRrec, respectively, formally replicate. In a full data set GWAS, eight further loci for ΔHRex and nine for ΔHRrec are genome-wide significant (P ≤ 5 × 10−8). In total, 30 loci are discovered, 8 being common across traits. Processes of neural development and modulation of adrenergic activity by the autonomic nervous system are enriched in these results. Our findings reinforce current understanding of HR response to exercise and recovery and could guide future studies evaluating its contribution to cardiovascular risk prediction.Genome-wide association studies have identified multiple loci for resting heart rate (HR) but the genetic factors associated with HR increase during and HR recovery after exercise are less well studied. Here, the authors examine both traits in a two-stage GWAS design in up to 67,257 individuals from UK Biobank.

The biological impact of blood pressure associated genetic variants in the natriuretic peptide receptor C gene on human vascular smooth muscle.

Abstract Elevated blood pressure (BP) is a major global risk factor for cardiovascular disease. Genome-wide association studies have identified several genetic variants at the NPR3 locus associated with BP, but the functional impact of these variants remains to be determined. Here we confirmed, by a genome-wide association study within UK Biobank, the existence of two independent BP-related signals within NPR3 locus. Using human primary vascular smooth muscle cells (VSMCs) and endothelial cells (ECs) from different individuals, we found that the BP-elevating alleles within one linkage disequilibrium block identified by the sentinel variant rs1173771 was associated with lower endogenous NPR3 mRNA and protein levels in VSMCs, together with reduced levels in open chromatin and nuclear protein binding. The BP-elevating alleles also increased VSMC proliferation, angiotensin II-induced calcium flux and cell contraction. However, an analogous genotype-dependent association was not observed in vascular ECs. Our study identifies novel, putative mechanisms for BP-associated variants at the NPR3 locus to elevate BP, further strengthening the case for targeting NPR-C as a therapeutic approach for hypertension and cardiovascular disease prevention.

Capturing genomic relationships that matter

There is a strong interrelationship within the cell nucleus between form and function of the genome. This connection is exhibited across multiple hierarchies, ranging from grand-scale positioning of chromosomes and their intersection with specific nuclear functional activities, the segregation of chromosome structure into distinct domains and long-range regulatory contacts that drive spatial and temporal expression patterns of genes. Fifteen years ago, the development of the chromosome conformation capture method placed the nature of specific, long-range regulatory interactions under scrutiny. However, its development and integration with next-generation sequencing technologies has greatly expanded the breadth and scope of what is detected. The sheer scale of data offered by these important advances has come with new and challenging bottlenecks that are both experimental and bioinformatical. Here, we discuss the recent and prospective development and implementation of new methodologies and analytical tools that are allowing an in-depth, yet focussed characterisation of genomic contacts that are associated with functional activities in the nucleus.

Genome-wide association study of response to methotrexate in early rheumatoid arthritis patients

Methotrexate (MTX) monotherapy is a common first treatment for rheumatoid arthritis (RA), but many patients do not respond adequately. In order to identify genetic predictors of response, we have combined data from two consortia to carry out a genome-wide study of response to MTX in 1424 early RA patients of European ancestry. Clinical endpoints were change from baseline to 6 months after starting treatment in swollen 28-joint count, tender 28-joint count, C-reactive protein and the overall 3-component disease activity score (DAS28). No single nucleotide polymorphism (SNP) reached genome-wide statistical significance for any outcome measure. The strongest evidence for association was with rs168201 in NRG3 (p = 10−7 for change in DAS28). Some support was also seen for association with ZMIZ1, previously highlighted in a study of response to MTX in juvenile idiopathic arthritis. Follow-up in two smaller cohorts of 429 and 177 RA patients did not support these findings, although these cohorts were more heterogeneous.

Divergent wiring of repressive and active chromatin interactions between mouse embryonic and trophoblast lineages

The establishment of the embryonic and trophoblast lineages is a developmental decision underpinned by dramatic differences in the epigenetic landscape of the two compartments. However, it remains unknown how epigenetic information and transcription factor networks map to the 3D arrangement of the genome, which in turn may mediate transcriptional divergence between the two cell lineages. Here, we perform promoter capture Hi-C experiments in mouse trophoblast (TSC) and embryonic (ESC) stem cells to understand how chromatin conformation relates to cell-specific transcriptional programmes. We find that key TSC genes that are kept repressed in ESCs exhibit interactions between H3K27me3-marked regions in ESCs that depend on Polycomb repressive complex 1. Interactions that are prominent in TSCs are enriched for enhancer–gene contacts involving key TSC transcription factors, as well as TET1, which helps to maintain the expression of TSC-relevant genes. Our work shows that the first developmental cell fate decision results in distinct chromatin conformation patterns establishing lineage-specific contexts involving both repressive and active interactions.The role of the genome structure in the establishment of the embryonic and trophoblast lineages is still not well understood. Here the authors perform promoter capture Hi-C in mouse trophoblast and embryonic stem cells and find divergent networks of repressive and active chromatin interactions between the two lineages.