Impact Factor 5.695

The #1 most cited and #1 largest open-access journal in Immunology

Review ARTICLE

Front. Immunol., 25 October 2012 | http://dx.doi.org/10.3389/fimmu.2012.00320

From biomarkers to a clue of biology: a computation-aided perspective of immune gene expression profiles in human type 1 diabetes

Dongmei Han1, Xiaodong Cai2, Ji Wen2, Norma S. Kenyon1,3,4 and Zhibin Chen1,5*
  • 1Diabetes Research Institute, University of Miami Miller School of Medicine, Miami, FL, USA
  • 2Department of Electrical and Computer Engineering, University of Miami, Coral Gables, FL, USA
  • 3Department of Medicine, University of Miami Miller School of Medicine, Miami, FL, USA
  • 4Department of Surgery, University of Miami Miller School of Medicine, Miami, FL, USA
  • 5Department of Microbiology and Immunology, University of Miami Miller School of Medicine, Miami, FL, USA

Dysregulated expression of key immune genes may cause breakdown of immunological tolerance and development of autoimmune disorders such as type 1 diabetes (T1D). General immune insufficiencies have also been implicated as a trigger of autoimmunity, due to their potential impact on immune homeostasis. Recent studies have detected evidence of systemic reduction in immune gene expression in long-term diabetic patients but the changes were not present before or at T1D onset. The changes could not be merely correlated with alteration in metabolic parameters. The studies also identified a dynamic expression pattern of several well-known as well as little-studied, immune-related genes during the course of T1D. An intriguing “ratio profile” of immune regulatory genes, such as CTLA4 and members of the S100 family, versus “baseline” immune genes, such as CD3G, prompted us to further examine immune gene expression relationships for a set of molecules representing T cells, B cells, and myeloid cells. No evidence was found to suggest an overall breach of tolerance equilibrium in T1D. Perplexingly, patients with long-term T1D presented a gene expression profile that was surprisingly more coordinated in analyses of “networking” relationship. Computational analyses of the “ratio profiles” or “relationship profiles” of immune gene expression might provide a clue for further studies of immunobiology in human T1D and other autoimmune diseases, as to how the profiles may be related to the pathogenic cause of the disease, to the effect of the diseases on immune homeostasis, or to an immunological process associated with the course of the diseases but is neither a direct cause nor a direct effect of the diseases.

Immunobiology and Immune Biomarkers in T1D

It is generally accepted that type 1 diabetes (T1D) develops as a result of breakdown in immunological tolerance induction mechanisms. The contribution of various tolerance mechanisms has been clearly delineated with the aid of animal models of T1D. It is postulated that T1D and other autoimmune diseases are likely caused by specific breaches in a network of immunological tolerance mechanisms (Mathis and Benoist, 2010). However, a general defect at a system level has also been implicated by clinical and experimental evidence (Gallegos and Bevan, 2004; King et al., 2004; Schuetz et al., 2010). Indeed, studies with animal models suggested that autoimmune diabetes could be caused by either specific defects in immune dysregulation (for review, see Lehuen et al., 2010), or a general immune insufficiency (lymphopenia; King et al., 2004). However, it is difficult to pinpoint the breakdown of immunological tolerance in human patients. In common cases of human T1D, the status of the immune system remains a vaguely delineated framework that awaits characterization at molecular and cellular levels.

In the mean time, studies of biomarkers that are associated with T1D are expected to provide important tools for clinical management of the disease. Biomarkers can be derived from a broad spectrum of factors (Purohit and She, 2008), for example, T1D-associated genetic loci, antibodies against endocrine pancreatic tissues and products, potentially pathogenic or protective cytokines, and expression profiles of immunological or metabolic genes. The best examples are perhaps the predicative values of anti-insulin and anti-islet-cell antibodies for T1D development (Eisenbarth et al., 1998). Another well-recognized example is the measurement of hemoglobin A1c (HbA1c), which is commonly used for clinical diabetes management. Since immunological processes are believed to play critical roles in all stages of T1D development, immune gene expression profiles can potentially be used as biomarkers for staging T1D development as well as gauging the impact of the disease.

Microarray gene expression profiles of peripheral blood provided new insights into pathogenesis of T1D (Reynier et al., 2010). Such studies could reveal differences in immune responsiveness between patients with T1D and healthy controls, and may identify changes in gene expression that associate with progression of T1D. Studies of the peripheral immune system of new-onset T1D patients have shown significantly higher levels of IL-1α, IFN-γ, and TNF-α as compared to normal controls (Hussain et al., 1996; Kallmann et al., 1997). Patients with childhood-onset T1D show abnormal monocyte gene expression levels with an altered gene expression network, implicating monocyte abnormalities in susceptibility to diabetes (Beyan et al., 2010). Two distinct gene expression profiles have been reported in monocytes from peripheral blood: a proinflammatory profile mainly associated to adult-onset and a chemotaxis, adhesion, and metabolism profile mainly associated to juvenile-onset diabetes (Padmos et al., 2008). Recently, a study showed that healthy first-degree relatives of patients with T1D exhibited significant differences in expression pattern of genes involved in the regulation of innate immune responses such as TLR signaling and CCR3 signaling in eosinophils, costimulation, and cytokine responses mediated by CD137, CD40, and CD28 signaling and IL-1 proinflammatory pathway, as compared to healthy controls (Stechova et al., 2011).

