Decreased expression of B cell related genes in leukocytes of women with Parkinson's disease
© Kedmi et al; licensee BioMed Central Ltd. 2011
Received: 24 January 2011
Accepted: 23 September 2011
Published: 23 September 2011
Parkinson's disease (PD) is a complex disorder caused by genetic, environmental and age-related factors, and it is more prevalent in men. We aimed to identify differentially expressed genes in peripheral blood leukocytes (PBLs) that might be involved in PD pathogenesis. Transcriptomes of 30 female PD-patients and 29 age- and sex-matched controls were profiled using GeneChip Human Exon 1.0 ST Arrays. Samples were from unrelated Ashkenazi individuals, non-carriers of LRRK2 G2019S or GBA founder mutations.
Differential expression was detected in 115 genes (206 exons), with over-representation of immune response annotations. Thirty genes were related to B cell functions, including the uniquely B cell-expressed IGHM and IGHD, the B cell surface molecules CD19, CD22 and CD79A, and the B cell gene regulator, PAX5. Quantitative-RT-PCR confirmation of these 6 genes in 79 individuals demonstrated decreased expression, mainly in women patients, independent of PD-pharmacotherapy status.
Our results suggest that the down regulation of genes related to B cell activity reflect the involvement of these cells in PD in Ashkenazi individuals and represents a molecular aspect of gender-specificity in PD.
Parkinson's disease (PD) is a complex disorder involving multiple affected genes and multiple environmental risk factors . It is the second most common neurodegenerative disorder, affecting about 1.8% of the population over the age of 65 years . Since the identification of mutations in the α-synuclein gene , 16 chromosomal loci, and mutations in 9 genes have been associated with familial and sporadic PD (reviewed by ). Abnormalities in multiple cellular pathways including the ubiquitin-proteasome, mitochondrial and apoptotic pathways and impaired protection from oxidative stress were suggested to be involved in the accumulation of synuclein and the selective loss of dopaminergic and other neurons (reviewed by ). Transcriptional profiling by microarray methodology of the substantia nigra (SN) from PD-patients [6–8], as well as from peripheral blood leukocytes (PBL) from PD-patients [9, 10] further demonstrated expression changes in genes belonging to these pathways. The identification of similar expression changes in SN and PBL encouraged our use of PBL as a surrogate tissue to search for novel genes and cellular pathways that might be involved in PD pathogenesis.
Genetic analysis of isolated populations derived from relatively few founders offers a powerful route for the identification of genetic risk factors . The Ashkenazi Jewish population has preserved its homogeneous genetic makeup, and has been valuable for the identification of genes associated with increased risk for many common complex diseases, relevant to the world population-at-large [12, 13]. In our previous studies, mutations in the LRRK2 and GBA genes were detected in a surprisingly high proportion (more than a third) of the Ashkenazi PD-patients tested [14, 15], allowing the sub-classification of these patients based on their LRRK2 or GBA carrier status.
Gender differences in PD are well characterized, with PD prevalence being about 1.5 times lower in women than in men in Western populations [16, 17]. Phenotypic changes between men and women have also been described, for example, age at PD onset is ~3 years later in women [16, 17], and men with PD have higher risk of developing cognitive impairment .
To increase the possibility of identifying PD-related expression changes in patients' PBL, we studied the relatively genetically homogenous population of Ashkenazi patients and controls. To further increase the homogeneity of the studied population in the initial step of the microarray expression profiling, we included only female patients and controls that do not carry the LRRK2 or GBA Ashkenazi founder mutations. The microarray methodology was designed to capture the expression intensity of each exon individually, allowing the identification of differentially expressed genes in PD-patients compared to healthy controls.
Quality control and batch effects removal
Principal component analysis (PCA) revealed two methodological factors that affected the expression levels, the Affymetrix whole-transcript target labeling kit batches (2 batches, Additional file 1, Figure S1, A) and the two researchers who extracted the RNAs (Additional file 1, Figure S1, B). Successful elimination of the batch effects and data homogeneity were demonstrated in PCA done following the batch-removed ANOVA function (Partek® Genomics Suite Version 6.4, Additional file 1, Figures S1, C and D).
Expression changes detected in PBL from Parkinson's disease patients
Genes associated with immune response, specifically with B cell-related functions, are differentially expressed in PD-patients
Enriched KEGG pathways among the 115 changed transcripts
Genes in pathway
B cell receptor signaling pathway
BLNK, CD19, CD22, CD72, CD79A, CD79B, CR2, RASGRP3
O = 8;E = 0.419;R = 19.111;P = 2.21e-8
Hematopoietic cell lineage
CD19, CD22, CD8B, CR2, FCER2, IL1R2
O = 6;E = 0.546;R = 10.995;P = 2.60e-5
Type II diabetes mellitus
CACNA1C, CACNA1D, IRS2
O = 3;E = 0.284;R = 10.560;P = 3.45e-3
Calcium signaling pathway
ADCY2, CACNA1C, CACNA1D, GNAL, P2RX5
O = 5;E = 1.114;R = 4.489;P = 5.78e-3
MAPK signaling pathway
CACNA1C, CACNA1D, CACNB4, IL1R2, RASGRF1, RASGRP3
O = 6;E = 1.779;R = 3.373;P = 9.64e-3
ABC transporters - General
O = 2;E = 0.254;R = 7.868;P = 2.92e-2
GnRH signaling pathway
ADCY2, CACNA1C, CACNA1D
O = 3;E = 0.650;R = 4.613;P = 2.93e-2
Cell adhesion molecules (CAMs)
CD22, CD40, CD8B
O = 3;E = 0.807;R = 3.716;P = 4.95e-2
B cell-related genes with significantly changed expression in PBL of PD-patients
GenBank Accession No
B cells involvement
A. Genes involved in B cell receptor signaling pathway
B cell receptor signaling pathway (a); B cell activation (b), (g); B cell differentiation (b)
B cell receptor signaling pathway (a); Expressed in B-cell (f), (g)
B cell receptor signaling pathway (a); Expressed in B-cell (f), (g)
B cell receptor signaling pathway (a); Expressed in B-cell (g)
CD79A antigen (immunoglobulin-associated alpha)
B cell receptor signaling pathway (a),(b); B cell differentiation (b); B cell proliferation (b); B cell receptor complex (c)
CD79B antigen (immunoglobulin-associated beta)
B cell receptor signaling pathway (a),(b); B cell receptor complex (c)
complement component (3d/epstein barr virus) receptor 2
B cell receptor signaling pathway (a); B Lymphocyte Cell Surface Molecules (d)
ras guanyl releasing protein 3 (calcium and dag-regulated)
B cell receptor signaling pathway (a)
B. Genes involved in other B cell-related annotations
B-cell scaffold protein with ankyrin repeats 1
B cell activation (b)
CD40 antigen (tnf receptor superfamily member 5)
Positive regulation of B cell activation (b); Positive regulation of B cell proliferation (b); B cell mediated immunity (b)
early b-cell factor
B Cell Development (e)
FC fragment of ige, low affinity ii, receptor for (cd23a)
Expressed in B-cell (f), (g)
Positive regulation of B cell activation (b); B cell mediated immunity (b)
immunoglobulin heavy constant delta
Expressed in primary B-Cells (f)
immunoglobulin heavy constant mu
Positive regulation of B cell proliferation B cell receptor signaling pathway (b); B cell receptor complex (c); B Cell Development; B cell activation (e); Expressed in primary B-Cells (f)
interleukin 1 receptor, type ii
Expressed in B-cell (f)
musculin (activated b-cell factor-1)
Expressed in B-cell (f)
spi-b transcription factor (spi-1/pu.1 related)
Expressed in B-cell (f)
T-cell leukemia/lymphoma 1a
Expressed in B-cell (f)
C. Additional B cell-related genes
Expressed mainly on mature B cells (h); May be involved in the life/death decision of B-cells (j)
cysteine-rich secretory protein 3
Specifically expressed in pre-B cell  (i)
B-cell novel protein 1
Specifically expressed in B-lymphocytes (k)
FC receptor-like 1
May serve as an activating coreceptor on B cells  (i)
FC receptor-like 2
Expression is limited to the mature B-cell lines (h); May have an regulatory role in normal and neoplastic B cell development (j)
FC receptor-like and mucin-like 1
Specifically expressed in B-cells(k)
nuclear factor (erythroid-derived 2)-like 3
Highly expressed in B-cell (h)
Might play a role in B cell lysosomal function  (i)
5'-nucleotidase, ecto (CD73)
Mediates B cell adhesion  (i)
paired box gene 5 (b-cell lineage specific activator)
Key regulator of the B-cell-restricted expression of the CD23a isoforms  (i)
tank-binding kinase 1
TBK1-mediated signaling in haematopoietic cells was critical for the induction of antigen-specific B and CD4(+) T cells  (i)
We further generated a combined list of genes included in all B cell-related annotations, using DAVID functional analysis tools (342 genes listed in Additional file 4, Table S3 from GO, KEGG, Biocarta, BBID, UniProt Tissue and UniProt Keywords databases; http://david.abcc.ncifcrf.gov/ [23, 24]). A significant enrichment of genes that carry B cell-related annotations was identified among the 115 changed genes (19/115, Tables 2A and 2B) compared to all genes in the core-level Exon 1.0 array (342/21,980, χ2 test, χ2 = 159.43, P < 0.0001), and to all genes expressed in our PBL samples (272/17,441, P < 0.0001, χ2 test, χ2 = 156.9, P < 0.0001). Furthermore, 11 additional genes that are associated with B cells were also detected via a manual search of the following databases: GeneCard, UniProt Tissue Specificity, UniProt Function and PubMed (Table 2C). In total, a high proportion (26%) of the significantly changed genes (30/115, Table 2) was related to B cell functions. Although some of these genes may carry additional roles, their common significantly enriched functions related to B cell biology. Notably, among these genes were IGHM and IGHD that are uniquely expressed in B cells, and the CD19 gene, which is frequently used as a marker for B cell identification.
Positive correlations between the differentially expressed genes
Of the 115 significantly changed genes, 10 genes were chosen for validation, representing both up- and down-regulated genes as well as different GO-annotations, including immune and defense responses (CD22, CD19, CD79A, IGHD and IGHM), regulation of cell proliferation and differentiation (EREG), regulation of transcription (PAX5), voltage-gated ion channel (KCNH8) and cell-cell adhesion (PCDH9). The SNCA gene, whose expression levels were not changed between patients and controls, was used as negative control. Validation of the expression results was done using RNA samples from 9 PD-patients and 6 controls that were also hybridized to the Human Exon 1.0 ST array. Highly significant positive correlations were demonstrated between the expression levels resulting from both methods in all tested genes [ranged between r = 0.864 (P < 0.0005) and r = 0.666 (P = 0.007); Additional file 6, Table S5].
Are the expression levels of B cell-related genes affected by disease state, disease duration, medication or gender?
We expanded the number of RNA samples analyzed by Real-Time PCR (total of 79) to include an additional 64 samples from males and females patients and controls who were not tested on the microarrays. Six B cell-related genes were analyzed: CD22, CD19 and CD79A that encode B cell surface molecules, IGHD and IGHM that encode B cell specific immunoglobulins, and PAX5 that encodes a regulator of B cell genes. In order to define whether disease state (PD), rather than PD pharmacological therapy, affected the expression levels, we included 20 samples from patients who were not treated at the time of enrollment ("naïve"-PD group, with disease duration of 1.65 ± 1.42 years). Two additional groups of medically treated patients were included in this extended analysis: 18 patients with short disease duration (PD-SDD, less than 2 years between onset of motor symptoms to enrollment), and 20 with long disease duration (PD-LDD, more than 5 years between onset of motor symptoms to enrollment). Importantly, all four experimental groups of patients and controls included men and women in ~1:1 ratio (total of 39 men and 40 women).
Studying the B cell related genes expression in other PD patient populations and in other neurological diseases
The expression levels of four of the B cell related genes (IGHD, IGHM, CD19 and CD22) were also examined in 10 non-Ashkenazi female PD patients and 11 non-Ashkenazi female controls. No changes in the expression levels of these genes were detected between PD-patients and controls (Additional file 7, Figure S2). Notably, the expression levels in non-Ashkenazi PD women (patients and controls) were similar to the expression levels in the Ashkenazi control women. These results might suggest that the observed changes in B cell related genes are specific to the Ashkenazi PD women population.
To examine the expression of the B cell related genes in other independent PD cohort and in patients with other neurological diseases we re-analyzed the previously published data by Scherzer at al. (; GSE6613). This dataset compares PD patients to healthy controls and to other neurological disease controls. Based on our criteria for detecting differentially expressed genes (P < 0.01 and FC > 1.5 or < -1.5), no significant changes in the expression of the B cell related genes were observed between any of the three experimental groups, in the whole cohort or when genders were analyzed separately (with the exception of IL1R2, which was up-regulated in males PD patients and was also one of the only up-regulated B cell related genes in our data). These results might suggest that the changes we observed in the expression of B cell related genes are not common to other neurological diseases. However, since we could not confirm the changes also in the PD-patients it is possible that other factors might influenced these results, including demographic properties such as ethnicity and age of the female patients and controls, as well as technical aspects, such as the array used and the normalization method in the dataset of Scherzer et al.. For example, in the original article the authors used the MAS5 algorithm to normalize the array data, while in our re-analysis we utilized the GCRMA algorithm. As a result, out of the 24 probesets defined by Scehrzer at al. as differentially expressed (Supplementalry Table 3 ), only 8 probesets showed P < 0.01 when analyzing the data with GCRMA. Taken together, these results might suggest that the B cell related changes are specific to Ashkenazi Jewish PD female patients, however, the results could also be referred to other causes.
While transcriptional events in brain tissues from deceased PD patients have shed light on disease pathogenesis, recent studies have demonstrated the potential of PBL as a surrogate tissue in PD research. Mammalian lymphocytes are capable of synthesizing dopamine (DA) and other catecholamines, as well as DA receptors and transporters (DAT) [26, 27]. Furthermore, when PBL from PD-patients were compared to controls, differences in DA signaling, including reduced DA content, impaired DAT immunoreactivity  and higher levels of dopamine D1-like and D2-like receptors, were observed . The hypothesis that PBL can echo some of the changes occurring in the substantia nigra of PD-patients is also supported by our results. Changes in expression levels, in the same direction (up- or down-regulation), of the ADCY2, CCDC92, CELSR1, HECTD2 and KIF1B genes, were demonstrated both here, in patients' PBL, and previously, in the substantia nigra of PD-patients [6, 19]. In addition, genes involved in PD-related pathways, such as the ubiquitin-proteasome and apoptosis, and mitochondrial function were differentially expressed in patients' PBL both in our study and in recent reports [9, 10]. These results suggest that PBL from PD-patients serve as an important, easily accessible tool, that might help in the study of mechanisms underlying Parkinson's disease pathogenesis. Our data also demonstrated that the selection of a relatively homogeneous group of RNA samples from Ashkenazi women that do not carry either GBA or LRRK2 founder mutations [14, 15] increased the ability to detect novel expression changes in PD patients' PBL.
The involvement of the immune system, and particularly T cells, in PD has been recognized (reviewed by ). Our expression analysis of PBL from PD-patients demonstrated that the most prominent group of differentially expressed genes were those involved in immune system processes. Notably, a significantly decreased expression of innate and humeral immune response genes was detected, mainly of genes related to B cell functions. The down-regulation of dozens of genes encoding the CD surface molecules, the B cell specific immunoglobulins IGHM and IGHD, and regulators of B cell differentiation and activation, suggest a decrease in the B cell population among women with PD. Such a decrease in the number and percentage of peripheral CD19+ B cells was previously demonstrated in PD-patients . Indeed, decreased mRNA levels of CD19 were confirmed in our study by both microarray and quantitative real-time PCR analyses. Changes in peripheral B and T lymphocytes were also described in other neurodegenerative disorders, such as Alzheimer's disease  and amyotrophic lateral sclerosis . Our findings suggest that B cells might in fact be an additional tissue involved in PD, which has been recently appreciated as multi-systemic, beyond the central nervous system, involving the enteric and autonomic nervous systems as well as the eye (reviewed by  and ).
The detection of expression changes in tissues taken from patients under pharmacotherapy raises the question of whether changes are related to the disease, therapy, or both. Therefore, the quantitative expression confirmation analysis included both treated and untreated (naïve) patients. Our data demonstrated that the reduced expression levels of B cell-related genes did not result from anti-parkinsonian medications, but were related to the presence of PD. This is in agreement with the decreased number and percentage of CD19+ B cells observed in both treated and untreated PD-patients .
While our initial microarray expression study tested only women, we greatly expanded the number of samples for the confirmation analysis, testing six B cell-related genes in male and female PD-patients and controls. Interestingly, the expression levels of all six genes were significantly under-expressed only in women patients compared to women controls, whereas in men, the expression levels were similarly lower in both patients and controls. Such a gender-dependent differential expression of B cell-related genes is reminiscent of the variations detected in the levels of different subsets of lymphocyte populations between healthy men and women [36, 37], specifically, the higher count of CD19+ B cells in healthy women compared to men . Additionally, since our analysis of Ashkenazi PD women samples demonstrated an expression changes in B cell related genes both in the initial stage and in the confirmation study, and given that we could not detect the same changes in expression patterns in non-Ashkenazi PD women (in Jewish women with PD from Eastern and North-African origin, as well as in Scherzer et al.  cohort), it is possible that these changes might be specific to Ashkenazi PD women.
Gender-associated phenotypes are well recognized in PD, particularly the greater prevalence and the younger age of motor symptom onset in men . Specific differences were also detected in peripheral blood, where serum uric acid levels have been inversely correlated with disease duration and daily levodopa dosage in male, but not female PD patients . Recently, gender-dependent gene expression changes were demonstrated in the SN of PD-patients, with only a small overlap of the differentially expressed genes between males and females . The expression changes detected in B cell-related genes in our study are yet another example of molecular variations between men and women with PD. It will be important to explore the relationship between the down regulation of B cell-related genes and the development of PD: Are women with down regulation of B cell-related genes more prone to develop PD, with a higher expression of B cell-related genes possibly carrying a protective effect? or is it the disease process itself that decreases the expression of B cell-related genes in women?
The results presented here demonstrated down-regulated expression of large number of B cell related genes, including genes that are specifically expressed in B cells and genes encoding B cell surface molecules. Moreover, we showed that their expression decreased specifically in women PD-patients. Our results provide compelling evidence for the involvement of B cells in PD even before the initiation of anti-Parkinsonian drugs. Further studies will be needed to elucidate the possible roles of the B cell down-regulation in PD pathogenesis.
Patients and controls
Age at enrollment (AAE) and age at onset (AAO) of the patients' and controls' samples used for the q-RT-PCR assays
61.45 ± 9.10
64.30 ± 8.07
62.30 ± 12.82
60.50 ± 12.79
64.50 ± 9.38
62.90 ± 10.08
65.70 ± 7.82
51.80 ± 10.79
67.10 ± 10.58
57.50 ± 12.83
63.00 ± 8.14
62.11 ± 8.37
64.89 ± 8.65
64.44 ± 8.85
RNA isolation, target labeling and hybridization
PBL were obtained from 2.5 mL of peripheral blood that was lyzed within 2 h from blood taking (buffer EL, QIAGEN, Germantown, MD). The PBL pellets were immediately frozen and kept at -80°C. Total RNA was extracted using QIAamp RNA blood Mini Kit in the presence of DNase (QIAGEN). RNA quality and quantity were measured using NanoDrop spectrophotometer ND-1000 (NanoDrop Technologies, Wilmington, DE). Only RNAs that passed the quality control cutoffs of 260/280 ratio of 1.8 and 260/230 ratio of 1.6 were used. GeneChip whole-transcript sense target labeling of 1 μg of total RNA was done according to the manufacturer's instructions (GeneChip WT cDNA Synthesis and Amplification and GeneChip WT Terminal Labeling Kits, Affymetrix Inc., Santa Clara, CA). The labeled targets were hybridized to the Affymetrix GeneChip Human Exon 1.0 ST Arrays. Hybridized arrays were washed and stained on a GeneChip Fluidics Station 450 and scanned on a GCS3000 Scanner (Affymetrix). Target labeling, hybridization and scanning were done in batches of 8 samples, each containing an equal number of samples from all experimental groups. Microarray experiments were designed to comply with minimum information about a microarray experiment (MIAME) guidelines .
Statistical and bioinformatic analyses
Normalization, filtering and statistical analyses were done using Partek® Genomics Suite Version 6.4 (Partek Inc., St. Louis, MO). Background adjustment, normalization, and probe-level summarization of the microarray data were done using the robust multichip average (RMA) algorithm  and summarized signals were log2 transformed. Gene-level summarization was done by calculating the mean signal for each gene, based on the meta-probeset file provided by Affymetrix. To eliminate the batch effects (scan date and personnel), the ANOVA-batch-removal tool implanted in Partek was applied. The first filtering step was done to include only the core-level data of the Human Exon 1.0 ST array as defined by Affymentrix, which includes 21,980 gene and 232,448 exon probesets, consisting of RefSeq and full-length GenBank mRNAs. To avoid analysis of non-expressed genes, only exons with signal values of ≥ 3.0 were selected for further analysis. To detect transcriptional changes, 3-way ANOVA was conducted on the exon-level data, with the disease status and the two batch-removed effects as the ANOVA factors. To eliminate false-positive results due to probes that site on SNPs , differentially expressed exon-level probesets that contained SNPs in more than 50% of their probes were identified using the SNPinProbe1.0 database  and removed. P-values for the ANOVA model were calculated using log-transformed data. Expression fold changes were calculated by least squares (LS) mean. Since Partek's alternative-splicing ANOVA algorithm is affected by exon-to-exon differences, Partek's gene view tool was used to determine alternatively-spliced genes.
Pearson correlation test was used to examine possible correlations between expression levels of the differentially expressed transcripts, which were employed for expression-based network (Expression correlation plug-in at Cytoscape software 2.6.1 ). To determine functional significance among the changed transcripts, over-represented gene ontology (GO) and KEGG pathways were identified, compared to their representation on the Human Exon 1.0 ST array, using the online software WebGestalt (http://bioinfo.vanderbilt.edu/webgestalt/ ).
Quantitative RT-PCR analysis
cDNA was synthesized from 0.5 μg total RNA using High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA) in the presence of an RNase inhibitor. The expression levels of the following genes were validated: CD19 (Hs01047409_g1), CD22 (Hs00233533_m1), CD79A (Hs00998120_g1), EREG (Hs00914313_m1), KCNH8 (Hs00976951_m1), IGHD (Hs00378878_m1), IGHM (Hs00378435_m1), PAX5 (Hs01045956_m1) and PCDH9 (Hs01009199_m1). The expression of SNCA (Hs01103386_m1) was used here as negative control, since its expression levels were not changed between the studied groups in the Affymetrix Exon 1.0 array experiment (data not shown). The expression levels of all tested genes were normalized relative to the expression levels of the GUSB (Hs99999908_m1) and TBP (Hs99999910_m1) genes. These reference genes were selected based on their microarray expression out of several candidate reference genes using NormFinder algorithm, which allows the selection of genes with the least expression changes between the different samples and the experimental groups (, data not shown). Quantification was done using TaqMan MGB probe and TaqMan Universal PCR Master Mix with 10 ng of cDNA on the StepOne Plus Real-Time PCR system (Applied Biosystems). Reactions were performed for 2 min at 50°C, 10 min at 95°C, and then 45 cycles of 15 s at 95°C and 1 min at 60°C. The expression levels were determined using the comparative threshold cycle (CT) quantification method . The geometric mean of CT values of the internal control genes, GUSB and TBP, was subtracted from CT values of each target genes (ΔCT) and signal values are expressed as 2 (-ΔCT). Pearson correlation test was applied to examine the correlation between results obtained from the Exon 1.0 ST Array and the Real-Time PCR. The expanded confirmation analysis (in 79 samples from men and women) included the following genes: CD19, CD22, CD79A, IGHD, IGHM, PAX5 and SNCA. Their expression signal values were compared between the study groups, and 2-way ANOVA with gender and group as the two independent variables, followed by 1-way ANOVA statistics were applied.
Analyzing published blood expression profile data of Parkinson's and other neurological diseases
To detect whether the expression changes in B cell related genes seen in our PD patients can be also detected in other neurological diseases' patients and in other PD cohort, we re-analyzed the published data of Scherzer et al. (; GSE6613). This dataset includes 50 PD patients, 22 healthy controls and 33 controls with other neurological diseases. The CEL files were downloaded to Partek® Genomics Suite Version 6.5. Background adjustment and normalization of the microarray data were done using the GCRMA algorithm. One sample of healthy control (GSM155358) was detected as an outlier by PCA, and was removed from further analyses. Since the gender of the subjects was not provided, we used the expression of XIST to distinguish between males and females. According to the XIST expression, the cohort included 41 females (12 PD patients, 11 healthy controls and 18 other neurological diseases controls) and 63 males (38 PD patients, 11 healthy controls and 15 other neurological diseases controls). ANOVAs were done to detect significant changes between PD patients, healthy controls and controls with other neurological diseases among: all subjects, only females and only males.
We thank Ziv Gan-Or, Dr. Anat Mirelman, Chen Shifrin, Hila Kobo, Dr. Mali Gana-Weisz, Shiran Levi, Ilana Elroi, Orna Moore, and Ariela Hilel for their assistance, and Dr. Boris Tartakovsky for scientific advice. This work was funded by grants from Kahn and Wolfson Foundations, Israel; Tel Aviv Sourasky Medical Center Grant of Excellence; National Parkinson Foundation, Miami, USA, and by the Legacy Heritage Bio-Medical Program of the Israel Science Foundation (grant No. 1922/08).
- Shulman JM, De Jager PL: Evidence for a common pathway linking neurodegenerative diseases. Nat Genet. 2009, 41: 1261-1262. 10.1038/ng1209-1261.PubMedView ArticleGoogle Scholar
- de Rijk MC, Launer LJ, Berger K, Breteler MM, Dartigues JF, Baldereschi M, Fratiglioni L, Lobo A, Martinez-Lage J, Trenkwalder C, Hofman A: Prevalence of Parkinson's disease in Europe: A collaborative study of population-based cohorts. Neurologic Diseases in the Elderly Research Group. Neurology. 2000, 54: S21-23.PubMedGoogle Scholar
- Polymeropoulos MH, Lavedan C, Leroy E, Ide SE, Dehejia A, Dutra A, Pike B, Root H, Rubenstein J, Boyer R, et al: Mutation in the alpha-synuclein gene identified in families with Parkinson's disease. Science. 1997, 276: 2045-2047. 10.1126/science.276.5321.2045.PubMedView ArticleGoogle Scholar
- Lesage S, Brice A: Parkinson's disease: from monogenic forms to genetic susceptibility factors. Hum Mol Genet. 2009, 18: R48-59. 10.1093/hmg/ddp012.PubMedView ArticleGoogle Scholar
- Shin JH, Dawson VL, Dawson TM: SnapShot: pathogenesis of Parkinson's disease. Cell. 2009, 139: 440-e441-442PubMedView ArticleGoogle Scholar
- Simunovic F, Yi M, Wang Y, Macey L, Brown LT, Krichevsky AM, Andersen SL, Stephens RM, Benes FM, Sonntag KC: Gene expression profiling of substantia nigra dopamine neurons: further insights into Parkinson's disease pathology. Brain. 2009, 132: 1795-1809. 10.1093/brain/awn323.PubMedPubMed CentralView ArticleGoogle Scholar
- Hauser MA, Li YJ, Xu H, Noureddine MA, Shao YS, Gullans SR, Scherzer CR, Jensen RV, McLaurin AC, Gibson JR, et al: Expression profiling of substantia nigra in Parkinson disease, progressive supranuclear palsy, and frontotemporal dementia with parkinsonism. Arch Neurol. 2005, 62: 917-921. 10.1001/archneur.62.6.917.PubMedView ArticleGoogle Scholar
- Grunblatt E, Mandel S, Jacob-Hirsch J, Zeligson S, Amariglo N, Rechavi G, Li J, Ravid R, Roggendorf W, Riederer P, Youdim MB: Gene expression profiling of parkinsonian substantia nigra pars compacta; alterations in ubiquitin-proteasome, heat shock protein, iron and oxidative stress regulated proteins, cell adhesion/cellular matrix and vesicle trafficking genes. J Neural Transm. 2004, 111: 1543-1573. 10.1007/s00702-004-0212-1.PubMedView ArticleGoogle Scholar
- Scherzer CR, Eklund AC, Morse LJ, Liao Z, Locascio JJ, Fefer D, Schwarzschild MA, Schlossmacher MG, Hauser MA, Vance JM, et al: Molecular markers of early Parkinson's disease based on gene expression in blood. Proc Natl Acad Sci USA. 2007, 104: 955-960. 10.1073/pnas.0610204104.PubMedPubMed CentralView ArticleGoogle Scholar
- Mutez E, Larvor L, Lepretre F, Mouroux V, Hamalek D, Kerckaert JP, Perez-Tur J, Waucquier N, Vanbesien-Mailliot C, Duflot A, et al: Transcriptional profile of Parkinson blood mononuclear cells with LRRK2 mutation. Neurobiol Aging. 2010Google Scholar
- Thorgeirsson TE, Geller F, Sulem P, Rafnar T, Wiste A, Magnusson KP, Manolescu A, Thorleifsson G, Stefansson H, Ingason A, et al: A variant associated with nicotine dependence, lung cancer and peripheral arterial disease. Nature. 2008, 452: 638-642. 10.1038/nature06846.PubMedPubMed CentralView ArticleGoogle Scholar
- Orr-Urtreger A, Bar-Shira A, Bercovich D, Matarasso N, Rozovsky U, Rosner S, Soloviov S, Rennert G, Kadouri L, Hubert A, et al: RNASEL mutation screening and association study in Ashkenazi and non-Ashkenazi prostate cancer patients. Cancer Epidemiol Biomarkers Prev. 2006, 15: 474-479. 10.1158/1055-9965.EPI-05-0606.PubMedView ArticleGoogle Scholar
- Ostrer H: A genetic profile of contemporary Jewish populations. Nat Rev Genet. 2001, 2: 891-898.PubMedView ArticleGoogle Scholar
- Orr-Urtreger A, Shifrin C, Rozovski U, Rosner S, Bercovich D, Gurevich T, Yagev-More H, Bar-Shira A, Giladi N: The LRRK2 G2019S mutation in Ashkenazi Jews with Parkinson disease: is there a gender effect?. Neurology. 2007, 69: 1595-1602. 10.1212/01.wnl.0000277637.33328.d8.PubMedView ArticleGoogle Scholar
- Gan-Or Z, Giladi N, Rozovski U, Shifrin C, Rosner S, Gurevich T, Bar-Shira A, Orr-Urtreger A: Genotype-phenotype correlations between GBA mutations and Parkinson disease risk and onset. Neurology. 2008, 70: 2277-2283. 10.1212/01.wnl.0000304039.11891.29.PubMedView ArticleGoogle Scholar
- Haaxma CA, Bloem BR, Borm GF, Oyen WJ, Leenders KL, Eshuis S, Booij J, Dluzen DE, Horstink MW: Gender differences in Parkinson's disease. J Neurol Neurosurg Psychiatry. 2007, 78: 819-824. 10.1136/jnnp.2006.103788.PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor KS, Cook JA, Counsell CE: Heterogeneity in male to female risk for Parkinson's disease. J Neurol Neurosurg Psychiatry. 2007, 78: 905-906. 10.1136/jnnp.2006.104695.PubMedPubMed CentralView ArticleGoogle Scholar
- Uc EY, McDermott MP, Marder KS, Anderson SW, Litvan I, Como PG, Auinger P, Chou KL, Growdon JC: Incidence of and risk factors for cognitive impairment in an early Parkinson disease clinical trial cohort. Neurology. 2009, 73: 1469-1477. 10.1212/WNL.0b013e3181bf992f.PubMedPubMed CentralView ArticleGoogle Scholar
- Cantuti-Castelvetri I, Keller-McGandy C, Bouzou B, Asteris G, Clark TW, Frosch MP, Standaert DG: Effects of gender on nigral gene expression and parkinson disease. Neurobiol Dis. 2007, 26: 606-614. 10.1016/j.nbd.2007.02.009.PubMedPubMed CentralView ArticleGoogle Scholar
- Duke DC, Moran LB, Kalaitzakis ME, Deprez M, Dexter DT, Pearce RK, Graeber MB: Transcriptome analysis reveals link between proteasomal and mitochondrial pathways in Parkinson's disease. Neurogenetics. 2006, 7: 139-148. 10.1007/s10048-006-0033-5.PubMedView ArticleGoogle Scholar
- Moran LB, Duke DC, Deprez M, Dexter DT, Pearce RK, Graeber MB: Whole genome expression profiling of the medial and lateral substantia nigra in Parkinson's disease. Neurogenetics. 2006, 7: 1-11. 10.1007/s10048-005-0020-2.PubMedView ArticleGoogle Scholar
- Zhang Y, James M, Middleton FA, Davis RL: Transcriptional analysis of multiple brain regions in Parkinson's disease supports the involvement of specific protein processing, energy metabolism, and signaling pathways, and suggests novel disease mechanisms. Am J Med Genet B Neuropsychiatr Genet. 2005, 137B: 5-16. 10.1002/ajmg.b.30195.PubMedView ArticleGoogle Scholar
- Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.PubMedView ArticleGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.PubMedView ArticleGoogle Scholar
- Song Y, Dayalu R, Matthews SA, Scharenberg AM: TRPML cation channels regulate the specialized lysosomal compartment of vertebrate B-lymphocytes. Eur J Cell Biol. 2006, 85: 1253-1264. 10.1016/j.ejcb.2006.08.004.PubMedView ArticleGoogle Scholar
- McKenna F, McLaughlin PJ, Lewis BJ, Sibbring GC, Cummerson JA, Bowen-Jones D, Moots RJ: Dopamine receptor expression on human T- and B-lymphocytes, monocytes, neutrophils, eosinophils and NK cells: a flow cytometric study. J Neuroimmunol. 2002, 132: 34-40. 10.1016/S0165-5728(02)00280-1.PubMedView ArticleGoogle Scholar
- Josefsson E, Bergquist J, Ekman R, Tarkowski A: Catecholamines are synthesized by mouse lymphocytes and regulate function of these cells by induction of apoptosis. Immunology. 1996, 88: 140-146. 10.1046/j.1365-2567.1996.d01-653.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Caronti B, Antonini G, Calderaro C, Ruggieri S, Palladini G, Pontieri FE, Colosimo C: Dopamine transporter immunoreactivity in peripheral blood lymphocytes in Parkinson's disease. J Neural Transm. 2001, 108: 803-807. 10.1007/s007020170030.PubMedView ArticleGoogle Scholar
- Barbanti P, Fabbrini G, Ricci A, Cerbo R, Bronzetti E, Caronti B, Calderaro C, Felici L, Stocchi F, Meco G, et al: Increased expression of dopamine receptors on lymphocytes in Parkinson's disease. Mov Disord. 1999, 14: 764-771. 10.1002/1531-8257(199909)14:5<764::AID-MDS1008>3.0.CO;2-W.PubMedView ArticleGoogle Scholar
- Hirsch EC, Hunot S: Neuroinflammation in Parkinson's disease: a target for neuroprotection?. Lancet Neurol. 2009, 8: 382-397. 10.1016/S1474-4422(09)70062-6.PubMedView ArticleGoogle Scholar
- Bas J, Calopa M, Mestre M, Mollevi DG, Cutillas B, Ambrosio S, Buendia E: Lymphocyte populations in Parkinson's disease and in rat models of parkinsonism. J Neuroimmunol. 2001, 113: 146-152. 10.1016/S0165-5728(00)00422-7.PubMedView ArticleGoogle Scholar
- Richartz-Salzburger E, Batra A, Stransky E, Laske C, Kohler N, Bartels M, Buchkremer G, Schott K: Altered lymphocyte distribution in Alzheimer's disease. J Psychiatr Res. 2007, 41: 174-178. 10.1016/j.jpsychires.2006.01.010.PubMedView ArticleGoogle Scholar
- Mantovani S, Garbelli S, Pasini A, Alimonti D, Perotti C, Melazzini M, Bendotti C, Mora G: Immune system alterations in sporadic amyotrophic lateral sclerosis patients suggest an ongoing neuroinflammatory process. J Neuroimmunol. 2009, 210: 73-79. 10.1016/j.jneuroim.2009.02.012.PubMedView ArticleGoogle Scholar
- Wakabayashi K, Mori F, Tanji K, Orimo S, Takahashi H: Involvement of the peripheral nervous system in synucleinopathies, tauopathies and other neurodegenerative proteinopathies of the brain. Acta Neuropathol. 120: 1-12.Google Scholar
- Djaldetti R, Lev N, Melamed E: Lesions outside the CNS in Parkinson's disease. Mov Disord. 2009, 24: 793-800. 10.1002/mds.22172.PubMedView ArticleGoogle Scholar
- Jentsch-Ullrich K, Koenigsmann M, Mohren M, Franke A: Lymphocyte subsets' reference ranges in an age- and gender-balanced population of 100 healthy adults--a monocentric German study. Clin Immunol. 2005, 116: 192-197. 10.1016/j.clim.2005.03.020.PubMedView ArticleGoogle Scholar
- Reichert T, DeBruyere M, Deneys V, Totterman T, Lydyard P, Yuksel F, Chapel H, Jewell D, Van Hove L, Linden J, et al: Lymphocyte subset reference ranges in adult Caucasians. Clin Immunol Immunopathol. 1991, 60: 190-208. 10.1016/0090-1229(91)90063-G.PubMedView ArticleGoogle Scholar
- Huppert FA, Solomou W, O'Connor S, Morgan K, Sussams P, Brayne C: Aging and lymphocyte subpopulations: whole-blood analysis of immune markers in a large population sample of healthy elderly individuals. Exp Gerontol. 1998, 33: 593-600. 10.1016/S0531-5565(98)00033-3.PubMedView ArticleGoogle Scholar
- Andreadou E, Nikolaou C, Gournaras F, Rentzos M, Boufidou F, Tsoutsou A, Zournas C, Zissimopoulos V, Vassilopoulos D: Serum uric acid levels in patients with Parkinson's disease: their relationship to treatment and disease duration. Clin Neurol Neurosurg. 2009, 111: 724-728. 10.1016/j.clineuro.2009.06.012.PubMedView ArticleGoogle Scholar
- Simunovic F, Yi M, Wang Y, Stephens R, Sonntag KC: Evidence for gender-specific transcriptional profiles of nigral dopamine neurons in Parkinson disease. PLoS One. 2010, 5: e8856-10.1371/journal.pone.0008856.PubMedPubMed CentralView ArticleGoogle Scholar
- Gan-Or Z, Bar-Shira A, Mirelman A, Gurevich T, Kedmi M, Giladi N, Orr-Urtreger A: LRRK2 and GBA mutations differentially affect the initial presentation of Parkinson disease. Neurogenetics. 2010, 11: 121-125. 10.1007/s10048-009-0198-9.PubMedView ArticleGoogle Scholar
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, et al: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.PubMedView ArticleGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.PubMedView ArticleGoogle Scholar
- Benovoy D, Kwan T, Majewski J: Effect of polymorphisms within probe-target sequences on olignonucleotide microarray experiments. Nucleic Acids Res. 2008, 36: 4417-4423. 10.1093/nar/gkn409.PubMedPubMed CentralView ArticleGoogle Scholar
- Duan S, Zhang W, Bleibel WK, Cox NJ, Dolan ME: SNPinProbe_1.0: A database for filtering out probes in the Affymetrix GeneChip(R) Human Exon 1.0 ST array potentially affected by SNPs. Bioinformation. 2008, 2: 469-470.PubMedPubMed CentralView ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang B, Kirov S, Snoddy J: WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 2005, 33: W741-748. 10.1093/nar/gki475.PubMedPubMed CentralView ArticleGoogle Scholar
- Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64: 5245-5250. 10.1158/0008-5472.CAN-04-0496.PubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
- Pfisterer P, Konig H, Hess J, Lipowsky G, Haendler B, Schleuning WD, Wirth T: CRISP-3, a protein with homology to plant defense proteins, is expressed in mouse B cells under the control of Oct2. Mol Cell Biol. 1996, 16: 6160-6168.PubMedPubMed CentralView ArticleGoogle Scholar
- Leu CM, Davis RS, Gartland LA, Fine WD, Cooper MD: FcRH1: an activation coreceptor on human B cells. Blood. 2005, 105: 1121-1126.PubMedView ArticleGoogle Scholar
- Airas L, Jalkanen S: CD73 mediates adhesion of B cells to follicular dendritic cells. Blood. 1996, 88: 1755-1764.PubMedGoogle Scholar
- Visan I, Goller M, Berberich I, Kneitz C, Tony HP: Pax-5 is a key regulator of the B cell-restricted expression of the CD23a isoform. Eur J Immunol. 2003, 33: 1163-1173. 10.1002/eji.200323747.PubMedView ArticleGoogle Scholar
- Ishii KJ, Kawagoe T, Koyama S, Matsui K, Kumar H, Kawai T, Uematsu S, Takeuchi O, Takeshita F, Coban C, Akira S: TANK-binding kinase-1 delineates innate and adaptive immune responses to DNA vaccines. Nature. 2008, 451: 725-729. 10.1038/nature06537.PubMedView ArticleGoogle 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.