Genetic and epigenetic variability contributes to the susceptibility and pathogenesis of autoimmune diseases. T cells play an important role in several autoimmune conditions, including lupus, which is more common and more severe in people of African descent. To investigate inherent epigenetic differences in T cells between ethnicities, we characterized genome-wide DNA methylation patterns in naïve CD4+ T cells in healthy African-Americans and European-Americans, and then confirmed our findings in lupus patients. Impressive ethnicity-specific clustering of DNA methylation profiling in naïve CD4+ T cells was revealed. Hypomethylated loci in healthy African-Americans were significantly enriched in pro-apoptotic and pro-inflammatory genes. We also found hypomethylated genes in African-Americans to be disproportionately related to autoimmune diseases including lupus. We then confirmed that these genes, such as IL32, CD226, CDKN1A, and PTPRN2 were similarly hypomethylated in lupus patients of African-American compared to European-American descent. Using patch DNA methylation and luciferase reporter constructs, we showed that methylation of the IL32 promoter region reduces gene expression in vitro. Importantly, bisulfite DNA sequencing demonstrated that cis-acting genetic variants within and directly disrupting CpG sites account for some ethnicity-specific variability in DNA methylation. Ethnicity-specific inherited epigenetic susceptibility loci in CD4+ T cells provide clues to explain differences in the susceptibility to autoimmunity and possibly other T cell-related diseases between populations.

Objective We examined genome-wide DNA methylation changes in CD8+ T cells from patients with lupus and controls and investigated the functional relevance of some of these changes in lupus. Methods Genome-wide DNA methylation of lupus and age, sex and ethnicity-matched control CD8+ T cells was measured using the Infinium MethylationEPIC arrays. Measurement of relevant cell subsets was performed via flow cytometry. Gene expression was quantified by qPCR. Inhibiting STAT1 and CIITA was performed using fludarabine and CIITA siRNA, respectively. Results Lupus CD8+ T cells had 188 hypomethylated CpG sites compared with healthy matched controls. Among the most hypomethylated were sites associated with HLA-DRB1. Genes involved in the type-I interferon response, including STAT1, were also found to be hypomethylated. IFNα upregulated HLA-DRB1 expression on lupus but not control CD8+ T cells. Lupus and control CD8+ T cells significantly increased STAT1 mRNA levels after treatment with IFNα. The expression of CIITA, a key interferon/STAT1 dependent MHC-class II regulator, is induced by IFNα in lupus CD8+ T cells, but not healthy controls. CIITA knockdown and STAT1 inhibition experiments revealed that HLA-DRB1 expression in lupus CD8+ T cells is dependent on CIITA and STAT1 signalling. Coincubation of naïve CD4+ T cells with IFNα-treated CD8+ T cells led to CD4+ T cell activation, determined by increased expression of CD69 and cytokine production, in patients with lupus but not in healthy controls. This can be blocked by neutralising antibodies targeting HLA-DR. Conclusions Lupus CD8+ T cells are epigenetically primed to respond to type-I interferon. We describe an HLA-DRB1+ CD8+ T cell subset that can be induced by IFNα in patients with lupus. A possible pathogenic role for CD8+ T cells in lupus that is dependent on a high type-I interferon environment and epigenetic priming warrants further characterisation.

As immunotherapies targeting the PDL1 checkpoint have become a mainstay of treatment for a subset of head and neck squamous cell carcinoma (HNSCC) patients, a detailed understanding of the mechanisms underlying PDL1-mediated immune evasion is needed. To elucidate factors regulating expression of PDL1 in HNSCC cells, a genome-wide CRISPR profiling approach was implemented to identify genes and pathways conferring altered PDL1 expression in an HNSCC cell line model. Our screen nominated several candidate PDL1 drivers, including Toll-like Receptor 2 (TLR2). Depletion of TLR2 blocks interferon-γ-induced PDL1 expression, and stimulation of TLR2 with either Staphylococcus aureus or a bacterial lipopeptide mimetic, Pam3CSK4, enhanced PDL1 expression in multiple models. The data herein demonstrate a role for TLR2 in modulating the expression of PDL1 in HNSCC models and suggest that microbiota may directly modulate immunosuppression in cancer cells. Our study represents a step toward disentangling the diverse pathways and stimuli regulating PDL1 expression in HNSCC and underscores a need for future work to characterize the complex microbiome in HNSCC patients treated with immunotherapy.

The MHC region encodes HLA genes and is the most complex region in the human genome. The extensively polymorphic nature of the HLA hinders accurate localization and functional assessment of disease risk loci within this region. Using targeted capture sequencing and constructing individualized genomes for transcriptome alignment, we identified 908 novel transcripts within the human MHC region. These include 593 novel isoforms of known genes, 137 antisense strand RNAs, 119 novel long intergenic noncoding RNAs, and 5 transcripts of 3 novel putative protein-coding human endogenous retrovirus genes. We revealed allele-dependent expression imbalance involving 88% of all heterozygous transcribed single nucleotide polymorphisms throughout the MHC transcriptome. Among these variants, the genetic variant associated with Behçet’s disease in the HLA-B/MICA region, which tags HLA-B*51, is within novel long intergenic noncoding RNA transcripts that are exclusively expressed from the haplotype with the protective but not the disease risk allele. Further, the transcriptome within the MHC region can be defined by 14 distinct coexpression clusters, with evidence of coregulation by unique transcription factors in at least 9 of these clusters. Our data suggest a very complex regulatory map of the human MHC, and can help uncover functional consequences of disease risk loci in this region.

Objective: Behçet's disease (BD) is a systemic inflammatory disease with an incompletely understood etiology. Despite the identification of multiple common genetic variants associated with BD, rare genetic variants have been less explored. We undertook this study to investigate the role of rare variants in BD by performing whole exome sequencing in BD patients of European descent.Methods: Whole exome sequencing was performed in a discovery set comprising 14 German BD patients of European descent. For replication and validation, Sanger sequencing and Sequenom genotyping were performed in the discovery set and in 2 additional independent sets of 49 German BD patients and 129 Italian BD patients of European descent. Genetic association analysis was then performed in BD patients and 503 controls of European descent. Functional effects of associated genetic variants were assessed using bioinformatic approaches.Results: Using whole exome sequencing, we identified 77 rare variants (in 74 genes) with predicted protein‐damaging effects in BD. These variants were genotyped in 2 additional patient sets and then analyzed to reveal significant associations with BD at 2 genetic variants detected in all 3 patient sets that remained significant after Bonferroni correction. We detected genetic association between BD and LIMK2 (rs149034313), involved in regulating cytoskeletal reorganization, and between BD and NEIL1 (rs5745908), involved in base excision DNA repair (P = 3.22 × 10⁻⁴ and P = 5.16 × 10⁻⁴, respectively). The LIMK2 association is a missense variant with predicted protein damage that may influence functional interactions with proteins involved in cytoskeletal regulation by Rho GTPase, inflammation mediated by chemokine and cytokine signaling pathways, T cell activation, and angiogenesis (Bonferroni‐corrected P = 5.63 × 10⁻¹⁴, P = 7.29 × 10⁻⁶, P = 1.15 × 10⁻⁵, and P = 6.40 × 10⁻³, respectively). The genetic association in NEIL1 is a predicted splice donor variant that may introduce a deleterious intron retention and result in a noncoding transcript variant.Conclusion: We used whole exome sequencing in BD for the first time and identified 2 rare putative protein‐damaging genetic variants associated with this disease. These genetic variants might influence cytoskeletal regulation and DNA repair mechanisms in BD and might provide further insight into increased leukocyte tissue infiltration and the role of oxidative stress in BD.