Ricardo C. Ferreira

Widespread seasonal gene expression reveals annual differences in human immunity and physiology.

Seasonal variations are rarely considered a contributing component to human tissue function or health, although many diseases and physiological process display annual periodicities. Here we find more than 4,000 protein-coding mRNAs in white blood cells and adipose tissue to have seasonal expression profiles, with inverted patterns observed between Europe and Oceania. We also find the cellular composition of blood to vary by season, and these changes, which differ between the United Kingdom and The Gambia, could explain the gene expression periodicity. With regards to tissue function, the immune system has a profound pro-inflammatory transcriptomic profile during European winter, with increased levels of soluble IL-6 receptor and C-reactive protein, risk biomarkers for cardiovascular, psychiatric and autoimmune diseases that have peak incidences in winter. Circannual rhythms thus require further exploration as contributors to various aspects of human physiology and disease.

Evidence for CTLA4 as a susceptibility gene for systemic lupus erythematosus.

Several lines of evidence implicate the Cytotoxic T Lymphocyte Antigen 4 (CTLA4) gene in susceptibility to autoimmune disease. We have examined the association of systemic lupus erythematosus (SLE) with polymorhisms within the CTLA4 gene that were previously proposed to regulate CTLA-4 function: a single nucleotide polymorphism (SNP) in position +49 of exon 1 and a dinucleotide repeat in the 3′ untranslated region (3’UTR). The 3′UTR repeat showed a significant association with SLE, with one allele conferring susceptibility and another conferring protection to the disease. The associated alleles do not support previous suggestions of an allele size-dependent effect of the 3’ UTR polymorphism in autoimmunity development and instead suggest that it is in linkage disequilibrium with a true causative locus. No association of the exon 1 SNP with SLE was found in our population. Given the conflicting results obtained in different studies on the association of SLE with this polymorphism, we performed a meta-analysis including seven previously published studies and the present one. Significantly increased and decreased risks for SLE were found for carriers of the G allele and the A allele, respectively. The functional characterization of disease-associated CTLA4 gene variants is now required to elucidate their role in the pathogenesis of SLE and other autoimmune diseases.

A type I interferon transcriptional signature precedes autoimmunity in children genetically at risk for type 1 diabetes.

Diagnosis of the autoimmune disease type 1 diabetes (T1D) is preceded by the appearance of circulating autoantibodies to pancreatic islets. However, almost nothing is known about events leading to this islet autoimmunity. Previous epidemiological and genetic data have associated viral infections and antiviral type I interferon (IFN) immune response genes with T1D. Here, we first used DNA microarray analysis to identify IFN-β–inducible genes in vitro and then used this set of genes to define an IFN-inducible transcriptional signature in peripheral blood mononuclear cells from a group of active systemic lupus erythematosus patients (n = 25). Using this predefined set of 225 IFN signature genes, we investigated the expression of the signature in cohorts of healthy controls (n = 87), patients with T1D (n = 64), and a large longitudinal birth cohort of children genetically predisposed to T1D (n = 109; 454 microarrayed samples). Expression of the IFN signature was increased in genetically predisposed children before the development of autoantibodies (P = 0.0012) but not in patients with established T1D. Upregulation of IFN-inducible genes was transient, temporally associated with a recent history of upper respiratory tract infections (P = 0.0064), and marked by increased expression of SIGLEC-1 (CD169), a lectin-like receptor expressed on CD14+ monocytes. DNA variation in IFN-inducible genes altered T1D risk (P = 0.007), as exemplified by IFIH1, one of the genes in our IFN signature for which increased expression is a known risk factor for disease. These findings identify transient increased expression of type I IFN genes in preclinical diabetes as a risk factor for autoimmunity in children with a genetic predisposition to T1D.

Functional IL6R 358Ala Allele Impairs Classical IL-6 Receptor Signaling and Influences Risk of Diverse Inflammatory Diseases

