Transcriptional Changes in Regulatory T Cells From Patients With Autoimmune Polyendocrine Syndrome Type 1 Suggest Functional Impairment of Lipid Metabolism and Gut Homing

Autoimmune polyendocrine syndrome type I (APS-1) is a monogenic model disorder of organ-specific autoimmunity caused by mutations in the Autoimmune regulator (AIRE) gene. AIRE facilitates the expression of organ-specific transcripts in the thymus, which is essential for efficient removal of dangerous self-reacting T cells and for inducing regulatory T cells (Tregs). Although reduced numbers and function of Tregs have been reported in APS-I patients, the impact of AIRE deficiency on gene expression in these cells is unknown. Here, we report for the first time on global transcriptional patterns of isolated Tregs from APS-1 patients compared to healthy subjects. Overall, we found few differences between the groups, although deviant expression was observed for the genes TMEM39B, SKIDA1, TLN2, GPR15, FASN, BCAR1, HLA-DQA1, HLA-DQB1, HLA-DRA, GPSM3 and AKR1C3. Of significant interest, the consistent downregulation of GPR15 may indicate failure of Treg gut homing which could be of relevance for the gastrointestinal manifestations commonly seen in APS-1. Upregulated FASN expression in APS-1 Tregs points to increased metabolic activity suggesting a putative link to faulty Treg function. Functional studies are needed to determine the significance of these findings for the immunopathogenesis of APS-1 and for Treg immunobiology in general.

Introduction

Immunological homeostasis and maintenance of health rely on a complex network of structured processes. These are organized in several layers to ensure that the immune system will efficiently stop invaders without harming the host. T cells are central in initiating and maintaining specific immune responses. T cells with suppressive capacity, i.e. regulatory T cells (Tregs), constitute about 5-7% of the CD4⁺ T cell repertoire, and have a critical role in maintaining homeostasis by dampening inexpedient immune responses in the peripheral blood and tissues (1, 2). They are characterized by expression of the linage-defining transcription factor forkhead box P3 (FOXP3), and are equipped with a distinct epigenomic signature, especially at the FOXP3 locus (1–3).

The primary location for development and maturation of Tregs is in the thymus through positive and negative selection in response to self-antigen stimulation (thymic or natural Tregs). Lately, it has become apparent that Tregs are also generated outside the thymus (peripheral Tregs or pTregs) and, at least in mice, can be converted from Cd4⁺Cd62l⁺ central memory T cells in blood (4). Mature Tregs are present in peripheral blood, tissues and especially in draining lymph nodes, where they control antigen presentation. They may also recirculate back into the adult thymus in order to control thymic output of de novo Tregs (5). Another important aspect of Tregs is their ability to mimic phenotypes of other T cells, e.g. Th1 and Th17 through expression of the relevant signature transcription factors (T-bet and RORγt, respectively) (6, 7), and modulate the environment with their cytokine profile. This has challenged the previous view that the CD4⁺ cells are static with either a T effector cell or a Treg phenotype.

Autoimmune polyendocrine syndrome type I (APS-1) is a severe monogenic autoimmune disorder characterized by autoimmune manifestations in endocrine glands and ectodermal tissues. Mutations in the Autoimmune Regulator (AIRE) gene (8–10) cause impaired immunological tolerance. Interestingly, AIRE is involved in all layers of tolerance induction: It orchestrates the expression of tissue-specific antigens in the thymus to be presented to developing T cells during negative selection (11) and has a further distinct role in the generation of Tregs (12). Aire specifically induces fetal generation of distinct subtypes of Cd4⁺FoxP3⁺ nTregs in mice (13), is involved in Ccr6-dependent recirculation of n-/p-Tregs into the thymus (5), and in generating pTregs in the periphery (14). We and others have previously found either decreased frequencies and/or impaired function of Tregs in patients with APS-1 (15–17), however no in-depth studies of the Tregs transcriptomes in these patients have been conducted. Correct function of Tregs is essential for a healthy immune system, exemplified by FoxP3 deficient mice and patients with FOXP3 mutations, both suffering from numerous autoimmune manifestations (3, 18). Also, single nucleotide polymorphisms (SNPs) linked to common autoimmune diseases are enriched in DNA regions of CpG demethylation crucial for Treg cell development and function (19).

With an increase in available immune therapies aiming to utilize the specific capacity and plasticity of Tregs to downregulate immune responses, it is crucial to understand suppressive and counter-acting activities of Tregs, and their contribution to autoimmune disease.

In order to investigate whether Treg transcriptional activity and not just selection is affected by AIRE dysfunction, we have undertaken a descriptive study of isolated Tregs from APS-I patients and compared both protein expression and transcriptomes to healthy controls. We hypothesize that transcriptome analyses of purified Tregs in APS-1 patients, with their dual failures of central and peripheral tolerance, could shed light on the pleiotropic nature of AIRE in Tregs and immune function in general.

Methods

Patients and Controls

Nine patients with APS-1 (confirmed with AIRE mutational analysis and autoantibody screening against IFN-ω) [44% males, mean age 49.2 years (range 30-63)] and 21 healthy subjects were included in this study (Table 1). The patients have been reported previously (20). Seventeen of the healthy controls were included in the RNA sequencing of Tregs (HC1-17), 59% males, mean age 49.2 years (range 20-76), while eight healthy subjects were included in the flow cytometry analyses of Tregs (HC1-4+18-21), 87% males, mean age 44.4 years (range 25-66).

TABLE 1

Table 1 Demographics of included APS-1 patients.

Flow Cytometry Analysis

The Treg flow cytometry panel was inspired by a clinical Tregs workshop (21). Blood was sampled in heparin-containing tubes and peripheral mononuclear cells (PBMCs) were isolated by Ficoll density gradient separation following standard procedures before freezing down at -150°C in AB serum with 10% DMSO. For flow analysis, thawed PBMCs from 5 patients and 8 controls (see Table 1) were added two microliter of Fc block for 15 min before addition of the extracellular antibodies anti-CTLA4 BV421 (BN13, Biolegend, San Diego, CA), anti-CD39 PE (ebioA1, Invitrogen, Carlsbad, California, USA) and anti-CD3 V500 (UCHT1); anti-CD4 Alexa Fluor 700 (RPA-T4); anti-CD8 PerCP-Cy5.5 (SKI); anti-CD25 PE-Cy7 (2A3); anti-CD45RA APC-H7 (HI100) and anti-CD31 BV650 (L133.1) from BD (Franklin Lakes, New Jersey, USA). The cells were resuspended in 1 mL PBS and live/dead Fixable Yellow Dead Cell stain kit (Invitrogen) was added in a 1:1000 dilution. Staining of intracellular markers [anti-FOXP3 PE-CF594 (259D/C7, BD), anti-HELIOS APC (22F6, Biolegend), and anti-Ki-67 (20Raj1, Invitrogen)] was subsequently performed using the eBioscience Anti-human FOXP3 Staining Set as explained by the manufacturer with the change that the cells were fixed overnight. Acquisition was executed by a LSR Fortessa flow cytometer and data was analyzed by FlowJo version 10.4.

