Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance.

BACKGROUND: Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states. METHODS: We exposed human primary monocytes from healthy donors (n = 6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2 h) and two models of endotoxin tolerance: repeated stimulations (6 + 6 h) and prolonged exposure to endotoxin (24 h). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing. RESULTS: We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6,884 of 34,616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease. CONCLUSION: This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

Variation at 10p12.2 and 10p14 influences risk of childhood B-cell acute lymphoblastic leukemia and phenotype.

Acute lymphoblastic leukemia (ALL) is the major pediatric cancer diagnosed in economically developed countries with B-cell precursor (BCP)-ALL, accounting for approximately 70% of ALL. Recent genome-wide association studies (GWAS) have provided the first unambiguous evidence for common inherited susceptibility to BCP-ALL, identifying susceptibility loci at 7p12.2, 9p21.3, 10q21.2, and 14q11.2. To identify additional BCP-ALL susceptibility loci, we conducted a GWAS and performed a meta-analysis with a published GWAS totaling 1658 cases and 4723 controls, with validation in 1449 cases and 1488 controls. Combined analysis identified novel loci mapping to 10p12.2 (rs10828317, odds ratio [OR] = 1.23; P = 2.30 × 10(-9)) and 10p14 marked by rs3824662 (OR = 1.31; P = 8.62 × 10(-12)). The single nucleotide polymorphism rs10828317 is responsible for the N215S polymorphism in exon 7 of PIP4K2A, and rs3824662 localizes to intron 3 of the transcription factor and putative tumor suppressor gene GATA3. The rs10828317 association was shown to be specifically associated with hyperdiploid ALL, whereas the rs3824662-associated risk was confined to nonhyperdiploid non-TEL-AML1 + ALL. The risk allele of rs3824662 was correlated with older age at diagnosis (P < .001) and significantly worse event-free survivorship (P < .0001). These findings provide further insights into the genetic and biological basis of inherited genetic susceptibility to BCP-ALL and the influence of constitutional genotype on disease development.

Comprehensive epigenomic profiling reveals the extent of disease-specific chromatin states and informs target discovery in ankylosing spondylitis.

Ankylosing spondylitis (AS) is a common, highly heritable inflammatory arthritis characterized by enthesitis of the spine and sacroiliac joints. Genome-wide association studies (GWASs) have revealed more than 100 genetic associations whose functional effects remain largely unresolved. Here, we present a comprehensive transcriptomic and epigenomic map of disease-relevant blood immune cell subsets from AS patients and healthy controls. We find that, while CD14+ monocytes and CD4+ and CD8+ T cells show disease-specific differences at the RNA level, epigenomic differences are only apparent upon multi-omics integration. The latter reveals enrichment at disease-associated loci in monocytes. We link putative functional SNPs to genes using high-resolution Capture-C at 10 loci, including PTGER4 and ETS1, and show how disease-specific functional genomic data can be integrated with GWASs to enhance therapeutic target discovery. This study combines epigenetic and transcriptional analysis with GWASs to identify disease-relevant cell types and gene regulation of likely pathogenic relevance and prioritize drug targets.

Epigenomic analysis reveals a dynamic and context-specific macrophage enhancer landscape associated with innate immune activation and tolerance.

BACKGROUND: Chromatin states and enhancers associate gene expression, cell identity and disease. Here, we systematically delineate the acute innate immune response to endotoxin in terms of human macrophage enhancer activity and contrast with endotoxin tolerance, profiling the coding and non-coding transcriptome, chromatin accessibility and epigenetic modifications. RESULTS: We describe the spectrum of enhancers under acute and tolerance conditions and the regulatory networks between these enhancers and biological processes including gene expression, splicing regulation, transcription factor binding and enhancer RNA signatures. We demonstrate that the vast majority of differentially regulated enhancers on acute stimulation are subject to tolerance and that expression quantitative trait loci, disease-risk variants and eRNAs are enriched in these regulatory regions and related to context-specific gene expression. We find enrichment for context-specific eQTL involving endotoxin response and specific infections and delineate specific differential regions informative for GWAS variants in inflammatory bowel disease and multiple sclerosis, together with a context-specific enhancer involving a bacterial infection eQTL for KLF4. We show enrichment in differential enhancers for tolerance involving transcription factors NFκB-p65, STATs and IRFs and prioritize putative causal genes directly linking genetic variants and disease risk enhancers. We further delineate similarities and differences in epigenetic landscape between stem cell-derived macrophages and primary cells and characterize the context-specific enhancer activities for key innate immune response genes KLF4, SLAMF1 and IL2RA. CONCLUSIONS: Our study demonstrates the importance of context-specific macrophage enhancers in gene regulation and utility for interpreting disease associations, providing a roadmap to link genetic variants with molecular and cellular functions.

Author Correction: Genetic and regulatory mechanism of susceptibility to high-hyperdiploid acute lymphoblastic leukaemia at 10q21.2.

This corrects the article DOI: 10.1038/ncomms14616.

Genetic and regulatory mechanism of susceptibility to high-hyperdiploid acute lymphoblastic leukaemia at 10p21.2.

