The development of inflammatory diseases depends on complex interactions between several genes and various environmental factors. Discovering new genetic risk factors and understanding the mechanisms whereby they influence disease development is of paramount importance. We previously reported that deficiency in Themis1, a new actor of TCR signaling, impairs regulatory T cell (Treg) function and predisposes Brown–Norway (BN) rats to spontaneous inflammatory bowel disease (IBD). In this study, we reveal that the epistasis between Themis1 and Vav1 controls the occurrence of these phenotypes. Indeed, by contrast with BN rats, Themis1 deficiency in Lewis rats neither impairs Treg suppressive functions nor induces pathological manifestations. By using congenic lines on the BN genomic background, we show that the impact of Themis1 deficiency on Treg suppressive functions depends on a 117-kb interval coding for a R63W polymorphism that impacts Vav1 expression and functions. Indeed, the introduction of a 117-kb interval containing the Lewis Vav1-R63 variant restores Treg function and protects Themis1-deficient BN rats from spontaneous IBD development. We further show that Themis1 binds more efficiently to the BN Vav1-W63 variant and is required to stabilize its recruitment to the transmembrane adaptor LAT and to fully promote the activation of Erk kinases. Together, these results highlight the importance of the signaling pathway involving epistasis between Themis1 and Vav1 in the control of Treg suppressive function and susceptibility to IBD development.

Chronic inflammatory disorders are multifactorial diseases whose development depends on a combination of environmental triggers and genetic risk factors (1). To date, most genetic studies performed in human populations have used a single-locus analysis strategy, allowing the association of gene variants with a specific phenotype. This has led to the identification of a great number of gene variants, most of them conferring only a small risk increment. In most cases, the identified gene variants explain <20% of disease heritability (2). One potential source of this low heritability rate may result from interactions between polymorphic genes, also known as epistasis, that are not taken into account by the current methods of analysis (24). Indeed, many susceptibility genes may exhibit effects that are totally or partially dependent on interactions with other genes (5). This highlights the importance of integrating gene–gene interactions in the genetic studies. In this regard, mouse and rat models of human diseases represent invaluable tools not only for the identification of gene–gene interactions involved in susceptibility to immune-mediated diseases (6) but also for deciphering the mechanisms whereby these interactions affect the disease pathophysiology.

Signaling through the TCR directs the functional outcome of T cell activation. Thus, because TCR signaling integrity controls T cell functions, any polymorphism or deficiency for TCR signaling genes could predispose to immune-mediated diseases. Indeed, our previous genetic studies using animal models have independently identified the role of the TCR signaling genes Vav1 (79) and Themis1 (10) in the susceptibility to autoimmune and inflammatory diseases. Vav1 is a guanine nucleotide exchange factor (GEF) for Rho GTPases, which plays a nonredundant role in transducing TCR signals important for T cell development and functions (1114). Vav1 phosphorylation relieves its catalytic domain and causes Vav1 to promote GDP-GTP exchange on the Rho, Rac1, and Cdc42 small GTPases. Besides this well-documented GEF activity, Vav1 also plays a role as an adaptor protein. Indeed, the generation of mice that express a mutant Vav1 protein lacking GEF activity clearly demonstrated that Vav1 enzymatic activity is not required for TCR-induced calcium flux and activation of the ERK MAPK pathway (15). Our previous genetic studies in Lewis (LEW) and Brown–Norway (BN) rat strains identified a major locus of 117-kb on rat chromosome 9 that controls the opposite susceptibility of those two strains to Th1 and Th2 immune-mediated diseases (7, 9, 16). Fine mapping of this locus identified a nonsynonymous exonic single nucleotide polymorphism in Vav1 leading to the substitution of an arginine by a tryptophan (R63W) in BN rats (7). Functional studies showed that the BN Vav1-W63 variant is constitutively active and exhibits normal GEF activity but reduced adaptor functions (7). Vav1-W63 variant was associated with a reduction of inflammatory cytokine production by T cells and decreased susceptibility to CNS inflammation (8).

Themis1 has been recently described as a new component of the TCR signaling machinery that controls thymic development of T cells. The ablation of Themis1 in C57BL/6 mice results in impaired T cell development, particularly of the CD4+ lineage (1721). After TCR engagement, Themis1 is phosphorylated by Lck and Zap-70, and recruited to the transmembrane adaptor protein LAT (22, 23). Themis1 has been shown to interact with Vav1 and positively regulates its phosphorylation (23). However, it remains unclear whether Themis1 acts as a positive or negative regulator of TCR signaling (17, 2325). We previously discovered a spontaneous mutation in the BN rat strain leading to the disruption of the Themis1 gene (10). This mutation leads to a defect in thymic development of CD4 T cells similar to that observed in Themis1-deficient C57BL/6 mice. Strikingly, the Themis1-deficient BN rats develop spontaneously inflammatory bowel disease (IBD) (10). The functional studies revealed that Themis1-deficient BN regulatory T cells (Tregs) migrate normally into the intestine but exhibit a defect in their suppressive functions. In contrast, no Treg dysfunction or spontaneous pathological manifestations have been reported in mouse models of Themis1 deficiency. We investigated in this study whether these differences reflect an epistasis phenomenon on the BN genetic background or species discrepancies.

We show in this article that the impact of Themis1 deficiency on Treg suppressive function and IBD development depends on the genetic background. In contrast with the BN rat strain, the disruption of the Themis1 gene in LEW rats does not lead to the spontaneous development of IBD or functional Treg defect, suggesting that an epistasis phenomenon in the BN genetic background controls these phenotypes. In addition, we show that the gene(s) responsible for this epistasis is (are) contained within the previously described 117-kb locus on rat chromosome 9 that controls susceptibility to immune-mediated diseases in LEW and BN rats (7, 9). Finally, a biochemical analysis of TCR signaling pointed out an epistasis between Themis1 and Vav1-W63 variant in the BN rats.

All breeding and experimental procedures were carried out in accordance with the European Union guidelines and were approved by our local ethical committee (license number: 31259). Rats were maintained under specific pathogen-free conditions in our facility. Specific pathogen-free BN and Lewis (LEW) rats were obtained from the Centre d’Elevage R. Janvier (Le Genest St. Isle, France). The BNThemis1−/− rats were bred from an autosomal recessive mutation that occurred spontaneously in a BN rat colony from our laboratory (10). The generation of the LEWThemis1−/− rat strain was performed according to the speed congenics strategy (26) as described previously (27). In brief, repeated backcrosses from (LEW×BNThemis1−/−) F1 hybrids to LEW rats were performed. In rats issued from the backcrosses, genomic composition was assessed using 200 polymorphic microsatellite markers selected to cover the whole genome with ∼10-cM intervals. At each backcross, selected recombinants carrying both the higher percentage of LEW genomic origin and the genomic interval from BN origin containing the mutated Themis1 gene were backcrossed to LEW rats. At the end of four successive backcrosses, all microsatellite markers showed homozygosis for LEW genome, but three markers on top of chromosome 1 extending on a 21-Mb interval, in the region bearing the Themis1 gene. Recombinants rats were then intercrossed to generate the LEWThemis1−/− rats that were homozygous for the Themis1 knockdown mutation from BN origin. The absence of Themis1 protein expression was confirmed by Western blot analysis (Supplemental Fig. 1A, 1B). The BN.LEWc9-Bf congenic line (Bf) was generated as described previously (7, 28). The BfThemis1−/− line was raised by intercrossing (BNThemis−/−×Bf) F1 rats and selecting recombinants carrying homozygosis for both the Themis1 mutation and the 117-kb interval from LEW origin that characterizes the Bf line. In this BfThemis1−/− line, Western blotting for Themis1 and Vav1 confirmed the absence of Themis1 and the increase in Vav1 protein amount that characterizes the LEW Vav1-R63 variant as compared with the BN Vav1-WR63 variant (Supplemental Fig. 1C, 1D). Themis1-deficient mice were kindly provided by Dr. Paul E. Love (National Institutes of Health, Bethesda, MD) and have been reported previously (17).

