A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification.

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

Massively parallel reporter assay confirms regulatory potential of hQTLs and reveals important variants in lupus and other autoimmune diseases.

We designed a massively parallel reporter assay (MPRA) in an Epstein-Barr virus transformed B cell line to directly characterize the potential for histone post-translational modifications, i.e., histone quantitative trait loci (hQTLs), expression QTLs (eQTLs), and variants on systemic lupus erythematosus (SLE) and autoimmune (AI) disease risk haplotypes to modulate regulatory activity in an allele-dependent manner. Our study demonstrates that hQTLs, as a group, are more likely to modulate regulatory activity in an MPRA compared with other variant classes tested, including a set of eQTLs previously shown to interact with hQTLs and tested AI risk variants. In addition, we nominate 17 variants (including 11 previously unreported) as putative causal variants for SLE and another 14 for various other AI diseases, prioritizing these variants for future functional studies in primary and immortalized B cells. Thus, we uncover important insights into the mechanistic relationships among genotype, epigenetics, and gene expression in SLE and AI disease phenotypes.

Massively Parallel Reporter Assay Confirms Regulatory Potential of hQTLs and Reveals Important Variants in Lupus and Other Autoimmune Diseases.

Objective: To systematically characterize the potential for histone post-translational modifications, i.e., histone quantitative trait loci (hQTLs), expression QTLs (eQTLs), and variants on systemic lupus erythematosus (SLE) and autoimmune (AI) disease risk haplotypes to modulate gene expression in an allele dependent manner. Methods: We designed a massively parallel reporter assay (MPRA) containing ~32K variants and transfected it into an Epstein-Barr virus transformed B cell line generated from an SLE case. Results: Our study expands our understanding of hQTLs, illustrating that epigenetic QTLs are more likely to contribute to functional mechanisms than eQTLs and other variant types, and a large proportion of hQTLs overlap transcription start sites (TSS) of noncoding RNAs. In addition, we nominate 17 variants (including 11 novel) as putative causal variants for SLE and another 14 for various other AI diseases, prioritizing these variants for future functional studies primary and immortalized B cells. Conclusion: We uncover important insights into the mechanistic relationships between genotype, epigenetics, gene expression, and SLE and AI disease phenotypes.

Single-cell multiomic analysis reveals the involvement of Type I interferon-responsive CD8+ T cells in amyloid beta-associated memory loss.

Alzheimer's disease (AD) is the leading cause of dementia worldwide, but there are limited therapeutic options and no current cure. While the involvement of microglia in AD has been highly appreciated, the role of other innate and adaptive immune cells remains largely unknown, partly due to their scarcity and heterogeneity. This study aimed to study non-microglial immune cells in wild type and AD-transgenic mouse brains across different ages. Our results uncovered the presence of a unique CD8+ T cell population that were selectively increased in aging AD mouse brains, here referred to as "disease-associated T cells (DATs)". These DATs were found to express an elevated tissue-resident memory and Type I interferon-responsive gene signature. Further analysis of aged AD mouse brains showed that these CD8+ T cells were not present in peripheral or meningeal tissues. Preventing CD8+ T cell development in AD-transgenic mice via genetic deletion of beta-2 microglobulin ( B2m ) led to a reduction of amyloid-ß plaque formation in aged mice, and improved memory in AD-transgenic mice as early as four months of age. The integration of transcriptomic and epigenomic profiles at the single-cell level revealed potential transcription factors that reshape the regulomes of CD8+ T cells. These findings highlight a critical role for DATs in the progression of AD and provide a new avenue for treatment.

Systemic lupus erythematosus variants modulate the function of an enhancer upstream of TNFAIP3.

TNFAIP3/A20 is a prominent autoimmune disease risk locus that is correlated with hypomorphic TNFAIP3 expression and exhibits complex chromatin architecture with over 30 predicted enhancers. This study aimed to functionally characterize an enhancer ∼55 kb upstream of the TNFAIP3 promoter marked by the systemic lupus erythematosus (SLE) risk haplotype index SNP, rs10499197. Allele effects of rs10499197, rs58905141, and rs9494868 were tested by EMSA and/or luciferase reporter assays in immune cell types. Co-immunoprecipitation, ChIP-qPCR, and 3C-qPCR were performed on patient-derived EBV B cells homozygous for the non-risk or SLE risk TNFAIP3 haplotype to assess haplotype-specific effects on transcription factor binding and chromatin regulation at the TNFAIP3 locus. This study found that the TNFAIP3 locus has a complex chromatin regulatory network that spans ∼1M bp from the promoter region of IL20RA to the 3' untranslated region of TNFAIP3. Functional dissection of the enhancer demonstrated co-dependency of the RelA/p65 and CEBPB binding motifs that, together, increase IL20RA and IFNGR1 expression and decreased TNFAIP3 expression in the context of the TNFAIP3 SLE risk haplotype through dynamic long-range interactions up- and downstream. Examination of SNPs in linkage disequilibrium (D' = 1.0) with rs10499197 identified rs9494868 as a functional SNP with risk allele-specific increase in nuclear factor binding and enhancer activation in vitro. In summary, this study demonstrates that SNPs carried on the ∼109 kb SLE risk haplotype facilitate hypermorphic IL20RA and IFNGR1 expression, while suppressing TNFAIP3 expression, adding to the mechanistic potency of this critically important locus in autoimmune disease pathology.

