Systematic identification of genotype-dependent enhancer variants in eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is a rare atopic disorder associated with esophageal dysfunction, including difficulty swallowing, food impaction, and inflammation, that develops in a small subset of people with food allergies. Genome-wide association studies (GWASs) have identified 9 independent EoE risk loci reaching genome-wide significance (p < 5 × 10-8) and 27 additional loci of suggestive significance (5 × 10-8 < p < 1 × 10-5). In the current study, we perform linkage disequilibrium (LD) expansion of these loci to nominate a set of 531 variants that are potentially causal. To systematically interrogate the gene regulatory activity of these variants, we designed a massively parallel reporter assay (MPRA) containing the alleles of each variant within their genomic sequence context cloned into a GFP reporter library. Analysis of reporter gene expression in TE-7, HaCaT, and Jurkat cells revealed cell-type-specific gene regulation. We identify 32 allelic enhancer variants, representing 6 genome-wide significant EoE loci and 7 suggestive EoE loci, that regulate reporter gene expression in a genotype-dependent manner in at least one cellular context. By annotating these variants with expression quantitative trait loci (eQTL) and chromatin looping data in related tissues and cell types, we identify putative target genes affected by genetic variation in individuals with EoE. Transcription factor enrichment analyses reveal possible roles for cell-type-specific regulators, including GATA3. Our approach reduces the large set of EoE-associated variants to a set of 32 with allelic regulatory activity, providing functional insights into the effects of genetic variation in this disease.

Epigenetic and transcriptional dysregulation in CD4+ T cells in patients with atopic dermatitis.

Atopic dermatitis (AD) is one of the most common skin disorders among children. Disease etiology involves genetic and environmental factors, with 29 independent AD risk loci enriched for risk allele-dependent gene expression in the skin and CD4+ T cell compartments. We investigated the potential epigenetic mechanisms responsible for the genetic susceptibility of CD4+ T cells. To understand the differences in gene regulatory activity in peripheral blood T cells in AD, we measured chromatin accessibility (an assay based on transposase-accessible chromatin sequencing, ATAC-seq), nuclear factor kappa B subunit 1 (NFKB1) binding (chromatin immunoprecipitation with sequencing, ChIP-seq), and gene expression levels (RNA-seq) in stimulated CD4+ T cells from subjects with active moderate-to-severe AD, as well as in age-matched non-allergic controls. Open chromatin regions in stimulated CD4+ T cells were highly enriched for AD genetic risk variants, with almost half of the AD risk loci overlapping AD-dependent ATAC-seq peaks. AD-specific open chromatin regions were strongly enriched for NF-κB DNA-binding motifs. ChIP-seq identified hundreds of NFKB1-occupied genomic loci that were AD- or control-specific. As expected, the AD-specific ChIP-seq peaks were strongly enriched for NF-κB DNA-binding motifs. Surprisingly, control-specific NFKB1 ChIP-seq peaks were not enriched for NFKB1 motifs, but instead contained motifs for other classes of human transcription factors, suggesting a mechanism involving altered indirect NFKB1 binding. Using DNA sequencing data, we identified 63 instances of altered genotype-dependent chromatin accessibility at 36 AD risk variant loci (30% of AD risk loci) that might lead to genotype-dependent gene expression. Based on these findings, we propose that CD4+ T cells respond to stimulation in an AD-specific manner, resulting in disease- and genotype-dependent chromatin accessibility alterations involving NFKB1 binding.

Shared and distinct interactions of type 1 and type 2 Epstein-Barr Nuclear Antigen 2 with the human genome.