The mAbs used for flow cytometry and purification were W3/25 (anti-rat CD4), OX6 (anti-rat MHC class II), OX33 (anti-rat B220), OX8 (anti-rat CD8α), R73 (anti-rat TCR αβ), V65 (anti-rat TCR γδ), JJ319 (anti-rat CD28), 341 (anti-rat CD8β), 3.2.3 (anti-rat NKR-P1), OX39 (anti-rat CD25), OX40 (anti-rat CD134), FJK-16s (anti-Foxp3), and HRL1 (anti-CD62L). The fluorescent-conjugated Abs were purchased from eBioscience, BD Biosciences, and Biolegend. The Abs used for biochemical studies were anti-Vav1 (C-14), anti–P-Vav1 Tyr174, anti-PLCγ1 (1249) from Santa Cruz, anti-Themis1 from Millipore, anti-LAT (1D-1) from Thermo Scientific, anti-Flag Ab (M2) from Sigma, anti-Akt, anti–P-Akt Ser473 (D9E), anti–P-ERK1/2 Thr202/Tyr204 (D13.14.4E), anti-ERK (3A7), anti–Myc-Tag Ab (71D10), HRP-linked anti-rabbit IgG, HRP-linked anti-mouse IgG, and HRP-linked anti-rabbit IgG (conformation specific) (L27A9) from Cell Signaling.

Spleen cells were teased apart in RPMI 1640 (Sigma-Aldrich, St. Louis, MO). The erythrocytes were lysed with NH4Cl-type ACK buffer and the cells were washed twice. CD4 T cells were negatively selected from spleen cells using anti-mouse IgG magnetic microbeads (Dynabeads, Life Technologies). In brief, cells were washed and incubated 20 min on ice with a mixture of the following mAbs: OX8, OX6, OX33, 3.2.3, and V65. After washing and incubating with anti-mouse IgG-coupled microbeads under agitation, CD4 T cells were purified by magnetic depletion. TCR+ CD4+ CD8 CD62L+ CD25 and TCR+ CD4+ CD8 CD25bright populations were purified using a FACSAria II-Sorp (BD Biosciences) cell sorter. Foxp3 expression was assessed using the eBioscience Fix-Perm Kit and anti-Foxp3 Ab. Purity of purified populations was tested before culture.

The culture medium was RPMI 1640 (Life Technologies Life Technologies, Cergy Pontoise, France) containing 10% FCS, 1% pyruvate, 1% nonessential amino acids, 1% l-glutamine, 1% penicillin-streptomycin, and 2 × 10−5 M 2-ME. For cytokine production assays, total CD4+ T cells or naive CD4+ CD62L+ CD25 T cells (105/well) were purified and stimulated with plate-bound anti-TCR (R73, 1 μg/ml) and soluble anti-CD28 (JJ319, 0.2 μg/ml) for 48 h. Cytokine production was examined in the cell supernatants by Luminex multiplex kit (Millipore). Induced Tregs were generated by stimulating BN and BNThemis1−/− CD4+ CD62L+ CD25 T cells with anti-CD3/anti-CD28 beads in the presence of IL-2 (100 U/ml) and TGF-β (5 ng/ml). Four days later, the CD4+ CD25bright were sorted using a FACSAria II-Sorp (BD Biosciences) cell sorter to test their suppressive functions in vitro. Coculture experiments were performed using naive CD4+ T cells (105/well) stained with CellTrace violet (CellTrace, Invitrogen) stimulated with irradiated syngeneic APC (0.5 × 105/well; T cell–depleted splenocytes) and plate-bound anti-CD3 mAb (G4.18; 0.05 μg/ml) in the presence of various amounts of natural Tregs purified from the spleen or induced Tregs generated in vitro. CellTrace violet dilution was assessed after 60 h of culture.

Salivary glands, stomach, kidneys, pancreas, skin, liver, pancreas, lungs, heart, brain, ovaries, testis, thyroid, intestines, and muscle were dissected and fixed for 24 h in 4% paraformaldehyde, then 24 h in ethanol, before embedding in paraffin and processing for conventional H&E staining. Those organs were blindly analyzed in male and female LEW and LEWThemis1−/− rats at the ages of 12–15 and 32–37 wk old. Intestinal lesions were analyzed in BN, BNThemis1−/−, LEW, LEWThemis1−/− Bf, and BfThemis1−/− rats and scored blindly by a pathologist using Wallace score. We also analyzed myeloperoxidase (MPO) activity as an index of granulocyte infiltration. In brief, tissue samples were homogenized in a solution of 0.5% hexadecyltrimethylammonium bromide (Sigma) with Precellys beads (Ozyme). Supernatants were added to O-dianisidine dihydrochloride solution supplemented with 1% hydrogen peroxide (Sigma). OD was read at 450 nm.

Themis1/Vav1 protein interaction was analyzed on lysates from HEK293 cells 24 h after transfection with C-myc tagged human wild type Vav1 or W63 mutated Vav1 in combination with Flag-tagged human Themis1. Transfections were performed using Effectene transfection reagent (Qiagen), and cells were stimulated with pervanadate. Total cellular proteins were extracted with ice-cold lysis buffer containing a mixture of protease and phosphatase inhibitors (complete Mini, EDTA-free; Roche), 10 mM Tris HCl, 1% Triton X-100, 50 mM NaF, 1 mM Na3VO4, and 1 mM DTT. For immunoprecipitation, clarified homogenates were incubated overnight at 4°C with beads coated with anti-flag M2 Ab (Sigma). After washes, protein complexes were incubated in buffer (150 mM NaCl, 50 mM Tris pH 7.4) completed with 200 μg/ml flag peptide for competitive elution. Eluates were completed with Laemmli buffer for migration followed by Western blotting on Immobilon-P membranes (Millipore) with appropriate Abs (Flag Ab, clone M2 [Sigma]; Myc-Tag Ab, clone 71D10 [Cell Signaling]). Immunoreactive bands were detected by chemiluminescence with the SuperSignal detection system (Pierce Chemical, Rockford, IL). For LAT immunoprecipitation, thymocyte suspensions were stimulated using preformed complexes of biotinylated R73 (anti-TCR, 10 μg/10 × 106 cells), W325 (anti-CD4, 10 μg/10 × 106 cells) and streptavidin (10 μg/10 × 106 cells). Stimulation was stopped by the addition of twice-concentrated lysis buffer (100 mM Tris, pH 7.5, 270 mM NaCl, 1 mM EDTA, 20% glycerol, and 0.2% n-dodecyl-β-maltoside) containing protease and phosphatase inhibitors. After 10 min of incubation on ice, cell lysates were centrifuged at 20,000 × g for 5 min at 4°C. Immunoprecipitations were performed using the LAT-1D1 mAb (Thermo Scientific) and Prot-G Sepharose beads (Santa Cruz). After washes, proteins were eluted with Laemmli buffer and analyzed by SDS-PAGE followed by Western blotting on PVDF membranes (Immobilon). ECL Prime (Amersham) was used as revelation substrate. Signal intensity quantification was performed using ImageJ software (1.47v).

Data were expressed as mean ± SEM. Statistical analyses were carried out using the GraphPad Instat statistical package (GraphPad Software, La Jolla, CA). Results were analyzed by the Mann–Whitney U test. Results were considered statistically significant when the p value was <0.05: *p < 0.05, **p < 0.01, ***p < 0.001.

To investigate the impact of the genomic background on the pathophysiological consequences of Themis1 deficiency, we transferred the recessive autosomal Themis1 mutation from BNThemis1−/− rats to the LEW genomic background (Supplemental Fig. 1A) using the speed congenic strategy (27). The Themis1 disruption was confirmed by the absence of Themis1 protein expression in thymocytes (Supplemental Fig. 1B). We next investigated the thymic development of T cells in LEWThemis1−/− rats. Although there were no changes in CD4CD8 double negative or CD4+CD8+ double positive thymocytes, the proportion (Fig. 1A) and absolute numbers (Fig. 1B) of CD4+ single-positive cells were greatly reduced in LEWThemis1−/− rats, together with a milder decrease in CD8+ single-positive thymocytes. Accordingly, the proportion (Fig. 1C) and absolute numbers (Fig. 1D) of CD4 T cells were markedly reduced in spleen (Fig. 1C, 1D) and lymph nodes of LEWThemis1−/− rats (data not shown). These abnormalities affected to a lesser extent the CD8 T cell compartment (Fig. 1C, 1D). Cohorts of LEWThemis1−/− rats were followed for 38 wk to analyze whether Themis1 deficiency in the LEW genetic background is accompanied by the development of pathological disorders. We did not observe any weight loss or clinical signs of spontaneous pathological manifestations. Furthermore, blinded histological analysis of numerous organs (salivary glands, stomach, kidneys, pancreas, skin, liver, pancreas, lungs, heart, brain, ovaries, testis, thyroid, intestines, and muscle) collected from 12- or 38-wk-old rats did not reveal any inflammatory cell infiltration in LEWThemis1−/− rats (data not shown). Thus, Themis1 disruption in the LEW genomic background reproduces the cardinal CD4 T cell lymphopenia phenotype described in BNThemis1−/− rats and in C57BL/6Themis1−/− mice (1721). However, in contrast with what was observed in BN rats, this deficiency does not predispose to the development of spontaneous pathological diseases.