Isolation of Tregs and Purification of RNA

Fresh blood samples (15-25 ml) were collected in EDTA tubes and Tregs were enriched using the MacsExpress Treg isolation kit according to the manufacturer’s protocol (Miltenyi Biotech, Bergisch Gladbach, Germany), resulting in samples constituting 43-77% Tregs (mean controls 65%; mean patients 60%), validated by flow cytometry (CD4⁺CD25^hiCD127^-). Positively selected cells were prepared for flow cytometry sorting (FACS) by staining with antibodies obtained from Miltenyi Biotech: anti-CD4 FITC (M-T466), anti-CD25 PE (4E3) and anti-CD127 APC (cMB15-18C9). After 30 min incubation in the dark at 4°C, the cells were washed by addition of 2 ml rinsing buffer (Miltenyi Biotech) and centrifugation for 10 min at 350g at 4°C. The resulting pellet was resuspended in PBS with 0.5% BSA and CD4⁺CD25^hiCD127^- cells (Tregs) were sorted into aliquots of 5000 cells by the use of BD FACSAria SORP or FACSymphony S6 directly into the buffer containing 75 μl RLT plus 10% v/v β-mercaptoethanol. RNA from the lysate was immediately purified using the RNeasy Micro Plus Kit, following the recommendations from the supplier (Qiagen, Venlo, Netherlands). The RNA quality was assessed by the Agilent Bioanalyzer, using the Agilent 6000 Pico kit or by TapeStation with the High Sensitivity D5000 ScreenTape kit (Agilent, Santa Clara, CA, USA). Concentrations were measured by Qubit RNA HS Assay kit (Invitrogen). Samples were stored at -80°C until use.

RNA Sequencing

cDNA was made from RNA and amplified using the SMART-Seq v4 Ultra Low Input RNA Kit according to the manufacturer’s protocol (Takara Bio, Madison, Wisconsin, USA). Undiluted RNA samples from the 5000 cells aliquots were used as input. cDNA amplification followed a program of 10 (for samples P1-P4 and HC1-HC4) and 14 (for samples P5-P9, R1 and HC 5-17) cycles (see Table 1 for subject information). The cDNA was purified by the Agencourt AMPure XP Kit (Beckman Coulter, Pasadena, CA, USA) and assessed by an Agilent Bioanalyzer (Agilent High Sensitivity DNA Kit) or TapeStation (High Sensitivity D5000 ScreenTape kit) according to the manufacturer’s protocols (Agilent). The samples were then stored at -20°C before use.

Library preparation was performed using the Nextera XT DNA Library Preparation kit and the Nextera XT Index kit according to the manufacturer’s protocol (Illumina, San Diego, CA, USA). The input cDNA was diluted according to the results from the Agilent Bioanalyzer/TapeStation and Qubit™ dsDNA HS Assay Kit (concentrations between 11-63 ng/microliter and addition of 3.5 to 5 microliter for library preparations). Unique combinations of i5 (5 μl) and i7 (5 μl) index adapters were added to each respective sample. Constructed libraries were then assessed with the High sensitivity D5000 ScreenTape kit, by TapeStation (Agilent). The final libraries consisted of fragments mostly between 250-2500 bp, with a peak of 500-1000. Sequencing of the products was performed by the Bergen Genomics core facility, using a HiSeq4000 sequencer according to the Illumina TruSeq Stranded mRNA protocol.

Bioinformatics Analysis of RNA Sequencing Data

Quality control of the sequencing data was performed using the FastQC software (v.0.11.9) (22), and reads were aligned to the GRCh38.p12 (build 100) reference transcriptome by the Kallisto pseudoaligner [v 0.46.2 (23)] in quant pair-end mode using a k-mer length of 31. The pseudoalignments were subsequently aggregated to the genome using the tximport (v.1.16.1) (24) package in R (v.4.0.2) and annotated using biomaRt (v.2.44.4) (25). Differential expression analysis was performed using the DESeq2 package (v.1.28.1) (26) using the default significant level of 10% FDR. In order to maximize power in this analysis the Independent Hypothesis Weighting (IHW R package) (27) multiple testing method was used to control for false discovery rate, and the Approximate posterior estimation for GLM coefficients (apeglm R package) (28) method for log2 fold change shrinkage was used to remove noise. BAM files for gene visualization were generated using STAR aligner (v. 2.7.5a) (29) and indexed using Sambamba (v. 0.8.0) (30).

Statistical Analysis and Plot Generation

All datasets were imported into RStudio [v. 1.3.1073 (31)], and were treated using tidy data principles with the tidyverse R packages (v. 1.3.1) (32). Tidyverse package ggplot2 was used to visualize the data in various volcano plots, violin plots, dot plots, box plots, and heatmaps. The cowplot package (v. 1.1.1) was utilized to combine different ggplots, tidyverse package scales was used to make sensible axis scales, ggrepel (v. 0.9.1) was used to annotate plots without overlapping labels, and ggpubr (v. 0.4.0) was used to make some publication ready plots. Colorblind friendly color scales were generated using RColorBrewer (v. 1.1-2) and Viridis [v. 0.6.0 (33)]. A non-parametric, two-tailed Mann-Whitney U test by the rstatix package (v. 0.7.0) was applied to analyze the differences between patients and controls regarding the Tregs subpopulations in flow cytometry, and results were adjusted for multiple testing using Benjamini-Hochberg FDR correction. Results were found statistically significant if FDR < 0.05. Two-sided Mann-whitney U tests were also performed on the age distribution data of the sequenced control and patient groups, with a significance threshold level of P <= 0.05.

Gene visualizations of reads were performed using the Gviz R package [v. 1.34.1 (34)] where EnsDb.Hsapiens.v100 were used as the reference genome supplied by the AnnotationHub R package (v. 2.22.1).

Pathway Analysis