One of the major subsets in the adoptive immune system, B lymphocytes, have been implicated to play a role in T1D development and disease progression (Serreze et al., 1998; Silveira and Grey, 2006), but the actual role that B cells play remains to be elucidated (Serreze et al., 1998; Wong et al., 2004). Several studies from animal models of autoimmune disease have shown that B cells can serve as antigen-presenting cells and prime β cell-specific T cells (Falcone et al., 1998; Serreze et al., 1998; Bouaziz et al., 2007). Data from a recent clinical trial demonstrated that treatment of new-onset T1D patients with anti-CD20 (Rituxan) preserved residual insulin secretion for at least 1 year (O’Neill et al., 2009; Pescovitz et al., 2009). B cells may also contribute to immune responses through the secretion of effector cytokines (Harris et al., 2005a,b; Duddy et al., 2007). Gene expression levels of two key B cell markers, activation-induced cytidine deaminase and immunoglobulin G gamma chain, were found significantly lower in long-term T1D patients as compared to healthy controls or new-onset patients (Han et al., 2011). The cytotoxic lymphocyte gene products granzyme B (GB), perforin, and fas ligand (FasL) have been shown to play a vital part in the T1D development (Kreuwel et al., 1999; Moriwaki et al., 1999; Kreuwel and Sherman, 2001; Yoon and Jun, 2001). mRNA levels of perforin and FasL genes were found significantly lower in patients with long-term T1D as compared to healthy controls (Giordano et al., 1995; Han et al., 2005).

It can perhaps be agreed that there are a plethora of examples of conflicting reports in immune-related changes in the course of T1D. For example, conflicting evidence exists on whether an active Th1-like immune response destroys beta cells, followed by presentation of autoantigens during the prediabetic phase (Karlsson Faresjo et al., 2004; Hedman et al., 2006). Th17 cells have been shown to play a critical role in the induction of autoimmune tissue injury and inflammation, and might be involved in exacerbation of diabetes, but its role in T1D remains to be clarified (Cooke, 2006; Bending et al., 2009; Emamaullee et al., 2009; Martin-Orozco et al., 2009; Honkanen et al., 2010). Increased levels of IL-17 secreting T cells were found in children with new-onset T1D (Marwaha et al., 2010). There are also conflicting reports about the levels of cytokine secretion and cytokine mRNA expression in at-risk, new-onset, and long-term T1D patients (Halminen et al., 2001; Nicoletti et al., 2002). These discrepancies might be contributed by different stages of the disease even within the same group, or sampling variations.

Another example is CTLA4, an immunoregulatory molecule that plays a key role in negatively regulating T cell responses (for review, see Teft et al., 2006). Genetic variations in the CTLA4 locus is associated with a number of autoimmune diseases (for review, see Kristiansen et al., 2000). Many groups studied potential dysregulation of CTLA4 expression in T1D. The mRNA expression level for CTLA4, as well as ICOS and GITR, was found lower in regulatory T (Treg) cells of children with newly diagnosed diabetes as compared to the healthy controls (Luczynski et al., 2009). CTLA4 protein expression was reported lower in the patients with diabetes (Haseda et al., 2011; Ryden et al., 2012) and with autoimmune thyroiditis as compared in controls (Kucharska et al., 2009). The expression of CTLA4 is influenced by genetic polymorphisms, although it remains debated how exactly each of CTLA4 polymorphisms impact human T1D genetics and how they contribute to CTLA4 expression variation (Ueda et al., 2003; Anjos et al., 2005). It should be noted that the reported changes in CTLA4 expression were subtle. A subtle reduction of CTLA4 could indeed impact T1D development (Chen et al., 2006), with mechanism remain to be elucidated. On the other hand, induction of CTLA4 is associated with tolerogenic effect of a therapeutic agent, the murine analog of anti-thymocyte globulin (ATG), in a transgenic mouse model of T1D (Lu et al., 2011). In experimental settings of CTLA4-deficiencies, i.e., “all-or-nothing” modeling, or anti-CTLA4-antibody-mediated blockade, it has been shown that CTLA4 acts on both Treg and effector T (Teff) cells (Wing et al., 2008; Peggs et al., 2009; Ise et al., 2010; Jain et al., 2010; Pruitt et al., 2011; Miska et al., 2012). However, in a recent study of antitumor immunity, Teff cells exhibited unique sensitivity to subtle reduction of CTLA4 (Miska et al., 2012).

In our studies to identify potential biomarkers in association with T1D progression, gene expression analyses were performed with quantitative RT-PCR for the mRNA levels of a set of immune-related genes using commercially tested primer and probe sets (Han et al., 2011, 2012). In this initial stage of studies, whole blood samples from at-risk, new-onset, and long-term T1D patients, as well as healthy controls were preserved for analyses after collection, without further processing or fractioning, to avoid any potential loss and changes caused by processing. It was found that IFN-γ, IL-4, and IL-10 mRNA levels were significantly higher in new-onset as compared to at-risk and long-term T1D patients (Han et al., 2011). The gene expression levels of most cytokines and effector molecules were suppressed in long-term T1D patients as compared to healthy controls (Han et al., 2011). Surprisingly, CTLA4 expression levels per se were not changed in the at-risk or new-onset stages, but reduced in long-term diabetic patients (Han et al., 2012). Unexpectedly, we found a significant difference between healthy controls and T1D groups in mRNA levels for “baseline” immune gene such as CD3G (representing T cells), CD20 (representing B cells), and CD11b (representing myeloid cells). When CTLA4 expression was examined in reference to CD3G expression, as a ratio of CTLA4/CD3G, we detected a decreased ratio in the sample from at-risk and new-onset patients but an increased ratio in samples from long-term T1D patients. The distinct “ratio profiles” for various immune regulatory genes (Han et al., 2012) prompted us to further examine the relationship of immune gene expression in different stages of T1D.