FIGURE 1.

Themis1 deficiency in LEW rats affects thymic development of T cells. (A) Representative dot plots of CD4 and CD8 expression in LEW and LEWThemis1−/− thymocytes. Numbers represent cell percentages in the boxed area. (B) Absolute numbers of LEW (n = 4) and LEWThemis1−/− (n = 4) thymocyte subsets. (C) Representative dot plots of CD4 and CD8 expression in LEW and LEWThemis1−/− spleen T cells. Numbers represent cell percentages in the boxed area. (D) Absolute numbers of splenocytes, T cells, CD4+ T cells, CD8+ T cells, and B cells. LEW, white columns; LEWThemis1−/−: black columns. Results are representative of one of three independent experiments, and data are presented as mean ± SEM. *p ≤ 0.05. DN, double negative; DP, double positive; ns, not significant; SP, single positive.

FIGURE 1.

Themis1 deficiency in LEW rats affects thymic development of T cells. (A) Representative dot plots of CD4 and CD8 expression in LEW and LEWThemis1−/− thymocytes. Numbers represent cell percentages in the boxed area. (B) Absolute numbers of LEW (n = 4) and LEWThemis1−/− (n = 4) thymocyte subsets. (C) Representative dot plots of CD4 and CD8 expression in LEW and LEWThemis1−/− spleen T cells. Numbers represent cell percentages in the boxed area. (D) Absolute numbers of splenocytes, T cells, CD4+ T cells, CD8+ T cells, and B cells. LEW, white columns; LEWThemis1−/−: black columns. Results are representative of one of three independent experiments, and data are presented as mean ± SEM. *p ≤ 0.05. DN, double negative; DP, double positive; ns, not significant; SP, single positive.

Close modal

We first analyzed whether Themis1 deficiency affects splenic CD4 T cell function in a LEW genetic background. Consistent with our findings in BNThemis1−/− rats (10), we observed an overexpression of OX40 and a downregulation of CD62L expression in LEWThemis1−/− CD4 T cells (Fig. 2A), indicating an increased proportion of activated CD4 T cells. The cytokine production was then measured in both total and naive (CD62L+ CD25) CD4 T cells after anti-TCR and anti-CD28 stimulation. As compared with CD4 T cells from LEW rats, both total (Fig. 2B) and naive CD4 T cells (Fig. 2C) from LEWThemis1−/− rats displayed increased production of IL-4, IL-10, and IL-17, thus confirming the skewed cytokine production previously associated with Themis1 deficiency in the BN genomic background. Of note, the IFN-γ production was not decreased by Themis1 deficiency in the LEW background, contrary to its effect in the BN background. Next, we investigated the impact of Themis1 deficiency on LEW Treg development and function. Although the proportion of Tregs among spleen CD4 T cells was increased in LEWThemis1−/− rats, the total number of Tregs was actually decreased, as a result of CD4 lymphopenia (Fig. 3A, 3B). We next studied Treg function in a coculture assay using CD4+ CD8 CD25bright T cells. Sorted Tregs from LEW and LEWThemis1−/− rats showed similar purity and expression levels of the Treg transcription factor Foxp3 (Fig. 3C, 3D). Interestingly, the Treg suppressive function was unaltered in LEWThemis1−/− rats (Fig. 3E, 3F). These observations concord with the absence of spontaneous disease because the development of IBD in BNThemis1−/− rats is linked to a functional defect of Tregs. Similarly, the suppressive function of Tregs purified from Themis1-deficient C57BL/6 mice was similar to that of Tregs originating from wild-type mice (Supplemental Fig. 2). Thus, the differential consequence of Themis1 deficiency observed between BNThemis1−/− rats and mice rather reflects an epistatic phenomenon occurring in the BN genetic background than species discrepancies.

FIGURE 2.

Themis1 deficiency in LEW rats affects effector CD4 T cell functions. (A) Proportions of OX40+ and CD62L+ cells among CD4+ T cells from spleens of LEW (n = 4; white columns) and LEWThemis1−/− (n = 4; black columns) rats. (B and C) Cytokine production by total CD4+ T cells (B) or naive CD62L+ CD25 CD4+ T cells (C) purified from LEW or LEWThemis1−/− rats and stimulated with anti-TCR and anti-CD28 mAb for 48 h. Results are representative of one of three independent experiments. Data are presented as mean ± SEM. *p ≤ 0.05. ns, not significant.

FIGURE 2.

Themis1 deficiency in LEW rats affects effector CD4 T cell functions. (A) Proportions of OX40+ and CD62L+ cells among CD4+ T cells from spleens of LEW (n = 4; white columns) and LEWThemis1−/− (n = 4; black columns) rats. (B and C) Cytokine production by total CD4+ T cells (B) or naive CD62L+ CD25 CD4+ T cells (C) purified from LEW or LEWThemis1−/− rats and stimulated with anti-TCR and anti-CD28 mAb for 48 h. Results are representative of one of three independent experiments. Data are presented as mean ± SEM. *p ≤ 0.05. ns, not significant.

Close modal
FIGURE 3.

Themis1 deficiency in LEW rats does not alter Treg suppressive functions. (A) Proportions of CD25high and Foxp3+ cells among CD4 T cells from spleens of LEW (n = 4) and LEWThemis1−/− (n = 4) rats. (B) Absolute numbers of Foxp3+ CD4+ T cells in spleen of LEW and LEWThemis1−/− rats. (C) Percentages of Foxp3+ cells and (D) Foxp3 expression levels among sorted CD25bright CD4+ CD8 TCR+ Tregs used in coculture experiments. Data are presented as mean ± SEM. (E and F) Suppressive activity of CD25bright CD4+ T cells from LEW and LEWThemis1−/− rats assessed in coculture experiments with CellTrace-labeled naive LEW CD4+ T cells as effector cells stimulated with plate-bound anti-CD3 and irradiated syngeneic APC. For each experiment, we pooled the purified Tregs from three rats per group because of the low numbers of Tregs present in Themis1-deficient rats. (E) Representative histograms of the proliferation of effector cells alone or in the presence of LEW and LEWThemis1−/− Tregs (Treg/Teff ratio 1:2). Percentages indicate the proportion of CellTrace low effector cells. Results are representative of one of three independent experiments. (F) Suppressive activity of CD25bright CD4+ T cells from LEW and LEWThemis1−/− rats at different Treg/Teff ratios. Data are expressed as percentage of inhibition and presented as mean values ± SEM obtained from three independent experiments. *p ≤ 0.05. ns, not significant.

FIGURE 3.

Themis1 deficiency in LEW rats does not alter Treg suppressive functions. (A) Proportions of CD25high and Foxp3+ cells among CD4 T cells from spleens of LEW (n = 4) and LEWThemis1−/− (n = 4) rats. (B) Absolute numbers of Foxp3+ CD4+ T cells in spleen of LEW and LEWThemis1−/− rats. (C) Percentages of Foxp3+ cells and (D) Foxp3 expression levels among sorted CD25bright CD4+ CD8 TCR+ Tregs used in coculture experiments. Data are presented as mean ± SEM. (E and F) Suppressive activity of CD25bright CD4+ T cells from LEW and LEWThemis1−/− rats assessed in coculture experiments with CellTrace-labeled naive LEW CD4+ T cells as effector cells stimulated with plate-bound anti-CD3 and irradiated syngeneic APC. For each experiment, we pooled the purified Tregs from three rats per group because of the low numbers of Tregs present in Themis1-deficient rats. (E) Representative histograms of the proliferation of effector cells alone or in the presence of LEW and LEWThemis1−/− Tregs (Treg/Teff ratio 1:2). Percentages indicate the proportion of CellTrace low effector cells. Results are representative of one of three independent experiments. (F) Suppressive activity of CD25bright CD4+ T cells from LEW and LEWThemis1−/− rats at different Treg/Teff ratios. Data are expressed as percentage of inhibition and presented as mean values ± SEM obtained from three independent experiments. *p ≤ 0.05. ns, not significant.

