Plasma proteomic profiles from disease-discordant monozygotic twins suggest that molecular pathways are shared in multiple systemic autoimmune diseases*
© O'Hanlon et al.; licensee BioMed Central Ltd. 2011
Received: 5 January 2011
Accepted: 1 November 2011
Published: 1 November 2011
Although systemic autoimmune diseases (SAID) share many clinical and laboratory features, whether they also share some common features of pathogenesis remains unclear. We assessed plasma proteomic profiles among different SAID for evidence of common molecular pathways that could provide insights into pathogenic mechanisms shared by these diseases.
Differential quantitative proteomic analyses (one-dimensional reverse-phase liquid chromatography-mass spectrometry) were performed to assess patterns of plasma protein expression. Monozygotic twins (four pairs discordant for systemic lupus erythematosus, four pairs discordant for juvenile idiopathic arthritis and two pairs discordant for juvenile dermatomyositis) were studied to minimize polymorphic gene effects. Comparisons were also made to 10 unrelated, matched controls.
Multiple plasma proteins, including acute phase reactants, structural proteins, immune response proteins, coagulation and transcriptional factors, were differentially expressed similarly among the different SAID studied. Multivariate Random Forest modeling identified seven proteins whose combined altered expression levels effectively segregated affected vs. unaffected twins. Among these seven proteins, four were also identified in univariate analyses of proteomic data (syntaxin 17, α-glucosidase, paraoxonase 1, and the sixth component of complement). Molecular pathway modeling indicated that these factors may be integrated through interactions with a candidate plasma biomarker, PON1 and the pro-inflammatory cytokine IL-6.
Together, these data suggest that different SAID may share common alterations of plasma protein expression and molecular pathways. An understanding of the mechanisms leading to the altered plasma proteomes common among these SAID may provide useful insights into their pathogeneses.
Keywordsautoimmunity inflammation human proteome
Systemic autoimmune diseases (SAID) (for example, systemic lupus erythematosus (SLE), rheumatoid arthritis, scleroderma, and dermatomyositis) result in significant morbidity and mortality and a large socioeconomic burden in the United States, where they are estimated to afflict more than five percent of the population . Evidence for immune-mediated pathologies associated with these heterogeneous syndromes comes from the frequent finding of autoantibodies, chronic inflammation of multiple organ systems, and clinical improvement with immunosuppressive therapy. Familial disease associations but limited disease concordance between monozygotic (MZ) twins, ethnogeographic and seasonal clustering of disease onset, and the identification of shared genetic risk factors support the hypothesis that chronic immune activation in SAID is triggered by specific environmental exposures in genetically susceptible individuals .
Proteomic analyses of human biological fluids (for example, plasma, urine, saliva, cerebral spinal and synovial fluids) have enabled the differential quantitation of large numbers of protein molecules between healthy and diseased subjects. Studies utilizing bio-fluid proteomics have identified multiple, pathologic markers and molecular pathways associated with different disease phenotypes, severities, and therapeutic responses [3, 4]. Yet, despite these in-roads, considerable variability in the published SAID literature exists and likely results from multiple factors including different proteomic methodologies (for example, 2-D electrophoresis, mass spectrometry, antibody array), choice of bio-fluids or tissues analyzed, and the inherent heterogeneity of SAID phenotypes, patient histories, and human genetic variations. Nevertheless, some consensus has emerged in multiple, independent lines of proteomic research in the rheumatic diseases . These common findings in multiple rheumatic diseases to date include Type I interferon inducible proteins, autoantibodies, numerous inflammatory cytokines/chemokines, and markers of molecular pathways associated with chronic immune activation (for example, NF-kB, TNFα, and complement fixation), oxidative stress, coagulation, protein degradation and lipid metabolism [3–8].
Proteomic analysis of blood plasma has several useful research advantages despite its technical complexity. Blood plasma has an exceedingly complex proteome consisting of approximately 1,000 distinct polypeptides, whose concentrations vary over several orders of magnitude . The vast majority of total plasma protein, however, is comprised of a smaller number of more abundant proteins (for example, albumin, immunoglobulins and haptoglobin), which necessitate their pre-depletion to enhance the detection of other minor protein constituents present at much lower concentrations. Despite these methodologic challenges, the plasma proteome is one of the most extensively characterized bio-fluids in humans [10, 11]. Moreover, plasma samples are more easily obtained using a minimally invasive procedure, and are an ideal source of circulating disease-associated markers as well as those derived from dead or leaking cells from pathologic tissues throughout the body [3, 4, 9].
In human proteomic studies, statistically significant differences in protein levels among experimental and control subjects are often subtle (that is, 1.5- to 4-fold variations) and influenced potenlially by the degree of genetic variation that exists among human study subjects [3, 4, 11]. To help mitigate the potentially confounding effects of human genetic polymorphisms in our study population, we utilized liquid chromatography electrospray ionization mass spectrometry (LC-ESI-MS) to measure quantitative differences in the plasma proteome of SAID-discordant MZ twins and unrelated, matched controls. In a hypothesis-generating study, we sought to compare plasma proteomes with the expectation of identifying putative disease-associated markers among study subjects with greater genetic similarity, but possibly different environmental and/or epigenetic influences. To this end, we have identified multiple molecular pathways and possible biomarkers common among different SAID.
Materials and methods
MZ twin pairs discordant for SAID and unrelated, matched, healthy controls (n = 10) were identified for this study. These subjects were selected among those enrolled and providing informed consent between 2001 and 2006 in the NIH investigational review board-approved Twins-Sib study assessing the pathogenesis of SAID. Ethical approval for this proteomics study was obtained from the NIH investigational review board and all human subjects provided informed consent. Study subjects included nine Caucasian twin pairs and one twin pair of Hispanic descent. Patients were defined as those meeting American College of Rheumatology (ACR) criteria for systemic lupus erythematosus (SLE), juvenile idiopathic arthritis (JIA), or juvenile dermatomyositis (JDM) and required the exclusion of inherited, metabolic, infectious diseases or other mimics of SAID; patients were within four years of diagnosis. Twin monozygosity was confirmed by short tandem repeat analysis of genomic DNAs (Proactive Genetics Inc., Martinez, GA, USA). Study subjects comprised three groups: (1) 10 SAID probands (4 SLE, 4 JIA, and 2 JDM); (2) probands' 10 autoimmune disease unaffected MZ twins; and (3) 10 unrelated, matched controls who were also free of SAID. The 10 sets of twin pairs included 6 juveniles (mean age 12.2 years) and 4 adult cases (mean age 25.8 years). The mean ages of juvenile and adult unrelated, healthy controls were 9.8 and 27 years, respectively. Each study group had seven females and three males. Physical global disease activity assessments were determined on a visual analogue scale (0 to 100 mm): SLE (mean 13.2, range 4 to 30); JIA (mean 25.2, range 0 to 40); JDM (mean 4.5, range 2 to 7). To minimize potential confounders, plasma samples were collected in the morning with immunosuppressive therapy held at least 24 hours prior to collection. Unrelated controls were age- (within six years), gender- and ethnically-matched to twins, were free of infections, trauma, vaccines and surgeries for eight weeks and had no first degree family members with SAID.
Proteomic differential expression analysis
Plasma samples were collected and frozen within one hour at -80°C. All samples were shipped on dry ice to PPD Inc. Biomarker Discovery Sciences (Menlo Park, CA, USA). Upon processing, thawed samples were stabilized with a sodium azide and a protease inhibitor cocktail containing 100 ug/mL aprotinin and 5% (w/v) sodium azide, which were added to the plasma at a volume ratio of 1:100. Experimental run order was prepared within a block randomization scheme consisting of matched twin and control samples (10 blocks of 3 subjects each). The order of processing and analyzing samples was separately randomized within each block. Plasma proteins were analyzed by mass spectrometric analysis using a one-dimensional (1-D) separation approach as described below .
For the proteomic analysis, plasma was pre-depleted for the six most abundant proteins (albumin, IgG, IgA, haptoglobin, transferrin and α1 anti-trypsin) by an antibody-based affinity column. The remaining proteins were denatured, reduced, and sulfhydryl groups carboxymethylated prior to trypsin digestion. Also, prior to the trypsin digestion, low molecular weight molecules were excluded during a buffer exchange step with a 5 kDa cut-off filter. Tryptic peptides were then profiled by liquid chromatography-electro-spray ionization-mass spectrometry (LC-ESI-MS) on a high-resolution (R > 5,000) time-of-flight (TOF) instrument (Waters Corp., model LCT) Milford, MA, USA using a capillary chromatography column. The on-line chromatography pump (Agilent, model capillary 1100) Santa Clara, CA, USA was used for reverse-phase (RP) separation with a water/acetonitrile gradient and 0.1% formic acid added to aid in ionization efficiency and chromatographic behavior. A total of 9,549 molecular components were tracked and quantified in the 1-D analysis.
Quality control samples from a large human plasma pool were chemically processed and analyzed along with the clinical samples with an average frequency ratio of one QC sample per eight clinical samples. Process quality control samples were required to maintain coefficients of variation (CV) for many endogenous biomolecules of less than 20%.
Peptides of interest (significantly different in plasma levels) were linked by accurate mass and chromatographic retenion time to separate tandem mass spectrometry (MS/MS) experiments on an ion-trap mass spectrometer (Thermo, model LTQ, West Palm Beach, FL, USA). The resulting MS/MS spectra contained fragmentation patterns with characteristic peptide backbone cleavages. Each MS/MS raw spectrum from an isolated precursor ion was compared with in silico protein digestion and fragmentation data using the NCBI RefSeq sequence database to find a match-quality score and subsequent identification. Mascot software from Matrix Science (Boston, MA, USA) was used for peptide identification. To help separate correct from incorrect database search results, probabilities of correct identification were computed by unsupervised machine learning with an expectation-maximization (EM) algorithm . Here, the probabilities are based both on Mascot scores and on the differences between observed and predicted retention time or retention index. The retention time is predicted using amino acid composition throughout the peptide and specifically at the amino-terminus, as well as peptide length, following the approach previously published , but trained on a data set similar to that acquired here. In this study the probability minimum threshold was set to 0.8.
A label-free differential quantification method was employed that relies on changes in analyte signal intensities directly reflecting their concentrations in one sample relative to another [12, 15]. This quantification technology employs overall spectral intensity normalization by employing signals of molecules that do not significantly change concentration from sample to sample. A simple correction can be applied for any differences in sample concentrations and/or drift over time in LC-MS instrument response. The computation performs normalization by determining the median of the ratios for a large number of molecular ions (spectral components). Analysis of the data included spectral smoothing, baseline subtraction, noise evaluation, peak identification, intensity evaluation, inter-scan evaluation to construct chromatographic peaks and to establish molecular components, and final signal quantification [12, 15]. All processed, primary data are provided as a supplementary submission to this article (Additional files 1, 2, 3).
If the data of the different study groups were approximately normally distributed as determined by the Shapiro-Wilk test, then a two-sided t-test was used; if not, the nonparametric rank test (Wilcoxon or Kruskal-Wallis test) was applied. These comparisons are paired for the two draw times from each individual. Fold-changes in quantitative expression and P-values were determined. All tests of hypotheses in this exploratory study were two-sided and a P-value of < 0.05 was considered significant.
As an alternative means of data interpretation, we determined the relative importance that combined sets of protein components confer upon the accurate classification of the individual study groups (affected twins, unaffected twins, and unrelated, matched controls) using the Random Forests (RF) algorithm developed by Breiman and Cutler [16, 17]. The quantitative expression levels of all factors identified in the 1-D differential expression analysis of disease-discordant twin pairs were classified using RF models (decision trees = 500, node size = 3). Individual decision trees were constructed from combined, unmatched cases and control training data sets utilizing bootstrap sampling with replacement and random variable selection. Classification was performed by a majority vote across the separate trees using test cases and controls omitted from the modeling data set from each of the respective decision trees. In this approach, training and test data are randomly re-utilized in the construction of individual decision trees with an "out-of-bag" (oob) estimate of error rates equalling 20%. All factors in test populations were ranked by their relative importance (RI) in accurately classifying case and control study subjects.
Data were analyzed using the Ingenuity Pathways Analysis (IPA) informatics platform (Ingenuity® Systems, Redwood City, CA, USA). For univariate component analysis, the complete data set, including protein identifiers, corresponding quantitative expression and P-values was utilized. Each protein identifier was mapped to its corresponding gene object and overlaid onto a global molecular network developed from information contained in the IPA Knowledge Base. Networks of genes were then generated algorithmically based on their connectivity as established in the published literature. Fischer's exact test was used to calculate a P-value determining the probability that each biologic function and/or pathway assigned to the data set is due to chance alone.
In a separate analysis, plasma protein components identified as having high relative importance values in the RF multivariate analysis were used to explore putative biologic interactions using IPA Grow, Connect, and Path Explorer applications.
Protein blot analysis
Plasma protein samples (30 μg each) from discordant twins and unrelated, matched controls were resolved by SDS-PAGE (10% precast Criterion gels, Bio-Rad, Hercules, CA, USA) and subsequently dry-blotted to PVDF membranes (iBlot system, Invitrogen, Carlsbad, CA, USA). Protein blots were blocked and incubated with rabbit polyclonal, primary antibodies recognizing human plasma PON1, RBP1, or LRG1 and transferrin (TF) as an internal control for 1 to 24 hours in TBS/0.05% Tween-20 (Abcam, Cambridge, MA, USA). Blots were washed and incubated for 30 minutes with a secondary antibody-HRP conjugate (goat anti-rabbit heavy and light chain IgG (Abcam)). Washed blots were incubated for one minute with chemiluminescent substrate and visualized using a GBOX HR50 molecular imaging system (Syngene, Frederick, MD, USA). Syngene GeneSnap imaging and analysis software was used to quantify and normalize replicate analyses of plasma protein levels.
Plasma proteomic differential expression analysis
Summary of significant differences detected (P < 0.05) among MZ twins discordant for SAID and unrelated, matched controls*
Affected vs. unaffected twins
Keratin type II
Complement C8 precursor
Acute phase reaction
Retinoblastoma-binding protein 1
Complement C4 precursor
Plasma protease C1 inhibitor precursor
Hematopoietic-specific transmembrane-4 protein
Ubiquitin carboxyl-terminal hydrolase 29
Complement C6 precursor
Complement C9 precursor
Leucine-rich α-2-glycoprotein precursor
Serum paraoxonase/arylesterase 1
Apolipoprotein A-II precursor
Minor histocompatibility antigen H13
Affected Twins vs. Unrelated, Matched Controls
Leucine-rich α-2-glycoprotein precursor
Apolipoprotein C-IV precursor
Apolipoprotein E precursor
Nuclear receptor co-activator 6
Plasma retinol-binding protein precursor
Retinoblastoma-binding protein 1
Apolipoprotein C-III precursor
Coagulation factor XII precursor
Peroxisomal carnitine octanoyl transferase
Metabolism fatty acid
Ficolin 3 precursor
Probable ATP-dependent helicase DHX37
Polycystin 1 precursor
Complement C1q A chain precursor
NK cell receptor 3DL1 precursor
Serum paraoxonase/arylesterase 1
Unaffected Twins vs. Unrelated, Matched Controls
Apolipoprotein C-IV precursor
Apolipoprotein C-III precursor
Histone deacetylase 3
Ficolin 3 precursor
Keratin type II
Serum paraoxonase/arylesterase 1
In cases involving comparisons of affected twins to either unaffected twins or unrelated controls, multiple acute phase reactants and markers of immune activation are apparent. The PON1 gene product, paraoxonase 1, was the only marker exhibiting significant differences in expression levels in each of the three two-group comparisons (Figure 1). PON1 levels were reduced in the plasma of affected cases compared to either unaffected twins or unrelated, matched controls. Two additional markers, RBP1 and LRG1, were detected at modestly increased levels (approximately 1.2- to 1.5-fold) in affected twins compared to either unaffected twins or unrelated controls.
Random Forest (RF) multivariate analyses
Protein blot analysis
While certain autoimmune diseases share selected genetic, clinical and laboratory features, it is not clear if shared pathogenic mechanisms might link a number of SAID. One approach to the study of disease pathogenesis is the use of MZ twins as a means of controlling for the inherent genetic variability of study subjects in order to better assess the contribution of genetic, epigenetic and environmental factors . MZ twins, however, are not genetically identical owing to various post-meiotic and age-related epigenetic modifications. Despite these differences, microarray analyses suggest that RNA expression levels of polymorphic genes are more tightly controlled in MZ twins than other first degree family members or unrelated controls [19, 20].
In the present study, we have evaluated biologic pathways altered among multiple SAID by studying levels of plasma proteins using LC-ESI-MS from MZ twins discordant for SAID and unrelated, matched controls. Blood plasma is well-suited to the study of systemic or multi-organ diseases given its capacity to sample proteins from damaged tissues and detect changes in other physiologic pathways associated with complex host responses to disease processes and infectious and/or other environmental agents . The human plasma proteome is one of the most complex and better characterized human bio-fluids wherein the identity and expression levels of its approximately 1,000 distinct protein constituents are currently cataloged .
Previous studies have examined human tissue and bio-fluid proteomes in autoimmune conditions with the goal of identifying disease-specific biomarkers to aid in improved disease diagnosis and understanding of underlying pathogeneses [4, 21–29]. These findings point consistently to coordinated changes in the levels of multiple proteins involved in such canonical pathways as immune activation, signal transduction, cell adhesion, apoptosis, and acute phase responses, in addition to various transcription factors, structural and transport proteins. In fact, composite phenotypic profiles of coordinated changes in multiple protein factors and physiologic pathways rather than solitary biomarkers may prove more reliable in differentiating complex and sometimes overlapping autoimmune syndromes.
We examined multiple SAID in an attempt to uncover shared biomarkers or proteomic profiles, with the understanding that these otherwise heterogeneous disorders often share many clinical features, immunologic abnormalities, genetic risk factors and serum autoantibodies [30, 31]. We hypothesized that certain proteomic profiles may be similar among patients with different SAID and that those profiles will differ from those of unaffected MZ twins. Moreover, we asked whether the proteomic profiles of unaffected twins more closely resembled that of unrelated, matched controls or possibly shared some features with their affected twins as a consequence of their genetic similarity and/or shared environmental exposures.
Collectively, our proteomics data from affected MZ twins was consisten with that from other published studies of human autoimmune diseases. Namely, the apparent coordinated regulation of multiple proteins from several canonical pathways (for example, immune regulation, acute phase response, protein and lipid homeostasis, apoptosis and signal transduction) appears to be associated with these chronic inflammatory conditions. In univariate analyses, we observed multiple proteins whose plasma levels were statistically different in affected twins compared to either unaffected twins or unrelated controls. Some of these proteins (for example, α1-microglobulin, fibrinogen, apolipoproteins A and E, complement C3 and C4B, and retinol binding protein) may exhibit altered plasma levels as a consequence of chronic inflammation as they were also reported as up-regulated in synovial fluid from osteoarthritis (OA) patients . Increased levels of apolipoprotein A were also observed in isolated peripheral blood mononuclear cells from SLE patients and muscle biopsies of patients with inclusion body myositis [21, 24, 28]. Similarly, the leucine-rich α2 glycoprotein marker (LRG1) - a molecule involved in signal transduction, cell adhesion, and granulocyte differentiation - was elevated in plasma from our affected twins and was also found elevated in both the cerebrospinal fluid and serum proteomes from multiple sclerosis patients . More recently, LRG1 was identified as a novel, serum pro-inflammatory biomarker for RA and Crohn's disease . Molecular Pathways analysis of our total proteomics data set comparing SAID discordant MZ twins, helped us identify numerous acute phase reactants, immune complement components, coagulation factors, and retinol binding proteins as potentially important mediators of disease. Together, these data suggest that many of the physiological pathways altered in these patients are not necessarily disease-specific but rather may contribute to inflammatory processes shared by multiple SAID.
Proteomic data sets with large and complex arrays of candidate markers mapping across multiple biologic pathways present limits to the interpretation of univariate data by disregarding potential protein-protein interactions as a basis for accurate disease profiling. Investigators have employed machine learning algorithms for the multivariate analysis of large proteomic data sets derived from cancer prevention trials and human autoimmune disease studies [33, 34]. Liu et al. described the use of a support vector machine algorithm to effectively classify RA patients and controls using serum proteomic component peaks . Among the several decision tree ensemble methods available, we utilized the Random Forests algorithm to create a model which accurately classified affected vs. unaffected twin pairs. Putative interactions among seven proteins (STX17, MGAM, PON1, C6, SYNE1, PLEKHG5 and AZGP1) accounted for the majority of this effect. Several of these proteins were likewise identified in our univariate analyses (STX17, MGAM, PON1 and C6). The STX17 marker was one of three proteins whose altered plasma levels was unique to the comparison of discordant MZ twins, while PON1 was the only marker identified with statistically different levels in each of the three two-group comparisons.
The PON1 gene product, paraoxonase 1, is an arylesterase that serves an important role in several physiological pathways including the detoxification of xenobiotics - most notably organophosphorus metabolites associated with pesticide exposures - as well as reducing oxidative damage when associated with circulating high and low density lipoproteins [35–37]. Interestingly, functional polymorphisms in the PON1 gene influence expression levels and activity of the enzyme and have been associated with several immune-mediated conditions, atherosclerotic risk, and possibly influence responses to anti-TNF-α therapy in RA [38–41].
Several independent lines of evidence implicate reduced plasma PON1 levels as a potential biomarker for a subset of SAID [39, 42, 43]. In our present study, we observed an apparent gradient of decreasing PON1 levels among our three study groups in univariate analyses whereby PON1 levels were lowest in SAID-affected twins and highest in unrelated controls. Also, PON1 was identified as an informative marker in a multivariate RF model, which effectively segregated SAID affected vs. unaffected twins. In molecular pathway modeling, PON1 mapped as a central node in interactions predicted among all the relevant factors in the RF analysis. More recently, certain PON1 polymorphic variants were implicated as risk factors for other chronic inflammatory diseases, including RA and types 1 and 2 diabetes [44, 45]. Plasma protein blot analysis of our twin pairs and matched, unrelated controls demonstrated reduced plasma PON1 levels in 50% of the twin cases independent of disease phenotype. We speculate that shared or similar environmental factors, such as pesticide exposures, might influence the development of different SAID by a common mechanism .
There are several limitations to our plasma proteomics study design. Most importantly, small sample sizes and the resulting decrease in statistical power owing to the difficulties associated with the identification and recruitment of SAID-discordant MZ twins with recent disease onset. Also, the heterogeneity of human study subjects, including variations in environmental exposures, clinical phenotypes, disease activity and duration and immunosuppressive therapies may influence plasma protein composition and present potential confounders. Additionally, given the capacity of mass spectrometric techniques to detect several thousand component peaks from individual plasma samples, higher false discovery rates (FDR) are anticipated in the absence of corrections for multiple statistical comparisons. Despite these limitations, most of the candidate markers and molecular pathways identified in our study are consistent with those identified in other studies of individual human autoimmune disease [21–28, 44, 47].
We have described proteomic profiles common to multiple, different SAID. We analyzed SAID-discordant MZ twins to minimize polymorphic gene effects and found that, in comparison to affected twins, plasma proteomes of unaffected twins more closely resemble those of unrelated, matched controls. These data suggest that in addition to genetic predispositions, disease pathogenesis in MZ twins who develop SAID are likely influenced by post-meiotic genetic events (for example, copy number variations between MZ twins), different epigenetic modifications, epistatic protein interactions, and/or environmental exposures that promote pro-inflammatory biologic pathways. Moreover, the use of complex proteomic profiles - rather than individual biomarkers - may provide a more highly integrated description of immune dysfunction and disease pathogeneses. Our hope is that such studies might lead to earlier and more accurate diagnostics, and more effective, targeted therapeutics.
sixth component of complement
coefficients of variation
false discovery rates
Ingenuity Pathways Analysis
liquid chromatography electrospray ionization mass spectrometry
tandem mass spectrometry
systemic autoimmune disease
systemic lupus erythematosus
The authors thank Drs. Ejaz Shamim, Michael Fessler, Howard Schulman, Christopher Becker, and Aaron Kantor for their critical reviews of the manuscript. The authors also thank the following individuals who referred discordant twin pairs to the Twin-Sib proteomics study: Drs. Linda Clark, Harry Gewanter, Yukiko Kimura, Deborah Rothman, Bracha Shaham, and Fernando Silva.
The funding sponsors had no role in the design, conduct or interpretation of the study.
The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does the mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.
This work was supported by the intramural research program project Z01 ES101074 of the National Institute of Environmental Health Sciences, NIH.
- NIH Autoimmune Diseases Coordinating Committee Report. [http://www.niaid.nih.gov/topics/autoimmune/Documents/adccfinal.pdf]
- Gourley M, Miller FW: Mechanisms of disease: Environmental factors in the pathogenesis of rheumatic disease. Nat Clin Pract Rheumatol. 2007, 3: 172-180. 10.1038/ncprheum0435.View ArticlePubMedGoogle Scholar
- Vanarsa K, Mohan C: Proteomics in rheumatology: the dawn of a new era. F1000 Med Rep. 2010, 2: 87-PubMed CentralPubMedGoogle Scholar
- De Franceschi L, Bosello S, Scambi C, Biasi D, De Santis M, Caramaschi P, Peluso G, La Verde V, Bambara LM, Ferraccioli G: Proteome analysis of biological fluids from autoimmune-rheumatological disorders. Proteomics Clin Appl. 2011, 5: 78-89. 10.1002/prca.201000069.View ArticlePubMedGoogle Scholar
- Bilgic H, Ytterberg SR, Amin S, McNallan KT, Wilson JC, Koeuth T, Ellingson S, Newman B, Bauer JW, Peterson EJ, Baechler EC, Reed AM: Interleukin-6 and type I interferon-regulated genes and chemokines mark disease activity in dermatomyositis. Arthritis Rheum. 2009, 60: 3436-3446. 10.1002/art.24936.View ArticlePubMedGoogle Scholar
- Carlsson A, Wuttge DM, Ingvarsson J, Bengtsson AA, Sturfelt G, Borrebaeck CA, Wingren C: Serum protein profiling of systemic lupus erythematosus and systemic sclerosis using recombinant antibody microarrays. Mol Cell Proteomics. 2011, 10: M110.005033Google Scholar
- Greenberg SA: Dermatomyositis and type 1 interferons. Curr Rheumatol Rep. 2010, 12: 198-203. 10.1007/s11926-010-0101-6.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang L, Dai Y, Qi S, Sun B, Wen J, Zhang L, Tu Z: Comparative proteome analysis of peripheral blood mononuclear cells in systemic lupus erythematosus with iTRAQ quantitative proteomics. Rheumatol Int. 2010,Google Scholar
- Zhang Q, Faca V, Hanash S: Mining the plasma proteome for disease applications across seven logs of protein abundance. J Proteome Res. 2011, 10: 46-50. 10.1021/pr101052y.View ArticlePubMedGoogle Scholar
- Anderson NL, Polanski M, Pieper R, Gatlin T, Tirumalai RS, Conrads TP, Veenstra TD, Adkins JN, Pounds JG, Fagan R, Lobley A: The human plasma proteome: a nonredundant list developed by combination of four separate sources. Mol Cell Proteomics. 2004, 3: 311-326. 10.1074/mcp.M300127-MCP200.View ArticlePubMedGoogle Scholar
- Anderson NL, Anderson NG: The human plasma proteome: history, character, and diagnostic prospects. Mol Cell Proteomics. 2002, 1: 845-867. 10.1074/mcp.R200007-MCP200.View ArticlePubMedGoogle Scholar
- Roy SM, Becker CH: Quantification of proteins and metabolites by mass spectrometry without isotopic labeling. Methods Mol Biol. 2007, 359: 87-105. 10.1007/978-1-59745-255-7_6.View ArticlePubMedGoogle Scholar
- Keller A, Nesvizhskii AI, Kolker E, Aebersold R: Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem. 2002, 74: 5383-5392. 10.1021/ac025747h.View ArticlePubMedGoogle Scholar
- Krokhin OV, Craig R, Spicer V, Ens W, Standing KG, Beavis RC, Wilkins JA: An improved model for prediction of retention times of tryptic peptides in ion pair reversed-phase HPLC: its application to protein peptide mapping by off-line HPLC-MALDI MS. Mol Cell Proteomics. 2004, 3: 908-919. 10.1074/mcp.M400031-MCP200.View ArticlePubMedGoogle Scholar
- Wang W, Zhou H, Lin H, Roy S, Shaler TA, Hill LR, Norton S, Kumar P, Anderle M, Becker CH: Quantification of proteins and metabolites by mass spectrometry without isotopic labeling or spiked standards. Anal Chem. 2003, 75: 4818-4826. 10.1021/ac026468x.View ArticlePubMedGoogle Scholar
- Random Forests - Leo Breiman and Adele Cutler. [http://stat-www.berkeley.edu/users/breiman/RandomForests/]
- Svetnik V, Liaw A, Tong C, Culberson JC, Sheridan RP, Feuston BP: Random Forest: a classification and regression tool for compound classification and QSAR modeling. J Chem Inf Comput Sci. 2003, 43: 1947-1958. 10.1021/ci034160g.View ArticlePubMedGoogle Scholar
- Petronis A: Epigenetics and twins: three variations on the theme. Trends Genet. 2006, 22: 347-350. 10.1016/j.tig.2006.04.010.View ArticlePubMedGoogle Scholar
- Sharma A, Sharma VK, Horn-Saban S, Lancet D, Ramachandran S, Brahmachari SK: Assessing natural variations in gene expression in humans by comparing with monozygotic twins using microarrays. Physiol Genomics. 2005, 21: 117-123. 10.1152/physiolgenomics.00228.2003.View ArticlePubMedGoogle Scholar
- Cheung VG, Bruzel A, Burdick JT, Morley M, Devlin JL, Spielman RS: Monozygotic twins reveal germline contribution to allelic expression differences. Am J Hum Genet. 2008, 82: 1357-1360. 10.1016/j.ajhg.2008.05.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Dai Y, Hu C, Huang Y, Huang H, Liu J, Lv T: A proteomic study of peripheral blood mononuclear cells in systemic lupus erythematosus. Lupus. 2008, 17: 799-804. 10.1177/0961203308089444.View ArticlePubMedGoogle Scholar
- Liu W, Li X, Ding F, Li Y: Using SELDI-TOF MS to identify serum biomarkers of rheumatoid arthritis. Scand J Rheumatol. 2008, 37: 94-102. 10.1080/03009740701747152.View ArticlePubMedGoogle Scholar
- Gobezie R, Kho A, Krastins B, Sarracino DA, Thornhill TS, Chase M, Millett PJ, Lee DM: High abundance synovial fluid proteome: distinct profiles in health and osteoarthritis. Arthritis Res Ther. 2007, 9: R36-10.1186/ar2172.PubMed CentralView ArticlePubMedGoogle Scholar
- Liao H, Wu J, Kuhn E, Chin W, Chang B, Jones MD, O'Neil S, Clauser KR, Karl J, Hasler F, Roubenoff R, Zolg W, Guild BC: Use of mass spectrometry to identify protein biomarkers of disease severity in the synovial fluid and serum of patients with rheumatoid arthritis. Arthritis Rheum. 2004, 50: 3792-3803. 10.1002/art.20720.View ArticlePubMedGoogle Scholar
- Giusti L, Baldini C, Bazzichi L, Bombardieri S, Lucacchini A: Proteomic diagnosis of Sjogren's syndrome. Expert Rev Proteomics. 2007, 4: 757-767. 10.1586/14789418.104.22.1687.View ArticlePubMedGoogle Scholar
- O'Connor KC, Roy SM, Becker CH, Hafler DA, Kantor AB: Comprehensive phenotyping in multiple sclerosis: discovery based proteomics and the current understanding of putative biomarkers. Dis Markers. 2006, 22: 213-225.PubMed CentralView ArticlePubMedGoogle Scholar
- de Seny D, Fillet M, Meuwis MA, Geurts P, Lutteri L, Ribbens C, Bours V, Wehenkel L, Piette J, Malaise M, Merville MP: Discovery of new rheumatoid arthritis biomarkers using the surface-enhanced laser desorption/ionization time-of-flight mass spectrometry ProteinChip approach. Arthritis Rheum. 2005, 52: 3801-3812. 10.1002/art.21607.View ArticlePubMedGoogle Scholar
- Li J, Yin C, Okamoto H, Jaffe H, Oldfield EH, Zhuang Z, Vortmeyer AO, Rushing EJ: Proteomic analysis of inclusion body myositis. J Neuropathol Exp Neurol. 2006, 65: 826-833. 10.1097/01.jnen.0000228204.19915.69.View ArticlePubMedGoogle Scholar
- Miyamae T, Malehorn DE, Lemster B, Mori M, Imagawa T, Yokota S, Bigbee WL, Welsh M, Klarskov K, Nishomoto N, Vallejo AN, Hirsch R: Serum protein profile in systemic-onset juvenile idiopathic arthritis differentiates response versus nonresponse to therapy. Arthritis Res Ther. 2005, 7: R746-755. 10.1186/ar1723.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller FW: Inflammatory myopathies: polymyositis, dermatomyositis, and related conditions. Arthritis and Allied Conditions, A Textbook of Rheumatology. Edited by: Koopman W, Moreland L. 2004, Philadelphia: Lippincott, Williams and Wilkins, 15: 1593-1620.Google Scholar
- Fernando MM, Stevens CR, Walsh EC, De Jager PL, Goyette P, Plenge RM, Vyse TJ, Rioux JD: Defining the role of the MHC in autoimmunity: a review and pooled analysis. PLoS Genet. 2008, 4: e1000024-10.1371/journal.pgen.1000024.PubMed CentralView ArticlePubMedGoogle Scholar
- Serada S, Fujimoto M, Ogata A, Terabe F, Hirano T, Iijima H, Shinzaki S, Nishikawa T, Ohkawara T, Iwahori K, Ohguro N, Kishimoto T, Naka T: iTRAQ-based proteomic identification of leucine-rich alpha-2 glycoprotein as a novel inflammatory biomarker in autoimmune diseases. Ann Rheum Dis. 2010, 69: 770-774. 10.1136/ard.2009.118919.View ArticlePubMedGoogle Scholar
- Izmirlian G: Application of the random forest classification algorithm to a SELDI-TOF proteomics study in the setting of a cancer prevention trial. Ann N Y Acad Sci. 2004, 1020: 154-174. 10.1196/annals.1310.015.View ArticlePubMedGoogle Scholar
- Geurts P, Fillet M, de Seny D, Meuwis MA, Malaise M, Merville MP, Wehenkel L: Proteomic mass spectra classification using decision tree based ensemble methods. Bioinformatics. 2005, 21: 3138-3145. 10.1093/bioinformatics/bti494.View ArticlePubMedGoogle Scholar
- Costa LG, Cole TB, Jarvik GP, Furlong CE: Functional genomic of the paraoxonase (PON1) polymorphisms: effects on pesticide sensitivity, cardiovascular disease, and drug metabolism. Annu Rev Med. 2003, 54: 371-392. 10.1146/annurev.med.54.101601.152421.View ArticlePubMedGoogle Scholar
- Lacasana M, Lopez-Flores I, Rodriguez-Barranco M, Aguilar-Garduno C, Blanco-Munoz J, Perez-Mendez O, Gamboa R, Gonzalez-Alzaga B, Bassol S, Cebrian ME: Interaction between organophosphate pesticide exposure and PON1 activity on thyroid function. Toxicol Appl Pharmacol. 2010, 249: 16-24. 10.1016/j.taap.2010.07.024.View ArticlePubMedGoogle Scholar
- Furlong CE, Suzuki SM, Stevens RC, Marsillach J, Richter RJ, Jarvik GP, Checkoway H, Samii A, Costa LG, Griffith A, Roberts JW, Yearout D, Zabetian CP: Human PON1, a biomarker of risk of disease and exposure. Chem Biol Interact. 2010, 187: 355-361. 10.1016/j.cbi.2010.03.033.PubMed CentralView ArticlePubMedGoogle Scholar
- Franceschi C, Olivieri F, Marchegiani F, Cardelli M, Cavallone L, Capri M, Salvioli S, Valensin S, De Benedictis G, Di Iorio A, Caruso C, Paolisso G, Monti D: Genes involved in immune response/inflammation, IGF1/insulin pathway and response to oxidative stress play a major role in the genetics of human longevity: the lesson of centenarians. Mech Ageing Dev. 2005, 126: 351-361. 10.1016/j.mad.2004.08.028.View ArticlePubMedGoogle Scholar
- Liu C, Batliwalla F, Li W, Lee A, Roubenoff R, Beckman E, Khalili H, Damle A, Kern M, Furie R, Dupuis J, Plenge RM, Coenen MJ, Behrens TW, Carulli JP, Gregersen PK: Genome-wide association scan identifies candidate polymorphisms associated with differential response to anti-TNF treatment in rheumatoid arthritis. Mol Med. 2008, 14: 575-581.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsatsakis AM, Zafiropoulos A, Tzatzarakis MN, Tzanakakis GN, Kafatos A: Relation of PON1 and CYP1A1 genetic polymorphisms to clinical findings in a cross-sectional study of a Greek rural population professionally exposed to pesticides. Toxicol Lett. 2009, 186: 66-72. 10.1016/j.toxlet.2008.10.018.View ArticlePubMedGoogle Scholar
- Ginsberg G, Neafsey P, Hattis D, Guyton KZ, Johns DO, Sonawane B: Genetic polymorphism in paraoxonase 1 (PON1): Population distribution of PON1 activity. J Toxicol Environ Health B Crit Rev. 2009, 12: 473-507. 10.1080/10937400903158409.View ArticlePubMedGoogle Scholar
- Srivastava R, Yu S, Parks BW, Black LL, Kabarowski JH: Autoimmune-mediated reduction of high-density lipoprotein-cholesterol and paraoxonase 1 activity in systemic lupus erythematosus-prone gld mice. Arthritis Rheum. 2011, 63: 201-211. 10.1002/art.27764.PubMed CentralView ArticlePubMedGoogle Scholar
- Jakubowski H: The role of paraoxonase 1 in the detoxification of homocysteine thiolactone. Adv Exp Med Biol. 2010, 660: 113-127. 10.1007/978-1-60761-350-3_11.View ArticlePubMedGoogle Scholar
- Kerekes G, Szekanecz Z, Dér H, Sándor Z, Lakos G, Muszbek L, Csipö I, Sipka S, Seres I, Paragh G, Kappelmayer J, Szomják E, Veres K, Szegedi G, Shoenfeld Y, Soltész P: Endothelial dysfunction and atherosclerosis in rheumatoid arthritis: a multiparametric analysis using imaging techniques and laboratory markers of inflammation and autoimmunity. J Rheumatol. 2008, 35: 398-406.PubMedGoogle Scholar
- Precourt LP, Amre D, Denis MC, Lavoie JC, Delvin E, Seidman E, Levy E: The three-gene paraoxonase family: physiologic roles, actions and regulation. Atherosclerosis. 2011, 214: 20-36. 10.1016/j.atherosclerosis.2010.08.076.View ArticlePubMedGoogle Scholar
- Parks CG, Walitt BT, Pettinger M, Chen JC, De Roos AJ, Hunt J, Sarto G, Howard BV: Insecticide use and risk of rheumatoid arthritis and systemic lupus erythematosus in the women's health initiative observational study. Arthritis Care Res (Hoboken). 2011, 63: 184-194. 10.1002/acr.20335.View ArticleGoogle Scholar
- Wu T, Mohan C: Proteomic toolbox for autoimmunity research. Autoimmun Rev. 2009, 8: 595-598. 10.1016/j.autrev.2009.01.019.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.