“Relationship Profile” of Innate and Adaptive Immune Gene Expression

It is thought that “cross-talk” between different subsets of immune cells is crucial for immunological tolerance (Lehuen et al., 2010; Suzuki et al., 2010), but it remains a challenge to study immune cell interaction in human patients. The gene expression values were obtained for T lymphocyte genes (CD3G and CTLA4), B lymphocyte genes (CD19 and CD20), and myeloid cell-related genes (CD11b, TLR9, and ARG1), as well as a subset of members of the S100 family that has been implicated in immune regulation (Han et al., 2012). We probed the relationship of these values, by a computational approach, to gain a clue to the systemic and regulatory relationships among T-, B-, and myeloid cells. Whereas the well-established immunological relationships between CD3G and CTLA4, as well as CD19 and CD20, should validate the experimental and computational methods of gene expression relationship and network analyses, the inclusion of the S100 family members in the expression analyses (Han et al., 2012) attests the approaches’ utility to explore novel relationships between a gene with a well-defined immunological function and another gene with a yet-to-be characterized role in immunity and tolerance.

Gene expression networks were analyzed by calculating Pearson’s correlation coefficient among the expression values of each pair of genes, measured in a previously reported study (Han et al., 2012). Correlation between any two genes with an absolute coefficient value greater than 0.8 was presented in gene expression network graphs (Figure 1). A robust correlation between four pairs of genes was expected in the healthy samples, based on the established function of these genes: (1) CTLA4 and CD3G; (2) CD19 and CD20; (3) CD11b and TLR9; and (4) S100A8 and S100A9. The anticipated relationships were designed to be a “positive” control for the computational method. Indeed, strong correlations of expression between all four pairs of genes were validated in the healthy control samples (Figure 1A). The correlation between the T cell pair and the B cell pair appeared to be relatively weak (<0.8) in the healthy group (Figure 1A). The sample size in the at-risk group (n = 19; Han et al., 2012) may be too small to estimate correlation coefficients reliably. Therefore, the at-risk group was not included in this computational analysis. In the healthy control group, one sample was discarded from this analysis because it contained statistical outliers (n = 69; Han et al., 2012).

FIGURE 1
www.frontiersin.org

FIGURE 1. Gene expression network analysis to examine potential coordination among innate and adaptive immune gene expression. Quantitative RT-PCR was used to assess the levels of mRNA expression for thirteen innate (CD11b, TLR9, and S100 family) and adaptive (CD3, CTLA4, CD19, and CD20) immune genes. Gene expression relationships were analyzed for healthy controls (A) and compared to that of new-onset (B) and long-term T1D patients (C). Solid red lines, absolute value of correlation coefficient ≥0.9; dashed blue lines, absolute value of correlation coefficient ≥0.8 but <0.9.

The gene expression network was constructed using a standard formula for Pearson’s correlation coefficient. The computation analysis with this formula does not predict how the correlation coefficients of pairs of genes in one subject group would change when the average expression levels of a set of genes are all reduced, for example, whether the absolute values of the correlation coefficients would increase, decrease, or remain similar. Intuitively, we expected that the onset and duration of T1D might be associated with a disruption of a gene expression correlation network that exists in the healthy controls. Contrary to our expectation, we did not detect an apparent disruption of the gene expression correlations in the samples from the new-onset diabetes group (n = 33; Figure 1B). Instead, it was clear that the basic correlation pattern was preserved in the new-onset group as compared to the HT controls. Surprisingly, when compared to healthy controls (Figure 1A), the gene expression correlation network was much strengthened overall in long-term diabetic patients (n = 59; Figure 1C).

As expected, the expression of several S100 genes correlated with the mRNA level of CD11b, an indicator of myeloid cells. Novel relationships of gene expressions were identified in the samples from the healthy group between S100 family genes and T cell genes (CTLA4 and CD3G), but not with B cell genes (CD20 and CD19). Surprisingly, although healthy controls and new-onset T1D patients exhibited a similar pattern of gene expression relationships, more correlative relationships of gene expression were identified in the long-term diabetic group than in healthy controls, such that the expression of CD3G, CD11b, S100A6, S100A9, S100A10, S100A11, or S100A13 correlated with a majority of the other genes analyzed in the study, regardless of their primary functional denotation in T-, B-, or myeloid cells. Distinct gene expression relationships also emerged in the samples from the new-onset and long-term T1D groups between the B cell gene CD19 and TLR9, or between CD19 and members of S100 gene family (Figures 1B,C).

The correlation coefficients were further used as a distance metric in a hierarchical clustering analysis that yielded a dendrogram, a commonly used method for depiction of relationship analyses. The healthy, new-onset and long-term T1D groups generally exhibited similar clusters (Figure 2). Thus, the results from the hierarchical clustering analyses support the finding from network analyses, suggesting that the expression relationships of the selected immune genes were well-preserved during progression of T1D.

FIGURE 2
www.frontiersin.org

FIGURE 2. Hierarchical clustering analysis of expression relationship among innate and adaptive immune genes. A dendrogram was presented based on correlation coefficient of gene expression levels. The HT, NO, and LT groups exhibited a similar pattern of gene expression correlation. Most of the 13 genes appeared in two clusters: a B cell cluster and a cluster consisting of both T cell and myeloid cell genes. The sample size in the AT group (n = 19) may not be sufficient for the clustering technique.

Cross-Group Analysis of Gene–Gene Interactions