BACKGROUND: There are two major genetic types of Epstein-Barr Virus (EBV): type 1 (EBV-1) and type 2 (EBV-2). EBV functions by manipulating gene expression in host B cells, using virus-encoded gene regulatory proteins including Epstein-Barr Nuclear Antigen 2 (EBNA2). While type 1 EBNA2 is known to interact with human transcription factors (hTFs) such as RBPJ, EBF1, and SPI1 (PU.1), type 2 EBNA2 shares only ~ 50% amino acid identity with type 1 and thus may have distinct binding partners, human genome binding locations, and functions. RESULTS: In this study, we examined genome-wide EBNA2 binding in EBV-1 and EBV-2 transformed human B cells to identify shared and unique EBNA2 interactions with the human genome, revealing thousands of type-specific EBNA2 ChIP-seq peaks. Computational predictions based on hTF motifs and subsequent ChIP-seq experiments revealed that both type 1 and 2 EBNA2 co-occupy the genome with SPI1 and AP-1 (BATF and JUNB) hTFs. However, type 1 EBNA2 showed preferential co-occupancy with EBF1, and type 2 EBNA2 preferred RBPJ. These differences in hTF co-occupancy revealed possible mechanisms underlying type-specific gene expression of known EBNA2 human target genes: MYC (shared), CXCR7 (type 1 specific), and CD21 (type 2 specific). Both type 1 and 2 EBNA2 binding events were enriched at systemic lupus erythematosus (SLE) and multiple sclerosis (MS) risk loci, while primary biliary cholangitis (PBC) risk loci were specifically enriched for type 2 peaks. CONCLUSIONS: This study reveals extensive type-specific EBNA2 interactions with the human genome, possible differences in EBNA2 interaction partners, and a possible new role for type 2 EBNA2 in autoimmune disorders. Our results highlight the importance of considering EBV type in the control of human gene expression and disease-related investigations.

Epstein-Barr virus nuclear antigen 2 extensively rewires the human chromatin landscape at autoimmune risk loci.

The interplay between environmental and genetic factors plays a key role in the development of many autoimmune diseases. In particular, the Epstein-Barr virus (EBV) is an established contributor to multiple sclerosis, lupus, and other disorders. Previously, we showed that the EBV nuclear antigen 2 (EBNA2) transactivating protein occupies up to half of the risk loci for a set of seven autoimmune disorders. To further examine the mechanistic roles played by EBNA2 at these loci on a genome-wide scale, we globally examined gene expression, chromatin accessibility, chromatin looping, and EBNA2 binding in a B cell line that was (1) uninfected, (2) infected with a strain of EBV lacking EBNA2, or (3) infected with a strain that expresses EBNA2. We identified more than 400 EBNA2-dependent differentially expressed human genes and more than 5000 EBNA2 binding events in the human genome. ATAC-seq analysis revealed more than 2000 regions in the human genome with EBNA2-dependent chromatin accessibility, and HiChIP data revealed more than 1700 regions where EBNA2 altered chromatin looping interactions. Autoimmune genetic risk loci were highly enriched at the sites of these EBNA2-dependent chromatin-altering events. We present examples of autoimmune risk genotype-dependent EBNA2 events, nominating genetic risk mechanisms for autoimmune risk loci such as ZMIZ1 Taken together, our results reveal important interactions between host genetic variation and EBNA2-driven disease mechanisms. Further, our study highlights a critical role for EBNA2 in rewiring human gene regulatory programs through rearrangement of the chromatin landscape and nominates these interactions as components of genetic mechanisms that influence the risk of multiple autoimmune diseases.

A functional mechanism for a non-coding variant near AGTR2 associated with risk for preterm birth.

BACKGROUND: Preterm birth (PTB), defined as delivery before 37 gestational weeks, imposes significant public health burdens. A recent maternal genome-wide association study of spontaneous PTB identified a noncoding locus near the angiotensin II receptor type 2 (AGTR2) gene. Genotype-Tissue Expression data revealed that alleles associated with decreased AGTR2 expression in the uterus were linked to an increased risk of PTB and shortened gestational duration. We hypothesized that a causative variant in this locus modifies AGTR2 expression by altering transcription factor (TF) binding. METHODS: To investigate this hypothesis, we performed bioinformatics analyses and functional characterizations at the implicated locus. Potential causal single nucleotide polymorphisms (SNPs) were prioritized, and allele-dependent binding of TFs was predicted. Reporter assays were employed to assess the enhancer activity of the top PTB-associated non-coding variant, rs7889204, and its impact on TF binding. RESULTS: Our analyses revealed that rs7889204, a top PTB-associated non-coding genetic variant is one of the strongest eQTLs for the AGTR2 gene in uterine tissue samples. We observed differential binding of CEBPB (CCAAT enhancer binding protein beta) and HOXA10 (homeobox A10) to the alleles of rs7889204. Reporter assays demonstrated decreased enhancer activity for the rs7889204 risk "C" allele. CONCLUSION: Collectively, these results demonstrate that decreased AGTR2 expression caused by reduced transcription factor binding increases the risk for PTB and suggest that enhancing AGTR2 activity may be a preventative measure in reducing PTB risk.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus.

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Genetic Associations with Gestational Duration and Spontaneous Preterm Birth.