Close modal

The striking differences in the impact of Themis1 deficiency in BN or LEW genomic backgrounds raised the hypothesis that an epistasis phenomenon occurs in the BN rats. Thus, a peculiar BN gene(s) variant(s) could be required to confer the Treg defect and IBD development in the absence of Themis1. The TCR signaling molecule Vav1 could be the gene responsible for this epistatic phenomenon for the following reasons: 1) BN rats harbor a peculiar Vav1 variant (Vav1-W63) with reduced adaptor function and signaling capacities (7), and 2) Themis1 has been shown to interact and to positively regulate Vav1 phosphorylation in thymocytes (23). Therefore, because Themis1 is also expressed on rat Tregs (data not shown), the effect of Themis1 deficiency on BN Treg function could result from a cumulative impairment of the TCR signaling due to the association of Themis1 deficiency with the presence of the Vav1-W63 variant. To test this hypothesis, we crossed BNThemis1−/− rats (bearing the Vav1-W63 variant) with BN congenic rats bearing a 117-kb interval from LEW origin containing the wild type Vav1-R63 variant (congenic line named Bf) to generate the BfThemis1−/− rats. Those BfThemis1−/− rats of BN genomic background carry both the Themis1 deficiency and the Vav1-R63 variant (Supplemental Fig. 1C, 1D).

As observed in BNThemis1−/− and LEWThemis1−/−, BfThemis1−/− rats also showed an impaired thymic T cell development leading to a peripheral T cell lymphopenia, an increased frequency of activated peripheral CD4 T cells, and a biased cytokine production by total or naive CD4 T cells (Supplemental Fig. 3). We then investigated whether the differential impact of Themis1 deficiency on Treg functions observed between the BN and LEW strains depend on the 117-kb interval containing the LEW Vav1-R63 variant by analyzing Treg functions in the BfThemis1−/− strain. First, we controlled that Themis1 deficiency had no impact on the percentage of Foxp3-expressing cells among TCR+ CD4+ CD8 CD25bright sorted Tregs, nor on the level of Foxp3 expression (Fig. 4A, 4B), irrespective of the Vav1 variant. Then, using coculture experiments, we compared the suppressive activity of Tregs from BN and BNThemis1−/− (Vav1-W63 variant), on the one hand, and of Bf and Bf Themis1−/− (Vav1-R63 variant), on the other hand. As previously described, a clear defect in the suppressive activity of BNThemis1−/− Tregs was found. By contrast, under the same experimental conditions, BfThemis1−/− Tregs exhibited an intact in vitro suppressive activity as compared with Bf Tregs (Fig. 4C). These results reveal that the defect in Treg suppressive activity associated with Themis1 deficiency depends on the 117-kb interval containing the Vav1-W63 polymorphism. This defect concerns natural Tregs but not induced Tregs because BN and BN Themis1−/− Tregs generated in the presence of TGF-β exhibit similar suppressive activity (Supplemental Fig. 4).

FIGURE 4.

Effect of Themis1 deficiency on Treg functions depends on the 117-kb interval containing Vav1 variants. (A) Percentage of Foxp3+ cells and (B) Foxp3 expression levels in sorted CD25bright CD4+ CD8 TCR+ Tregs used in coculture experiments. (C) Suppressive activity of CD25bright CD4+ T cells from BN and BNThemis1−/− rats (left) or Bf and BfThemis1−/− rats (right) assessed in coculture experiments with CellTrace-labeled naive effector cells. Representative histograms of one experiment of four independent experiments are shown. The suppressive activity of CD25bright CD4+ T cells at different Treg/Teff ratios is expressed as percentage of inhibition of the response of effector cells alone. Results are presented as mean values ± SEM obtained from four independent experiments, each experiment including pools of four rats per group. *p ≤ 0.05. ns, not significant.

FIGURE 4.

Effect of Themis1 deficiency on Treg functions depends on the 117-kb interval containing Vav1 variants. (A) Percentage of Foxp3+ cells and (B) Foxp3 expression levels in sorted CD25bright CD4+ CD8 TCR+ Tregs used in coculture experiments. (C) Suppressive activity of CD25bright CD4+ T cells from BN and BNThemis1−/− rats (left) or Bf and BfThemis1−/− rats (right) assessed in coculture experiments with CellTrace-labeled naive effector cells. Representative histograms of one experiment of four independent experiments are shown. The suppressive activity of CD25bright CD4+ T cells at different Treg/Teff ratios is expressed as percentage of inhibition of the response of effector cells alone. Results are presented as mean values ± SEM obtained from four independent experiments, each experiment including pools of four rats per group. *p ≤ 0.05. ns, not significant.

Close modal

Because the Treg defect observed in the BNThemis1−/− is linked to the development of an IBD, we reasoned that BfThemis1−/− rats that harbor Tregs with normal suppressive function would not develop signs of inflammation. To test this hypothesis, we examined organs from BfThemis1−/−, BNThemis1−/−, and age-matched control Bf and BN rats. BfThemis1−/− rats were devoid of tissue inflammation, whereas BNThemis1−/− rats showed inflammation of the intestinal tract as revealed by the macroscopic lesions including edema, erythema, and increased mucous secretion. The disease prevalence rate reached 51% in BNThemis1−/− rats, whereas only 1 BfThemis1−/− rat out of 46 showed minor signs of inflammation (Fig. 5A). Histopathological studies (Fig. 5B, 5C) and quantitation of MPO activity in the various parts of the intestinal tract (Fig. 5D) confirmed that the introduction of the 117-kb interval from LEW origin protected BNThemis1−/− rats from IBD development. Thus, as anticipated, the restoration of the suppressive function of Tregs in Themis1-deficient rats bearing the 117 kb of LEW origin is accompanied by the absence of IBD.

FIGURE 5.

Effect of Themis1 deficiency on intestinal homeostasis depends on the 117-kb interval containing Vav1 variants. (A) Macroscopical scores of IBD in 12- to 16-wk-old BN (n = 31), Bf (n = 37), BNThemis1−/− (n = 37), and BfThemis1−/− rats (n = 46). Numbers represent the percentage of rats with scores >1. (B) H&E staining of duodenum, jejunum, ileum, and colon from BNThemis1−/− (Vav1-W63) and BfThemis1−/− (Vav1-R63) rats. The immune cell infiltrates, granulomas, and thickening of the intestinal wall are representative of microscopic lesions observed in all affected BNThemis1−/− rats and were not observed in BfThemis1−/−, BN, and Bf rats (original magnification ×100). (C) Microscopical scores of IBD in BNThemis1−/− with macroscopic lesions (n = 5) and BfThemis1−/− rats (n = 11). (D) MPO activity in the duodenum, jejunum, ileum, and colon tissue samples from BNThemis1−/− rats with macroscopic lesions (n = 11) and BfThemis1−/− rats (n = 16). Results are presented as mean ± SEM. *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001, ****p ≤ 0.0001.

FIGURE 5.

Effect of Themis1 deficiency on intestinal homeostasis depends on the 117-kb interval containing Vav1 variants. (A) Macroscopical scores of IBD in 12- to 16-wk-old BN (n = 31), Bf (n = 37), BNThemis1−/− (n = 37), and BfThemis1−/− rats (n = 46). Numbers represent the percentage of rats with scores >1. (B) H&E staining of duodenum, jejunum, ileum, and colon from BNThemis1−/− (Vav1-W63) and BfThemis1−/− (Vav1-R63) rats. The immune cell infiltrates, granulomas, and thickening of the intestinal wall are representative of microscopic lesions observed in all affected BNThemis1−/− rats and were not observed in BfThemis1−/−, BN, and Bf rats (original magnification ×100). (C) Microscopical scores of IBD in BNThemis1−/− with macroscopic lesions (n = 5) and BfThemis1−/− rats (n = 11). (D) MPO activity in the duodenum, jejunum, ileum, and colon tissue samples from BNThemis1−/− rats with macroscopic lesions (n = 11) and BfThemis1−/− rats (n = 16). Results are presented as mean ± SEM. *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001, ****p ≤ 0.0001.