For experimental biologists to identify gene expression relationship, the most straightforward method is perhaps to simply analyze the ratio of gene expression levels for genes that are known to be related (Han et al., 2012), as described above. Statistical analyses of correlation coefficients, as described in Figures 1 and 2, are also intuitively comprehensible and desirable by experimental biologists to detect potential relationship. However, other types of interactive relationships could exist and may necessitate more sophisticated computation methods. In genetic association studies, interactions of several genetic loci were found, in addition to the other well-established loci such as HLA, to be associated with T1D (Cordell et al., 1995; Ide et al., 2004). Gene–gene interactions are believed to be an important factor that may account for missing heritability for T1D (Pociot et al., 2010) and other complex human diseases (Manolio et al., 2009). We therefore employed a computational approach to identify gene–gene interactions that were associated with a T1D group vis-à-vis the healthy controls, using the gene expression values.

When we compared a T1D group with the healthy group, we used an indicator variable Y to indicate whether a subject is from the T1D group (Y = 1) or the healthy group (Y = 0), and then employed the following logistic regression model to investigate the effect of gene expression and the age on the status of the subject,

log(P(Y=1)P(Y=0))=β0+Σi=113Xiβi+Σi=112Σj=i+113XiXjβij+age×βa,

where P(Y = 1) or P(Y = 0) is the probability that the subject belongs to the T1D or healthy group, Xi is the expression level of gene i, and βi, βij, βa are regression coefficients to be determined. This logistic regression model is similar to the one used to detect gene–gene interactions in genome-wide association analysis (Cordell, 2009), except the variables here are gene expression values and the age instead of genotypes. Since the number of variables in the model is greater than the number of samples available, the traditional maximum likelihood approach cannot be used to determine the regression coefficients. We thus employed a method named elastic net (Zou and Trevor, 2005), which is capable of selecting relevant variables from a set of a large number of variables. The elastic net was originally used to predict if a leukemia patient has type 1 or 2 leukemia using the expression levels of 7129 genes (Zou and Trevor, 2005). Specifically, the elastic net determined the values of regression coefficients by maximizing the following penalized likelihood function,

L(y|β)λ(1α)2(Σi=113βi2+Σi=112Σj=i+113βij2)λα(Σi=113|βi|+Σi=112Σj=i+113|βij|),

where L(y|β) is the log likelihood function of the data with β standing for all regression coefficients, α (0 < α ≤ 1) and λ (>0) are two constants that could be determined with cross-validation. We used an efficient program named glmnet (Friedman et al., 2010) that implemented the elastic net to fit the data to the logistic regression model. The constants α and λ were chosen from leave-one-out cross-validation that yielded the smallest deviance. The standard errors of non-zero coefficients obtained from glmnet were then calculated from the sandwich formula for the penalized regression models (Fan and Li, 2001).

In the analysis for LT versus HT group, cross-validation yielded α = 0.71 and λ = 0.01393. Hosmer–Lemeshow test for the goodness of fit of the model gave a p-value of 0.567 which implies that the data fitted the model properly. The area under the receiver operating characteristic (ROC) curve is 0.831, which also indicates that the data fitted the model well. The regression coefficients obtained from the glmnet were presented in Table 1. Three pair-wise gene interactions were identified associating with long-term T1D vis-à-vis healthy controls: ARG1 versus S100A6, ARG1 versus S100A9, and CD3G versus CD20, since their corresponding regression coefficient is ≠0 with a statistical significance <0.05. To illustrate the effect of gene-gene interactions, let us denote the probability of a subject belonging to the long-term T1D group as p and correspondingly the probability of subject being in the healthy group as 1 - p. The odds for the subject being associated with the long-term T1D group is then o = p/(1 - p). Taking the analysis of CD3G versus CD20 for illustration, the regression coefficient is 0.21 (Table 1).

TABLE 1
www.frontiersin.org

TABLE 1. Elastic net regression analyses of gene–gene interaction between the HT and LT-T1D groups.

In a hypothesized situation where the product of expression levels of CD3G and CD20 is increased by onefold but all other expression levels do not change, the odds of the subject being associated with long-term T1D is denoted in this situation as o1 and the odds ratio is defined as o1/o. The odds ratio o1/o = e0.21 = 1.23 implies that the relative likelihood of association with the long-term T1D group is increased by 1.23 times if the product of expression levels of CD3G and CD20 is increased by onefold under this hypothesized situation. Of course, this kind of hypothetical situation unlikely occurs, because (1) if the product of expression levels of CD3G and CD20 changes, most likely expression levels of CD3G and CD20 themselves change too, and (2) as we have seen expression levels of some genes are highly correlated; meaning, if expression levels of CD3G and CD20 change, expression levels of some other genes may change too. Therefore, the number does not really specify exact relative risk, but rather suggests that a characteristic CD3G and CD20 interaction indicates a biased association with the long-term T1D versus healthy condition. The immunological significance of this finding, as to how it is related to T1D, remains to be tested. It might suggest, although purely speculative at this point, an altered T–B collaboration in an established diabetes condition. Collaboration of T- and B-lymphocytes is critical for a productive immune response and its regulation.

Notably, this penalized regression model identified only one main effect of single genes between long-term T1D patients and healthy individuals (S100A13, p < 0.0005; Table 1), although both Student’s t test and linear regression analyses involving single genes and age indicated significant reduction of expression in 10 of the 13 genes (Han et al., 2012). This was due to the fact that the correlations among genes were high and the multiple penalized regression analysis only picked out the gene whose difference between the two groups was most significant. This approach did not identify any gene–gene interaction in cross-group analyses between HT versus AT, or HT versus NO.

Summary