Inflammation, which is directly regulated by interleukin-6 (IL-6) signaling, is implicated in the etiology of several chronic diseases. Although a common, non-synonymous variant in the IL-6 receptor gene (IL6R Asp358Ala; rs2228145 A>C) is associated with the risk of several common diseases, with the 358Ala allele conferring protection from coronary heart disease (CHD), rheumatoid arthritis (RA), atrial fibrillation (AF), abdominal aortic aneurysm (AAA), and increased susceptibility to asthma, the variants effect on IL-6 signaling is not known. Here we provide evidence for the association of this non-synonymous variant with the risk of type 1 diabetes (T1D) in two independent populations and confirm that rs2228145 is the major determinant of the concentration of circulating soluble IL-6R (sIL-6R) levels (34.6% increase in sIL-6R per copy of the minor allele 358Ala; rs2228145 [C]). To further investigate the molecular mechanism of this variant, we analyzed expression of IL-6R in peripheral blood mononuclear cells (PBMCs) in 128 volunteers from the Cambridge BioResource. We demonstrate that, although 358Ala increases transcription of the soluble IL6R isoform (P = 8.3×10−22) and not the membrane-bound isoform, 358Ala reduces surface expression of IL-6R on CD4+ T cells and monocytes (up to 28% reduction per allele; P≤5.6×10−22). Importantly, reduced expression of membrane-bound IL-6R resulted in impaired IL-6 responsiveness, as measured by decreased phosphorylation of the transcription factors STAT3 and STAT1 following stimulation with IL-6 (P≤5.2×10−7). Our findings elucidate the regulation of IL-6 signaling by IL-6R, which is causally relevant to several complex diseases, identify mechanisms for new approaches to target the IL-6/IL-6R axis, and anticipate differences in treatment response to IL-6 therapies based on this common IL6R variant.

IL-21 production by CD4+ effector T cells and frequency of circulating follicular helper T cells are increased in type 1 diabetes patients

Aims/hypothesisType 1 diabetes results from the autoimmune destruction of insulin-secreting pancreatic beta cells by T cells. Despite the established role of T cells in the pathogenesis of the disease, to date, with the exception of the identification of islet-specific T effector (Teff) cells, studies have mostly failed to identify reproducible alterations in the frequency or function of T cell subsets in peripheral blood from patients with type 1 diabetes.MethodsWe assessed the production of the proinflammatory cytokines IL-21, IFN-γ and IL-17 in peripheral blood mononuclear cells from 69 patients with type 1 diabetes and 61 healthy donors. In an additional cohort of 30 patients with type 1 diabetes and 32 healthy donors, we assessed the frequency of circulating T follicular helper (Tfh) cells in whole blood. IL-21 and IL-17 production was also measured in peripheral blood mononuclear cells (PBMCs) from a subset of 46 of the 62 donors immunophenotyped for Tfh.ResultsWe found a 21.9% (95% CI 5.8, 40.2; p = 3.9 × 10−3) higher frequency of IL-21+ CD45RA− memory CD4+ Teffs in patients with type 1 diabetes (geometric mean 5.92% [95% CI 5.44, 6.44]) compared with healthy donors (geometric mean 4.88% [95% CI 4.33, 5.50]). Consistent with this finding, we found a 14.9% increase in circulating Tfh cells in the patients (95% CI 2.9, 26.9; p = 0.016).Conclusions/interpretationThese results indicate that increased IL-21 production is likely to be an aetiological factor in the pathogenesis of type 1 diabetes that could be considered as a potential therapeutic target.

Compensatory T-Cell Regulation in Unaffected Relatives of SLE Patients, and Opposite IL-2/CD25-Mediated Effects Suggested by Coreferentiality Modeling