Despite high-hyperdiploid acute lymphoblastic leukaemia (HD-ALL) being the most common subgroup of paediatric ALL, its aetiology remains unknown. Genome-wide association studies have demonstrated association at 10q21.2. Here, we sought to determine how this region influences HD-ALL risk. We impute genotypes across the locus, finding the single nucleotide polymorphism rs7090445 highly associated with HD-ALL (P=1.54 × 10-38), and residing in a predicted enhancer element. We show this region physically interacts with the transcription start site of ARID5B, that alleles of rs7090445 have differential enhancer activity and influence RUNX3 binding. RUNX3 knock-down reduces ARID5B expression and rs7090445 enhancer activity. Individuals carrying the rs7090445-C risk allele also have reduced ARID5B expression. Finally, the rs7090445-C risk allele is preferentially retained in HD-ALL blasts consistent with inherited genetic variation contributing to arrest of normal lymphocyte development, facilitating leukaemic clonal expansion. These data provide evidence for a biological mechanism underlying hereditary risk of HD-ALL at 10q21.2.

Identification of 19 new risk loci and potential regulatory mechanisms influencing susceptibility to testicular germ cell tumor.

Genome-wide association studies (GWAS) have transformed understanding of susceptibility to testicular germ cell tumors (TGCTs), but much of the heritability remains unexplained. Here we report a new GWAS, a meta-analysis with previous GWAS and a replication series, totaling 7,319 TGCT cases and 23,082 controls. We identify 19 new TGCT risk loci, roughly doubling the number of known TGCT risk loci to 44. By performing in situ Hi-C in TGCT cells, we provide evidence for a network of physical interactions among all 44 TGCT risk SNPs and candidate causal genes. Our findings implicate widespread disruption of developmental transcriptional regulators as a basis of TGCT susceptibility, consistent with failed primordial germ cell differentiation as an initiating step in oncogenesis. Defective microtubule assembly and dysregulation of KIT-MAPK signaling also feature as recurrently disrupted pathways. Our findings support a polygenic model of risk and provide insight into the biological basis of TGCT.

Finite-connectivity spin-glass phase diagrams and low-density parity check codes.

We obtain phase diagrams of regular and irregular finite-connectivity spin glasses. Contact is first established between properties of the phase diagram and the performance of low-density parity check (LDPC) codes within the replica symmetric (RS) ansatz. We then study the location of the dynamical and critical transition points of these systems within the one step replica symmetry breaking theory (RSB), extending similar calculations that have been performed in the past for the Bethe spin-glass problem. We observe that the location of the dynamical transition line does change within the RSB theory, in comparison with the results obtained in the RS case. For LDPC decoding of messages transmitted over the binary erasure channel we find, at zero temperature and rate , an RS critical transition point at while the critical RSB transition point is located at , to be compared with the corresponding Shannon bound . For the binary symmetric channel we show that the low temperature reentrant behavior of the dynamical transition line, observed within the RS ansatz, changes its location when the RSB ansatz is employed; the dynamical transition point occurs at higher values of the channel noise. Possible practical implications to improve the performance of the state-of-the-art error correcting codes are discussed.

Cytochrome P450 Allele CYP3A7*1C Associates with Adverse Outcomes in Chronic Lymphocytic Leukemia, Breast, and Lung Cancer.

CYP3A enzymes metabolize endogenous hormones and chemotherapeutic agents used to treat cancer, thereby potentially affecting drug effectiveness. Here, we refined the genetic basis underlying the functional effects of a CYP3A haplotype on urinary estrone glucuronide (E1G) levels and tested for an association between CYP3A genotype and outcome in patients with chronic lymphocytic leukemia (CLL), breast, or lung cancers. The most significantly associated SNP was rs45446698, an SNP that tags the CYP3A7*1C allele; this SNP was associated with a 54% decrease in urinary E1G levels. Genotyping this SNP in 1,008 breast cancer, 1,128 lung cancer, and 347 CLL patients, we found that rs45446698 was associated with breast cancer mortality (HR, 1.74; P = 0.03), all-cause mortality in lung cancer patients (HR, 1.43; P = 0.009), and CLL progression (HR, 1.62; P = 0.03). We also found borderline evidence of a statistical interaction between the CYP3A7*1C allele, treatment of patients with a cytotoxic agent that is a CYP3A substrate, and clinical outcome (Pinteraction = 0.06). The CYP3A7*1C allele, which results in adult expression of the fetal CYP3A7 gene, is likely to be the functional allele influencing levels of circulating endogenous sex hormones and outcome in these various malignancies. Further studies confirming these associations and determining the mechanism by which CYP3A7*1C influences outcome are required. One possibility is that standard chemotherapy regimens that include CYP3A substrates may not be optimal for the approximately 8% of cancer patients who are CYP3A7*1C carriers.

Genetic Predisposition to Chronic Lymphocytic Leukemia Is Mediated by a BMF Super-Enhancer Polymorphism.