BACKGROUND: Despite evidence that genetic factors contribute to the duration of gestation and the risk of preterm birth, robust associations with genetic variants have not been identified. We used large data sets that included the gestational duration to determine possible genetic associations. METHODS: We performed a genomewide association study in a discovery set of samples obtained from 43,568 women of European ancestry using gestational duration as a continuous trait and term or preterm (<37 weeks) birth as a dichotomous outcome. We used samples from three Nordic data sets (involving a total of 8643 women) to test for replication of genomic loci that had significant genomewide association (P<5.0×10-8) or an association with suggestive significance (P<1.0×10-6) in the discovery set. RESULTS: In the discovery and replication data sets, four loci (EBF1, EEFSEC, AGTR2, and WNT4) were significantly associated with gestational duration. Functional analysis showed that an implicated variant in WNT4 alters the binding of the estrogen receptor. The association between variants in ADCY5 and RAP2C and gestational duration had suggestive significance in the discovery set and significant evidence of association in the replication sets; these variants also showed genomewide significance in a joint analysis. Common variants in EBF1, EEFSEC, and AGTR2 showed association with preterm birth with genomewide significance. An analysis of mother-infant dyads suggested that these variants act at the level of the maternal genome. CONCLUSIONS: In this genomewide association study, we found that variants at the EBF1, EEFSEC, AGTR2, WNT4, ADCY5, and RAP2C loci were associated with gestational duration and variants at the EBF1, EEFSEC, and AGTR2 loci with preterm birth. Previously established roles of these genes in uterine development, maternal nutrition, and vascular control support their mechanistic involvement. (Funded by the March of Dimes and others.).

Macrophage memories of early-life injury drive neonatal nociceptive priming.

The developing peripheral nervous and immune systems are functionally distinct from those of adults. These systems are vulnerable to early-life injury, which influences outcomes related to nociception following subsequent injury later in life (i.e., "neonatal nociceptive priming"). The underpinnings of this phenomenon are unclear, although previous work indicates that macrophages are trained by inflammation and injury. Our findings show that macrophages are both necessary and partially sufficient to drive neonatal nociceptive priming, possibly due to a long-lasting remodeling in chromatin structure. The p75 neurotrophic factor receptor is an important effector in regulating neonatal nociceptive priming through modulation of the inflammatory profile of rodent and human macrophages. This "pain memory" is long lasting in females and can be transferred to a naive host to alter sex-specific pain-related behaviors. This study reveals a mechanism by which acute, neonatal post-surgical pain drives a peripheral immune-related predisposition to persistent pain following a subsequent injury.

Filaggrin loss-of-function variants are associated with atopic dermatitis phenotypes in a diverse, early-life prospective cohort.

Loss-of-function (LoF) variants in the filaggrin (FLG) gene are the strongest known genetic risk factor for atopic dermatitis (AD), but the impact of these variants on AD outcomes is poorly understood. We comprehensively identified genetic variants through targeted region sequencing of FLG in children participating in the Mechanisms of Progression of Atopic Dermatitis to Asthma in Children cohort. Twenty FLG LoF variants were identified, including 1 novel variant and 9 variants not previously associated with AD. FLG LoF variants were found in the cohort. Among these children, the presence of 1 or more FLG LoF variants was associated with moderate/severe AD compared with those with mild AD. Children with FLG LoF variants had a higher SCORing for Atopic Dermatitis (SCORAD) and higher likelihood of food allergy within the first 2.5 years of life. LoF variants were associated with higher transepidermal water loss (TEWL) in both lesional and nonlesional skin. Collectively, our study identifies established and potentially novel AD-associated FLG LoF variants and associates FLG LoF variants with higher TEWL in lesional and nonlesional skin.