In human systemic lupus erythematosus (SLE), diverse autoantibodies accumulate over years before disease manifestation. Unaffected relatives of SLE patients frequently share a sustained production of autoantibodies with indiscriminable specificity, usually without ever acquiring the disease. We studied relations of IgG autoantibody profiles and peripheral blood activated regulatory T-cells (aTregs), represented by CD4+CD25bright T-cells that were regularly 70–90% Foxp3+. We found consistent positive correlations of broad-range as well as specific SLE-associated IgG with aTreg frequencies within unaffected relatives, but not patients or unrelated controls. Our interpretation: unaffected relatives with shared genetic factors compensated pathogenic effects by aTregs engaged in parallel with the individual autoantibody production. To study this further, we applied a novel analytic approach named coreferentiality that tests the indirect relatedness of parameters in respect to multivariate phenotype data. Results show that independently of their direct correlation, aTreg frequencies and specific SLE-associated IgG were likely functionally related in unaffected relatives: they significantly parallelled each other in their relations to broad-range immunoblot autoantibody profiles. In unaffected relatives, we also found coreferential effects of genetic variation in the loci encoding IL-2 and CD25. A model of CD25 functional genetic effects constructed by coreferentiality maximization suggests that IL-2-CD25 interaction, likely stimulating aTregs in unaffected relatives, had an opposed effect in SLE patients, presumably triggering primarily T-effector cells in this group. Coreferentiality modeling as we do it here could also be useful in other contexts, particularly to explore combined functional genetic effects.

Regulatory T Cell Responses in Participants with Type 1 Diabetes after a Single Dose of Interleukin-2: A Non-Randomised, Open Label, Adaptive Dose-Finding Trial

Background Interleukin-2 (IL-2) has an essential role in the expansion and function of CD4+ regulatory T cells (Tregs). Tregs reduce tissue damage by limiting the immune response following infection and regulate autoreactive CD4+ effector T cells (Teffs) to prevent autoimmune diseases, such as type 1 diabetes (T1D). Genetic susceptibility to T1D causes alterations in the IL-2 pathway, a finding that supports Tregs as a cellular therapeutic target. Aldesleukin (Proleukin; recombinant human IL-2), which is administered at high doses to activate the immune system in cancer immunotherapy, is now being repositioned to treat inflammatory and autoimmune disorders at lower doses by targeting Tregs. Methods and Findings To define the aldesleukin dose response for Tregs and to find doses that increase Tregs physiologically for treatment of T1D, a statistical and systematic approach was taken by analysing the pharmacokinetics and pharmacodynamics of single doses of subcutaneous aldesleukin in the Adaptive Study of IL-2 Dose on Regulatory T Cells in Type 1 Diabetes (DILT1D), a single centre, non-randomised, open label, adaptive dose-finding trial with 40 adult participants with recently diagnosed T1D. The primary endpoint was the maximum percentage increase in Tregs (defined as CD3+CD4+CD25highCD127low) from the baseline frequency in each participant measured over the 7 d following treatment. There was an initial learning phase with five pairs of participants, each pair receiving one of five pre-assigned single doses from 0.04 × 106 to 1.5 × 106 IU/m2, in order to model the dose-response curve. Results from each participant were then incorporated into interim statistical modelling to target the two doses most likely to induce 10% and 20% increases in Treg frequencies. Primary analysis of the evaluable population (n = 39) found that the optimal doses of aldesleukin to induce 10% and 20% increases in Tregs were 0.101 × 106 IU/m2 (standard error [SE] = 0.078, 95% CI = −0.052, 0.254) and 0.497 × 106 IU/m2 (SE = 0.092, 95% CI = 0.316, 0.678), respectively. On analysis of secondary outcomes, using a highly sensitive IL-2 assay, the observed plasma concentrations of the drug at 90 min exceeded the hypothetical Treg-specific therapeutic window determined in vitro (0.015–0.24 IU/ml), even at the lowest doses (0.040 × 106 and 0.045 × 106 IU/m2) administered. A rapid decrease in Treg frequency in the circulation was observed at 90 min and at day 1, which was dose dependent (mean decrease 11.6%, SE = 2.3%, range 10.0%–48.2%, n = 37), rebounding at day 2 and increasing to frequencies above baseline over 7 d. Teffs, natural killer cells, and eosinophils also responded, with their frequencies rapidly and dose-dependently decreased in the blood, then returning to, or exceeding, pretreatment levels. Furthermore, there was a dose-dependent down modulation of one of the two signalling subunits of the IL-2 receptor, the β chain (CD122) (mean decrease = 58.0%, SE = 2.8%, range 9.8%–85.5%, n = 33), on Tregs and a reduction in their sensitivity to aldesleukin at 90 min and day 1 and 2 post-treatment. Due to blood volume requirements as well as ethical and practical considerations, the study was limited to adults and to analysis of peripheral blood only. Conclusions The DILT1D trial results, most notably the early altered trafficking and desensitisation of Tregs induced by a single ultra-low dose of aldesleukin that resolves within 2–3 d, inform the design of the next trial to determine a repeat dosing regimen aimed at establishing a steady-state Treg frequency increase of 20%–50%, with the eventual goal of preventing T1D. Trial Registration ISRCTN Registry ISRCTN27852285; ClinicalTrials.gov NCT01827735

