RESUMO
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.
Assuntos
Epigênese Genética , Monócitos , Humanos , Monócitos/metabolismo , Tolerância à Endotoxina , Epigenômica , Cromatina/genética , Imunidade Inata/genética , Transcriptoma , Endotoxinas/toxicidade , Proteínas de Membrana/genéticaRESUMO
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.
RESUMO
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.
Assuntos
Elementos Facilitadores Genéticos , Epigenômica , Cromatina , Endotoxinas , Humanos , Imunidade Inata/genética , MacrófagosRESUMO
This corrects the article DOI: 10.1038/ncomms14616.
RESUMO
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.
Assuntos
Estudo de Associação Genômica Ampla , Neoplasias Embrionárias de Células Germinativas/genética , Neoplasias Testiculares/genética , Adulto , Cromatina/genética , Perfilação da Expressão Gênica , Predisposição Genética para Doença , Genótipo , Humanos , Masculino , Pessoa de Meia-Idade , Anotação de Sequência Molecular , Neoplasias Embrionárias de Células Germinativas/epidemiologia , Polimorfismo de Nucleotídeo Único , Risco , Neoplasias Testiculares/epidemiologia , Adulto JovemRESUMO
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.
Assuntos
Cromossomos Humanos Par 10/genética , Regulação Leucêmica da Expressão Gênica , Loci Gênicos/genética , Predisposição Genética para Doença , Leucemia-Linfoma Linfoblástico de Células Precursoras/genética , Alelos , Linhagem Celular Tumoral , Criança , Pré-Escolar , Subunidade alfa 3 de Fator de Ligação ao Core/genética , Subunidade alfa 3 de Fator de Ligação ao Core/metabolismo , Proteínas de Ligação a DNA/genética , Proteínas de Ligação a DNA/metabolismo , Diploide , Elementos Facilitadores Genéticos/genética , Epigênese Genética , Variação Genética , Estudo de Associação Genômica Ampla , Humanos , Polimorfismo de Nucleotídeo Único , Análise de Sequência de RNA , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismoRESUMO
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.
Assuntos
Proteínas Adaptadoras de Transdução de Sinal/genética , Elementos Facilitadores Genéticos , Predisposição Genética para Doença , Leucemia Linfocítica Crônica de Células B/genética , Polimorfismo de Nucleotídeo Único , Proteínas Proto-Oncogênicas c-bcl-2/genética , Fator de Transcrição RelA/genética , Proteínas Adaptadoras de Transdução de Sinal/metabolismo , Alelos , Linfócitos B/metabolismo , Linfócitos B/patologia , Sítios de Ligação , Linhagem Celular Tumoral , Cromatina/química , Cromatina/metabolismo , Mapeamento Cromossômico , Cromossomos Humanos Par 15 , Loci Gênicos , Estudo de Associação Genômica Ampla , Histonas/genética , Histonas/metabolismo , Humanos , Leucemia Linfocítica Crônica de Células B/metabolismo , Leucemia Linfocítica Crônica de Células B/patologia , Razão de Chances , Ligação Proteica , Proteínas Proto-Oncogênicas c-bcl-2/metabolismo , Risco , Fator de Transcrição RelA/metabolismoRESUMO
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.
Assuntos
Mieloma Múltiplo/genética , Proteínas Adaptadoras de Transdução de Sinal/genética , Proteína 5 Relacionada à Autofagia/genética , Estudos de Casos e Controles , Proteínas Cromossômicas não Histona , Inibidor p16 de Quinase Dependente de Ciclina , Inibidor de Quinase Dependente de Ciclina p18/genética , Predisposição Genética para Doença , Estudo de Associação Genômica Ampla , Fatores de Troca do Nucleotídeo Guanina/genética , Humanos , Complexo Repressor Polycomb 2/genética , RNA Longo não Codificante/genética , Fatores de Transcrição/genética , Ubiquitina-Proteína Ligases/genéticaRESUMO
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.
Assuntos
Hidrocarboneto de Aril Hidroxilases/genética , Neoplasias da Mama/genética , Leucemia Linfocítica Crônica de Células B/genética , Neoplasias Pulmonares/genética , Polimorfismo de Nucleotídeo Único/genética , Adulto , Idoso , Idoso de 80 Anos ou mais , Alelos , Neoplasias da Mama/urina , Citocromo P-450 CYP3A , Estrona/urina , Feminino , Genótipo , Glucuronídeos/urina , Humanos , Leucemia Linfocítica Crônica de Células B/urina , Neoplasias Pulmonares/urina , Masculino , Pessoa de Meia-Idade , Adulto JovemRESUMO
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.
Assuntos
Cromatina/metabolismo , Neoplasias Colorretais/genética , Loci Gênicos , Predisposição Genética para Doença , Pareamento de Bases/genética , Linhagem Celular Tumoral , Cromossomos Humanos Par 8/genética , Humanos , Hibridização in Situ Fluorescente , Anotação de Sequência Molecular , Motivos de Nucleotídeos/genética , Fatores de Risco , Estatística como AssuntoRESUMO
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).
Assuntos
Genoma Humano , Genômica/métodos , Software , Mapeamento Cromossômico , Bases de Dados Genéticas , Humanos , InternetRESUMO
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.
Assuntos
Cromossomos Humanos Par 10 , Fator de Transcrição GATA3/genética , Fenótipo , Fosfotransferases (Aceptor do Grupo Álcool)/genética , Polimorfismo de Nucleotídeo Único , Leucemia-Linfoma Linfoblástico de Células Precursoras B/genética , Leucemia-Linfoma Linfoblástico de Células Precursoras/genética , Fatores Etários , Alelos , Linfócitos B/metabolismo , Linfócitos B/patologia , Criança , Éxons , Feminino , Loci Gênicos , Predisposição Genética para Doença , Estudo de Associação Genômica Ampla , Humanos , Íntrons , Masculino , Pessoa de Meia-Idade , Leucemia-Linfoma Linfoblástico de Células Precursoras B/mortalidade , Leucemia-Linfoma Linfoblástico de Células Precursoras B/patologia , Leucemia-Linfoma Linfoblástico de Células Precursoras/mortalidade , Leucemia-Linfoma Linfoblástico de Células Precursoras/patologia , Análise de SobrevidaRESUMO
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.
Assuntos
Aberrações Cromossômicas , Cromossomos Humanos , Predisposição Genética para Doença , Mieloma Múltiplo/genética , Estudos de Casos e Controles , HumanosRESUMO
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.