Variants on the UBE2L3/YDJC Autoimmune Disease Risk Haplotype Increase UBE2L3 Expression by Modulating CCCTC-Binding Factor and YY1 Binding.

OBJECTIVE: Genetic variants spanning UBE2L3 are associated with increased expression of the UBE2L3-encoded E2 ubiquitin-conjugating enzyme H7 (UbcH7), which facilitates activation of proinflammatory NF-κB signaling and susceptibility to autoimmune diseases. We undertook this study to delineate how genetic variants carried on the UBE2L3/YDJC autoimmune risk haplotype function to drive hypermorphic UBE2L3 expression. METHODS: We used bioinformatic analyses, electrophoretic mobility shift assays, and luciferase reporter assays to identify and functionally characterize allele-specific effects of risk variants positioned in chromatin accessible regions of immune cells. Chromatin conformation capture with quantitative polymerase chain reaction (3C-qPCR), chromatin immunoprecipitation (ChIP)-qPCR, and small interfering RNA (siRNA) knockdown assays were performed on patient-derived Epstein-Barr virus-transformed B cells homozygous for the UBE2L3/YDJC nonrisk or risk haplotype to determine if the risk haplotype increases UBE2L3 expression by altering the regulatory chromatin architecture in the region. RESULTS: Of the 7 prioritized variants, 5 demonstrated allele-specific increases in nuclear protein binding affinity and regulatory activity. High-throughput sequencing of chromosome conformation capture coupled with ChIP (HiChIP) and 3C-qPCR uncovered a long-range interaction between the UBE2L3 promoter (rs140490, rs140491, rs11089620) and the downstream YDJC promoter (rs3747093) that was strengthened in the presence of the UBE2L3/YDJC risk haplotype, and correlated with the loss of CCCTC-binding factor (CTCF) and gain of YY1 binding at the risk alleles. Depleting YY1 by siRNA disrupted the long-range interaction between the 2 promoters and reduced UBE2L3 expression. CONCLUSION: The UBE2L3/YDJC autoimmune risk haplotype increases UBE2L3 expression through strengthening a YY1-mediated interaction between the UBE2L3 and YDJC promoters.

Role of Systemic Lupus Erythematosus Risk Variants With Opposing Functional Effects as a Driver of Hypomorphic Expression of TNIP1 and Other Genes Within a Three-Dimensional Chromatin Network.

OBJECTIVE: Genetic variants in the region of tumor necrosis factor-induced protein 3-interacting protein 1 (TNIP1) are associated with autoimmune disease and reduced TNIP1 gene expression. The aim of this study was to define the functional genetic mechanisms driving TNIP1 hypomorphic expression imparted by the systemic lupus erythematosus-associated TNIP1 H1 risk haplotype. METHODS: Dual luciferase expression and electrophoretic mobility shift assays were used to evaluate the allelic effects of 11 risk variants on enhancer function and nuclear protein binding in immune cell line models (Epstein-Barr virus [EBV]-transformed human B cells, Jurkat cells, and THP-1 cells), left in a resting state or stimulated with phorbol 12-myristate 13-acetate/ionomycin. HiChIP was used to define the regulatory 3-dimensional (3-D) chromatin network of the TNIP1 haplotype by detecting in situ long-range DNA contacts associated with H3K27ac-marked chromatin in EBV B cells. Then, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used to determine the expression of genes within the 3-D chromatin network. RESULTS: Bioinformatics analyses of 50 single-nucleotide polymorphisms on the TNIP1 H1 risk haplotype identified 11 non-protein-coding variants with a high likelihood of influencing TNIP1 gene expression. Eight variants in EBV B cells, 5 in THP-1 cells, and 2 in Jurkat cells exhibited various allelic effects on enhancer activation, resulting in a cumulative suppressive effect on TNIP1 expression (net effect of risk variants -7.14 fold, -6.80 fold, and -2.44 fold, respectively; n > 3). Specifically, in EBV B cells, only 2 variants (rs10057690 and rs13180950) exhibited allele-specific loss of both enhancer activity and nuclear protein binding (each P < 0.01 relative to nonrisk alleles). In contrast, the rs10036748 risk allele reduced binding affinities of the transcriptional repressors basic helix-loop-helix family member 40/differentially expressed in chondrocytes 1 (bHLHe40/DEC1) (P < 0.05 relative to nonrisk alleles) and CREB-1 (P not significant) in EBV B cells, resulting in a gain of enhancer activity (P < 0.05). HiChIP and qRT-PCR analyses revealed that overall transcriptional repression of the TNIP1 haplotype extended to the neighboring genes DCTN4 and GMA2, both of which also showed decreased expression in the presence of the TNIP1 risk haplotype (P < 0.001 and P < 0.01, respectively, relative to the nonrisk haplotype); notably, it was found that these genes share a 3-D chromatin network. CONCLUSION: Hypomorphic TNIP1 expression results from the combined concordant and opposing effects of multiple risk variants carried on the TNIP1 risk haplotype, with the strongest regulatory effect in B lymphoid lineage cells. Furthermore, the TNIP1 risk haplotype effect extends to neighboring genes within a shared chromatin network.