String (https://string-db.org/cgi/network) was utilized to perform pathway analysis for the differentially expressed genes identified. This web tool can mediate identification of affected pathways by finding associations between selected proteins based on text-mining, co-mentions in PubMed abstracts, known interactions (experimentally determined or from curated databases) or predicted interactions.

Ethics

This project was conducted in compliance with the Declaration of Helsinki and approved by the Regional Ethics Committee of Western Norway (approval numbers 2009/2555 and 2018/1417). All patients were recruited from the Registry for organ specific autoimmune disorders (ROAS), Haukeland University Hospital, Norway, and gave written informed consent for participation. Samples from healthy subjects were obtained from the Haukeland University Hospital blood bank. Data was stored and analyzed on secure servers throughout the whole project (SAFE, University of Bergen; TSD, University of Oslo and HUNT Cloud, Norwegian University of Science and Technology, all Norway).

Results

Flow Cytometry Profiling of Tregs From APS-1 Patients

In order to investigate the impact of AIRE mutations on Treg distribution and protein expression within Tregs, PBMCs from five of the nine included APS-1 patients and eight healthy subjects were profiled for common Tregs surface markers using flow cytometry (Figure 1A). Tregs were identified by the markers CD4⁺CD25^hiFOXP3⁺, and subtypes hinting at functional properties were then assessed within this cohort of cells according to the gating strategy outlined in Figure 1B and shown in Supplementary Figure 1. There were no major differences between patients and controls in the frequency of CD4⁺ cells within CD3⁺ PBMCs, nor in Tregs within the CD4⁺ parent cohort [Tregs mean patients 5.5% (range 4.4-6.7%) versus mean controls 6.1% (range 3.0-9.9%)]. Furthermore, there were no notable variation in the thymic recent emigrant/epithelial marker CD31, the ectoenzyme for functional capacity of Tregs CD39, the naïve marker CD45RA, nor the activation marker Ki67 although one outlier in each of CD45RA and Ki-67 in the patient group exhibited increased expression compared to controls. An overall trend towards lower amounts of HELIOS-expression and higher levels of CTLA-4- signals could be observed within the APS-1 cohort, but these observations did not reach statistical significance after multiple testing corrections (Figure 1A).

FIGURE 1

Figure 1 Violin plots of flow cytometry data from a subset of APS-1 patients and healthy subjects. (A) Plots show different cell populations positive for markers and their percentages relative to the gated parent population. Statistical analysis was performed by using Mann-Whitney U test with Benjamini-Hochberg correction (FDR) in R. Gating strategy is displayed in (B).

Sample Characterization and Quality Controls of the cDNA Libraries Made From Tregs-RNA

Tregs were isolated from whole blood by magnetic beads and then FACS from a total of 9 APS-1 patients and 17 age and sex matched healthy subjects (Supplementary Figure 2A). The largest group of patients (n=5) shared the 13bp-deletion mutation in AIRE. Although two of these were homozygous for c.976-979del13, P1 had this mutation together with a large deletion spanning at least exon 1 to intron 8 (35) while two were compound heterozygous for the 13 bp deletion together with two different missense mutations (classified in Table 1). There were no differences in age or sex distribution between patients and controls (Supplementary Figures 2B, C).

Because of limitations regarding the collection of Tregs from APS-1 patients, the sample collection, library preparation and sequencing were performed in two batches, visualized in Supplementary Figure 2D. There were no statistical differences in the age distribution nor genders within each batch between patients and controls (Supplementary Figures 2E, F).

Five-thousand cells were FACS sorted from up-concentrated Tregs of each patient and control, RNA was isolated and then sequenced. The first batch was sequenced to an average read depth of 70.3 million bases, while the second batch was sequenced to an average read depth of 42 million bases (Supplementary Figure 3).

Principal component analysis (PCA) plots were compiled for the sequencing data with the two largest PCs explaining 17% and 12% of the variance, respectively. The metadata variable closest to explaining this variance was batch 1 or 2, indicating a possible issue with batch effects (Supplementary Figure 4). Hence, batch-variance were included as a covariate in the differential expression model design to analyze gene expression in Tregs from APS-1 patients compared to controls. However, as a quality control measure, two individuals were sampled in both sequencing batches representing two different timepoints for Treg-isolation. PCA plots of these patients displayed very little difference between batches, indicating that the impact of any batch effect may be small. We included only the initial samples for downstream analyses.

Treg Characterization by Transcriptomic Profiling

Successful Treg isolation and sequencing was verified by exploring the gene expression of five genes characteristic of Tregs (36). Interrogating the gene expression of the master regulator of the Treg pathway FOXP3, HELIOS (IKZF2), Interleukin-2 receptor alpha and beta chains (IL2RA and IL2RB), as well as TNF receptor superfamily member 1B (TNFRSF1B also known as CD120b) showed relatively high expression in transcripts per million (TPM) (Figure 2A). They were further expressed at similar levels with no significant difference in distribution between patients and control. This expression signature confirmed that the isolated cells were Tregs, and that there is little difference in these signature Treg genes between the healthy controls and APS-1 patients within the dataset. Analyzing whether the distribution of reads conforms with the exons within the gene or whether these reads may be misaligned can reveal false differential gene expression findings caused by low expression levels. In this case the RNA sequencing read distributions across these genes showed stable exon expression (data not shown).

FIGURE 2

Figure 2 Violin plots of RNA expression levels in transcripts per million (TPM) of important regulatory T cell markers in Tregs of APS-1 patients (cyan) and healthy subjects (red). (A) APS-1 mutations are indicated by the color of the dots within the violin plot. AIRE expression in TPM for the different groups shown in (B) Histogram of distribution of patients and controls according to genotype in (C) FDR values taken from differential expression analysis using DESeq2 in R.

As the AIRE gene, which is impaired in APS-1 patients, is implicated in Tregs development and function (14–17, 37, 38), we performed a detailed analysis of its’ expression. The results show low level (sub 2 TPM) AIRE expression in most patients and controls (evenly distributed between the different mutations). However, there were no significant difference in gene expression between the groups, and no consistent difference in expression of AIRE or the 5 genes characteristic of Tregs between the AIRE mutations encountered in this dataset (Figures 2B, C). While expression levels were relatively low, read distribution conformed mainly to within exons in AIRE, confirming AIRE mRNA expression.

Differential Expression of Genes in APS-1 Patients Versus Controls

Differential expression analysis was performed using DESeq2, taking advantage of an experimental design consisting of batch and condition, where condition indicated APS-1 or control and batch was treated as a covariate in the statistical modelling. Overall, the results indicated very small differences between APS-1 patients and healthy subjects as most genes were not differentially expressed (Tables 2A, B and Figure 3A). The most differentially expressed genes were Transmembrane protein 39 (TMEM39B), Ski/Dach domain-containing protein 1 (SKIDA1), and Talin 2 (TLN2), expressed at a 5% FDR level at log2 fold changes of 0.61, 2.12, and -3.82e-06 respectively. Using a less stringent 10% FDR level, the genes G Protein-Coupled Receptor 15 (GPR15), Fatty Acid Synthase (FASN), and Cas Scaffolding Protein Family Member 1 (BCAR1) were also found to be differentially expressed, with log2 fold changes of -1.46e-05, 0.98, and 3.03. Of these, SKIDA1, BCAR1 and TLN2 were lowly expressed at below 2 TPM (Figure 3B). Gene visualization of read distribution using the Gviz r package shows reads aligning to exons in BCAR1, FASN, GPAR15 and TMEM39B, while reads in SKIDA1 and TLN2 seem more spurious (Supplementary Figure 5).

TABLE 2A

Table 2A Gene expression significantly up or down-regulated in APS-1 patients’ Tregs’ compared with healthy subjects (Adjusted p<0.1).

TABLE 2B

Table 2B Gene expression significantly up or down-regulated in APS-1 patients’ Tregs’ compared with healthy controls (FDR<0.1).

FIGURE 3

Figure 3 Differential expression analysis of regulatory T cells from APS-1 patients and healthy subjects. (A). Volcano plot of the most significantly differentially expressed genes with the magnitude of differential expression in log2 fold change and the significance level in -log10 FDR. (B) Violin plots showing the absolute expression level in transcripts per million (TPM) in all samples for the genes identified as significantly differentially expressed at a 10% FDR level in (A). FDR values on violin plots taken from differential expression analysis using DESeq2 (v. 1.30.1)in R (v. 4.0.2.).

Because clinical and experimental findings suggest that different AIRE mutations have different impact on AIRE function and activity (e.g. from complete loss of function to residual activity) (39, 40), differential expression analysis was also performed using only the subset of APS-1 patients (n=5) with the 13bp deletion mutation, the most common mutation in the dataset. Results of this subgroup analysis showed that the MHC-II complex alleles HLA-DQA1, HLA-DQB1, HLA-DRA, and G Protein Signaling Modulator 3 (GPSM3) were significantly differentially expressed at a 5% FDR significance level between patients with the 13bp deletion and healthy controls (Figure 4A), with log2 fold changes of -2.14e-07, -5.26e-07, -1.21e-06, and -1.90e-07 respectively. Reducing the significance level to 10% FDR, TMEM39B again appeared significant in this subgroup and Aldo-Keto Reductase Family 1 Member C (AKR1C3) showed differential expression with log2 fold changes of 0.56 and -6.91e-06. Notably, for GPSM3 and the HLA-DR and HLA-DQ hits, this might be driven by high expression in a few of the controls (Figure 4B).

FIGURE 4

Figure 4 Differential expression analysis of regulatory T cells from APS-1 patients with the 13bp deletion and controls. (A) Volcano plot of the most significantly differentially expressed genes with the magnitude of differential expression in log2 fold change and the significance level in -log10 FDR. (B) Violin plots showing the absolute expression level in transcripts per million (TPM) in all samples for the genes identified as significantly differentially expressed at a 10% FDR level in (A) FDR values on violin plots taken from differential expression analysis using DESeq2 in R.

To investigate whether subtle changes in selected Treg related genes could be missed by a stringent per gene significant level as reported by Ferraro et al. (41), volcano plots were generated on the subset of genes reported as “193 genes that are most differentially expressed between Tregs and effector T cells”. However, none of these genes were found differentially expressed at a sub significant level in the dataset when comparing Tregs from APS-1 patients to healthy controls (data not shown).

Pathway Analyses

Analysis using String was performed to investigate whether any functional link could be identified between the differentially expressed genes in this study, and additional well-established Treg associated genes (Supplementary Figure 6). The five highly relevant Tregs-genes FOXP3, IKZF2, IL2RA, IL2RB, and TRNRSF1B were included in the model. Furthermore, String gene ontology enrichment analysis was performed in order to investigate any links between gene ontology and biological processes for the differentially expressed genes (Supplementary Table 1). As antigen presenting cells and their HLA-proteins which present peptides to T cells have pivotal roles in immune reactivity and tolerance, it is to no surprise that String reported an association between FOXP3 and HLA-DQB1 as “being co-mentioned” in PubMed-abstracts in numerous articles, e.g (42). The gene ontology enrichment strengthened this association by pointing out the relevance of HLA-molecules for proper immune function. String further highlighted a link between SKIDA1 and IKZF2, based on experimental evidence and text-mining on interaction of putative homologues of these proteins in other organisms (C. elegans and D. melanogaster). TLN2 and BCAR1 were suggested to occasionally be co-expressed in homo sapiens, which could imply functional associations in one or several of the following functional pathways; membrane “ruffle” as part of cell motility, Rap1-signalling and/or in cellular focal adhesion.

Discussion

We have shown, for the first time, that transcriptional changes within Tregs of APS-1 patients could imply impairment in Treg lipid metabolism and gut homing, respectively. Although our flow cytometry and transcriptomic data show little overall difference between APS-1 patients and healthy controls, some important identified deteriorations could interfere with proper Treg function. For the flow cytometry profiling, there were some trends among the APS-1 patients of increased frequencies of the inhibitory receptor CTLA-4, and decreased frequencies of HELIOS and CD31 expressing cells. However, with few patients the data is difficult to interpret. These findings could however corroborate previously suggested impairment of thymic generation of Tregs in APS-1 (15, 37).

Using differential expression analysis based on RNAseq data, we found that the six genes GPR15, TMEM39B, FASN, TLN2, SKIDA1 and BCAR1 differ significantly in their expression level between patients with APS-1 and controls. Although this low number of differentially expressed genes supports that there are no extensive impairments of Tregs in APS-1 patients, several of the identified putative dysregulated Tregs-genes have important functional implications for Tregs.

One of our more interesting findings concerns GPR15, encoding a G protein-coupled receptor that acts as a homing receptor for the gut and the lamina propria (43). The protein is regulated by the Ahr-RORγt-FOXP3 axis (44) and has been found to be upregulated in immune cells in blood and affected tissues in several autoimmune disorders; thus it could be important in modulating inflammation (45–48). Interestingly, the gut-homing property seems to include FOXP3⁺ Tregs in particular, indicating that GPR15 plays a role in mucosal immune tolerance largely by regulating the influx of Tregs. Recently, it was demonstrated that patients with colorectal cancer have increased frequencies of GPR15⁺ Tregs in peripheral blood, correlating with high numbers of GPR15⁺ Tregs infiltrating the colonic cancerous lesions (49). These observations indicate that GPR15 directs Tregs into the colon, thereby modifying the tumor microenvironment and promoting intestinal tumorigenesis. Considering this phenomenon, it is important to bear in mind that a substantial fraction of APS-1 patients have gastrointestinal manifestations and highly dysregulated immune responses against intestinal commensal microbiota (50, 51). The latter has also been shown to be associated with Tregs dysfunction. In addition, gut biopsies of APS-1 patients have revealed reduced numbers of FOXP3+ cells in these tissues compared to healthy controls (51). It may be possible that the consistently downregulated GPR15 expression in APS-1 patients’ Tregs could reflect an opposite pattern than colorectal cancer: Lower expression of GPR15 in peripheral blood Tregs, alternatively lower numbers of GPR15⁺ Tregs in general, and consequently lower numbers of Tregs in the gastrointestinal compartments potentially unleashing a proinflammatory and autoimmunity-prone microenvironment. GPR15⁺ Tregs may also be indirectly linked to AIRE deficiency as tumor associated GPR15⁺ Tregs in colorectal cancer has been shown to have a thymic origin (49).

The most statistically significant finding in this study was the identification of higher expression of TMEM39B in APS-1 patients. Both TMEM39B and its paralogue TMEM39A has been shown to be important for autophagy, possibly by controlling the spatial localization and synthesis of phosphatidylinositol 4-phosphate (52). Although autophagy mediates Tregs lineage stability and survival fitness (53), we cannot confirm any connection between TMEM39B, Tregs and autophagy. TMEM39B is, according to the Human Protein Atlas, expressed in many tissues, but particularly in the thymus and lymphoid tissues. It has also been shown to be overexpressed in lymphoid malignancies, indicating involvement in lymphocytic cellular growth and survival (54). Additional studies are needed to investigate the consequence of increased TMEM39B expression, and autophagy in APS-1 patients.

FASN is another differentially expressed gene which showed upregulation in Tregs from APS-1 patients compared to healthy controls. This enzyme is important in fatty acid synthesis and metabolism by catalyzing the condensation of acetyl-CoA and malonyl-CoA to generate longer acyl chains, which can subsequently be used as a precursor to generate fatty acids and phospholipids. Of particular relevance for adaptive immunity is the requirement of metabolic reprogramming in T cell activation, and fatty acid metabolism involvement in controlling the balance between an effector and a regulatory phenotype (55). Tregs require more active glycolysis, mitochondrial respiration, and fatty acid oxidation, compared to effector T cells indicated by elevated expression levels of FASN amongst other genes encoding metabolic enzymes (56). Specifically, it has been found that de novo fatty-acid synthesis mediated by FASN determines Treg maturation and if FASN is deleted from Tregs, tumor growth is inhibited (57). Indeed, in vivo inhibition of FASN in Th17 cells leads to reduction of autoimmune disease in mice (58), perhaps reflecting the crucial balance between Tregs and Th17 in preventing autoimmunity. We found an upregulation of FASN, indicating that Tregs of APS-1 patients may have more active fatty acid metabolism than counterparts from healthy controls. However, further functional or metabolic studies are needed to confirm this connection.

BCAR1 and SKIDA1 were identified as differentially expressed in our study, although the absolute expression of these were below 1 TPM. Read distributions were consistently mapped to exons for BCAR1 however, which may indicate that while low, the expression is reliable. SKIDA1 on the other hand had no consistent read distributions between samples and few reads conforming to exons suggesting that this may be a spurious finding and its biologically relevance questionable. BCAR1 belongs to the CAS (Crk-associated substrate) adaptor protein family and is involved in Rap1 and chemokine signaling pathways. The protein is expressed in nearly all kinds of cells and serves as an important mediator for several infectious diseases although not much is known about its role in tolerance and autoimmunity (59). TLN2, being downregulated in APS-1s’ Tregs, is also involved in Rap1 signaling, but little is known about the biological function of this protein. Its paralogue, TLN1, encodes Talin-1 which is a cytoskeletal protein essential in mediating integrin activation. Talin-1 has been reported to have a role in Treg cell–mediated maintenance of immune homeostasis since T cell–specific deletion of Talin-1 in mice yielded lymphocyte activation due to Tregs deficiency (60, 61).

When only considering patients with the common 13bp deletion in exon 8 of AIRE, the genes TMEM39B, AKR1C, GPSM3 and several members of the HLA-DR/DQ gene family emerged as differentially expressed compared to healthy controls. While it reduces the power of the differential expression analysis to include only a sub-cohort of patients, this is the largest subgroup in this study (N=5/9), and the 13bp deletion is also the most common Norwegian APS-1 mutation (20, 62). Although driven by a subgroup of healthy controls with high HLA-DR/DQ-expression, significantly lower levels of these genes were obvious in our analysis in patients with the 13bp AIRE mutation. HLA-DR/DQ-molecules are often expressed at high levels in Tregs and are considered as general activation markers for T cells. In fact, human CD4⁺CD25^hi T cells expressing HLA-DR/DQ have long been recognized as a functionally distinct population of Tregs that induces early contact-dependent suppression associated with high FOXP3 expression (63). Low levels of these could potentially indicate that the patients’ Tregs are less active than healthy controls, but functional studies are needed to confirm these findings. In addition to TMEM39B, the two final genes that were differentially expressed between APS-1 patients with the 13bp AIRE deletion and controls were GPSM3 and AKR1C3. AKR1C3, also known as 17β-hydroxysteroid dehydrogenase type 5 is a steroidogenic enzyme (17β-HSD5) capable of converting prostaglandins to metabolizing estrogen and progesterone. GPSM3 in turn has a selective immune cell expression distribution and SNPs in this gene have been associated with autoimmune diseases. GPSM3 is upregulated when the level of myeloid-derived suppressor cells (MDSCs) expands, which again is correlated to increased levels of autoimmunity or inflammation. These MDSCs can also promote Treg development (64).

For our transcriptomic analyses, we have isolated Tregs based on the markers CD4⁺CD25^hiCD127^- which should correspond to CD4⁺CD25^hiFOXP3⁺ (65). Notably, even though APS-1 patients have fewer Tregs than control subjects in their blood, we have here included equal cell counts for subsequent analysis. Hence, it is not surprising that the classic five Tregs-markers were similar between patients and controls in our study. It has been previously reported that FOXP3 expression levels are reduced in Tregs from APS-1 patients compared to healthy controls (16, 17), however we did not find a significant difference in FOXP3 expression between patients and controls in this study. This may be explained by the lack of power due to a low number of patients included here or a selection bias as CD127^-, which is correlated to FOXP3+ expression, was part of the cell sorting criteria for the RNAseq analysis. Nevertheless, sorting on CD4⁺D25^hiCD127^- cells defines the nTregs coming from the thymus where AIRE has its most profound role, which is the reason why we chose to study this cell cohort. We do recognize that this leads to neglection of other blood borne suppressor cells which might have notable contributions to the total T suppressor cell pool as even patients with IPEX (hallmarked by FOXP3 deleterious mutations) have heterogeneous Treg-like cells in their blood (66).

In conclusion, we report here the first detailed transcriptomic analysis of Tregs from APS-1 patients. Although differences between Tregs of APS-1 patients and healthy controls were relatively subtle, we believe our findings provide important and additional clues to the dysfunctional Treg responses previously reported for APS-1 patients. Specifically, our findings implicate deficient gut homing, increased autophagy and cellular turnover, and increased metabolic activity in Tregs of APS-1 patients compared to healthy controls. Additional functional or metabolic studies are however needed to confirm the consequences of these transcriptional dysfunctions.

Data Availability Statement

The datasets presented in this article are not readily available because of restrictions related to ethical matters regarding sharing of human transcriptomic data. Requests to access the datasets should be directed to the corresponding author (QW5ldHRlLmJvZUB1aWIubm8=). Access to the dataset requires an ethical approval from the requesting party.

Ethics Statement

The studies involving human participants were reviewed and approved by the Regional Ethics Committee of Western Norway (approval numbers 2009/2555 and 2018/1417). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

AB, EB, SJ, EH, BO, and AW conceived and designed the study. ØB and EH recruited patients. AB, EB, TS, MH, TT, ØB, SJ, EH, BO, and AW acquired, processed, analyzed, and interpreted the data. AB, EB, SJ, EH, and AW drafted the manuscript. All authors critically revised the manuscript for important intellectual content. All authors contributed to the article and approved the submitted version.

Funding

KG Jebsen Center for Autoimmune Disorders, Western Norway Health Authorities, and the Novo Nordisk Foundation (grant NNF17OC0027492).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The technical assistance from Elisabeth Halvorsen, Elin Theodorsen and Hajirah Muneer is greatly appreciated. We will further thank the core facility for genomics, University of Bergen v/Rita Holdhus for excellent service. The study would not be possible without contributions from patients and physicians of the Norwegian Patient Registry for organ specific autoimmune disorders (ROAS).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.722860/full#supplementary-material

References

1. Polansky JK, Kretschmer K, Freyer J, Floess S, Garbe A, Baron U, et al. DNA Methylation Controls Foxp3 Gene Expression. Eur J Immunol (2008) 38(6):1654–63. doi: 10.1002/eji.200838105

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sakaguchi S. Naturally Arising Foxp3-Expressing CD25+CD4+ Regulatory T Cells in Immunological Tolerance to Self and Non-Self. Nat Immunol (2005) 6(4):345–52. doi: 10.1038/ni1178

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Sakaguchi S, Sakaguchi N, Asano M, Itoh M, Toda M. Immunologic Self-Tolerance Maintained by Activated T Cells Expressing IL-2 Receptor Alpha-Chains (CD25). Breakdown of a Single Mechanism of Self-Tolerance Causes Various Autoimmune Diseases. J Immunol (1995) 155(3):1151–64.

PubMed Abstract | Google Scholar

4. Zhang X, Chang Li X, Xiao X, Sun R, Tian Z, Wei H. CD4(+)CD62L(+) Central Memory T Cells can be Converted to Foxp3(+) T Cells. PLoS One (2013) 8(10):e77322. doi: 10.1371/journal.pone.0077322

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Cowan JE, Baik S, McCarthy NI, Parnell SM, White AJ, Jenkinson WE, et al. Aire Controls the Recirculation of Murine Foxp3(+) Regulatory T-Cells Back to the Thymus. Eur J Immunol (2018) 48(5):844–54. doi: 10.1002/eji.201747375

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Kim BS, Lu H, Ichiyama K, Chen X, Zhang YB, Mistry NA, et al. Generation of RORgammat(+) Antigen-Specific T Regulatory 17 Cells From Foxp3(+) Precursors in Autoimmunity. Cell Rep (2017) 21(1):195–207. doi: 10.1016/j.celrep.2017.09.021

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Levine AG, Mendoza A, Hemmers S, Moltedo B, Niec RE, Schizas M, et al. Stability and Function of Regulatory T Cells Expressing the Transcription Factor T-Bet. Nature (2017) 546(7658):421–5. doi: 10.1038/nature22360

PubMed Abstract | CrossRef Full Text | Google Scholar

8. The Finnish-German APECED Consortium. An Autoimmune Disease, APECED, Caused by Mutations in a Novel Gene Featuring Two PHD-Type Zinc-Finger Domains. Nat Genet (1997) 17(4):399–403. Autoimmune Polyendocrinopathy-Candidiasis-Ectodermal Dystrophy. doi: 10.1038/ng1297-399

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Husebye ES, Anderson MS, Kampe O. Autoimmune Polyendocrine Syndromes. N Engl J Med (2018) 378(12):1132–41. doi: 10.1056/NEJMra1713301

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Nagamine K, Peterson P, Scott HS, Kudoh J, Minoshima S, Heino M, et al. Positional Cloning of the APECED Gene. Nat Genet (1997) 17(4):393–8. doi: 10.1038/ng1297-393

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Anderson MS, Venanzi ES, Klein L, Chen Z, Berzins SP, Turley SJ, et al. Projection of an Immunological Self Shadow Within the Thymus by the Aire Protein. Science (2002) 298(5597):1395–401. doi: 10.1126/science.1075958

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Aschenbrenner K, D’Cruz LM, Vollmann EH, Hinterberger M, Emmerich J, Swee LK, et al. Selection of Foxp3+ Regulatory T Cells Specific for Self Antigen Expressed and Presented by Aire+ Medullary Thymic Epithelial Cells. Nat Immunol (2007) 8(4):351–8. doi: 10.1038/ni1444

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Yang S, Fujikado N, Kolodin D, Benoist C, Mathis D. Immune Tolerance. Regulatory T Cells Generated Early in Life Play a Distinct Role in Maintaining Self-Tolerance. Science (2015) 348(6234):589–94. doi: 10.1126/science.aaa7017

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Malchow S, Leventhal DS, Lee V, Nishi S, Socci ND, Savage PA. Aire Enforces Immune Tolerance by Directing Autoreactive T Cells Into the Regulatory T Cell Lineage. Immunity (2016) 44(5):1102–13. doi: 10.1016/j.immuni.2016.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kekalainen E, Tuovinen H, Joensuu J, Gylling M, Franssila R, Pontynen N, et al. A Defect of Regulatory T Cells in Patients With Autoimmune Polyendocrinopathy-Candidiasis-Ectodermal Dystrophy. J Immunol (2007) 178(2):1208–15. doi: 10.4049/jimmunol.178.2.1208

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Laakso SM, Laurinolli TT, Rossi LH, Lehtoviita A, Sairanen H, Perheentupa J, et al. Regulatory T Cell Defect in APECED Patients Is Associated With Loss of Naive FOXP3(+) Precursors and Impaired Activated Population. J Autoimmun (2010) 35(4):351–7. doi: 10.1016/j.jaut.2010.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wolff AS, Oftedal BE, Kisand K, Ersvaer E, Lima K, Husebye ES. Flow Cytometry Study of Blood Cell Subtypes Reflects Autoimmune and Inflammatory Processes in Autoimmune Polyendocrine Syndrome Type I. Scand J Immunol (2010) 71(6):459–67. doi: 10.1111/j.1365-3083.2010.02397.x

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Bennett CL, Christie J, Ramsdell F, Brunkow ME, Ferguson PJ, Whitesell L, et al. The Immune Dysregulation, Polyendocrinopathy, Enteropathy, X-Linked Syndrome (IPEX) Is Caused by Mutations of FOXP3. Nat Genet (2001) 27(1):20–1. doi: 10.1038/83713

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ohkura N, Yasumizu Y, Kitagawa Y, Tanaka A, Nakamura Y, Motooka D, et al. Regulatory T Cell-Specific Epigenomic Region Variants Are a Key Determinant of Susceptibility to Common Autoimmune Diseases. Immunity (2020) 52(6):1119–32.e4. doi: 10.1016/j.immuni.2020.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Bruserud O, Oftedal BE, Landegren N, Erichsen MM, Bratland E, Lima K, et al. A Longitudinal Follow-Up of Autoimmune Polyendocrine Syndrome Type 1. J Clin Endocrinol Metab (2016) 101(8):2975–83. doi: 10.1210/jc.2016-1821

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Santegoets SJ, Dijkgraaf EM, Battaglia A, Beckhove P, Britten CM, Gallimore A, et al. Monitoring Regulatory T Cells in Clinical Samples: Consensus on an Essential Marker Set and Gating Strategy for Regulatory T Cell Analysis by Flow Cytometry. Cancer Immunol Immunother (2015) 64(10):1271–86. doi: 10.1007/s00262-015-1729-x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. FastQC. A Quality Control Tool for High Throughput Sequence Data (2015). Available at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.

Google Scholar

23. Bray NL, Pimentel H, Melsted P, Pachter L. Near-Optimal Probabilistic RNA-Seq Quantification. Nat Biotechnol (2016) 34(5):525–7. doi: 10.1038/nbt.3519

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Soneson C, Love MI, Robinson MD. Differential Analyses for RNA-Seq: Transcript-Level Estimates Improve Gene-Level Inferences. F1000Res (2015) 4:1521. doi: 10.12688/f1000research.7563.1

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Durinck S, Spellman PT, Birney E, Huber W. Mapping Identifiers for the Integration of Genomic Datasets With the R/Bioconductor Package biomaRt. Nat Protoc (2009) 4(8):1184–91. doi: 10.1038/nprot.2009.97

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Love MI, Huber W, Anders S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2. Genome Biol (2014) 15(12):550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Ignatiadis N, Klaus B, Zaugg JB, Huber W. Data-Driven Hypothesis Weighting Increases Detection Power in Genome-Scale Multiple Testing. Nat Methods (2016) 13(7):577–80. doi: 10.1038/nmeth.3885

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Zhu A, Ibrahim JG, Love MI. Heavy-Tailed Prior Distributions for Sequence Count Data: Removing the Noise and Preserving Large Differences. Bioinformatics (2019) 35(12):2084–92. doi: 10.1093/bioinformatics/bty895

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics (2013) 29(1):15–21. doi: 10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Tarasov A, Vilella AJ, Cuppen E, Nijman IJ, Prins P. Sambamba: Fast Processing of NGS Alignment Formats. Bioinformatics (2015) 31(12):2032–4. doi: 10.1093/bioinformatics/btv098

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Team R. Studio: Integrated Development for R (2020). Available at: http://www.rstudio.com/.

Google Scholar

32. Wickham H, Averick M, Bryan J, Chang W, D’Agostino McGowan L, Francois R, et al. Welcome to the Tidyverse Vol. 4. Journal of Open Source Software (2019). p. 1686. doi: 10.21105/joss.01686

CrossRef Full Text | Google Scholar

33. Garnier S, Ross N, Rudis R, Camargo PA, Sciaini M, Scherer C. Viridis - Colorblind-Friendly Color Maps for R (2021). Available at: https://sjmgarnier.github.io/viridis/.

Google Scholar

34. Hahne F, Ivanek R. “Visualizing Genomic Data Using Gviz and Bioconductor”. In: Mathé E, Davis S, editors. Statistical Genomics Methods in Molecular Biology, vol. 1418. New York, NY: Humana Press (2016). doi: 10.1007/978-1-4939-3578-9_16

CrossRef Full Text | Google Scholar

35. Boe Wolff AS, Oftedal B, Johansson S, Bruland O, Lovas K, Meager A, et al. AIRE Variations in Addison’s Disease and Autoimmune Polyendocrine Syndromes (APS): Partial Gene Deletions Contribute to APS I. Genes Immun (2008) 9(2):130–6. doi: 10.1038/sj.gene.6364457

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Zemmour D, Zilionis R, Kiner E, Klein AM, Mathis D, Benoist C. Single-Cell Gene Expression Reveals a Landscape of Regulatory T Cell Phenotypes Shaped by the TCR. Nat Immunol (2018) 19(3):291–301. doi: 10.1038/s41590-018-0051-0

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Koivula TT, Laakso SM, Niemi HJ, Kekalainen E, Laine P, Paulin L, et al. Clonal Analysis of Regulatory T Cell Defect in Patients With Autoimmune Polyendocrine Syndrome Type 1 Suggests Intrathymic Impairment. Scand J Immunol (2017) 86(4):221–8. doi: 10.1111/sji.12586

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Yano M, Kuroda N, Han H, Meguro-Horike M, Nishikawa Y, Kiyonari H, et al. Aire Controls the Differentiation Program of Thymic Epithelial Cells in the Medulla for the Establishment of Self-Tolerance. J Exp Med (2008) 205(12):2827–38. doi: 10.1084/jem.20080046

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Halonen M, Kangas H, Ruppell T, Ilmarinen T, Ollila J, Kolmer M, et al. APECED-Causing Mutations in AIRE Reveal the Functional Domains of the Protein. Hum Mutat (2004) 23(3):245–57. doi: 10.1002/humu.20003

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Oftedal BE, Hellesen A, Erichsen MM, Bratland E, Vardi A, Perheentupa J, et al. Dominant Mutations in the Autoimmune Regulator AIRE Are Associated With Common Organ-Specific Autoimmune Diseases. Immunity (2015) 42(6):1185–96. doi: 10.1016/j.immuni.2015.04.021

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Ferraro A, D’Alise AM, Raj T, Asinovski N, Phillips R, Ergun A, et al. Interindividual Variation in Human T Regulatory Cells. Proc Natl Acad Sci USA (2014) 111(12):E1111–20. doi: 10.1073/pnas.1401343111

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Woo HJ, Yu C, Reifman J. Collective Genetic Interaction Effects and the Role of Antigen-Presenting Cells in Autoimmune Diseases. PLoS One (2017) 12(1):e0169918. doi: 10.1371/journal.pone.0169918

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kim SV, Xiang WV, Kwak C, Yang Y, Lin XW, Ota M, et al. GPR15-Mediated Homing Controls Immune Homeostasis in the Large Intestine Mucosa. Science (2013) 340(6139):1456–9. doi: 10.1126/science.1237013

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Xiong L, Dean JW, Fu Z, Oliff KN, Bostick JW, Ye J, et al. Ahr-Foxp3-RORgammat Axis Controls Gut Homing of CD4(+) T Cells by Regulating GPR15. Sci Immunol (2020) 5(48):eaaz7277. doi: 10.1126/sciimmunol.aaz7277

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Cartwright A, Schmutz C, Askari A, Kuiper JH, Middleton J. Orphan Receptor GPR15/BOB Is Up-Regulated in Rheumatoid Arthritis. Cytokine (2014) 67(2):53–9. doi: 10.1016/j.cyto.2014.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Schmutz C, Hulme A, Burman A, Salmon M, Ashton B, Buckley C, et al. Chemokine Receptors in the Rheumatoid Synovium: Upregulation of CXCR5. Arthritis Res Ther (2005) 7(2):R217–29. doi: 10.1186/ar1475

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Ammitzboll C, von Essen MR, Bornsen L, Petersen ER, McWilliam O, Ratzer R, et al. GPR15(+) T Cells Are Th17 Like, Increased in Smokers and Associated With Multiple Sclerosis. J Autoimmun (2019) 97:114–21. doi: 10.1016/j.jaut.2018.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Adamczyk A, Gageik D, Frede A, Pastille E, Hansen W, Rueffer A, et al. Differential Expression of GPR15 on T Cells During Ulcerative Colitis. JCI Insight (2017) 2(8):e90585. doi: 10.1172/jci.insight.90585

CrossRef Full Text | Google Scholar

49. Adamczyk A, Pastille E, Kehrmann J, Vu VP, Geffers R, Wasmer MH, et al. GPR15 Facilitates Recruitment of Regulatory T Cells to Promote Colorectal Cancer. Cancer Res (2021) 81(11):2970–82. doi: 10.1158/0008-5472.CAN-20-2133

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Kluger N, Jokinen M, Krohn K, Ranki A. Gastrointestinal Manifestations in APECED Syndrome. J Clin Gastroenterol (2013) 47(2):112–20. doi: 10.1097/MCG.0b013e31827356e1

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Hetemaki I, Jarva H, Kluger N, Baldauf HM, Laakso S, Bratland E, et al. Anticommensal Responses Are Associated With Regulatory T Cell Defect in Autoimmune Polyendocrinopathy-Candidiasis-Ectodermal Dystrophy Patients. J Immunol (2016) 196(7):2955–64. doi: 10.4049/jimmunol.1500301

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Orvedahl A, Sumpter R Jr, Xiao G, Ng A, Zou Z, Tang Y, et al. Image-Based Genome-Wide siRNA Screen Identifies Selective Autophagy Factors. Nature (2011) 480(7375):113–7. doi: 10.1038/nature10546

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Wei J, Long L, Yang K, Guy C, Shrestha S, Chen Z, et al. Autophagy Enforces Functional Integrity of Regulatory T Cells by Coupling Environmental Cues and Metabolic Homeostasis. Nat Immunol (2016) 17(3):277–85. doi: 10.1038/ni.3365

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Kim SJ, Sohn I, Do IG, Jung SH, Ko YH, Yoo HY, et al. Gene Expression Profiles for the Prediction of Progression-Free Survival in Diffuse Large B Cell Lymphoma: Results of a DASL Assay. Ann Hematol (2014) 93(3):437–47. doi: 10.1007/s00277-013-1884-0

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Cluxton D, Petrasca A, Moran B, Fletcher JM. Differential Regulation of Human Treg and Th17 Cells by Fatty Acid Synthesis and Glycolysis. Front Immunol (2019) 10:115. doi: 10.3389/fimmu.2019.00115

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Almeida L, Lochner M, Berod L, Sparwasser T. Metabolic Pathways in T Cell Activation and Lineage Differentiation. Semin Immunol (2016) 28(5):514–24. doi: 10.1016/j.smim.2016.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Lim SA, Wei J, Nguyen TM, Shi H, Su W, Palacios G, et al. Lipid Signalling Enforces Functional Specialization of Treg Cells in Tumours. Nature (2021) 591(7849):306–11. doi: 10.1038/s41586-021-03235-6

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Young KE, Flaherty S, Woodman KM, Sharma-Walia N, Reynolds JM. Fatty Acid Synthase Regulates the Pathogenicity of Th17 Cells. J Leukoc Biol (2017) 102(5):1229–35. doi: 10.1189/jlb.3AB0417-159RR

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Deneka A, Korobeynikov V, Golemis EA. Embryonal Fyn-Associated Substrate (EFS) and CASS4: The Lesser-Known CAS Protein Family Members. Gene (2015) 570(1):25–35. doi: 10.1016/j.gene.2015.06.062

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Klann JE, Remedios KA, Kim SH, Metz PJ, Lopez J, Mack LA, et al. Talin Plays a Critical Role in the Maintenance of the Regulatory T Cell Pool. J Immunol (2017) 198(12):4639–51. doi: 10.4049/jimmunol.1601165

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Klann JE, Kim SH, Remedios KA, He Z, Metz PJ, Lopez J, et al. Integrin Activation Controls Regulatory T Cell-Mediated Peripheral Tolerance. J Immunol (2018) 200(12):4012–23. doi: 10.4049/jimmunol.1800112

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Wolff AS, Erichsen MM, Meager A, Magitta NF, Myhre AG, Bollerslev J, et al. Autoimmune Polyendocrine Syndrome Type 1 in Norway: Phenotypic Variation, Autoantibodies, and Novel Mutations in the Autoimmune Regulator Gene. J Clin Endocrinol Metab (2007) 92(2):595–603. doi: 10.1210/jc.2006-1873

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Baecher-Allan C, Wolf E, Hafler DA. MHC Class II Expression Identifies Functionally Distinct Human Regulatory T Cells. J Immunol (2006) 176(8):4622–31. doi: 10.4049/jimmunol.176.8.4622

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Billard MJ, Gall BJ, Richards KL, Siderovski DP, Tarrant TK. G Protein Signaling Modulator-3: A Leukocyte Regulator of Inflammation in Health and Disease. Am J Clin Exp Immunol (2014) 3(2):97–106.

PubMed Abstract | Google Scholar

65. Liu W, Putnam AL, Xu-Yu Z, Szot GL, Lee MR, Zhu S, et al. CD127 Expression Inversely Correlates With FoxP3 and Suppressive Function of Human CD4+ T Reg Cells. J Exp Med (2006) 203(7):1701–11. doi: 10.1084/jem.20060772

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Zemmour D, Charbonnier LM, Leon J, Six E, Keles S, Delville M, et al. Single-Cell Analysis of FOXP3 Deficiencies in Humans and Mice Unmasks Intrinsic and Extrinsic CD4(+) T Cell Perturbations. Nat Immunol (2021) 22(5):607–19. doi: 10.1038/s41590-021-00910-8

PubMed Abstract | CrossRef Full Text | Google Scholar