Chronic lymphocytic leukemia (CLL) is an adult B cell malignancy. Genome-wide association studies show that variation at 15q15.1 influences CLL risk. We deciphered the causal variant at 15q15.1 and the mechanism by which it influences tumorigenesis. We imputed all possible genotypes across the locus and then mapped highly associated SNPs to areas of chromatin accessibility, evolutionary conservation, and transcription factor binding. SNP rs539846 C>A, the most highly associated variant (p = 1.42 × 10(-13), odds ratio = 1.35), localizes to a super-enhancer defined by extensive histone H3 lysine 27 acetylation in intron 3 of B cell lymphoma 2 (BCL2)-modifying factor (BMF). The rs539846-A risk allele alters a conserved RELA-binding motif, disrupts RELA binding, and is associated with decreased BMF expression in CLL. These findings are consistent with rs539846 influencing CLL susceptibility through differential RELA binding, with direct modulation of BMF expression impacting on anti-apoptotic BCL2, a hallmark of oncogenic dependency in CLL.

Genome-wide association study identifies multiple susceptibility loci for multiple myeloma.

Multiple myeloma (MM) is a plasma cell malignancy with a significant heritable basis. Genome-wide association studies have transformed our understanding of MM predisposition, but individual studies have had limited power to discover risk loci. Here we perform a meta-analysis of these GWAS, add a new GWAS and perform replication analyses resulting in 9,866 cases and 239,188 controls. We confirm all nine known risk loci and discover eight new loci at 6p22.3 (rs34229995, P=1.31 × 10(-8)), 6q21 (rs9372120, P=9.09 × 10(-15)), 7q36.1 (rs7781265, P=9.71 × 10(-9)), 8q24.21 (rs1948915, P=4.20 × 10(-11)), 9p21.3 (rs2811710, P=1.72 × 10(-13)), 10p12.1 (rs2790457, P=1.77 × 10(-8)), 16q23.1 (rs7193541, P=5.00 × 10(-12)) and 20q13.13 (rs6066835, P=1.36 × 10(-13)), which localize in or near to JARID2, ATG5, SMARCD3, CCAT1, CDKN2A, WAC, RFWD3 and PREX1. These findings provide additional support for a polygenic model of MM and insight into the biological basis of tumour development.

Capture Hi-C identifies the chromatin interactome of colorectal cancer risk loci.

Multiple regulatory elements distant from their targets on the linear genome can influence the expression of a single gene through chromatin looping. Chromosome conformation capture implemented in Hi-C allows for genome-wide agnostic characterization of chromatin contacts. However, detection of functional enhancer-promoter interactions is precluded by its effective resolution that is determined by both restriction fragmentation and sensitivity of the experiment. Here we develop a capture Hi-C (cHi-C) approach to allow an agnostic characterization of these physical interactions on a genome-wide scale. Single-nucleotide polymorphisms associated with complex diseases often reside within regulatory elements and exert effects through long-range regulation of gene expression. Applying this cHi-C approach to 14 colorectal cancer risk loci allows us to identify key long-range chromatin interactions in cis and trans involving these loci.

visPIG--a web tool for producing multi-region, multi-track, multi-scale plots of genetic data.

We present VISual Plotting Interface for Genetics (visPIG; http://vispig.icr.ac.uk), a web application to produce multi-track, multi-scale, multi-region plots of genetic data. visPIG has been designed to allow users not well versed with mathematical software packages and/or programming languages such as R, Matlab®, Python, etc., to integrate data from multiple sources for interpretation and to easily create publication-ready figures. While web tools such as the UCSC Genome Browser or the WashU Epigenome Browser allow custom data uploads, such tools are primarily designed for data exploration. This is also true for the desktop-run Integrative Genomics Viewer (IGV). Other locally run data visualisation software such as Circos require significant computer skills of the user. The visPIG web application is a menu-based interface that allows users to upload custom data tracks and set track-specific parameters. Figures can be downloaded as PDF or PNG files. For sensitive data, the underlying R code can also be downloaded and run locally. visPIG is multi-track: it can display many different data types (e.g association, functional annotation, intensity, interaction, heat map data,…). It also allows annotation of genes and other custom features in the plotted region(s). Data tracks can be plotted individually or on a single figure. visPIG is multi-region: it supports plotting multiple regions, be they kilo- or megabases apart or even on different chromosomes. Finally, visPIG is multi-scale: a sub-region of particular interest can be 'zoomed' in. We describe the various features of visPIG and illustrate its utility with examples. visPIG is freely available through http://vispig.icr.ac.uk under a GNU General Public License (GPLv3).

Common variation at 3q26.2, 6p21.33, 17p11.2 and 22q13.1 influences multiple myeloma risk.

To identify variants for multiple myeloma risk, we conducted a genome-wide association study with validation in additional series totaling 4,692 individuals with multiple myeloma (cases) and 10,990 controls. We identified four risk loci at 3q26.2 (rs10936599, P = 8.70 × 10(-14)), 6p21.33 (rs2285803, PSORS1C2, P = 9.67 × 10(-11)), 17p11.2 (rs4273077, TNFRSF13B, P = 7.67 × 10(-9)) and 22q13.1 (rs877529, CBX7, P = 7.63 × 10(-16)). These data provide further evidence for genetic susceptibility to this B-cell hematological malignancy, as well as insight into the biological basis of predisposition.