Human cytomegalovirus infection coopts chromatin organization to diminish TEAD1 transcription factor activity.

Human cytomegalovirus (HCMV) infects up to 80% of the world's population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hours after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon-6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV's established roles in these processes.

Macrophage epigenetic memories of early life injury drive neonatal nociceptive priming.

The developing peripheral nervous and immune systems are functionally distinct from adults. These systems are vulnerable to early life injury, which influences outcomes related to nociception following subsequent injury later in life (neonatal nociceptive priming). The underpinnings of this phenomenon are largely unknown, although previous work indicates that macrophages are epigenetically trained by inflammation and injury. We found that macrophages are both necessary and partially sufficient to drive neonatal nociceptive priming possibly due to a long-lasting epigenetic remodeling. The p75 neurotrophic factor receptor (NTR) was an important effector in regulating neonatal nociceptive priming through modulation of the inflammatory profile of rodent and human macrophages. This pain memory was long lasting in females and could be transferred to a naive host to alter sex-specific pain-related behaviors. This study reveals a novel mechanism by which acute, neonatal post-surgical pain drives a peripheral immune-related predisposition to persistent pain following a subsequent injury.

A DNA tumor virus globally reprograms host 3D genome architecture to achieve immortal growth.

Epstein-Barr virus (EBV) immortalization of resting B lymphocytes (RBLs) to lymphoblastoid cell lines (LCLs) models human DNA tumor virus oncogenesis. RBL and LCL chromatin interaction maps are compared to identify the spatial and temporal genome architectural changes during EBV B cell transformation. EBV induces global genome reorganization where contact domains frequently merge or subdivide during transformation. Repressed B compartments in RBLs frequently switch to active A compartments in LCLs. LCLs gain 40% new contact domain boundaries. Newly gained LCL boundaries have strong CTCF binding at their borders while in RBLs, the same sites have much less CTCF binding. Some LCL CTCF sites also have EBV nuclear antigen (EBNA) leader protein EBNALP binding. LCLs have more local interactions than RBLs at LCL dependency factors and super-enhancer targets. RNA Pol II HiChIP and FISH of RBL and LCL further validate the Hi-C results. EBNA3A inactivation globally alters LCL genome interactions. EBNA3A inactivation reduces CTCF and RAD21 DNA binding. EBNA3C inactivation rewires the looping at the CDKN2A/B and AICDA loci. Disruption of a CTCF site at AICDA locus increases AICDA expression. These data suggest that EBV controls lymphocyte growth by globally reorganizing host genome architecture to facilitate the expression of key oncogenes.

Global discovery of lupus genetic risk variant allelic enhancer activity.

Genome-wide association studies of Systemic Lupus Erythematosus (SLE) nominate 3073 genetic variants at 91 risk loci. To systematically screen these variants for allelic transcriptional enhancer activity, we construct a massively parallel reporter assay (MPRA) library comprising 12,396 DNA oligonucleotides containing the genomic context around every allele of each SLE variant. Transfection into the Epstein-Barr virus-transformed B cell line GM12878 reveals 482 variants with enhancer activity, with 51 variants showing genotype-dependent (allelic) enhancer activity at 27 risk loci. Comparison of MPRA results in GM12878 and Jurkat T cell lines highlights shared and unique allelic transcriptional regulatory mechanisms at SLE risk loci. In-depth analysis of allelic transcription factor (TF) binding at and around allelic variants identifies one class of TFs whose DNA-binding motif tends to be directly altered by the risk variant and a second class of TFs that bind allelically without direct alteration of their motif by the variant. Collectively, our approach provides a blueprint for the discovery of allelic gene regulation at risk loci for any disease and offers insight into the transcriptional regulatory mechanisms underlying SLE.

Genetic, Inflammatory, and Epithelial Cell Differentiation Factors Control Expression of Human Calpain-14.

