Search

The developing peripheral nervous and immune systems are functionally distinct from adults. These systems are vulnerable to effects of early life injury which can influence outcomes related to nociception following subsequent injury later in life (i.e. “neonatal nociceptive priming”). The underpinnings of this phenomenon are largely unknown, although macrophages can be epigenetically trained by injury. We found that macrophages are both necessary and partially sufficient to drive neonatal nociceptive priming possibly due to a long-lasting epigenetic remodeling of peripheral macrophages. The p75 neurotrophic factor receptor (NTR) was observed to be an important effector in regulating neonatal nociceptive priming. p75NTR modulates the inflammatory profile and responses of rodent and human macrophages. This “pain memory” was able to be transferred to a naive host to alter sex-specific pain-related behaviors. This study reveals a novel mechanism by which acute post-surgical pain may transition to chronic pain in children.Graphical Abstract

STUDY OBJECTIVES: Obstructive sleep apnea (OSA) is a chronic and widely prevalent disease associated with multiple health disorders. Current diagnostic strategies for OSA are limited because of cost, time, and access. Epigenetic signatures offer insight into the relationships between disease and environment and could play a significant role in developing both diagnostic and therapeutic tools for OSA. In the current study, a systematic literature search was conducted to investigate the existing evidence of OSA-associated epigenetic modifications. METHODS: A systematic literature search was performed using electronic academic databases including PubMed, CINAHL, Scopus, Embase, EBM Reviews, and Web of Science. However, the current study focused on screening for original, English-language articles pertaining to OSA and associated epigenetic mechanisms. To produce unbiased results, screening was performed independently by authors. RESULTS: We identified 2,944 publications in our systematic search. Among them, 65 research articles were related to OS A-associated differential gene expression, genetic variation, and epigenetic modifications. Although these 65 articles were considered for full manuscript review, only 12 articles met the criteria of OSA-associated epigenetic modifications in human and animal models. Human patients with OSA had unique epigenetic changes compared to healthy control patients and, interestingly, epigenetic signatures were commonly identified in genes associated with metabolic and inflammatory pathways. CONCLUSIONS: Although the available studies are limited, this research provides novel insights for the development of epigenetic markers for the diagnosis and treatment of OSA. Thorough genome-wide investigations will be required to develop cost-effective, robust biomarkers for the identification of OSA among children and adults. Here, we offer a study design for such efforts. CITATION: Leader BA, Koritala BSC, Moore CA, Dean EG, Kottyan LC, Smith DF. Epigenetics of obstructive sleep apnea syndrome: a systematic review. J Clin Sleep Med. 2021;17(12):2533–2541.

The molecular processes underlying human health and disease are highly complex. Often, genetic and environmental factors contribute to a given disease or phenotype in a non-additive manner, yielding a gene–environment (G × E) interaction. In this work, we broadly review current knowledge on the impact of gene–environment interactions on human health. We first explain the independent impact of genetic variation and the environment. We next detail well-established G × E interactions that impact human health involving environmental toxicants, pollution, viruses, and sex chromosome composition. We conclude with possibilities and challenges for studying G × E interactions.

