A number of publications have reported the use of microarray technology to identify gene expression signatures to infer mechanisms and pathways associated with systemic lupus erythematosus (SLE) in human peripheral blood mononuclear cells. However, meta-analysis approaches with microarray data have not been well-explored in SLE.
In this study, a pathway-based meta-analysis was applied to four independent gene expression oligonucleotide microarray data sets to identify gene expression signatures for SLE, and these data sets were confirmed by a fifth independent data set.
Differentially expressed genes (DEGs) were identified in each data set by comparing expression microarray data from control samples and SLE samples. Using Ingenuity Pathway Analysis software, pathways associated with the DEGs were identified in each of the four data sets. Using the leave one data set out pathway-based meta-analysis approach, a 37-gene metasignature was identified. This SLE metasignature clearly distinguished SLE patients from controls as observed by unsupervised learning methods. The final confirmation of the metasignature was achieved by applying the metasignature to a fifth independent data set.
The novel pathway-based meta-analysis approach proved to be a useful technique for grouping disparate microarray data sets. This technique allowed for validated conclusions to be drawn across four different data sets and confirmed by an independent fifth data set. The metasignature and pathways identified by using this approach may serve as a source for identifying therapeutic targets for SLE and may possibly be used for diagnostic and monitoring purposes. Moreover, the meta-analysis approach provides a simple, intuitive solution for combining disparate microarray data sets to identify a strong metasignature.
Microarrays are powerful tools with capability of measuring the transcript abundance of tens of thousands of genes simultaneously in biological samples. Microarray technology has matured over the past 15 years and is now employed for the study of gene expression signatures associated with disease [1–3]. The clinical utility of microarrays as prognostic tools can be evidenced by the approval of the US Food and Drug Administration (FDA) of a customized microarray, MammaPrint™ (Agendia, Amsterdam, The Netherlands) for predicting the outcomes in breast cancer patients on the basis of a 70-gene expression signature .
Some of the challenges associated with identification of gene expression signatures that differentiate the disease state from healthy controls are the availability of samples, sample size, heterogeneous data sets, and reproducibility. Thus, robustness of the gene expression signature derived from one study needs to be validated by other independent studies, preferably with large sample sizes. In practice, however, several studies with relatively small sample sizes are often used to identify gene expression signatures. In these circumstances, it is beneficial to combine the results of several individual studies using meta-analysis. This process enhances statistical power in identifying more robust and reliable gene signatures.
Several meta-analysis approaches have been proposed specifically for handling heterogeneous data sets. For example, Rhodes et al.  used the approach of utilizing P values of genes across studies to identify gene expression signatures that differentiate cancer tissues from normal tissues and to predict poor or good patient outcomes. Choi et al.  used an effect size estimate approach in a meta-analysis of two cDNA microarray data sets, human hepatocellular carcinoma and prostate cancer, to identify a transcriptional signature for cancer. A Bayesian approach was used by Wang et al. , who performed microarray studies on three different platforms and combined them to study differences in gene expression between B-cell chronic lymphocytic leukemia and normal B cells. Shen et al.  suggested a Bayesian mixture model incorporating the probability of expression measure.
Most of the currently used meta-analysis approaches first identify a set of commonly probed genes across studies and then derive a gene expression signature from these. A shortcoming of this approach is a potential loss of valuable information from individual data sets during the combining process. Thus, we propose a pathway-based meta-analysis approach whereby differentially expressed genes (DEGs) from individual studies are selected using a combination of P value and fold change and the results are combined at the pathway level instead of at the gene level (see Figure 1 and Methods). Additionally, while most other methods perform very little validation or rely solely on the biological plausibility of the obtained results to serve as validation, the approach proposed here includes statistical validation through the leave one data set out permutation method. The results are further confirmed using an independent data set.
A number of authors have reported the use of microarray technology to identify gene expression signatures in systemic lupus erythematosus (SLE) [9–14], mechanisms underpinning SLE [15–17], and pathways related to SLE [18–20]. However, meta-analysis approaches have not been explored sufficiently in the study of SLE microarray data . In the present study, the pathway-based meta-analysis method was applied to four independent gene expression oligonucleotide microarray data sets to identify gene expression signatures for SLE. These data sets were generated from peripheral blood mononuclear cell (PBMC) samples from SLE patients and healthy controls. The resulting signatures were then confirmed by testing on a fifth independent data set.
Data sets from five independent microarray studies comparing PBMC samples from SLE patients with those from healthy individuals were obtained from prominent SLE researchers. These data sets are referred to as data sets 1, 2, 3, 4, and 5. Data sets 1, 2, 4, and 5 are associated with peer-reviewed publications (11-14). Data set 3 is composed of unpublished data. Three of the studies (studies 1, 4, and 5) included only pediatric patients, while the remaining two included only adults. All studies employed the Affymetrix GeneChipmicroarray platform (Affymetrix, Inc., Santa Clara, CA, USA) but the versions of the array type varied (Table 1). In the case of two different array types used for the same study (that is, data sets 1 and 2), we treated them as separate data sets (data sets 1a, 1b, 2a, and 2b) during the meta-analysis. Raw data in the form of Affymetrix CEL files were provided for studies 1, 2, 3, and 5. For data set 4, however, expression values for a short-listed set of genes were provided. While data sets 1 to 4 were used in the meta-analysis workflow, data set 5 served as an independent data set to validate the gene signature derived from the meta-analysis.
Workflow of the pathway-based meta-analysis approach
The overall workflow of the pathway-based meta-analysis is summarized in Figure 1. The meta-analysis used a leave one data set out validation process. Both principal component analysis (PCA) and hierarchical cluster analysis (HCA) were used to visually inspect the leave one data set out cross-validation results. Last, the combined meta-signature obtained from the 4 data sets was validated against an independent fifth data set (data set 5).
For the individual quality control and data analysis steps mentioned below, each data set was considered separately. Additionally, since data sets 1 and 2 used two chip types each, they were considered as four different data sets (1a, 1b, 2a, and 2b) for the initial analysis.
Quality assessment was done for each data set using the Genedata Expressionist (Genedata, San Francisco, CA, USA)  (Figure 1, step 1). Only one sample in data set 2a was discarded from further analysis, because it had too high a value for defective area percentage.
Individual data processing and analysis
Following quality control assessment, each data set was analyzed individually using the ArrayTrack™ tool (US Food and Drug Administration's National Center for Toxicological Research, Jefferson, AR, USA) . ArrayTrack is a comprehensive tool for microarray data storage, analysis, and interpretation that has been developed at the FDA's National Center for Toxicological Research. To maintain consistency during the individual analysis of data sets, similar normalization methods, statistical tests, and parameters were used with all data sets. First, all data sets except data set 4 were normalized using Robust Multi-array Analysis. Then Welch's t test was performed on each data set individually. The P value and fold change filters (0.01 and 1.5, respectively) were used to identify a unique list of DEGs from each data set (Figure 1, step 2). This list represented genes that were either notably upregulated or downregulated in the PBMCs of SLE patients when compared to the PBMCs of healthy controls. Each DEG list was then used to identify biological pathways significantly represented in SLE samples compared to the healthy controls (P < 0.01) in each data set (Figure 1, step 3). This pathway analysis was done using Ingenuity Pathway Analysis (IPA) software (Ingenuity Systems Inc., Redwood City, CA, USA).
Pathways common to all of the data sets were identified from the individual lists of pathways enriched in SLE patients compared to healthy controls (one for each data set) (Figure 1, step 4). The resulting list of pathways was indicative of processes significantly affected in all of the SLE data sets and comprised a pathway signature representative of all data sets and of the disease. From this common pathway signature, gene markers that met all of the following criteria were selected: (1) exhibited a fold change greater than 2 in at least one of the data sets (stringency increased from 1.5-fold to 2-fold to obtain a robust signature), (2) present in the DEG list in at least one of the data sets, and (3) involved in at least one of the commonly enriched pathways (Figure 1, step 5). These DEGs composed the collective signature (Figure 1, step 6).
Validation with the leave one data set out permutation method
To validate this technique, a leave one data set out permutation approach was employed (Figure 1, step 7). The meta-analysis technique described above was reiterated four times, each time leaving out one of the four data sets (data sets 1 to 4) and performing the analysis using the remaining three data sets. This gave rise to four different scenarios (Table 2). The gene signature obtained using the three data sets (for example, data sets 1 to 3) was then applied to the data set left out (for example, data set 4). Unsupervised visualization techniques such as PCA and HCA were performed to examine how well the signature could differentiate SLE patients from healthy controls (Figures 2 and 3).
Gene markers were generated for each of the four scenarios as described in the Meta-analysis section. Gene markers present in at least three of the four scenarios were grouped to comprise a 37-gene metasignature.
Confirmation using a fifth independent data set
Confirmation of the final 37-gene metasignature was done using an independent fifth data set, data set 5 (Figure 1, step 8). Again, PCA and HCA were carried out to evaluate the ability of the metasignature to differentiate SLE patients from healthy controls in this independent data set (Figure 4).
Results and discussion
Individual data sets of SLE and healthy control data sets derived from Affymetrix microarrays were analyzed using ArrayTrack following quality control (Figure 1, step 1) and normalization procedures. DEGs for individual data sets were identified using a P value cutoff of 0.01 and a fold change cutoff of 1.5 (Figure 1, step 2).
Biological pathways identified in SLE patients through the leave one data set out permutation method
After applying the leave one data set out approach for each of the four scenarios (Table 2), commonly enriched biological pathways were identified using IPA software (Table 3). Three biological pathways were consistently enriched in SLE patients in all four scenarios: interferon (IFN) signaling, interleukin (IL)-10 signaling, and glucocorticoid receptor signaling. An additional pathway, LXR/RXR signaling, was identified only in scenario IV.
Previous studies have provided evidence of increased autoimmunity in patients undergoing IFN treatment . More specifically, there is evidence of women developing SLE during IFN-α treatment . Several studies have shown upregulation of the IFN signaling pathway in SLE patients [9, 12, 17, 26–30]. Therefore, it is understandable that IFN signaling appears to be affected across all data sets.
IL-10 signaling appears to be dysregulated and may be indicative of the inflammatory processes involved in SLE. IL-10 binds to IL-10 receptor 1 on immune cells and activates the JAK-STAT signaling pathway, which is the key IFN signaling mechanism . In support of the hypothesis that IL-10 is involved in SLE, IL-10 has been identified as one of risk loci for SLE in a large genome-wide association study .
Glucocorticoid receptors are also believed to influence cytokine signaling and may be indirectly involved in the pathways underpinning SLE . In fact, glucocorticoids are routinely used in the treatment of SLE patients.
Genes differentially expressed in SLE
Each of the four scenarios produced a gene signature: scenario I produced a signature comprising 51 genes, scenario II produced a signature with 31 genes, scenario III produced a signature with 34 genes, and scenario IV produced a signature with 28 genes. These DEGs represent the three main SLE disease pathways (IFN signaling, IL-10 signaling pathway, and glucocorticoid signaling pathways) as discussed in the section above.
A separate analysis of the pediatric and adult data sets used to identify DEGs and pathways in the two populations was conducted. Similar gene expression patterns were observed in adult and pediatric populations, although the extent of upregulation of some of the genes was higher in the pediatric data sets (unpublished results).
Validation of the meta-analysis approach
For each scenario, the signature obtained using the three data sets was applied to the fourth data set to observe how effectively the expression of the signature genes could distinguish between the SLE and healthy populations.
The PCAs and HCAs obtained for each scenario are presented in Figures 2 and 3, respectively. The PCA and HCA produced similar and consistent results. Grouping of samples based on the expression of signature genes alone produced a clear distinction between SLE patients and healthy controls. The results suggest that the DEG signatures derived by using the leave one data set out permutation approach in the four scenarios (Table 2) can potentially identify a robust gene expression signature for SLE.
Gene expression signatures for SLE and Systemic Lupus Erythematosus Disease Activity Index scores
The Systemic Lupus Erythematosus Disease Activity Index (SLEDAI) is a validated scoring system that can be used to describe the range of disease activity and comprises a weighted score calculated by the presence or absence of 24 symptoms. The association of SLEDAI scores to expression profiles of SLE patients was evaluated. While the majority of the samples were grouped into their respective classes (SLE or control; see Figure 2), 12 SLE patients exhibited expression profiles similar to the control samples. On closer examination of these samples, it was found that the scores for nine of the patients indicated that they were either in remission (SLEDAI score 0) or had mild activity of the disease (SLEDAI score 2 or 3). These findings lend further credence to the ability of the pathway-based meta-analysis approach used here in distinguishing SLE patients from healthy controls. Correlation between SLEDAI scores and gene expression signatures has also been reported in the literature [9, 26, 34].
Metasignature for SLE
A 37-gene signature was generated by the meta-analysis workflow (Table 4). Many IFN-induced genes involved in the IFN signaling pathway (Figure 5), such as IFIT1, IFIT3, IFITM1, IFIT35, MX1, and OAS1, were present in the signature. Overexpression of IFN-regulated genes in PBMCs of SLE patients has been reported in several publications [9, 15, 26, 35–37]. In addition to genes involved in the IFN signaling pathway, genes in cytokine signal transduction (SOCS1 and SOCS3) were also among the DEGs in SLE patients. Differential expression of many other biomarkers associated with inflammatory and/or immune responses and with cellular proliferation was also observed, as shown in Table 4.
Confirmation of the metasignature using an independent data set
This signature was applied to an independent fifth data set (data set 5) to evaluate its ability to distinguish the SLE samples from the control samples. Figure 4 shows that the signature demonstrated clear differentiation between SLE patients and healthy controls. In the HCA analysis, nine of ten healthy samples clustered together and were clearly separated from the cluster of SLE samples (Figure 4A). The PCA analysis also showed that the majority of the SLE samples and healthy samples were grouped separately (Figure 4B). The one SLE sample that was clustered with the healthy samples had a SLEDAI score of 2, confirming our earlier observations with different data sets (Figure 2).
The novel pathway-based meta-analysis approach proved to be a useful technique for grouping disparate microarray data sets. This technique allowed for validated conclusions to be drawn across four different data sets and confirmed by testing on an independent fifth data set. Since the metasignature was obtained from pathways that were enriched in SLE samples across all of the data sets, it is highly representative of biological pathways related to SLE. The metasignature produced here may serve as a source for identifying potential therapeutic targets for SLE. Also, with further refinement, it might be made clinically more intuitive, which may also prove to be useful for the diagnosis and monitoring of SLE. Moreover, the meta-analysis approach outlined here provides a simple and intuitive solution for combining disparate microarray data sets to identify a robust metasignature.
Santos ES, Perez CA, Raez LE: How is gene-expression profiling going to challenge the future management of lung cancer?. Future Oncol. 2009, 5: 827-835. 10.2217/fon.09.60.
Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, Barrette T, Pandey A, Chinnaiyan AM: Large-scale meta-analysis of cancer microarray data identifies common transcriptional profiles of neoplastic transformation and progression. Proc Natl Acad Sci USA. 2004, 101: 9309-9314. 10.1073/pnas.0401994101.
Wang J, Coombes KR, Highsmith WE, Keating MJ, Abruzzo LV: Differences in gene expression between B-cell chronic lymphocytic leukemia and normal B cells: a meta-analysis of three microarray studies. Bioinformatics. 2004, 20: 3166-3178. 10.1093/bioinformatics/bth381.
Nikpour M, Dempsey AA, Urowitz MB, Gladman DD, Barnes DA: Association of a gene expression profile from whole blood with disease activity in systemic lupus erythaematosus. Ann Rheum Dis. 2008, 67: 1069-1075. 10.1136/ard.2007.074765.
Lyons PA, McKinney EF, Rayner TF, Hatton A, Woffendin HB, Koukoulaki M, Freeman TC, Jayne DR, Chaudhry AN, Smith KG: Novel expression signatures identified by transcriptional analysis of separated leucocyte subsets in systemic lupus erythematosus and vasculitis. Ann Rheum Dis. 2010, 69: 1208-1213. 10.1136/ard.2009.108043.
Allantaz F, Chaussabel D, Stichweh D, Bennett L, Allman W, Mejias A, Ardura M, Chung W, Smith E, Wise C, Palucka K, Ramilo O, Punaro M, Banchereau J, Pascual V: Blood leukocyte microarrays to diagnose systemic onset juvenile idiopathic arthritis and follow the response to IL-1 blockade. J Exp Med. 2007, 204: 2131-2144. 10.1084/jem.20070070.
Chaussabel D, Quinn C, Shen J, Patel P, Glaser C, Baldwin N, Stichweh D, Blankenship D, Li L, Munagala I, Bennett L, Allantaz F, Mejias A, Ardura M, Kaizer E, Monnet L, Allman W, Randall H, Johnson D, Lanier A, Punaro M, Wittkowski KM, White P, Fay J, Klintmalm G, Ramilo O, Palucka AK, Banchereau J, Pascual V: A modular analysis framework for blood genomics studies: application to systemic lupus erythematosus. Immunity. 2008, 29: 150-164. 10.1016/j.immuni.2008.05.012.
Sawalha AH, Jeffries M, Webb R, Lu Q, Gorelik G, Ray D, Osban J, Knowlton N, Johnson K, Richardson B: Defective T-cell ERK signaling induces interferon-regulated gene expression and overexpression of methylation-sensitive genes similar to lupus patients. Genes Immun. 2008, 9: 368-378. 10.1038/gene.2008.29.
Thibault DL, Graham KL, Lee LY, Balboni I, Hertzog PJ, Utz PJ: Type I interferon receptor controls B-cell expression of nucleic acid-sensing Toll-like receptors and autoantibody production in a murine model of lupus. Arthritis Res Ther. 2009, 11: R112-10.1186/ar2771.
Liu K, Li QZ, Delgado-Vega AM, Abelson AK, Sánchez E, Kelly JA, Li L, Liu Y, Zhou J, Yan M, Ye Q, Liu S, Xie C, Zhou XJ, Chung SA, Pons-Estel B, Witte T, de Ramón E, Bae SC, Barizzone N, Sebastiani GD, Merrill JT, Gregersen PK, Gilkeson GG, Kimberly RP, Vyse TJ, Kim I, D'Alfonso S, Martin J, Harley JB, Criswell LA, Profile Study Group; Italian Collaborative Group; German Collaborative Group; Spanish Collaborative Group; Argentinian Collaborative Group; SLEGEN Consortium, Wakeland EK, Alarcón-Riquelme ME, Mohan C: Kallikrein genes are associated with lupus and glomerular basement membrane-specific antibody-induced nephritis in mice and humans. J Clin Invest. 2009, 119: 911-923. 10.1172/JCI36728.
Tong W, Cao X, Harris S, Sun H, Fang H, Fuscoe J, Harris A, Hong H, Xie Q, Perkins R, Shi L, Casciano D: ArrayTrack: supporting toxicogenomic research at the U.S. Food and Drug Administration National Center for Toxicological Research. Environ Health Perspect. 2003, 111: 1819-1826. 10.1289/ehp.6497.
Kirou KA, Lee C, George S, Louca K, Peterson MG, Crow MK: Activation of the interferon-α pathway identifies a subgroup of systemic lupus erythematosus patients with distinct serologic features and active disease. Arthritis Rheum. 2005, 52: 1491-1503. 10.1002/art.21031.
Wang S, Yang N, Zhang L, Huang B, Tan H, Liang Y, Li Y, Yu X: Jak/STAT signaling is involved in the inflammatory infiltration of the kidneys in MRL/lpr mice. Lupus. 2010, 19: 1171-1180. 10.1177/0961203310367660.
Gateva V, Sandling JK, Hom G, Taylor KE, Chung SA, Sun X, Ortmann W, Kosoy R, Ferreira RC, Nordmark G, Gunnarsson I, Svenungsson E, Padyukov L, Sturfelt G, Jönsen A, Bengtsson AA, Rantapää-Dahlqvist S, Baechler EC, Brown EE, Alarcón GS, Edberg JC, Ramsey-Goldman R, McGwin G, Reveille JD, Vilá LM, Kimberly RP, Manzi S, Petri MA, Lee A, Gregersen PK, et al: A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus. Nat Genet. 2009, 41: 1228-1233. 10.1038/ng.468.
Feng X, Wu H, Grossman JM, Hanvivadhanakul P, FitzGerald JD, Park GS, Dong X, Chen W, Kim MH, Weng HH, Furst DE, Gorn A, McMahon M, Taylor M, Brahn E, Hahn BH, Tsao BP: Association of increased interferon-inducible gene expression with disease activity and lupus nephritis in patients with systemic lupus erythematosus. Arthritis Rheum. 2006, 54: 2951-2962. 10.1002/art.22044.
Pascual V, Allantaz F, Patel P, Palucka AK, Chaussabel D, Banchereau J: How the study of children with rheumatic diseases identified interferon-alpha and interleukin-1 as novel therapeutic targets. Immunol Rev. 2008, 223: 39-59. 10.1111/j.1600-065X.2008.00643.x.
The authors thank The Office of Women's Health, US Food and Drug Administration, for funding the research project (OWH study 07-09-0001-SA). The authors also thank all the prominent research groups in the SLE field who contributed data for the meta-analysis: Dr Thomas Aune, Dr Emily Baechler Gillespie, Dr Mary Crow, Dr Bernard Lauwerys, GeneLogic, Dr Dafna Gladman, Dr Virginia Pascual, and Dr Violeta Rus. Of all data sets, only data generated with oligonucleotide microarrays were selected for the meta-analysis described in this article. The views presented in this article do not necessarily reflect those of the US Food and Drug Administration.
Authors and Affiliations
National Center for Toxicological Research, US Food and Drug Administration, 3900 NCTR Road, Jefferson, AR, 72079, USA
Dhivya Arasappan & Weida Tong
Genomics Group, Office of Clinical Pharmacology, Office of Translational Science, Center for Drug Evaluation and Research, US Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD, 20993, USA
Padmaja Mummaneni & Shashi Amur
Z-Tech Corporation, an ICF International Company at the National Center for Toxicological Research, US Food and Drug Administration, 3900 NCTR Road, HFT 230, Jefferson, AR, 72079, USA
University of Texas, Austin, Texas Institute for Drug and Diagnostic Development, 2500, Speedway, TX, 78712, USA
The authors declare that they have no competing interests.
DA created the first draft of the manuscript and performed the meta-analysis. PM and SA performed the analysis to confirm the metasignature in an independent data set. SA, PM, and HF helped significantly to draft the manuscript. SA and WT helped coordinate the project and finalize the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License (
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Arasappan, D., Tong, W., Mummaneni, P. et al. Meta-analysis of microarray data using a pathway-based approach identifies a 37-gene expression signature for systemic lupus erythematosus in human peripheral blood mononuclear cells.
BMC Med9, 65 (2011). https://doi.org/10.1186/1741-7015-9-65