Eosinophilic esophagitis (EoE) is a chronic, food-driven allergic disease resulting in eosinophilic esophageal inflammation. We recently found that EoE susceptibility is associated with genetic variants in the promoter of CAPN14, a gene with reported esophagus-specific expression. CAPN14 is dynamically up-regulated as a function of EoE disease activity and after exposure of epithelial cells to interleukin-13 (IL-13). Herein, we aimed to explore molecular modulation of CAPN14 expression. We identified three putative binding sites for the IL-13-activated transcription factor STAT6 in the promoter and first intron of CAPN14 Luciferase reporter assays revealed that the two most distal STAT6 elements were required for the ∼10-fold increase in promoter activity subsequent to stimulation with IL-13 or IL-4, and also for the genotype-dependent reduction in IL-13-induced promoter activity. One of the STAT6 elements in the promoter was necessary for IL-13-mediated induction of CAPN14 promoter activity while the other STAT6 promoter element was necessary for full induction. Chromatin immunoprecipitation in IL-13 stimulated esophageal epithelial cells was used to further support STAT6 binding to the promoter of CAPN14 at these STAT6 binding sites. The highest CAPN14 and calpain-14 expression occurred with IL-13 or IL-4 stimulation of esophageal epithelial cells under culture conditions that allow the cells to differentiate into a stratified epithelium. This work corroborates a candidate molecular mechanism for EoE disease etiology in which the risk variant at 2p23 dampens CAPN14 expression in differentiated esophageal epithelial cells following IL-13/STAT6 induction of CAPN14 promoter activity.

Transcription factors operate across disease loci, with EBNA2 implicated in autoimmunity.

Explaining the genetics of many diseases is challenging because most associations localize to incompletely characterized regulatory regions. Using new computational methods, we show that transcription factors (TFs) occupy multiple loci associated with individual complex genetic disorders. Application to 213 phenotypes and 1,544 TF binding datasets identified 2,264 relationships between hundreds of TFs and 94 phenotypes, including androgen receptor in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of systemic lupus erythematosus risk loci are occupied by the Epstein-Barr virus EBNA2 protein and many coclustering human TFs, showing gene-environment interaction. Similar EBNA2-anchored associations exist in multiple sclerosis, rheumatoid arthritis, inflammatory bowel disease, type 1 diabetes, juvenile idiopathic arthritis and celiac disease. Instances of allele-dependent DNA binding and downstream effects on gene expression at plausibly causal variants support genetic mechanisms dependent on EBNA2. Our results nominate mechanisms that operate across risk loci within disease phenotypes, suggesting new models for disease origins.

The antiprotease SPINK7 serves as an inhibitory checkpoint for esophageal epithelial inflammatory responses.

Loss of barrier integrity has an important role in eliciting type 2 immune responses, yet the molecular events that initiate and connect this with allergic inflammation remain unclear. We reveal an endogenous, homeostatic mechanism that controls barrier function and inflammatory responses in esophageal allergic inflammation. We show that a serine protease inhibitor, SPINK7 (serine peptidase inhibitor, kazal type 7), is part of the differentiation program of human esophageal epithelium and that SPINK7 depletion occurs in a human allergic, esophageal condition termed eosinophilic esophagitis. Experimental manipulation strategies reducing SPINK7 in an esophageal epithelial progenitor cell line and primary esophageal epithelial cells were sufficient to induce barrier dysfunction and transcriptional changes characterized by loss of cellular differentiation and altered gene expression known to stimulate allergic responses (for example, FLG and SPINK5). Epithelial silencing of SPINK7 promoted production of proinflammatory cytokines including thymic stromal lymphopoietin (TSLP). Loss of SPINK7 increased the activity of urokinase plasminogen-type activator (uPA), which in turn had the capacity to promote uPA receptor-dependent eosinophil activation. Treatment of epithelial cells with the broad-spectrum antiserine protease, α1 antitrypsin, reversed the pathologic features associated with SPINK7 silencing. The relevance of this pathway in vivo was supported by finding genetic epistasis between variants in TSLP and the uPA-encoding gene, PLAU We propose that the endogenous balance between SPINK7 and its target proteases is a key checkpoint in regulating mucosal differentiation, barrier function, and inflammatory responses and that protein replacement with antiproteases may be therapeutic for select allergic diseases.