Close modal

A previous study showed that Themis1 could interact with Vav1 (23). Because Vav1 appeared to be the strongest candidate gene responsible for the epistasis, we investigated whether Vav1-R63 and Vav1-W63 variants interact with Themis1 differently using cotransfection experiments in HEK cells. The obtained results reveal that the BN Vav1-W63 variant interacts more efficiently with Themis1 as compared with the Vav1-R63 variant (Fig. 6A, 6B). Next, we investigated the impact of Themis1 deficiency on the recruitment of Vav1 to the LAT platform after TCR stimulation using thymocytes. This analysis revealed that the absence of Themis1 reduced the recruitment of the Vav1-W63 variant from BN origin to LAT (Fig. 6C, 6D), whereas it had no impact on the recruitment of Vav1-R63 from LEW origin. However, despite this reduced recruitment, the phosphorylation of the Vav1-W63 variant was not affected by Themis1 deficiency (Fig. 6E). In cell expressing Vav1-W63 variant, Themis1 deficiency led to a decreased recruitment of PLC-γ and Grb2 to LAT, whereas it did not exert such an effect in Vav1-R63 expressing cells (Fig. 6C, 6D). We next analyzed the consequences of this modified proximal TCR signaling complex assembly on the downstream signaling pathways Erk and Akt, known to be respectively dependent on Vav1 adaptor function and enzymatic activity (15). Themis1 deficiency reduced Erk activation in cells expressing the Vav1-W63 variant but not in those expressing Vav1-R63 (Fig. 6E). In contrast, Akt activation was not affected by Themis1 deficiency. Together, these results show that Themis1 interacts more efficiently with Vav1-W63 variant and is required to stabilize its recruitment to the transmembrane adaptor LAT to fully promote the activation of Erk kinases that are likely necessary for the suppressive functions of Tregs (2931).

FIGURE 6.

Association of Themis1 deficiency and Vav1-W63 polymorphism impacts on LAT signalosome and impairs TCR-induced Erk signaling. (A) Western blot analysis of Vav1-W63 or Vav1-R63 interaction with Themis1. c-Myc–tagged Vav1 variants and flag-tagged Themis1 were overexpressed in HEK cells, and Themis1 were immunoprecipitated in untreated or pervanadate-stimulated HEK cells. (B) Densitometric quantification of Vav1 after Themis1 immunoprecipitation; relative intensities are normalized to the respective Themis1 signal. (C) Western blot analysis of Vav1, PLC-γ, Grb2, and LAT after LAT immunoprecipitation in unstimulated thymocytes (left) or after 1 min anti-TCR/anti-CD4 stimulation (right). Protein levels in whole cell lysates (WCL) are shown in the lower panel. (D) Densitometric quantification of Vav1, PLC-γ, and Grb2 after LAT immunoprecipitation; relative intensities normalized to the respective LAT signal are shown. (E) Western blot analysis of P-Erk, P-Vav1, P-Akt activity, and total Vav1, Erk, and Akt protein levels in anti-TCR/anti-CD4–stimulated thymocytes. Results shown are representative of two independent experiments.

FIGURE 6.

Association of Themis1 deficiency and Vav1-W63 polymorphism impacts on LAT signalosome and impairs TCR-induced Erk signaling. (A) Western blot analysis of Vav1-W63 or Vav1-R63 interaction with Themis1. c-Myc–tagged Vav1 variants and flag-tagged Themis1 were overexpressed in HEK cells, and Themis1 were immunoprecipitated in untreated or pervanadate-stimulated HEK cells. (B) Densitometric quantification of Vav1 after Themis1 immunoprecipitation; relative intensities are normalized to the respective Themis1 signal. (C) Western blot analysis of Vav1, PLC-γ, Grb2, and LAT after LAT immunoprecipitation in unstimulated thymocytes (left) or after 1 min anti-TCR/anti-CD4 stimulation (right). Protein levels in whole cell lysates (WCL) are shown in the lower panel. (D) Densitometric quantification of Vav1, PLC-γ, and Grb2 after LAT immunoprecipitation; relative intensities normalized to the respective LAT signal are shown. (E) Western blot analysis of P-Erk, P-Vav1, P-Akt activity, and total Vav1, Erk, and Akt protein levels in anti-TCR/anti-CD4–stimulated thymocytes. Results shown are representative of two independent experiments.

Close modal

Tregs play a pivotal role in maintaining immunological self-tolerance and immune system homeostasis. They are defined by the expression of Foxp3, a specific transcription factor necessary for their development and function (3236). However, there is substantial evidence that Foxp3 is not the sole cell-intrinsic regulator involved in the maintenance of Treg suppressive function (37, 38). In this context, we previously reported that deficiency for the TCR signaling molecule Themis1 in BN rats impairs the suppressive function of Tregs, despite normal levels of Foxp3 expression, and leads to the spontaneous development of IBD (10). In addition, we showed that these Tregs migrate normally in the intestine, thus excluding the involvement of Treg migration defect in this phenotype (Supplemental Fig. 5).

In this study, we show that the pathophysiological consequence of Themis1 deficiency depends on another gene that plays key roles in TCR-dependent signaling. Indeed, Themis1-deficient LEW rats, although lymphopenic, did not develop IBD or other spontaneous diseases and their Tregs were fully functional. These results exclude lymphopenia as the factor responsible for the Treg defect in Themis1-deficient BN rats and favor the hypothesis of an intrinsic role of Themis1 in controlling Treg function. The striking differences observed in Themis1-deficient rats of the BN and LEW genetic backgrounds suggest that an epistasis between Themis1 and other gene(s) controls Treg function. We showed that the 117-kb locus on rat chromosome 9 previously implicated in modulating the susceptibility of LEW and BN rats to immune-mediated diseases (7, 9) contains the gene involved in this epistasis. Indeed, the functional defect of Tregs and the spontaneous IBD development were observed only in rats that combine Themis1 deficiency with the presence of the 117-kb locus of BN origin. This 117-kb region contains only four genes: C3 (complement 3 precursor), Gpr108 (G protein–coupled receptor 108), Cip4/Trip10 (Cdc42-interacting protein 4), and exons 1–15 of Vav1. The coding regions of these four genes were previously sequenced and led to the identification of only two nonsynonymous polymorphisms in exon 32 of C3 and in exon 1 of Vav1 (7). Vav1 appeared to be a strong candidate gene responsible for the epistasis in the BN rats for the following reasons. First, Themis1 and Vav1 are components of the same SLP76-LAT signalosome (22, 23, 39, 40). Second, Vav1 phosphorylation is decreased in Themis1-deficient mice after TCR engagement (23), suggesting that Themis1 is a positive regulator of Vav1. Third, the BN Vav1 contains a polymorphism (R63W) responsible for the reduction of its adaptor function (7). Fourth, the BN Vav1-W63 variant interacts more efficiently with Themis1 as compared with the LEW Vav1-R63 variant. Finally, we show that Themis1 deficiency reduces the recruitment of Vav1-W63 to LAT signalosome and impairs TCR signaling as revealed by the defect in Erk activation. In contrast, Themis1 deficiency has no apparent effect on ERK activation in rat thymocytes containing the Vav1-R63 variant, consistent with findings obtained in Themis1-deficient C57BL/6 mice, which harbor the Vav1-R63 variant (19). Thus, this study suggests that Themis1 and Vav1 form a cooperative signaling complex important for the optimal transmission of TCR-mediated signals necessary for Tregs to mediate their suppressive functions.