A large on-going effort has been devoted by many groups to identify biomarkers for T1D and other autoimmune diseases. Undoubtedly, those biomarkers, including differential gene expression profiles, will be instrumental to improve the clinical management of the diseases. For some autoimmune diseases, T1D as a well-known example, even though diagnosis criteria are straightforward, biomarkers may prove to be a great aid in staging the disease as well as predicting disease development before and after onset. Given the immune origin of an autoimmune disease, immune gene expression profiles are logical candidates of biomarkers. Such profiles may or may not necessarily be directly related to the disease cause or effect, but regardless may be useful indicators for the status of the immune system, the “immunological wellness,” during the disease progression (Han et al., 2012).

Successful prevention or cure of T1D and other autoimmune diseases will require a better understanding of the mechanisms of the diseases. In this regard, substantial progress has been made. Multiple genetic risks and environmental factors are implicated in the immunological tolerance breakdown that leads to T1D. Rare cases of T1D caused by Mendelian genetic mutations at either the AIRE or FOXP3 locus highlight the critical role of central tolerance by thymic deletion or peripheral tolerance by regulatory T lymphocytes, respectively (Mathis and Benoist, 2010). On the other hand, a surprisingly limited spectrum of autoimmune damage in Aire/Foxp3 double-deficient mice suggests robustness of an immunological tolerance framework reinforced by a network of known and unknown elements (Chen et al., 2005). Clinical observations and experimental studies suggest that a general immune insufficiency may cause autoimmune diseases (Gallegos and Bevan, 2004; King et al., 2004; Schuetz et al., 2010). For a T1D population consisting of common cases that are not caused by monoallelic genetic mutations, however, the overall immunological status during T1D initiation and progression remains a challenge to characterize, due to a lack of experimental approaches for a system level assessment. A daunting task remains to study the immunobiology and pathophysiology behind human T1D and other autoimmune diseases, due to obvious clinical and logistic limitations and ethical concerns.

Perhaps, it is fair to critique that the computational analyses would be difficult to interpret in the absence of experimental evidence, and therefore may not be suitable for an experimental journal unless experimental evidence is provided. Currently available experimental biology approaches are not handy yet to definitively tackle human immunology at a system level in a complex disease setting such as T1D. Subtle and specific imbalances in innate and adaptive immune regulation, e.g., CTLA4/CD3 ratios, could be detected and may be associated with the T1D development. Systemic perturbations might occur, apparently not as a cause of T1D but as a part associated with the T1D course (Han et al., 2012). Of course, the computational analysis of the gene expression relationships could be regarded as just another approach by another “blind man” to the disease “elephant.” However, we argue that it may provide a clue to the biology behind the disease, suggesting that T1D pathogenesis in humans is not due to a gross encroachment of the integrity of immune regulatory network, nor to a general insufficiency in the immune system. For a reason(s) yet to be determined, a strengthened correlation of immune gene expression was associated with T1D progression. It remains to be understood whether and how the altered relationships may impact immune function and immune tolerance induction. A tightened immune regulatory network might facilitate tolerance induction. On the other hand, an increased rigidity of a network might also affect its responsiveness. In that vein, the general insufficiency of immune gene expression in long-term diabetes, together with a tightened network of immune gene expression, might pose a secondary risk of suboptimal immunity, and a perturbation in the homeostasis of the immune system.

Conflict of Interest Statement

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.

Acknowledgments

This work was supported by the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK, award DP3DK085696). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIDDK or NIH. Xiaodong Cai was partly supported by the National Science Foundation (award 0746882. Stochastic Modeling, Analysis and Simulation of Gene Networks). The authors gratefully acknowledge the expert support of Drs J. S. Skyler, A. Pugliese, R. Alejandro, D. Mineo, and D. M. Berman and Ms. D. Matheson.

References

Anjos, S. M., Shao, W., Marchand, L., and Polychronakos, C. (2005). Allelic effects on gene regulation at the autoimmunity-predisposing CTLA4 locus: a re-evaluation of the 3′ +6230G>A polymorphism. Genes Immun. 6, 305–311.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bending, D., De La Pena, H., Veldhoen, M., Phillips, J. M., Uyttenhove, C., Stockinger, B., et al. (2009). Highly purified Th17 cells from BDC2.5NOD mice convert into Th1-like cells in NOD/SCID recipient mice. J. Clin. Invest. 119, 565–572.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Beyan, H., Drexhage, R. C., van der Heul Nieuwenhuijsen, L., de Wit, H., Padmos, R. C., Schloot, N. C., et al. (2010). Monocyte gene-expression profiles associated with childhood-onset type 1 diabetes and disease risk: a study of identical twins. Diabetes 59, 1751–1755.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bouaziz, J. D., Yanaba, K., Venturi, G. M., Wang, Y., Tisch, R. M., Poe, J. C., et al. (2007). Therapeutic B cell depletion impairs adaptive and autoreactive CD4+ T cell activation in mice. Proc. Natl. Acad. Sci. U.S.A. 104, 20878–20883.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Chen, Z., Benoist, C., and Mathis, D. (2005). How defects in central tolerance impinge on a deficiency in regulatory T cells. Proc. Natl. Acad. Sci. U.S.A. 102, 14735–14740.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Chen, Z., Stockton, J., Mathis, D., and Benoist, C. (2006). Modeling CTLA4-linked autoimmunity with RNA interference in mice. Proc. Natl. Acad. Sci. U.S.A. 103, 16400–16405.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cooke, A. (2006). Th17 cells in inflammatory conditions. Rev. Diabet. Stud. 3, 72–75.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cordell, H. J. (2009). Detecting gene–gene interactions that underlie human diseases. Nat. Rev. Genet. 10, 392–404.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cordell, H. J., Todd, J. A., Bennett, S. T., Kawaguchi, Y., and Farrall, M. (1995). Two-locus maximum lod score analysis of a multifactorial trait: joint consideration of IDDM2 and IDDM4 with IDDM1 in type 1 diabetes. Am. J. Hum. Genet. 57, 920–934.