Multi-parametric flow cytometric and genetic investigation of the peripheral B cell compartment in human type 1 diabetes

The appearance of circulating islet‐specific autoantibodies before disease diagnosis is a hallmark of human type 1 diabetes (T1D), and suggests a role for B cells in the pathogenesis of the disease. Alterations in the peripheral B cell compartment have been reported in T1D patients; however, to date, such studies have produced conflicting results and have been limited by sample size. In this study, we have performed a detailed characterization of the B cell compartment in T1D patients (n = 45) and healthy controls (n = 46), and assessed the secretion of the anti‐inflammatory cytokine interleukin (IL)‐10 in purified B cells from the same donors. Overall, we found no evidence for a profound alteration of the B cell compartment or in the production of IL‐10 in peripheral blood of T1D patients. We also investigated age‐related changes in peripheral B cell subsets and confirmed the sharp decrease with age of transitional CD19+CD27−CD24hiCD38hi B cells, a subset that has recently been ascribed a putative regulatory function. Genetic analysis of the B cell compartment revealed evidence for association of the IL2–IL21 T1D locus with IL‐10 production by both memory B cells (P = 6·4 × 10−4) and islet‐specific CD4+ T cells (P = 2·9 × 10−3). In contrast to previous reports, we found no evidence for an alteration of the B cell compartment in healthy individuals homozygous for the non‐synonymous PTPN22 Trp620 T1D risk allele (rs2476601; Arg620Trp). The IL2–IL21 association we have identified, if confirmed, suggests a novel role for B cells in T1D pathogenesis through the production of IL‐10, and reinforces the importance of IL‐10 production by autoreactive CD4+ T cells.

Postthymic Expansion in Human CD4 Naive T Cells Defined by Expression of Functional High-Affinity IL-2 Receptors

As the thymus involutes with age, the maintenance of peripheral naive T cells in humans becomes strongly dependent on peripheral cell division. However, mechanisms that orchestrate homeostatic division remain unclear. In this study we present evidence that the frequency of naive CD4 T cells that express CD25 (IL-2 receptor α-chain) increases with age on subsets of both CD31+ and CD31− naive CD4 T cells. Analyses of TCR excision circles from sorted subsets indicate that CD25+ naive CD4 T cells have undergone more rounds of homeostatic proliferation than their CD25− counterparts in both the CD31+ and CD31− subsets, indicating that CD25 is a marker of naive CD4 T cells that have preferentially responded to survival signals from self-Ags or cytokines. CD25 expression on CD25− naive CD4 T cells can be induced by IL-7 in vitro in the absence of TCR activation. Although CD25+ naive T cells respond to lower concentrations of IL-2 as compared with their CD25− counterparts, IL-2 responsiveness is further increased in CD31− naive T cells by their expression of the signaling IL-2 receptor β-chain CD122, forming with common γ-chain functional high-affinity IL-2 receptors. CD25 plays a role during activation: CD25+ naive T cells stimulated in an APC-dependent manner were shown to produce increased levels of IL-2 as compared with their CD25− counterparts. This study establishes CD25+ naive CD4 T cells, which are further delineated by CD31 expression, as a major functionally distinct immune cell subset in humans that warrants further characterization in health and disease.

Dissection of a complex disease susceptibility region using a Bayesian stochastic search approach to fine mapping

Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r2 ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4+ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data.