Besides the functional Treg defect observed in Themis1-deficient BN rats, we previously reported that effector CD4 T cells secreted less IFN-γ and higher amounts of Th2/Th17 cytokines (10). This effect may be directly related to the role played by Themis1 in TCR signaling in effector T cells or to a secondary effect of lymphopenia induced by Themis1 deficiency. In this article, we show that Themis1-deficient LEW and Bf CD4 T cells produce exacerbated Th2/Th17 cytokines, whereas their production of IFN-γ was not impacted by Themis1 deficiency. These data suggest that the cooperation between Themis1 and Vav1-W63 is mandatory for IFN-γ production, whereas the exacerbated production of Th2/Th17 cytokines might be secondary to Themis1 deficiency–induced lymphopenia. In this regard, it has been shown that lymphopenia may influence the orientation of immune responses by affecting the differentiation of T cell subsets (41, 42). Thus, depending on the genetic background, Themis1 might influence the development of immune-mediated diseases by affecting both Treg suppressive functions and the production of IFN-γ by effector CD4 T cells.

The mechanism whereby the combined presence of Vav1-W63 variant and Themis1 deficiency lead to functional defect of Tregs and reduced IFN-γ production by effector CD4 T cells is still unknown. We previously showed that the Vav1-W63 variant is constitutively active and exhibits normal GEF activity but reduced adaptor functions associated with a marked reduction in Vav1 protein levels (7). In this study, we showed that Themis1 interacts more efficiently with Vav1-W63, and this may promote proximal TCR signaling assembly to compensate for the signaling defect related to Vav1-W63. Indeed, we showed that Themis1 favors Vav1-W63 recruitment to LAT and this was associated with an increased Grb2 and PLC-γ recruitment to this signaling platform. In Themis1-deficient BN rats, we observed a reduction of Vav1-W63 recruitment to LAT that is likely responsible for the defect in Erk activation, because Erk activity depends on Vav1 adaptor function (15). Of note, in the absence of Themis1 in BN rats, PLCγ and Grb2 recruitment to LAT remains in a similar range compared with the Bf strain. These results indicate that the compensatory effect of Themis1 in the presence of the W63-Vav1 might involve, besides PLCγ and Grb2, other molecules to stabilize the LAT signalosome complex. This study revealed an association among the reduced Vav1-W63 recruitment to the LAT signalosome, the defective Erk activity, the Treg functional defect, and the decreased production of IFN-γ by effector T cells. These data are consistent with previous reports showing the critical role of the integrity of LAT signaling platform as well as Erk activation for the Treg suppressive function and the Th1/Th2 polarization (2931, 43, 44). Indeed, the disruption of LAT/ PLC-γ signalosome assembly reduces Erk and Ca2+ activation (29, 30) and despite normal Foxp3 expression, LAT- or PLC-γ–deficient Tregs are ineffective in controlling the proliferation of effector T cells (2931). In addition, sustained Erk activity and calcium signaling impact on the Th1/Th2 polarization, promoting Th1 differentiation (4345). Finally, it has been shown that TCR-dependent signals promote epigenetic modifications (46) that are crucial for the acquisition of a functional and stable Treg phenotype by developing Foxp3+ cells (47). It remains to be determined whether the functional defect of Tregs in Themis1-deficient BN rats depends on epigenetic modifications induced by the reduction of TCR signaling during Treg development.

This study describes for the first time, to our knowledge, the importance of the cooperation between Themis1 and Vav1 for the control of Treg suppressive function and the development of a spontaneous inflammatory disease. The epistatic interaction between Themis1 and Vav1 could have important repercussions in genetics of human pathologies. Indeed, Vav1 has been pointed out as a susceptibility gene for multiple sclerosis (8). Themis1 also has been identified as a potential candidate for Crohn’s disease (48), celiac disease (4951), and multiple sclerosis (52). Given the strong impact of the Themis1/Vav1 cooperation revealed by our study, it remains to be tested whether the association of variants of those two genes could exert synergetic effects on susceptibility to immune-mediated pathologies.

In conclusion, our study highlights the importance of studying gene–gene interactions in complex diseases and provides insights into the underlying cellular and molecular mechanisms. Given the large number of genome-wide association studies that have recently been performed, it is clear that genome-wide interaction testing will be the logical next step after single-locus testing. This will require the development of new powerful bioinformatics tools that take into account gene–gene interactions in GWAS data. Such approaches should become possible in the coming years (5355) and are likely to generate novel advances in our understanding of the genetic component in complex diseases.

We thank P. Rousset for technical help with MPO assays; D. Dunia, R. Liblau, and L. Mars for comments on the manuscript; and the flow cytometry core facility (Centre de Physiopathologie de Toulouse Purpan), the animal house staff, and the histopathology core facility for assistance (Inserm US006 ANEXPLO/CREFRE). This work is dedicated to the memory of our mentor Prof. Philippe Druet, who recently died.

This work was supported by INSERM (to G.J.F., I.B., A.S.D., and D.L.), the Association Française contre les Myopathies, the Agence Nationale de la Recherche (Grant ANR-08-GENO-041-01), the Fondation pour la Recherche Médicale (Grant DEQ2000326531), the Association de Recherche sur la Sclérose en Plaques, the Région Midi-Pyrénées and a Fight-MG European grant (Grant FP7-Health-2009-242210), the Ministère de l’Education Nationale de la Recherche et de la Technologie (to C.P.), the Fondation pour la Recherche Médicale (to C.P.), and the Centre National de la Recherche Scientifique (to A.S.).

The online version of this article contains supplemental material.

Abbreviations used in this article:

Bf

BN.LEWc9-Bf congenic line

BN

Brown–Norway

GEF

guanine nucleotide exchange factor

IBD

inflammatory bowel disease

LEW

Lewis

MPO

myeloperoxidase

Treg

regulatory T cell.