Pubmed Abstract | Pubmed Full Text

Duddy, M., Niino, M., Adatia, F., Hebert, S., Freedman, M., Atkins, H., et al. (2007). Distinct effector cytokine profiles of memory and naive human B cell subsets and implication in multiple sclerosis. J. Immunol. 178, 6092–6099.

Pubmed Abstract | Pubmed Full Text

Eisenbarth, G. S., Gianani, R., Yu, L., Pietropaolo, M., Verge, C. F., Chase, H. P., et al. (1998). Dual-parameter model for prediction of type I diabetes mellitus. Proc. Assoc. Am. Physicians 110, 126–135.

Pubmed Abstract | Pubmed Full Text

Emamaullee, J. A., Davis, J., Merani, S., Toso, C., Elliott, J. F., Thiesen, A., et al. (2009). Inhibition of Th17 cells regulates autoimmune diabetes in NOD mice. Diabetes 58, 1302–1311.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Falcone, M., Lee, J., Patstone, G., Yeung, B., and Sarvetnick, N. (1998). B lymphocytes are crucial antigen-presenting cells in the pathogenic autoimmune response to GAD65 antigen in nonobese diabetic mice. J. Immunol. 161, 1163–1168.

Pubmed Abstract | Pubmed Full Text

Fan, J., and Li, R. (2001). Variable selection via nonconcave penalized likelihood and its oracle properties. J. Am. Stat. Assoc. 96, 1348–1360.

Friedman, J., Hastie, T., and Tibshirani, R. (2010). Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33, 1–22.

Pubmed Abstract | Pubmed Full Text

Gallegos, A. M., and Bevan, M. J. (2004). Driven to autoimmunity: the nod mouse. Cell 117, 149–151.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Giordano, C., De Maria, R., Stassi, G., Todaro, M., Richiusa, P., Giordano, M., et al. (1995). Defective expression of the apoptosis-inducing CD95 (Fas/APO-1) molecule on T and B cells in IDDM. Diabetologia 38, 1449–1454.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Halminen, M., Simell, O., Knip, M., and Ilonen, J. (2001). Cytokine expression in unstimulated PBMC of children with type 1 diabetes and subjects positive for diabetes-associated autoantibodies. Scand. J. Immunol. 53, 510–513.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Han, D., Cai, X., Wen, J., Matheson, D., Skyler, J. S., Kenyon, N. S., et al. (2012). Innate and adaptive immune gene expression profiles as biomarkers in human type 1 diabetes. Clin. Exp. Immunol. 170, 131–138.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Han, D., Leith, J., Alejandro, R., Bolton, W., Ricordi, C., and Kenyon, N. S. (2005). Peripheral blood cytotoxic lymphocyte gene transcript levels differ in patients with long-term type 1 diabetes compared to normal controls. Cell Transplant. 14, 403–409.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Han, D., Leyva, C. A., Matheson, D., Mineo, D., Messinger, S., Blomberg, B. B., et al. (2011). Immune profiling by multiple gene expression analysis in patients at-risk and with type 1 diabetes. Clin. Immunol. 139, 290–301.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Harris, D. P., Goodrich, S., Gerth, A. J., Peng, S. L., and Lund, F. E. (2005a). Regulation of IFN-gamma production by B effector 1 cells: essential roles for T-bet and the IFN-gamma receptor. J. Immunol. 174, 6781–6790.

Pubmed Abstract | Pubmed Full Text

Harris, D. P., Goodrich, S., Mohrs, K., Mohrs, M., and Lund, F. E. (2005b). Cutting edge: the development of IL-4-producing B cells (B effector 2 cells) is controlled by IL-4, IL-4 receptor alpha, and Th2 cells. J. Immunol. 175, 7103–7107.

Pubmed Abstract | Pubmed Full Text

Haseda, F., Imagawa, A., Murase-Mishiba, Y., Sano, H., Hirano-Kuwata, S., Ueda, H., et al. (2011). Low CTLA-4 expression in CD4+ helper T-cells in patients with fulminant type 1 diabetes. Immunol. Lett. 139, 80–86.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hedman, M., Ludvigsson, J., and Faresjo, M. K. (2006). Nicotinamide reduces high secretion of IFN-gamma in high-risk relatives even though it does not prevent type 1 diabetes. J. Interferon Cytokine Res. 26, 207–213.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Honkanen, J., Nieminen, J. K., Gao, R., Luopajarvi, K., Salo, H. M., Ilonen, J., et al. (2010). IL-17 immunity in human type 1 diabetes. J. Immunol. 185, 1959–1967.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hussain, M. J., Peakman, M., Gallati, H., Lo, S. S., Hawa, M., Viberti, G. C., et al. (1996). Elevated serum levels of macrophage-derived cytokines precede and accompany the onset of IDDM. Diabetologia 39, 60–69.

Pubmed Abstract | Pubmed Full Text