BACKGROUND: Eosinophilic esophagitis (EoE) is an emerging, chronic, rare allergic disease associated with marked eosinophil accumulation in the esophagus. Prior genome-wide association studies (GWAS) have provided strong evidence for three genome-wide susceptibility loci. OBJECTIVE: We aimed to replicate known and suggestive EoE genetic risk loci and conduct a meta-analysis of previously reported datasets. METHODS: An EoE-Custom SNP Chip (EoE-CSC) containing 956 candidate EoE risk single-nucleotide polymorphisms (SNPs) was employed to genotype 627 cases and 365 controls. Statistical power was enhanced by adding 1990 external controls and performing meta-analyses with two independent EoE GWAS. RESULTS: Meta-analysis identified replicated association and genome-wide significance at six loci: 2p23 (two independent genetic effects), 5q22, 10p14, 11q13, and 16p13. Seven additional loci were identified at suggestive significance (p < 10(−6)): 1q31, 5q23, 6q15, 6q21, 8p21, 17q12, and 22q13. From these risk loci, 13 protein-coding EoE candidate risk genes were expressed in a genotype-dependent manner. EoE risk genes were expressed in disease-relevant cell types, including esophageal epithelia, fibroblasts, and immune cells, with some expressed as a function of disease activity. The genetic risk burden of EoE-associated genetic variants was markedly larger in cases relative to controls (P

ObjectiveSystemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci. Nevertheless, these loci only partially explain SLE heritability and their putative causal variants are rarely prioritised, which make challenging to elucidate disease biology. To detect new SLE loci and causal variants, we performed the largest genome-wide meta-analysis for SLE in East Asian populations.MethodsWe newly genotyped 10 029 SLE cases and 180 167 controls and subsequently meta-analysed them jointly with 3348 SLE cases and 14 826 controls from published studies in East Asians. We further applied a Bayesian statistical approach to localise the putative causal variants for SLE associations.ResultsWe identified 113 genetic regions including 46 novel loci at genome-wide significance (p−8). Conditional analysis detected 233 association signals within these loci, which suggest widespread allelic heterogeneity. We detected genome-wide associations at six new missense variants. Bayesian statistical fine-mapping analysis prioritised the putative causal variants to a small set of variants (95% credible set size ≤10) for 28 association signals. We identified 110 putative causal variants with posterior probabilities ≥0.1 for 57 SLE loci, among which we prioritised 10 most likely putative causal variants (posterior probability ≥0.8). Linkage disequilibrium score regression detected genetic correlations for SLE with albumin/globulin ratio (rg=−0.242) and non-albumin protein (rg=0.238).ConclusionThis study reiterates the power of large-scale genome-wide meta-analysis for novel genetic discovery. These findings shed light on genetic and biological understandings of SLE.

This study tests the hypothesis that activation of MAPK by physiologically relevant concentrations of IL-33 contributes to enhanced cytokine expression by IL-12 stimulated human NK cells. While IL-33 canonically triggers type 2 cytokine responses, this cytokine can also synergize with type 1 cytokines like IL-12 to provoke IFN-γ. We show that picogram concentrations of IL-12 and IL-33 are sufficient to promote robust secretion of IFN-γ by human NK cells that greatly exceeds resposes to either cytokine alone. Nanogram doses of IL-33, potentially consistent with levels in tissue microenvironments, synergize with IL-12 to induce secretion of additional cytokines, including TNF and GM-CSF. IL-33-induced activation of the p38 MAPK pathway in human NK cells is crucial for enhanced release of IFN-γ and TNF in response to IL-12. Mechanistically, IL-33-induced p38 MAPK signaling enhances stability of IFNG transcripts and triggers A disintegrin and metalloproteinase domain 17 (ADAM17) mediated cleavage of TNF from the cell surface. These data support our hypothesis and suggest that altered sensitivity of NK cells to IL-12 in the presence of IL-33 may have important consequences in diseases associated with mixed cytokine milieus, like asthma and chronic obstructive pulmonary disease.

Transcription factors read the genome, fundamentally connecting DNA sequence to gene expression across diverse cell types. Determining how, where, and when TFs bind chromatin will advance our understanding of gene regulatory networks and cellular behavior. The 2017 ENCODE-DREAM in vivo Transcription-Factor Binding Site (TFBS) Prediction Challenge highlighted the value of chromatin accessibility data to TFBS prediction, establishing state-of-the- art methods for TFBS prediction from DNase-seq. However, the more recent Assay-for- Transposase-Accessible-Chromatin (ATAC)-seq has surpassed DNase-seq as the most widely- used chromatin accessibility profiling method. Furthermore, ATAC-seq is the only such technique available at single-cell resolution from standard commercial platforms. While ATAC-seq datasets grow exponentially, suboptimal motif scanning is unfortunately the most common method for TFBS prediction from ATAC-seq. To enable community access to state-of-the-art TFBS prediction from ATAC-seq, we (1) curated an extensive benchmark dataset (127 TFs) for ATAC-seq model training and (2) built “maxATAC”, a suite of user-friendly, deep neural network models for genome-wide TFBS prediction from ATAC-seq in any cell type. With models available for 127 human TFs, maxATAC is the first collection of high-performance TFBS prediction models for ATAC-seq. maxATAC performance extends to primary cells and single-cell ATAC-seq, enabling improved TFBS prediction in vivo. We demonstrate maxATAC’s capabilities by identifying TFBS associated with allele-dependent chromatin accessibility at atopic dermatitis genetic risk loci.Author SummaryProteins called transcription factors interpret the genome, reading both DNA sequence and chromatin state, to orchestrate gene expression across the diversity of human cell types. In any given cell type, most chromatin is “inaccessible”, and only those parts of the genetic code needed or likely to be needed soon are “accessible” for transcription factor binding to affect gene expression and cellular behavior. Hundreds of transcription factors are expressed in a given cell type and context (e.g., age, disease), and knowledge of their context-specific DNA binding sites is key to uncovering how transcription factors regulate cellular behaviors in health or disease. However, experimentally profiling the >1,600 human transcription factors across all cell types and contexts is infeasible. We built a suite of computational models “maxATAC” to predict transcription factor binding from a measurement of accessible chromatin, ATAC-seq. Importantly, ATAC-seq is feasible even at single-cell resolution. Thus, this data type, in combination with maxATAC, can be used to infer transcription factor binding sites in directly-relevant cell types isolated from physiological and disease settings, enabling insights into disease mechanisms, including how genetic variants and cellular context impact transcription factor binding, gene expression patterns and disease risk.

Nature Communications 13, 6859 (2022). doi:10.1038/s41467-022-34456-6, Published by Nature Publishing Group UK, [London]

Eosinophilic esophagitis (EoE) is a chronic allergic inflammatory disease with a complex underlying genetic etiology. Herein, we conduct whole-exome sequencing of a multigeneration EoE pedigree (discovery set) and 61 additional multiplex families with EoE (replication set). A series of rare, heterozygous, missense variants are identified in the genes encoding the desmosome-associated proteins DSP and PPL in 21% of the multiplex families. Esophageal biopsies from patients with these variants retain dilated intercellular spaces and decrease DSP and PPL expression even during disease remission. These variants affect barrier integrity, cell motility and RhoGTPase activity in esophageal epithelial cells and have increased susceptibility to calpain-14–mediated degradation. An acquired loss of esophageal DSP and PPL is present in non-familial EoE. Taken together, herein, we uncover a pathogenic role for desmosomal dysfunction in EoE, providing a deeper mechanistic understanding of tissue-specific allergic responses., Eosinophilic esophagitis (EoE) is a chronic allergic inflammatory disease with a complex underlying genetic etiology. Here, the authors identify a series of rare variants in DSP and PPL in multiplex families with EoE and uncover a pathogenic role for desmosomal dysfunction in EoE.

IgA nephropathy (IgAN) is a progressive form of kidney disease defined by glomerular deposition of IgA. We performed a genome-wide association study involving 10,146 kidney biopsy-diagnosed IgAN cases and 28,751 matched controls across 17 international cohorts. We defined 30 independent genome-wide significant risk loci jointly explaining 11% of disease risk. A total of 16 loci were novel, including TNFSF4, REL, CD28, CXCL8/PF4V1, LY86, LYN, ANXA3, TNFSF8/15, REEP3, ZMIZ1, RELA, ETS1, IGH, IRF8, TNFRSF13B and FCAR. The SNP-based heritability of IgAN was estimated at 23%. We observed a positive genetic correlation between IgAN and total serum IgA levels, allergy, tonsillectomy, and several infections, and a negative correlation with inflammatory bowel disease. All significant non-HLA loci shared with serum IgA levels had a concordant effect on the risk of IgAN. Moreover, IgAN loci were globally enriched in gene orthologs causing abnormal IgA levels when genetically manipulated in mice. The explained heritability was enriched in the regulatory elements of cells from the immune and hematopoietic systems and intestinal mucosa, providing support for the pathogenic role of extra-renal tissues. The polygenic risk of IgAN was associated with early disease onset, increased lifetime risk of kidney failure, as well as hematuria and several other traits in a phenome-wide association study of 590,515 individuals. In the comprehensive functional annotation analysis of candidate causal genes across genome-wide significant loci, we observed the convergence of biological candidates on a common set of inflammatory signaling pathways and cytokine ligand-receptor pairs, prioritizing potential new drug targets.

Background Genetic variants in complement genes have been associated with a wide range of human disease states, but well-powered genetic association studies of complement activation have not been performed in large multiethnic cohorts. Methods We performed medical records?based genome-wide and phenome-wide association studies for plasma C3 and C4 levels among participants of the Electronic Medical Records and Genomics (eMERGE) network. Results In a GWAS for C3 levels in 3949 individuals, we detected two genome-wide significant loci: chr.1q31.3 (CFH locus; rs3753396-A; ?=0.20; 95% CI, 0.14 to 0.25; P=1.52x10(-11)) and chr.19p13.3 (C3 locus; rs11569470-G; ?=0.19; 95% CI, 0.13 to 0.24; P=1.29x10(-8)). These two loci explained approximately 2% of variance in C3 levels. GWAS for C4 levels involved 3998 individuals and revealed a genome-wide significant locus at chr.6p21.32 (C4 locus; rs3135353-C; ?=0.40; 95% CI, 0.34 to 0.45; P=4.58x10(-35)). This locus explained approximately 13% of variance in C4 levels. The multiallelic copy number variant analysis defined two structural genomic C4 variants with large effect on blood C4 levels: C4-BS (?=?0.36; 95% CI, ?0.42 to ?0.30; P=2.98x10(-22)) and C4-AL-BS (?=0.25; 95% CI, 0.21 to 0.29; P=8.11x10(-23)). Overall, C4 levels were strongly correlated with copy numbers of C4A and C4B genes. In comprehensive phenome-wide association studies involving 102,138 eMERGE participants, we cataloged a full spectrum of autoimmune, cardiometabolic, and kidney diseases genetically related to systemic complement activation. Conclusions We discovered genetic determinants of plasma C3 and C4 levels using eMERGE genomic data linked to electronic medical records. Genetic variants regulating C3 and C4 levels have large effects and multiple clinical correlations across the spectrum of complement-related diseases in humans.Significance Statement The complement pathway represents one of the critical arms of the innate immune system. We combined genome-wide and phenome-wide association studies using medical records data for C3 and C4 levels to discover common genetic variants controlling systemic complement activation. Three genome-wide significant loci had large effects on complement levels. These loci encode three critical complement genes: CFH, C3, and C4. We performed detailed functional annotations of the significant loci, including multiallelic copy number variant analysis of the C4 locus to define two structural genomic variants with large effects on C4 levels. Blood C4 levels were strongly correlated with the copy number of C4A and C4B genes. Lastly, using genome-wide genetic correlations and electronic health records?based phenome-wide association studies in 102,138 participants, we catalogued a spectrum of human diseases genetically related to systemic complement activation, including inflammatory, autoimmune, cardiometabolic, and kidney diseases.

Dysregulated transcription factors represent a major class of drug targets that mediate the abnormal expression of many critical genes involved in SLE and other autoimmune diseases. Although strong evidence suggests that natural human genetic variation affects basal and inducible gene expression, it is still a considerable challenge to establish a biological link between GWAS-identified non-coding genetic risk variants and their regulated gene targets. Here, we combine genetic data, epigenomic data, and CRISPR activation (CRISPRa) assays to screen for functional variants regulating IRF8 expression. Using CRISPR-mediated deletion and 3D chromatin structure analysis, we demonstrate that the locus containing rs2280381 is a cell-type-specific distal enhancer for IRF8 that spatially interacts with the IRF8 promoter. Further, rs2280381 mediates IRF8 expression through enhancer RNA AC092723.1, which recruits TET1 to the IRF8 promoter to modulate IRF8 expression by affecting methylation levels. The alleles of rs2280381 modulate PU.1 binding and chromatin state to differentially regulate AC092723.1 and IRF8 expression. Our work illustrates a strategy to define the functional genetic variants modulating transcription factor gene expression levels and identifies the biologic mechanism by which autoimmune diseases risk genetic variants contribute to the pathogenesis of disease.

The advancement of precision medicine requires new methods to coordinate and deliver genetic data from heterogeneous sources to physicians and patients. The eMERGE III Network enrolled >25,000 participants from biobank and prospective cohorts of predominantly healthy individuals for clinical genetic testing to determine clinically actionable findings. The network developed protocols linking together the 11 participant collection sites and 2 clinical genetic testing laboratories. DNA capture panels targeting 109 genes were used for testing of DNA and sample collection, data generation, interpretation, reporting, delivery, and storage were each harmonized. A compliant and secure network enabled ongoing review and reconciliation of clinical interpretations, while maintaining communication and data sharing between clinicians and investigators. A total of 202 individuals had positive diagnostic findings relevant to the indication for testing and 1,294 had additional/secondary findings of medical significance deemed to be returnable, establishing data return rates for other testing endeavors. This study accomplished integration of structured genomic results into multiple electronic health record (EHR) systems, setting the stage for clinical decision support to enable genomic medicine. Further, the established processes enable different sequencing sites to harmonize technical and interpretive aspects of sequencing tests, a critical achievement toward global standardization of genomic testing. The eMERGE protocols and tools are available for widespread dissemination.

The duration of pregnancy is influenced by fetal and maternal genetic and non-genetic factors. Here we report a fetal genome-wide association meta-analysis of gestational duration, and early preterm, preterm, and postterm birth in 84,689 infants. One locus on chromosome 2q13 is associated with gestational duration; the association is replicated in 9,291 additional infants (combined P = 3.96 × 10−14). Analysis of 15,588 mother-child pairs shows that the association is driven by fetal rather than maternal genotype. Functional experiments show that the lead SNP, rs7594852, alters the binding of the HIC1 transcriptional repressor. Genes at the locus include several interleukin 1 family members with roles in pro-inflammatory pathways that are central to the process of parturition. Further understanding of the underlying mechanisms will be of great public health importance, since giving birth either before or after the window of term gestation is associated with increased morbidity and mortality., Gestational duration depends on both maternal and fetal genetic influences. Here, the authors perform a fetal genome-wide association meta-analysis and find that a locus on 2q13 is associated with pregnancy duration and further show that the lead SNP rs7594852 changes the binding properties of transcriptional repressor HIC1.

BACKGROUND: The incidence of eosinophilic esophagitis (EoE) is greater in male than female subjects, and the underlying molecular basis for this sex bias remains unclear. OBJECTIVE: We sought to delineate the contribution of the sex hormone estrogen to the EoE phenotype and esophageal epithelial barrier function and remodeling. METHODS: We performed demographic and incidence analyses of EoE in male and female subjects from a single-center pediatric cohort. Estrogen-responsive gene expression analyses and estrogen receptor (ESR) immunofluorescence staining of esophageal biopsy specimens from patients with EoE and control subjects were performed. The effect of 17β-estradiol (E2) on IL-13–induced signaling pathways, gene expression, and esophageal epithelial architecture and barrier function in a primary human esophageal keratinocyte cell (EPC2) culture system (EPC2–air-liquid interface) was examined. RESULTS: We observed a male predominance in patients with EoE. Analyses of RNA sequencing data sets revealed a significant dysregulation of the estrogen-responsive gene network and expression of ESR1 and ESR2 in esophageal biopsy specimens from patients with EoE compared with control subjects. IL-13 stimulation of EPC2–air-liquid interface cells led to altered cellular architecture with induced dilation of intercellular spaces and barrier dysfunction. Pretreatment of EPC2s with E2 prior to IL-13 exposure abrogated IL-13– induced architectural changes and esophageal barrier dysfunction. Mechanistically, E2-protective effects were dependent on ESR2 and associated with diminishing of IL-13– induced tyrosine kinase 2 and signal transducer and activator of transcription 6 phosphorylation and EoE-dysregulated gene expression. CONCLUSIONS: Estrogen-responsive genes are modified in patients with EoE compared with control subjects. E2 attenuated IL-13– induced architectural changes and esophageal epithelial barrier dysfunction through inhibition of the IL-13/tyrosine kinase 2/ signal transducer and activator of transcription 6 pathway via ESR2-dependent process. Estrogen hormone signaling may protect against development of EoE in female subjects.

Background Evidence suggests early oral introduction of peanut prevents peanut allergy, particularly in a high risk infant population. Despite early oral introduction guidelines, there will continue to be infants and young children who develop peanut allergy. Current testing modalities utilized for evaluation of peanut allergy have poor positive predictive value, particularly when applied to patients with no history of oral exposure to peanut. The oral food challenge is a procedure conducted to confirm a food allergy diagnosis. This procedure is not without risk and has particular unique challenges in the younger patient population. Objective We sought to identify key differentially expressed genes that may distinguish peanut allergic infants among those at high risk for peanut allergy. Methods Twenty infants aged 4-11 months with egg allergy and/or severe atopic dermatitis (i.e. high risk) were prospectively recruited from a tertiary care clinic. Infants were classified as peanut allergic, peanut sensitized, but tolerant, or peanut non-allergic based on clinical history, skin prick testing, specific IgE, and/or oral food challenge. RNA was isolated from whole blood and samples underwent RNA-seq library preparation and sequencing. Differential expression analysis was performed. Results Eight genes were differentially expressed in peanut sensitized, but tolerant subjects compared to the peanut non-allergic subjects. Four genes were differentially expressed in peanut allergic versus peanut non-allergic subjects and in peanut allergic versus peanut sensitized, but tolerant subjects. Conclusion Differential gene expression in whole peripheral blood was identified between these three groups of high risk infants, however, whole blood may not be an optimal biospecimen for detection of relevant, subtle gene expression changes.

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., Thousands of genetic variants have been associated with lupus, but causal variants and mechanisms are unknown. Here, the authors combine a massively parallel reporter assay with genome-wide ChIP experiments to identify risk variants with allelic enhancer activity mediated through transcription factor binding.

SUMMARY Myeloid lineage cells use TLRs to recognize and respond to diverse microbial ligands. Although unique transcription factors dictate the outcome of specific TLR signaling, whether lineage-specific differences exist to further modulate the quality of TLR-induced inflammation remains unclear. Comprehensive analysis of global gene transcription in human monocytes, monocyte-derived macrophages, and monocyte-derived dendritic cells stimulated with various TLR ligands identifies multiple lineage-specific, TLR-responsive gene programs. Monocytes are hyperresponsive to TLR7/8 stimulation that correlates with the higher expression of the receptors. While macrophages and monocytes express similar levels of TLR4, macrophages, but not monocytes, upregulate interferon-stimulated genes (ISGs) in response to TLR4 stimulation. We find that TLR4 signaling in macrophages uniquely engages transcription factor IRF1, which facilitates the opening of ISG loci for transcription. This study provides a critical mechanistic basis for lineage-specific TLR responses and uncovers IRF1 as a master regulator for the ISG transcriptional program in human macrophages., Graphical abstract, In brief Song et al. examine the responses of human myeloid populations to bacterial and viral ligands and discover the lineage-specific induction of proinflammatory gene expression. TLR4 signaling in macrophages, but not monocytes, enables nuclear translocation of IRF1, which acts as a pioneer transcription factor to increase the chromatin accessibility of IFNB1 and ISG loci.

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.

Runt-related transcription factor 1 (Runx1) can act as both an activator and a repressor. Here we show that CRISPR-mediated deletion of Runx1 in mouse metanephric mesenchyme-derived mK4 cells results in large-scale genome-wide changes to chromatin accessibility and gene expression. Open chromatin regions near down-regulated loci enriched for Runx sites in mK4 cells lose chromatin accessibility in Runx1 knockout cells, despite remaining Runx2-bound. Unexpectedly, regions near upregulated genes are depleted of Runx sites and are instead enriched for Zeb transcription factor binding sites. Re-expressing Zeb2 in Runx1 knockout cells restores suppression, and CRISPR mediated deletion of Zeb1 and Zeb2 phenocopies the gained expression and chromatin accessibility changes seen in Runx1KO due in part to subsequent activation of factors like Grhl2. These data confirm that Runx1 activity is uniquely needed to maintain open chromatin at many loci, and demonstrate that Zeb proteins are required and sufficient to maintain Runx1-dependent genome-scale repression., Author summary Runt-related transcription factor (Runx) 1 & 2 impact development and disease by activating or repressing transcription. In this manuscript we used genome editing tools to remove Runx1, and as expected, observed widespread changes in chromatin accessibility. Newly closed areas contained Runx1 binding sites and were enriched near genes whose expression depended on Runx1. Interestingly, this occurred despite continued binding of Runx2 to the same regions of DNA, which suggests that Runx2 is insufficient to maintain open chromatin and expression of Runx1 target genes in this cellular context. By contrast, newly opened chromatin regions, many near genes that were upregulated in Runx1 knockout cells, did not enrich for Runx1 binding sites. Instead, these regions were enriched for sites for the repressor Zeb proteins. We found that the loss of Zeb 1 & 2 expression, direct transcriptional targets of Runx1, resulted in the opening of chromatin and upregulation of genes residing near the newly open sites in Runx1 knockout cells. The same sites were also open and nearby genes expressed in edited Zeb1 and Zeb2 knockout cells. Among them were transcription factors, such as the Grhl2 gene, which in turn bind to and upregulate their target genes. Thus, the loss of a single transcription factor initiates a cascade of direct and indirect ramifications with likely negative effects on development and health.

Runt-related transcription factor 1 (Runx1) can act as both an activator and a repressor. Here we show that CRISPR-mediated deletion of Runx1 in an embryonic kidney-derived cell (mK4) results in large-scale genome-wide changes to chromatin accessibility and gene expression. Open chromatin regions near down-regulated loci are enriched for Runx sites, remain bound by Runx2, but lose chromatin accessibility and expression in Runx1 knockout cells. Unexpectedly, regions near upregulated genes are depleted of Runx sites and are instead enriched for Zeb transcription factor binding sites. Re-expressing Zeb2 in Runx1 knockout cells restores suppression. These data confirm that Runx1 activity is uniquely needed to maintain open chromatin at many loci, and demonstrate that genome-scale derepression is an indirect consequence of losing Runx1-dependent Zeb expression.

Systemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci1-8. Nevertheless, these loci only partially explain SLE heritability and provide limited biological insight. We report the largest study of SLE in East Asians (13,377 cases and 194,993 controls), identifying 233 association signals within 113 (46 novel) genetic loci. We detect six new lead missense variants and prioritize ten most likely putative causal variants, one of which we demonstrate exhibits allele-specific regulatory effect on ACAP1 in vitro. We suggest 677 effector genes with potential for drug repurposing, and provide evidence that two distinct association signals at a single locus act on different genes (NCF2 and SMG7). We demonstrate that SLE-risk variants overlap with cell-specific active regulatory elements, notably EBNA2-mediated super-enhancers in Epstein-Barr Virus-transformed B cells, and implicate the role for non-immune cells in SLE biology. These findings shed light on genetic and biological understandings of SLE.

Introduction The etiology of steroid-sensitive nephrotic syndrome (SSNS) is not well understood. Genetic studies have established common single nucleotide polymorphisms (SNPs) that are associated with increased SSNS disease risk. We review previous genetic association studies of SSNS and nominate particular transcriptional regulators and immune cells as potential key players in the etiology of this disease. Methods A list of SNPs associated with SSNS was compiled from published genome wide association and candidate gene studies. The Regulatory Element Locus Intersection (RELI) tool was used to calculate the enrichment of the overlap between disease risk SNPs and the genomic coordinates of data from a collection of >10,000 chromatin immunoprecipitation sequencing experiments. Results After linkage disequilibrium expansion of the previously reported tag associated SNPs, we identified 192 genetic variants at 8 independent risk loci. Using the Regulatory Element Locus Intersection algorithm, we identified transcriptional regulators with enriched binding at SSNS risk loci (10-05 < Pcorrected < 10-124), including ZNF530, CIITA, CD74, RFX5, and ZNF425. Many of these regulators have well-described roles in the immune response. RNA polymerase II binding in B cells also demonstrated enriched binding at SSNS risk loci (10-37, Graphical abstract

The human genome contains millions of putative regulatory elements, which regulate gene expression. We are just beginning to understand the functional consequences of genetic variation within these regulatory elements. Since the bulk of common genetic variation impacting polygenic disease phenotypes localizes to these non-coding regions of the genome, understanding the consequences will improve our understanding of the mechanisms mediating genetic risk in human disease. Here, we define the systemic lupus erythematosus (SLE) risk variant rs2431369 as likely causal for SLE and show that it is located in a functional regulatory element that modulates miR-146a expression. We use epigenomic analysis and genome-editing to show that the rs2431697-containing region is a distal enhancer that specifically regulates miR-146a expression in a cell-type dependent manner. 3D chromatin structure analysis demonstrates physical interaction between the rs2431697-containing region and the miR-146a promoter. Further, our data show that NF-kB binds the disease protective allele in a sequence-specific manner, leading to increased expression of this immunoregulatory microRNA. Our work provides a strategy for using disease-associated variants to define the functional regulatory elements of non-coding RNA molecules such as miR-146a and provides mechanistic links between autoimmune disease risk genetic variation and disease etiology.

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 demonstrated 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 >400 EBNA2-dependent differentially expressed human genes and >5,000 EBNA2 binding events in the human genome. ATAC-seq analysis revealed >2,000 regions in the human genome with EBNA2-dependent chromatin accessibility, and HiChIP-seq data revealed >1,700 regions where EBNA2 altered chromatin looping interactions. Importantly, 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.

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.

Genome-wide association studies of Systemic Lupus Erythematosus (SLE) nominate 3,073 genetic variants at 91 risk loci. To systematically screen these variants for allelic transcriptional enhancer activity, we constructed 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 revealed 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, larger 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.

Large meta-analyses of rheumatoid arthritis (RA) susceptibility in European (EUR) and East Asian (EAS) populations have identified >100 RA risk loci, but genome-wide studies of RA in African-Americans (AAs) are absent. To address this disparity, we performed an analysis of 916 AA RA patients and 1392 controls and aggregated our data with genotyping data from >100 000 EUR and Asian RA patients and controls. We identified two novel risk loci that appear to be specific to AAs: GPC5 and RBFOX1 (P(AA) 90 RA risk loci. Addition of AA data to those of EUR and EAS descent enabled identification of seven novel high-confidence candidate pathogenic variants (defined by posterior probability > 0.8). In summary, our trans-ethnic analyses are the first to include AAs, identified several new RA risk loci and point to candidate pathogenic variants that may underlie this common autoimmune disease. These findings may lead to better ways to diagnose or stratify treatment approaches in RA.

BACKGROUND: Eosinophilic esophagitis (EoE) is an esophageal inflammatory disease associated with atopic diseases. Thymic stromal lymphopoietin (TSLP) and calpain 14 (CAPN14) genetic variations contribute to EoE, but how this relates to atopy is unclear. OBJECTIVE: The purpose of this study was to explore the relationship between EoE, atopy, and genetic risk. METHODS: EoE-atopy enrichment was tested by using 700 patients with EoE and 801 community control subjects. Probing 372 single nucleotide polymorphisms (SNPs) in 63 atopy genes, we evaluated EoE associations using 412 nonatopic and 868 atopic disease control subjects. Interaction and stratified analyses of EoE-specific and atopy-related SNPs were performed. RESULTS: Atopic disease was enriched in patients with EoE (P < .0001). Comparing patients with EoE and nonatopic control subjects, EoE associated strongly with IL-4/kinesin family member 3A (IL4/KIF3A) (P = 2.8 × 10(−6); odds ratio [OR], 1.87), moderately with TSLP (P = 1.5 × 10(−4); OR, 1.43), and nominally with CAPN14 (P = .029; OR, 1.35). Comparing patients with EoE with atopic disease control subjects, EoE associated strongly with ST2 (P = 3.5 × 10(−6); OR, 1.77) and nominally with IL4/KIF3A (P = .019; OR, 1.25); TSLP’s association persisted (P = 4.7 × 10(−5); OR, 1.37), and CAPN14’s association strengthened (P = .0001; OR, 1.71). Notably, there was gene-gene interaction between TSLP and IL4 SNPs (P = .0074). Children with risk alleles for both genes were at higher risk for EoE (P = 2.0 × 10(−10); OR, 3.67). CONCLUSIONS: EoE genetic susceptibility is mediated by EoE-specific and general atopic disease loci, which can have synergistic effects. These results might aid in identifying potential therapeutics and predicting EoE susceptibility.