1
Hunter
D. J.
2005
.
Gene-environment interactions in human diseases.
Nat. Rev. Genet.
6
:
287
298
.
2
Zuk
O.
,
Hechter
E.
,
Sunyaev
S. R.
,
Lander
E. S.
.
2012
.
The mystery of missing heritability: genetic interactions create phantom heritability.
Proc. Natl. Acad. Sci. USA
109
:
1193
1198
.
3
Cordell
H. J.
2009
.
Detecting gene-gene interactions that underlie human diseases.
Nat. Rev. Genet.
10
:
392
404
.
4
McKinney
B. A.
,
Pajewski
N. M.
.
2011
.
Six degrees of epistasis: statistical network models for GWAS.
Front. Genet.
2
:
109
.
5
Rose
A. M.
,
Bell
L. C.
.
2012
.
Epistasis and immunity: the role of genetic interactions in autoimmune diseases.
Immunology
137
:
131
138
.
6
Mackay
T. F.
2014
.
Epistasis and quantitative traits: using model organisms to study gene-gene interactions.
Nat. Rev. Genet.
15
:
22
33
.
7
Colacios
C.
,
Casemayou
A.
,
Dejean
A. S.
,
Gaits-Iacovoni
F.
,
Pedros
C.
,
Bernard
I.
,
Lagrange
D.
,
Deckert
M.
,
Lamouroux
L.
,
Jagodic
M.
, et al
.
2011
.
The p.Arg63Trp polymorphism controls Vav1 functions and Foxp3 regulatory T cell development.
J. Exp. Med.
208
:
2183
2191
.
8
Jagodic
M.
,
Colacios
C.
,
Nohra
R.
,
Dejean
A. S.
,
Beyeen
A. D.
,
Khademi
M.
,
Casemayou
A.
,
Lamouroux
L.
,
Duthoit
C.
,
Papapietro
O.
, et al
.
2009
.
A role for VAV1 in experimental autoimmune encephalomyelitis and multiple sclerosis.
Sci. Transl. Med.
1
:
10ra21
.
9
Pedros
C.
,
Papapietro
O.
,
Colacios
C.
,
Casemayou
A.
,
Bernard
I.
,
Garcia
V.
,
Lagrange
D.
,
Mariamé
B.
,
Andreoletti
O.
,
Fournié
G. J.
,
Saoudi
A.
.
2013
.
Genetic control of HgCl2-induced IgE and autoimmunity by a 117-kb interval on rat chromosome 9 through CD4 CD45RChigh T cells.
Genes Immun.
14
:
258
267
.
10
Chabod
M.
,
Pedros
C.
,
Lamouroux
L.
,
Colacios
C.
,
Bernard
I.
,
Lagrange
D.
,
Balz-Hara
D.
,
Mosnier
J. F.
,
Laboisse
C.
,
Vergnolle
N.
, et al
.
2012
.
A spontaneous mutation of the rat Themis gene leads to impaired function of regulatory T cells linked to inflammatory bowel disease.
PLoS Genet.
8
:
e1002461
.
11
Fischer
K. D.
,
Zmuldzinas
A.
,
Gardner
S.
,
Barbacid
M.
,
Bernstein
A.
,
Guidos
C.
.
1995
.
Defective T-cell receptor signalling and positive selection of Vav-deficient CD4+ CD8+ thymocytes.
Nature
374
:
474
477
.
12
Tarakhovsky
A.
,
Turner
M.
,
Schaal
S.
,
Mee
P. J.
,
Duddy
L. P.
,
Rajewsky
K.
,
Tybulewicz
V. L.
.
1995
.
Defective antigen receptor-mediated proliferation of B and T cells in the absence of Vav.
Nature
374
:
467
470
.
13
Turner
M.
,
Mee
P. J.
,
Walters
A. E.
,
Quinn
M. E.
,
Mellor
A. L.
,
Zamoyska
R.
,
Tybulewicz
V. L.
.
1997
.
A requirement for the Rho-family GTP exchange factor Vav in positive and negative selection of thymocytes.
Immunity
7
:
451
460
.
14
Tybulewicz
V. L.
2005
.
Vav-family proteins in T-cell signalling.
Curr. Opin. Immunol.
17
:
267
274
.
15
Saveliev
A.
,
Vanes
L.
,
Ksionda
O.
,
Rapley
J.
,
Smerdon
S. J.
,
Rittinger
K.
,
Tybulewicz
V. L.
.
2009
.
Function of the nucleotide exchange activity of vav1 in T cell development and activation.
Sci. Signal.
2
:
ra83
.
16
Fournié
G. J.
,
Cautain
B.
,
Xystrakis
E.
,
Damoiseaux
J.
,
Mas
M.
,
Lagrange
D.
,
Bernard
I.
,
Subra
J. F.
,
Pelletier
L.
,
Druet
P.
,
Saoudi
A.
.
2001
.
Cellular and genetic factors involved in the difference between Brown Norway and Lewis rats to develop respectively type-2 and type-1 immune-mediated diseases.
Immunol. Rev.
184
:
145
160
.
17
Lesourne
R.
,
Uehara
S.
,
Lee
J.
,
Song
K. D.
,
Li
L.
,
Pinkhasov
J.
,
Zhang
Y.
,
Weng
N. P.
,
Wildt
K. F.
,
Wang
L.
, et al
.
2009
.
Themis, a T cell-specific protein important for late thymocyte development.
Nat. Immunol.
10
:
840
847
.
18
Fu
G.
,
Vallée
S.
,
Rybakin
V.
,
McGuire
M. V.
,
Ampudia
J.
,
Brockmeyer
C.
,
Salek
M.
,
Fallen
P. R.
,
Hoerter
J. A.
,
Munshi
A.
, et al
.
2009
.
Themis controls thymocyte selection through regulation of T cell antigen receptor-mediated signaling.
Nat. Immunol.
10
:
848
856
.
19
Johnson
A. L.
,
Aravind
L.
,
Shulzhenko
N.
,
Morgun
A.
,
Choi
S. Y.
,
Crockford
T. L.
,
Lambe
T.
,
Domaschenz
H.
,
Kucharska
E. M.
,
Zheng
L.
, et al
.
2009
.
Themis is a member of a new metazoan gene family and is required for the completion of thymocyte positive selection.
Nat. Immunol.
10
:
831
839
.
20
Patrick
M. S.
,
Oda
H.
,
Hayakawa
K.
,
Sato
Y.
,
Eshima
K.
,
Kirikae
T.
,
Iemura
S.
,
Shirai
M.
,
Abe
T.
,
Natsume
T.
, et al
.
2009
.
Gasp, a Grb2-associating protein, is critical for positive selection of thymocytes.
Proc. Natl. Acad. Sci. USA
106
:
16345
16350
.
21
Kakugawa
K.
,
Yasuda
T.
,
Miura
I.
,
Kobayashi
A.
,
Fukiage
H.
,
Satoh
R.
,
Matsuda
M.
,
Koseki
H.
,
Wakana
S.
,
Kawamoto
H.
,
Yoshida
H.
.
2009
.
A novel gene essential for the development of single positive thymocytes.
Mol. Cell. Biol.
29
:
5128
5135
.
22
Paster
W.
,
Brockmeyer
C.
,
Fu
G.
,
Simister
P. C.
,
de Wet
B.
,
Martinez-Riaño
A.
,
Hoerter
J. A.
,
Feller
S. M.
,
Wülfing
C.
,
Gascoigne
N. R.
,
Acuto
O.
.
2013
.
GRB2-mediated recruitment of THEMIS to LAT is essential for thymocyte development.
J. Immunol.
190
:
3749
3756
.
23
Lesourne
R.
,
Zvezdova
E.
,
Song
K. D.
,
El-Khoury
D.
,
Uehara
S.
,
Barr
V. A.
,
Samelson
L. E.
,
Love
P. E.
.
2012
.
Interchangeability of Themis1 and Themis2 in thymocyte development reveals two related proteins with conserved molecular function.
J. Immunol.
189
:
1154
1161
.
24
Brockmeyer
C.
,
Paster
W.
,
Pepper
D.
,
Tan
C. P.
,
Trudgian
D. C.
,
McGowan
S.
,
Fu
G.
,
Gascoigne
N. R.
,
Acuto
O.
,
Salek
M.
.
2011
.
T cell receptor (TCR)-induced tyrosine phosphorylation dynamics identifies THEMIS as a new TCR signalosome component.
J. Biol. Chem.
286
:
7535
7547
.
25
Fu
G.
,
Casas
J.
,
Rigaud
S.
,
Rybakin
V.
,
Lambolez
F.
,
Brzostek
J.
,
Hoerter
J. A.
,
Paster
W.
,
Acuto
O.
,
Cheroutre
H.
, et al
.
2013
.
Themis sets the signal threshold for positive and negative selection in T-cell development.
Nature
504
:
441
445
.
26
Wakeland
E.
,
Morel
L.
,
Achey
K.
,
Yui
M.
,
Longmate
J.
.
1997
.
Speed congenics: a classic technique in the fast lane (relatively speaking).
Immunol. Today
18
:
472
477
.
27
Lagrange
D.
,
Fournié
G. J.
.
2010
.
Generation of congenic and consomic rat strains.
Methods Mol. Biol.
597
:
243
266
.
28
Mas
M.
,
Subra
J. F.
,
Lagrange
D.
,
Pilipenko-Appolinaire
S.
,
Kermarrec
N.
,
Gauguier
D.
,
Druet
P.
,
Fournié
G. J.
.
2000
.
Rat chromosome 9 bears a major susceptibility locus for IgE response.
Eur. J. Immunol.
30
:
1698
1705
.
29
Chuck
M. I.
,
Zhu
M.
,
Shen
S.
,
Zhang
W.
.
2010
.
The role of the LAT-PLC-gamma1 interaction in T regulatory cell function.
J. Immunol.
184
:
2476
2486
.
30
Fu
G.
,
Chen
Y.
,
Yu
M.
,
Podd
A.
,
Schuman
J.
,
He
Y.
,
Di
L.
,
Yassai
M.
,
Haribhai
D.
,
North
P. E.
, et al
.
2010
.
Phospholipase Cgamma1 is essential for T cell development, activation, and tolerance.
J. Exp. Med.
207
:
309
318
.
31
Koonpaew
S.
,
Shen
S.
,
Flowers
L.
,
Zhang
W.
.
2006
.
LAT-mediated signaling in CD4+CD25+ regulatory T cell development.
J. Exp. Med.
203
:
119
129
.
32
Hori
S.
,
Nomura
T.
,
Sakaguchi
S.
.
2003
.
Control of regulatory T cell development by the transcription factor Foxp3.
Science
299
:
1057
1061
.
33
Fontenot
J. D.
,
Gavin
M. A.
,
Rudensky
A. Y.
.
2003
.
Foxp3 programs the development and function of CD4+CD25+ regulatory T cells.
Nat. Immunol.
4
:
330
336
.
34
Brunkow
M. E.
,
Jeffery
E. W.
,
Hjerrild
K. A.
,
Paeper
B.
,
Clark
L. B.
,
Yasayko
S. A.
,
Wilkinson
J. E.
,
Galas
D.
,
Ziegler
S. F.
,
Ramsdell
F.
.
2001
.
Disruption of a new forkhead/winged-helix protein, scurfin, results in the fatal lymphoproliferative disorder of the scurfy mouse.
Nat. Genet.
27
:
68
73
.
35
Bennett
C. L.
,
Christie
J.
,
Ramsdell
F.
,
Brunkow
M. E.
,
Ferguson
P. J.
,
Whitesell
L.
,
Kelly
T. E.
,
Saulsbury
F. T.
,
Chance
P. F.
,
Ochs
H. D.
.
2001
.
The immune dysregulation, polyendocrinopathy, enteropathy, X-linked syndrome (IPEX) is caused by mutations of FOXP3.
Nat. Genet.
27
:
20
21
.
36
Wildin
R. S.
,
Ramsdell
F.
,
Peake
J.
,
Faravelli
F.
,
Casanova
J. L.
,
Buist
N.
,
Levy-Lahad
E.
,
Mazzella
M.
,
Goulet
O.
,
Perroni
L.
, et al
.
2001
.
X-linked neonatal diabetes mellitus, enteropathy and endocrinopathy syndrome is the human equivalent of mouse scurfy.
Nat. Genet.
27
:
18
20
.
37
Hill
J. A.
,
Feuerer
M.
,
Tash
K.
,
Haxhinasto
S.
,
Perez
J.
,
Melamed
R.
,
Mathis
D.
,
Benoist
C.
.
2007
.
Foxp3 transcription-factor-dependent and -independent regulation of the regulatory T cell transcriptional signature.
Immunity
27
:
786
800
.
38
Sugimoto
N.
,
Oida
T.
,
Hirota
K.
,
Nakamura
K.
,
Nomura
T.
,
Uchiyama
T.
,
Sakaguchi
S.
.
2006
.
Foxp3-dependent and -independent molecules specific for CD25+CD4+ natural regulatory T cells revealed by DNA microarray analysis.
Int. Immunol.
18
:
1197
1209
.
39
Zhang
W.
,
Sloan-Lancaster
J.
,
Kitchen
J.
,
Trible
R. P.
,
Samelson
L. E.
.
1998
.
LAT: the ZAP-70 tyrosine kinase substrate that links T cell receptor to cellular activation.
Cell
92
:
83
92
.
40
Perez-Villar
J. J.
,
Whitney
G. S.
,
Sitnick
M. T.
,
Dunn
R. J.
,
Venkatesan
S.
,
O’Day
K.
,
Schieven
G. L.
,
Lin
T. A.
,
Kanner
S. B.
.
2002
.
Phosphorylation of the linker for activation of T-cells by Itk promotes recruitment of Vav.
Biochemistry
41
:
10732
10740
.
41
King
C.
,
Ilic
A.
,
Koelsch
K.
,
Sarvetnick
N.
.
2004
.
Homeostatic expansion of T cells during immune insufficiency generates autoimmunity.
Cell
117
:
265
277
.
42
Le Campion
A.
,
Gagnerault
M. C.
,
Auffray
C.
,
Bécourt
C.
,
Poitrasson-Rivière
M.
,
Lallemand
E.
,
Bienvenu
B.
,
Martin
B.
,
Lepault
F.
,
Lucas
B.
.
2009
.
Lymphopenia-induced spontaneous T-cell proliferation as a cofactor for autoimmune disease development.
Blood
114
:
1784
1793
.
43
Jorritsma
P. J.
,
Brogdon
J. L.
,
Bottomly
K.
.
2003
.
Role of TCR-induced extracellular signal-regulated kinase activation in the regulation of early IL-4 expression in naive CD4+ T cells.
J. Immunol.
170
:
2427
2434
.
44
Yamane
H.
,
Zhu
J.
,
Paul
W. E.
.
2005
.
Independent roles for IL-2 and GATA-3 in stimulating naive CD4+ T cells to generate a Th2-inducing cytokine environment.
J. Exp. Med.
202
:
793
804
.
45
Noble
A.
,
Truman
J. P.
,
Vyas
B.
,
Vukmanovic-Stejic
M.
,
Hirst
W. J.
,
Kemeny
D. M.
.
2000
.
The balance of protein kinase C and calcium signaling directs T cell subset development.
J. Immunol.
164
:
1807
1813
.
46
Ohkura
N.
,
Hamaguchi
M.
,
Morikawa
H.
,
Sugimura
K.
,
Tanaka
A.
,
Ito
Y.
,
Osaki
M.
,
Tanaka
Y.
,
Yamashita
R.
,
Nakano
N.
, et al
.
2012
.
T cell receptor stimulation-induced epigenetic changes and Foxp3 expression are independent and complementary events required for Treg cell development.
Immunity
37
:
785
799
.
47
Ohkura
N.
,
Kitagawa
Y.
,
Sakaguchi
S.
.
2013
.
Development and maintenance of regulatory T cells.
Immunity
38
:
414
423
.
48
Jostins
L.
,
Ripke
S.
,
Weersma
R. K.
,
Duerr
R. H.
,
McGovern
D. P.
,
Hui
K. Y.
,
Lee
J. C.
,
Schumm
L. P.
,
Sharma
Y.
,
Anderson
C. A.
, et al
International IBD Genetics Consortium (IIBDGC)
.
2012
.
Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease.
Nature
491
:
119
124
.
49
Dubois
P. C.
,
Trynka
G.
,
Franke
L.
,
Hunt
K. A.
,
Romanos
J.
,
Curtotti
A.
,
Zhernakova
A.
,
Heap
G. A.
,
Adány
R.
,
Aromaa
A.
, et al
.
2010
.
Multiple common variants for celiac disease influencing immune gene expression.
Nat. Genet.
42
:
295
302
.
50
Bondar
C.
,
Plaza-Izurieta
L.
,
Fernandez-Jimenez
N.
,
Irastorza
I.
,
Withoff
S.
,
Wijmenga
C.
,
Chirdo
F.
,
Bilbao
J. R.
CEGEC
.
2014
.
THEMIS and PTPRK in celiac intestinal mucosa: coexpression in disease and after in vitro gliadin challenge.
Eur. J. Hum. Genet.
22
:
358
362
.
51
Senapati
S.
,
Gutierrez-Achury
J.
,
Sood
A.
,
Midha
V.
,
Szperl
A.
,
Romanos
J.
,
Zhernakova
A.
,
Franke
L.
,
Alonso
S.
,
Thelma
B. K.
, et al
.
2015
.
Evaluation of European coeliac disease risk variants in a north Indian population.
Eur. J. Hum. Genet.
23
:
530
535
.
52
International Multiple Sclerosis Genetics Consortium, Wellcome Trust Case Control Consortium, S. Sawcer, G. Hellenthal, M. Pirinen, C. C. Spencer, N. A. Patsopoulos, L. Moutsianas, A. Dilthey, Z. Su, et al. 2011. Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis. Nature 476: 214–219
.
53
Zhu
Z.
,
Tong
X.
,
Zhu
Z.
,
Liang
M.
,
Cui
W.
,
Su
K.
,
Li
M. D.
,
Zhu
J.
.
2013
.
Development of GMDR-GPU for gene-gene interaction analysis and its application to WTCCC GWAS data for type 2 diabetes.
PLoS One
8
:
e61943
.
54
Liao
Z.
,
Zeng
Q.
,
Liao
B.
,
Li
X.
.
2014
.
A novel two-stage approach for epistasis detection in genome-wide case-control studies.
Biochem. Genet.
52
:
403
414
.
55
Hemani
G.
,
Shakhbazov
K.
,
Westra
H. J.
,
Esko
T.
,
Henders
A. K.
,
McRae
A. F.
,
Yang
J.
,
Gibson
G.
,
Martin
N. G.
,
Metspalu
A.
, et al
.
2014
.
Detection and replication of epistasis influencing transcription in humans.
Nature
508
:
249
253
.

The authors have no financial conflicts of interest.

Supplementary data