Ide, A., Kawasaki, E., Abiru, N., Sun, F., Kobayashi, M., Fukushima, T., et al. (2004). Association between IL-18 gene promoter polymorphisms and CTLA-4 gene 49A/G polymorphism in Japanese patients with type 1 diabetes. J. Autoimmun. 22, 73–78.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ise, W., Kohyama, M., Nutsch, K. M., Lee, H. M., Suri, A., Unanue, E. R., et al. (2010). CTLA-4 suppresses the pathogenicity of self antigen-specific T cells by cell-intrinsic and cell-extrinsic mechanisms. Nat. Immunol. 11, 129–135.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Jain, N., Nguyen, H., Chambers, C., and Kang, J. (2010). Dual function of CTLA-4 in regulatory T cells and conventional T cells to prevent multiorgan autoimmunity. Proc. Natl. Acad. Sci. U.S.A. 107, 1524–1528.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kallmann, B. A., Huther, M., Tubes, M., Feldkamp, J., Bertrams, J., Gries, F. A., et al. (1997). Systemic bias of cytokine production toward cell-mediated immune regulation in IDDM and toward humoral immunity in Graves’ disease. Diabetes 46, 237–243.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Karlsson Faresjo, M. G., Ernerudh, J., and Ludvigsson, J. (2004). Cytokine profile in children during the first 3 months after the diagnosis of type 1 diabetes. Scand. J. Immunol. 59, 517–526.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

King, C., Ilic, A., Koelsch, K., and Sarvetnick, N. (2004). Homeostatic expansion of T cells during immune insufficiency generates autoimmunity. Cell 117, 265–277.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kreuwel, H. T., Morgan, D. J., Krahl, T., Ko, A., Sarvetnick, N., and Sherman, L. A. (1999). Comparing the relative role of perforin/granzyme versus Fas/Fas ligand cytotoxic pathways in CD8+ T cell-mediated insulin-dependent diabetes mellitus. J. Immunol. 163, 4335–4341.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kreuwel, H. T., and Sherman, L. A. (2001). The role of Fas–FasL in CD8+ T-cell-mediated insulin-dependent diabetes mellitus (IDDM). J. Clin. Immunol. 21, 15–18.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kristiansen, O. P., Larsen, Z. M., and Pociot, F. (2000). CTLA-4 in autoimmune diseases – a general susceptibility gene to autoimmunity? Genes Immun. 1, 170–184.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kucharska, A. M., Gorska, E., Wasik, M., Pyrzak, B., and Demkow, U. (2009). Expression of CD152 (CTLA-4) in children with autoimmune thyroiditis and +49 A/G polymorphism of exon 1 of the CTLA-4 gene. J. Physiol. Pharmacol. 60(Suppl. 5), 77–80.

Pubmed Abstract | Pubmed Full Text

Lehuen, A., Diana, J., Zaccone, P., and Cooke, A. (2010). Immune cell crosstalk in type 1 diabetes. Nat. Rev. Immunol. 10, 501–513.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lu, Y., Suzuki, J., Guillioli, M., Umland, O., and Chen, Z. (2011). Induction of self-antigen-specific Foxp3(+ regulatory T cells in the periphery by lymphodepletion treatment with anti-mouse thymocyte globulin in mice. Immunology 134, 50–59.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Luczynski, W., Wawrusiewicz-Kurylonek, N., Stasiak-Barmuta, A., Urban, R., Ilendo, E., Urban, M., et al. (2009). Diminished expression of ICOS, GITR and CTLA-4 at the mRNA level in T regulatory cells of children with newly diagnosed type 1 diabetes. Acta Biochim. Pol. 56, 361–370.

Pubmed Abstract | Pubmed Full Text

Manolio, T. A., Collins, F. S., Cox, N. J., Goldstein, D. B., Hindorff, L. A., Hunter, D. J., et al. (2009). Finding the missing heritability of complex diseases. Nature 461, 747–753.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Martin-Orozco, N., Chung, Y., Chang, S. H., Wang, Y. H., and Dong, C. (2009). Th17 cells promote pancreatic inflammation but only induce diabetes efficiently in lymphopenic hosts after conversion into Th1 cells. Eur. J. Immunol. 39, 216–224.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Marwaha, A. K., Crome, S. Q., Panagiotopoulos, C., Berg, K. B., Qin, H., Ouyang, Q., et al. (2010). Cutting edge: Increased IL-17-secreting T cells in children with new-onset type 1 diabetes. J. Immunol. 185, 3814–3818.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mathis, D., and Benoist, C. (2010). Levees of immunological tolerance. Nat. Immunol. 11, 3–6.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Miska, J., Bas, E., Devarajan, P., and Chen, Z. (2012). Autoimmunity-mediated antitumor immunity: tumor as an immunoprivileged self. Eur. J. Immunol. 42, 2584–2596.

Pubmed Abstract | Pubmed Full Text

Moriwaki, M., Itoh, N., Miyagawa, J., Yamamoto, K., Imagawa, A., Yamagata, K., et al. (1999). Fas and Fas ligand expression in inflamed islets in pancreas sections of patients with recent-onset Type I diabetes mellitus. Diabetologia 42, 1332–1340.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Nicoletti, F., Conget, I., Di Mauro, M., Di Marco, R., Mazzarino, M. C., Bendtzen, K., et al. (2002). Serum concentrations of the interferon-gamma-inducible chemokine IP-10/CXCL10 are augmented in both newly diagnosed Type I diabetes mellitus patients and subjects at risk of developing the disease. Diabetologia 45, 1107–1110.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

O’Neill, S. K., Liu, E., and Cambier, J. C. (2009). Change you can B(cell)eive in: recent progress confirms a critical role for B cells in type 1 diabetes. Curr. Opin. Endocrinol. Diabetes Obes. 16, 293–298.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Padmos, R. C., Schloot, N. C., Beyan, H., Ruwhof, C., Staal, F. J., de Ridder, D., et al. (2008). Distinct monocyte gene-expression profiles in autoimmune diabetes. Diabetes 57, 2768–2773.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Peggs, K. S., Quezada, S. A., Chambers, C. A., Korman, A. J., and Allison, J. P. (2009). Blockade of CTLA-4 on both effector and regulatory T cell compartments contributes to the antitumor activity of anti-CTLA-4 antibodies. J. Exp. Med. 206, 1717–1725.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pescovitz, M. D., Greenbaum, C. J., Krause-Steinrauf, H., Becker, D. J., Gitelman, S. E., Goland, R., et al. (2009). Rituximab, B-lymphocyte depletion, and preservation of beta-cell function. N. Engl. J. Med. 361, 2143–2152.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pociot, F., Akolkar, B., Concannon, P., Erlich, H. A., Julier, C., Morahan, G., et al. (2010). Genetics of type 1 diabetes: what’s next? Diabetes 59, 1561–1571.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Pruitt, S. K., Boczkowski, D., de Rosa, N., Haley, N. R., Morse, M. A., Tyler, D. S., et al. (2011). Enhancement of anti-tumor immunity through local modulation of CTLA-4 and GITR by dendritic cells. Eur. J. Immunol. 41, 3553–3563.

Pubmed Abstract | Pubmed Full Text

Purohit, S., and She, J. X. (2008). Biomarkers for type 1 diabetes. Int. J. Clin. Exp. Med. 1, 98–116.

Pubmed Abstract | Pubmed Full Text

Reynier, F., Pachot, A., Paye, M., Xu, Q., Turrel-Davin, F., Petit, F., et al. (2010). Specific gene expression signature associated with development of autoimmune type-I diabetes using whole-blood microarray analysis. Genes Immun. 11, 269–278.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ryden, A., Bolmeson, C., Jonson, C. O., Cilio, C. M., and Faresjo, M. (2012). Low expression and secretion of circulating soluble CTLA-4 in peripheral blood mononuclear cells and sera from type 1 diabetic children. Diabetes Metab. Res. Rev. 28, 84–96.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schuetz, C., Niehues, T., Friedrich, W., and Schwarz, K. (2010). Autoimmunity, autoinflammation and lymphoma in combined immunodeficiency (CID). Autoimmun. Rev. 9, 477–482.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Serreze, D. V., Bridgett, M., Chapman, H. D., Chen, E., Richard, S. D., and Leiter, E. H. (1998). Subcongenic analysis of the Idd13 locus in NOD/Lt mice: evidence for several susceptibility genes including a possible diabetogenic role for beta 2-microglobulin. J. Immunol. 160, 1472–1478.

Pubmed Abstract | Pubmed Full Text

Silveira, P. A., and Grey, S. T. (2006). B cells in the spotlight: innocent bystanders or major players in the pathogenesis of type 1 diabetes. Trends Endocrinol. Metab. 17, 128–135.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Stechova, K., Kolar, M., Blatny, R., Halbhuber, Z., Vcelakova, J., Hubackova, M., et al. (2011). Healthy first degree relatives of patients with type 1 diabetes exhibit significant differences in basal gene expression pattern of immunocompetent cells compared to controls: expression pattern as predeterminant of autoimmune diabetes. Scand. J. Immunol. 75, 210–219.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Suzuki, J., Ricordi, C., and Chen, Z. (2010). Immune tolerance induction by integrating innate and adaptive immune regulators. Cell Transplant. 19, 253–268.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Teft, W. A., Kirchhof, M. G., and Madrenas, J. (2006). A molecular perspective of CTLA-4 function. Annu. Rev. Immunol. 24, 65–97.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ueda, H., Howson, J. M., Esposito, L., Heward, J., Snook, H., Chamberlain, G., et al. (2003). Association of the T-cell regulatory gene CTLA4 with susceptibility to autoimmune disease. Nature 423, 506–511.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wing, K., Onishi, Y., Prieto-Martin, P., Yamaguchi, T., Miyara, M., Fehervari, Z., et al. (2008). CTLA-4 control over Foxp3+ regulatory T cell function. Science 322, 271–275.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wong, F. S., Wen, L., Tang, M., Ramanathan, M., Visintin, I., Daugherty, J., et al. (2004). Investigation of the role of B-cells in type 1 diabetes in the NOD mouse. Diabetes 53, 2581–2587.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Yoon, J. W., and Jun, H. S. (2001). Cellular and molecular pathogenic mechanisms of insulin-dependent diabetes mellitus. Ann. N. Y. Acad. Sci. 928, 200–211.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zou, H., and Trevor, T. (2005). Regularization and variable selection via the elastic net. J. R. Stat. Soc. Ser. B 67, 301–320.

Keywords: autoimmunity, computation, diabetes, human, gene expression

Citation: Han D, Cai X, Wen J, Kenyon NS and Chen Z (2012) From biomarkers to a clue of biology: a computation-aided perspective of immune gene expression profiles in human type 1 diabetes. Front. Immun. 3:320. doi: 10.3389/fimmu.2012.00320

Received: 18 July 2012; Accepted: 02 October 2012;
Published online: 25 October 2012.

Edited by:

Lucienne Chatenoud, Université Paris Descartes, France

Reviewed by:

Ye Zheng, Salk Institute for Biological Studies, USA
Daniel Gray, The Walter and Eliza Hall Institute for Medical Research, Australia

Copyright: © 2012 Han, Cai, Wen, Kenyon and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.

*Correspondence: Zhibin Chen, Department of Microbiology and Immunology, University of Miami Miller School of Medicine, RMSB 3035, 1600 Northwest 10th Avenue, Miami, FL 33136, USA. e-mail: zchen@med.miami.edu