Abstract
Until now, large-scale genome-wide association studies have identified 94 genes associated with Alzheimer's disease, Parkinson's disease, and multiple sclerosis. Expression quantitative trait locus (eQTL) analysis showed that six genetic variants around six of these 94 genes could drive both disease susceptibility and altered expression of six nearby genes including CD33 (rs3865444), PILRB (rs1476679), NUP160 (rs10838725), LRRK2 (rs76904798), RGS1 (rs1323292), and METTL21B (rs701006). However, two of these six genetic variants rs1476679 and rs76904798 variants could regulate the expression of PILRB and LRRK2 only in the human monocyte-derived microglia-like (MDMi) cells, but not in human peripheral blood monocytes. Here, we aim to verify these findings using another two eQTL datasets in human peripheral blood immune cell CD14+ monocytes. The results that showed that rs1476679 and rs76904798 variants or their proxy variants could significantly regulate the expression of PILRB and LRRK2 in immune cell CD14+ monocytes and human peripheral blood. We believe that these findings provide important supplementary information about the regulatory mechanisms by which both variants influence PILRB and LRRK2 gene expression and neurodegenerative disease risk.
Introduction
It is known that Alzheimer's disease, Parkinson's disease and multiple sclerosis are three common neurodegenerative diseases (Liu et al., , , ,,,, ; Beecham et al., ; Chang et al., ; Jun et al., ). In recent years, large-scale genome-wide association studies (GWAS) have identified some neurodegenerative disease risk variants (Liu et al., , , ,,,, ; Beecham et al., ; Zhang et al., , ; Li et al., ; Chang et al., ; Jiang et al., ; Jun et al., ). In 2013, the International Genomics of Alzheimer's Project (IGAP) conducted a meta-analysis of Alzheimer's disease GWAS datasets including 17,008 Alzheimer's disease cases and 37,154 controls in stage 1, and 8,572 Alzheimer's disease cases and 11,312 controls in stage 2. The meta-analysis of 74,046 individuals in stage 1 and stage 2 identified 11 new susceptibility loci for Alzheimer's disease (Lambert et al., ). In 2014, Nalls et al. conducted a large-scale meta-analysis of Parkinson's disease GWAS data including up to 13,708 PD cases and 95,282 control from 15 independent GWAS datasets of European descent, and identified six new Parkinson's disease risk loci (Nalls et al., ). In 2017, Chang et al. conducted a meta-analysis of Parkinson's disease GWAS data including 26,035 cases and 403,190 controls, and identified 17 novel Parkinson's disease risk loci (Chang et al., ).
Evidence shows that genetic variants could modify gene expression and cause disease risk (Liu et al., , ,,). In a recent study, Ryan and colleagues developed a human monocyte-derived microglia-like (MDMi) cellular model to assess the effects of neurodegenerative disease variants (Ryan et al., ). Using this model system, Ryan and colleagues conducted an expression quantitative trait locus (eQTL) analysis to examine 94 genes associated with Alzheimer's disease, Parkinson's disease, and multiple sclerosis (Ryan et al., ). One eQTL analysis was performed using MDMi cells from 95 young, healthy subjects of European ancestry (Ryan et al., ). The other eQTL analysis was performed using human peripheral blood monocytes derived from 211 young, healthy subjects of European ancestry (Ryan et al., ). They identified that genetic variants could drive both disease susceptibility and altered expression of six nearby genes including CD33 (rs3865444), PILRB (rs1476679), NUP160 (rs10838725), LRRK2 (rs76904798), RGS1 (rs1323292), and METTL21B (rs701006) (Ryan et al., ). However, rs1476679 and rs76904798 variants could regulate the expression of PILRB and LRRK2 only in the MDMi cells, but not in human peripheral blood monocytes (Ryan et al., ). Here, we aim to further verify their findings using multiple eQTL datasets in human peripheral blood immune cell CD14+ monocytes.
Materials and Methods
Functional Analysis
It is described that enhancers are DNA regulatory sequences, and could regulate tissue-specific gene expression (Ong and Corces, ; Shlyueva et al., ). Hence, we first performed an enhancer analysis using HaploReg (version 4.1) to evaluate the overlap of these six variants with predicted enhancers in each reference epigenome data from Roadmap Epigenomics project (Ward and Kellis, ). HaploReg (version 4.1) defined the enhancers using four different methods including the 15-state core model, the 25-state model incorporating imputed epigenomes, the H3K4me1/H3K4me3 peaks and the H3K27ac/H3K9ac peaks (Ward and Kellis, ). In brief, both the core 15-state model and the 25-state model are based on the imputed marks, H3K4me1/H3K4me3 is based on peaks from H3K4me1 and H3K4me3, and H3K27ac/H3K9ac is based on peaks from H3K27ac and H3K9ac (Ward and Kellis, ). In addition to enhancer analysis, we also conducted a promoter analysis and DNAse analysis. More detailed information has been described in recent study (Liu et al., ).
eQTL Analysis
The first eQTL dataset in monocytes is from the International Human Epigenome Consortium (IHEC) (Chen et al., ). The IHEC consortium performed an eQTL analysis in three immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals (Chen et al., ). Data from this project was produced in different institutes. Peripheral blood mononuclear cells were isolated from donors at University of Cambridge (Chen et al., ). The second eQTL dataset in monocytes (CD14+ monocytes) is from a recent study conducted by Fairfax et al. (). In the original study, primary monocytes from 432 healthy Europeans were exposed to interferon-γ (IFN-γ) or differing durations of lipopolysaccharide (LPS) (Fairfax et al., ). Then the impact of immune stimulation upon regulatory variant activity was systematically evaluated by an eQTL analysis (Fairfax et al., ). In brief, this eQTL analysis tested 609,704 genetic variants and expression data for 15,421 probes for 414 individuals in the naive state, 367 individuals after exposure to IFN-γ, 322 individuals after 24-h LPS, and 261 individuals after 2-h LPS (Fairfax et al., ).
In both eQTL datasets, genetic variants mapping to within 1 Mb (on each side) of each tested gene, were tested their association with gene expression using a linear regression analysis (Fairfax et al., ; Chen et al., ). Here, we utilize the summary results from the corresponding studies to evaluate the potential association of six genetic variants with the expression of their nearby genes. If any variant is not available in above datasets, we used HaploReg (version 4) to identify it proxy variants based on the linkage disequilibrium (LD) information from the 1000 Genomes Project (EUR) (Ward and Kellis, ). We defined the genetic variants tagged by any genetic variant with r2> = 0.8 (Ward and Kellis, ). In addition to the monocytes, we also evaluated the association of six genetic variants with gene expression using two eQTL datasets in human peripheral blood including 2,116 and 5,257 samples, respectively (Joehanes et al., ; Zhernakova et al., ).
Results
Enhancer Analysis
The enhancer analysis using the 15-state core model, the 25-state model, the H3K4me1/H3K4me3 model and the H3K27ac/H3K9ac model showed that these 6 genetic variants are predicted to be mainly located in enhancers of blood cell types. These findings indicate that these 6 genetic variants are likely to regulate gene expression in the blood cell types, especially CD14+ monocytes. In brief, rs1476679 variant is predicted to be located in enhancers of monocytes-CD14+ RO01746 Primary Cells using the H3K4me1/H3K4me3 model and the H3K27ac/H3K9ac model (Table 1). The rs76904798 variant is predicted to be located in enhancers of monocytes-CD14+ RO01746 Primary Cells using the 15-state core model, the H3K4me1/H3K4me3 model and the H3K27ac/H3K9ac model. To validate the enhancer analysis, we then performed an eQTLs analysis using multiple datasets in CD14+ monocytes.
Table 1
| SNP | 15-state core model | 25-state model | H3K4me1/H3K4me3 model | H3K27ac/H3K9ac model |
|---|---|---|---|---|
| rs1323292 | Primary T cells from cord blood, Primary T helper memory cells from peripheral blood 2, Primary T helper naive cells from peripheral blood, Primary T helper naive cells from peripheral blood, Primary T helper memory cells from peripheral blood 1, Primary T helper cells PMA-I stimulated, Primary T helper cells from peripheral blood, Primary T regulatory cells from peripheral blood, Primary T CD8+ naive cells from peripheral blood, Primary T CD8+ memory cells from peripheral blood | Primary T cells from cord blood, Primary T cells from peripheral blood, Primary hematopoietic stem cells short term culture, Primary T helper memory cells from peripheral blood 2, Primary T helper naive cells from peripheral blood, Primary T helper naive cells from peripheral blood, Primary T helper memory cells from peripheral blood 1, Primary T helper cells PMA-I stimulated, Primary T helper 17 cells PMA-I stimulated, Primary T helper cells from peripheral blood, Primary T regulatory cells from peripheral blood, Primary T cells effector/memory enriched from peripheral blood, Primary Natural Killer cells from peripheral blood, Primary T CD8+ naive cells from peripheral blood, Primary T CD8+ memory cells from peripheral blood, Primary hematopoietic stem cells G-CSF-mobilized Female, Primary mononuclear cells from peripheral blood, Fetal Thymus, Dnd41 TCell Leukemia Cell Line, GM12878 Lymphoblastoid Cells, K562 Leukemia Cells | Adipose Derived Mesenchymal Stem Cell Cultured Cells, Primary T cells from cord blood, Primary T cells from peripheral blood, Primary T helper memory cells from peripheral blood 2, Primary T helper naive cells from peripheral blood, Primary T helper naive cells from peripheral blood, Primary T helper memory cells from peripheral blood 1, Primary T helper cells PMA-I stimulated, Primary T helper 17 cells PMA-I stimulated, Primary T helper cells from peripheral blood, Primary T regulatory cells from peripheral blood, Primary T cells effector/memory enriched from peripheral blood, Primary T CD8+ naive cells from peripheral blood, Primary T CD8+ memory cells from peripheral blood, Cortex derived primary cultured neurospheres, Foreskin Fibroblast Primary Cells skin02, Colon Smooth Muscle, Duodenum Mucosa, Duodenum Smooth Muscle, Rectal Mucosa Donor 29, Rectal Mucosa Donor 31, Small Intestine, Thymus, Dnd41 TCell Leukemia Cell Line, NHLF Lung Fibroblast Primary Cells | Primary T cells from peripheral blood, Primary T helper memory cells from peripheral blood 2, Primary T helper naive cells from peripheral blood, Primary T helper naive cells from peripheral blood, Primary T helper memory cells from peripheral blood 1, Primary T helper cells PMA-I stimulated, Primary T helper 17 cells PMA-I stimulated, Primary T helper cells from peripheral blood, Primary T regulatory cells from peripheral blood, Primary T cells effector/memory enriched from peripheral blood, Primary T CD8+ memory cells from peripheral blood, Primary mononuclear cells from peripheral blood, Colon Smooth Muscle, Duodenum Smooth Muscle, Esophagus, Fetal Thymus, Rectal Smooth Muscle, Stomach Smooth Muscle, Thymus, Dnd41 TCell Leukemia Cell Line. |
| rs701006 | Cortex derived primary cultured neurospheres, Foreskin Melanocyte Primary Cells skin03, Brain Cingulate Gyrus, Brain Hippocampus Middle, Brain Inferior Temporal Lobe | Cortex derived primary cultured neurospheres, Ganglion Eminence derived primary cultured neurospheres, Brain Angular Gyrus, Brain Anterior Caudate, Brain Cingulate Gyrus, Brain Hippocampus Middle, Brain Inferior Temporal Lobe, Brain_Dorsolateral_Prefrontal_Cortex, Brain Substantia Nigra | hESC Derived CD56+ Ectoderm Cultured Cells, Primary T helper memory cells from peripheral blood 2, Cortex derived primary cultured neurospheres, Foreskin Melanocyte Primary Cells skin03, Brain Angular Gyrus, Brain Anterior Caudate, Brain Cingulate Gyrus, Brain Hippocampus Middle, Brain Inferior Temporal Lobe, Brain Substantia Nigra, Placenta, Lung, Small Intestine | Primary T helper memory cells from peripheral blood 2, Primary T helper naive cells from peripheral blood, Foreskin Fibroblast Primary Cells skin02, Foreskin Melanocyte Primary Cells skin03, Brain Angular Gyrus, Brain Anterior Caudate, Brain Cingulate Gyrus, Brain Hippocampus Middle, Brain Inferior Temporal Lobe, Brain_Dorsolateral_Prefrontal_Cortex, Brain Substantia Nigra, Right Ventricle, Monocytes-CD14+ RO01746 Primary Cells |
| rs1476679 | hESC Derived CD184+ Endoderm Cultured Cells, Mesenchymal Stem Cell Derived Adipocyte Cultured Cells, Adipose Derived Mesenchymal Stem Cell Cultured Cells, Bone Marrow Derived Cultured Mesenchymal Stem Cells, Breast variant Human Mammary Epithelial Cells (vHMEC), Foreskin Fibroblast Primary Cells skin01, Foreskin Keratinocyte Primary Cells skin03, Fetal Thymus, Stomach Mucosa, GM12878 Lymphoblastoid Cells, HMEC Mammary Epithelial Primary Cells, HSMM Skeletal Muscle Myoblasts Cells, HSMM cell derived Skeletal Muscle Myotubes Cells, K562 Leukemia Cells, NH-A Astrocytes Primary Cells, Osteoblast Primary Cells | Mesenchymal Stem Cell Derived Adipocyte Cultured Cells, Adipose Derived Mesenchymal Stem Cell Cultured Cells, Breast variant Human Mammary Epithelial Cells (vHMEC), Mesenchymal Stem Cell Derived Chondrocyte Cultured Cells, Foreskin Fibroblast Primary Cells skin02, Foreskin Keratinocyte Primary Cells skin02, HSMM cell derived Skeletal Muscle Myotubes Cells, NHDF-Ad Adult Dermal Fibroblast Primary Cells, NHLF Lung Fibroblast Primary Cells | ES-I3 Cells, H1 Derived Mesenchymal Stem Cells, H9 Derived Neuronal Progenitor Cultured Cells, hESC Derived CD184+ Endoderm Cultured Cells, hESC Derived CD56+ Ectoderm Cultured Cells, hESC Derived CD56+ Mesoderm Cultured Cells, IMR90 fetal lung fibroblasts Cell Line, Mesenchymal Stem Cell Derived Adipocyte Cultured Cells, Adipose Derived Mesenchymal Stem Cell Cultured Cells, Bone Marrow Derived Cultured Mesenchymal Stem Cells, Breast variant Human Mammary Epithelial Cells (vHMEC), Primary monocytes from peripheral blood, Primary T cells from peripheral blood, Mesenchymal Stem Cell Derived Chondrocyte Cultured Cells, Foreskin Fibroblast Primary Cells skin01, Foreskin Fibroblast Primary Cells skin02, Foreskin Keratinocyte Primary Cells skin02, Foreskin Keratinocyte Primary Cells skin03, Brain Angular Gyrus, Colonic Mucosa, Fetal Thymus, Pancreas, Stomach Mucosa, GM12878 Lymphoblastoid Cells, HeLa-S3 Cervical Carcinoma Cell Line, HMEC Mammary Epithelial Primary Cells, HSMM Skeletal Muscle Myoblasts Cells, HSMM cell derived Skeletal Muscle Myotubes Cells, Monocytes-CD14+ RO01746 Primary Cells, NH-A Astrocytes Primary Cells, NHEK-Epidermal Keratinocyte Primary Cells, Osteoblast Primary Cells | hESC Derived CD56+ Mesoderm Cultured Cells, Monocytes-CD14+ RO01746 Primary Cells |
| rs10838725 | H1 Derived Mesenchymal Stem Cells, hESC Derived CD56+ Mesoderm Cultured Cells, IMR90 fetal lung fibroblasts Cell Line, Mesenchymal Stem Cell Derived Adipocyte Cultured Cells, ES-UCSF4 Cells, Adipose Derived Mesenchymal Stem Cell Cultured Cells, Bone Marrow Derived Cultured Mesenchymal Stem Cells, Breast Myoepithelial Primary Cells, Breast variant Human Mammary Epithelial Cells (vHMEC), Primary monocytes from peripheral blood, Primary B cells from peripheral blood, Primary T cells from peripheral blood, Primary T helper memory cells from peripheral blood 2, Primary T helper naive cells from peripheral blood, Primary T helper naive cells from peripheral blood, Primary T helper memory cells from peripheral blood 1, Primary T helper cells PMA-I stimulated, Primary T helper 17 cells PMA-I stimulated, Primary T helper cells from peripheral blood, Primary T regulatory cells from peripheral blood, Primary T cells effector/memory enriched from peripheral blood, Primary Natural Killer cells from peripheral blood, Primary T CD8+ naive cells from peripheral blood, Primary T CD8+ memory cells from peripheral blood, Primary hematopoietic stem cells G-CSF-mobilized Female, Primary hematopoietic stem cells G-CSF-mobilized Male, Cortex derived primary cultured neurospheres, Ganglion Eminence derived primary cultured neurospheres, Foreskin Fibroblast Primary Cells skin01, Foreskin Fibroblast Primary Cells skin02, Foreskin Keratinocyte Primary Cells skin02, | H1 Derived Mesenchymal Stem Cells, hESC Derived CD56+ Ectoderm Cultured Cells, hESC Derived CD56+ Mesoderm Cultured Cells, Primary T helper memory cells from peripheral blood 2, Primary T helper naive cells from peripheral blood, Primary T helper cells PMA-I stimulated, Primary T helper 17 cells PMA-I stimulated, Primary T regulatory cells from peripheral blood, Primary hematopoietic stem cells G-CSF-mobilized Female, Foreskin Fibroblast Primary Cells skin01, Foreskin Keratinocyte Primary Cells skin03, Foreskin Melanocyte Primary Cells skin03, Adipose Nuclei, Brain Angular Gyrus, Brain Anterior Caudate, Brain Inferior Temporal Lobe, Brain_Dorsolateral_Prefrontal_Cortex, Colon Smooth Muscle, Pancreatic Islets, Psoas Muscle, | ||
| Foreskin Keratinocyte Primary Cells skin03, Foreskin Melanocyte Primary Cells skin01, Foreskin Melanocyte Primary Cells skin03, Adipose Nuclei, Aorta, Liver, Brain Angular Gyrus, Brain Anterior Caudate, Brain Cingulate Gyrus, Brain Germinal Matrix, Brain Hippocampus Middle, Brain_Dorsolateral_Prefrontal_Cortex, Brain Substantia Nigra, Colonic Mucosa, Colon Smooth Muscle, Duodenum Mucosa, Duodenum Smooth Muscle, Fetal Adrenal Gland, Fetal Brain Male, Fetal Brain Female, Fetal Heart, Fetal Intestine Small, Pancreatic Islets, Fetal Lung, Fetal Muscle Trunk, Fetal Muscle Leg, Fetal Stomach, Fetal Thymus, Left Ventricle, Lung, Ovary, Pancreas, Psoas Muscle, Rectal Mucosa Donor 31, Rectal Smooth Muscle, Right Atrium, Right Ventricle, Sigmoid Colon, Skeletal Muscle Male, Skeletal Muscle Female, Small Intestine, Stomach Mucosa, Stomach Smooth Muscle, Spleen, A549 EtOH 0.02pct Lung Carcinoma Cell Line, HeLa-S3 Cervical Carcinoma Cell Line, HSMM cell derived Skeletal Muscle Myotubes Cells, Monocytes-CD14+ RO01746 Primary Cells, NH-A Astrocytes Primary Cells, NHEK-Epidermal Keratinocyte Primary Cells, NHLF Lung Fibroblast Primary Cells, Osteoblast Primary Cells | Monocytes-CD14+ RO01746 Primary Cells, NH-A Astrocytes Primary Cells, Osteoblast Primary Cells | |||
| rs76904798 | Primary monocytes from peripheral blood, Primary neutrophils from peripheral blood, Monocytes-CD14+ RO01746 Primary Cells | Primary monocytes from peripheral blood, Primary neutrophils from peripheral blood, Primary B cells from cord blood, Monocytes-CD14+ RO01746 Primary Cells | Monocytes-CD14+ RO01746 Primary Cells | |
| rs3865444 | Primary monocytes from peripheral blood, Primary neutrophils from peripheral blood, Primary hematopoietic stem cells short term culture, Primary T regulatory cells from peripheral blood, Primary hematopoietic stem cells G-CSF-mobilized Female, Primary hematopoietic stem cells G-CSF-mobilized Male | Primary B cells from cord blood, Primary B cells from peripheral blood, Primary T cells from cord blood, Primary hematopoietic stem cells, Primary T regulatory cells from peripheral blood, Primary Natural Killer cells from peripheral blood, Primary hematopoietic stem cells G-CSF-mobilized Female, Primary mononuclear cells from peripheral blood, Duodenum Mucosa, Duodenum Smooth Muscle, Fetal Thymus, Dnd41 TCell Leukemia Cell Line, K562 Leukemia Cells, Monocytes-CD14+ RO01746 Primary Cells | Primary monocytes from peripheral blood, Primary neutrophils from peripheral blood, Primary hematopoietic stem cells, Primary hematopoietic stem cells short term culture, Primary T regulatory cells from peripheral blood, Primary Natural Killer cells from peripheral blood, Primary hematopoietic stem cells G-CSF-mobilized Female, Primary hematopoietic stem cells G-CSF-mobilized Male, Adipose Nuclei, Liver, Fetal Brain Male, Fetal Thymus, Rectal Mucosa Donor 29, Spleen, K562 Leukemia Cells, Monocytes-CD14+ RO01746 Primary Cells | iPS DF 6.9 Cells, Primary monocytes from peripheral blood, Primary T helper memory cells from peripheral blood 2, Primary T CD8+ memory cells from peripheral blood, Primary hematopoietic stem cells G-CSF-mobilized Female, Brain Anterior Caudate, Brain Hippocampus Middle, Brain Substantia Nigra, Lung, Spleen, K562 Leukemia Cells, Monocytes-CD14+ RO01746 Primary Cells |
Enhancer analysis of six neurodegenerative disease variants using HaploReg (version 4.1).
CD14+ cells are bolded.
Promoter Analysis
All these four models including 15-state core model, the 25-state model, the H3K4me1/H3K4me3 model and the H3K27ac/H3K9ac model showed that the rs3865444 variant was predicted to be mainly located in promoter of blood cell types, especially in immune cell types. Three of these four models including the 15-state core model, the H3K4me1/H3K4me3 model and the H3K27ac/H3K9ac model indicated that the rs1323292 variant was predicted to be mainly located in promoter of blood cell types, especially in immune cell types. Table 2 provides the results from promoter analysis of six neurodegenerative disease variants using HaploReg (version 4.1).
Table 2
| SNP | 15-state core model | 25-state model | H3K4me1/H3K4me3 model | H3K27ac/H3K9ac model |
|---|---|---|---|---|
| rs1323292 | Dnd41 TCell Leukemia Cell Line | Primary T helper naive cells from peripheral blood, Primary T helper memory cells from peripheral blood 1, Primary T helper cells PMA-I stimulated, Primary T regulatory cells from peripheral blood, Colon Smooth Muscle, Duodenum Smooth Muscle, Dnd41 TCell Leukemia Cell Line | Primary mononuclear cells from peripheral blood, Liver, Duodenum Mucosa, Rectal Mucosa Donor 31, Dnd41 TCell Leukemia Cell Line | |
| rs701006 | Brain Cingulate Gyrus | Primary T helper naive cells from peripheral blood, Brain Angular Gyrus, Brain Cingulate Gyrus, Brain Inferior Temporal Lobe | ||
| rs1476679 | Foreskin Fibroblast Primary Cells skin02 | H1 Derived Mesenchymal Stem Cells, H9 Cells, iPS-15b Cells, Mesenchymal Stem Cell Derived Adipocyte Cultured Cells, Adipose Derived Mesenchymal Stem Cell Cultured Cells, Breast Myoepithelial Primary Cells, Mesenchymal Stem Cell Derived Chondrocyte Cultured Cells, Brain Anterior Caudate, Brain Inferior Temporal Lobe, Colonic Mucosa, Stomach Mucosa, NH-A Astrocytes Primary Cells | ||
| rs10838725 | Brain Angular Gyrus, Colon Smooth Muscle, Small Intestine | Adipose Nuclei, Brain Anterior Caudate, Fetal Heart, Skeletal Muscle Male | ||
| rs76904798 | NA | NA | NA | NA |
| rs3865444 | hESC Derived CD184+ Endoderm Cultured Cells, Primary hematopoietic stem cells, HeLa-S3 Cervical Carcinoma Cell Line, K562 Leukemia Cells, Monocytes-CD14+ RO01746 Primary Cells | hESC Derived CD184+ Endoderm Cultured Cells, Primary monocytes from peripheral blood, Primary neutrophils from peripheral blood, Primary hematopoietic stem cells short term culture, Primary T helper memory cells from peripheral blood 2, Primary T helper naive cells from peripheral blood, Primary T helper naive cells from peripheral blood, Primary T helper cells from peripheral blood, Primary T CD8+ naive cells from peripheral blood, Primary hematopoietic stem cells G-CSF-mobilized Male, Adipose Nuclei, Liver, Fetal Brain Male, HeLa-S3 Cervical Carcinoma Cell Line | hESC Derived CD184+ Endoderm Cultured Cells, Primary neutrophils from peripheral blood, Primary hematopoietic stem cells, Primary hematopoietic stem cells short term culture, Primary hematopoietic stem cells G-CSF-mobilized Female, Primary hematopoietic stem cells G-CSF-mobilized Male, Primary mononuclear cells from peripheral blood, Liver, Brain Cingulate Gyrus, Brain Hippocampus Middle, Brain_Dorsolateral_Prefrontal_Cortex, Duodenum Mucosa, Fetal Adrenal Gland, Placenta, Fetal Thymus, Gastric, Lung, Pancreas, HeLa-S3 Cervical Carcinoma Cell Line, K562 Leukemia Cells, Monocytes-CD14+ RO01746 Primary Cells | Primary mononuclear cells from peripheral blood, HeLa-S3 Cervical Carcinoma Cell Line, K562 Leukemia Cells, Monocytes-CD14+ RO01746 Primary Cells |
Promoter analysis of six neurodegenerative disease variants using HaploReg (version 4.1).
DNAse Analysis
The DNAse analysis showed that rs1323292 variant was predicted to be mainly located in DNAse of blood cell types, especially in primary T cells from cord blood, primary T cells from peripheral blood, primary Natural Killer cells from peripheral blood, fetal Thymus, and small Intestine. The rs701006 variant was predicted to be mainly located in DNAse of Foreskin Melanocyte Primary Cells skin01 and Fetal Brain Female. The rs1476679 variant was predicted to be mainly located in DNAse of iPS DF 19.11 Cells, Breast variant Human Mammary Epithelial Cells (vHMEC), HSMM cell derived Skeletal Muscle Myotubes Cells. However, we did identify any rs10838725, rs76904798, rs3865444.
eQTLs Analysis in CD14+ Monocytes
In the first dataset including up to 197 individuals, the results showed significant association of all these six genetic variants (rs3865444, rs1476679, rs10838725, rs76904798, rs1323292, and rs701006) with the expression of nearby genes including CD33, PILRB, NUP160, LRRK2, RGS1, and METTL21B. In addition, there are also some other genes. Here, we list the significant results with P < 0.01 and the corresponding false discovery rate (FDR) in Table 3.
Table 3
| SNP | Chr | Phenotype or probe ID | Gene | P-value | FDR | Dataset | Reference |
|---|---|---|---|---|---|---|---|
| rs1323292* | 1 | ENSG00000090104.7 | RGS1 | 3.07E-34 | 6.25E-31 | CD14+ monocytes | Chen et al., |
| rs701006 | 12 | ENSG00000123427.11 | METTL21B | 3.67E-26 | 1.55E-24 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000121716.12 | PILRB | 6.74E-20 | 5.33E-19 | CD14+ monocytes | Chen et al., |
| rs10838725 | 11 | ENSG00000030066.9 | NUP160 | 5.73E-18 | 6.15E-18 | CD14+ monocytes | Chen et al., |
| rs76904798 | 12 | ENSG00000188906.9 | LRRK2 | 9.42E-15 | 3.48E-11 | CD14+ monocytes | Chen et al., |
| rs3865444 | 19 | ENSG00000105383.10 | CD33 | 5.26E-14 | 6.09E-11 | CD14+ monocytes | Chen et al., |
| rs76904798 | 12 | ENSG00000260943.1 | RP11-476D10.1 | 3.84E-12 | 4.71E-11 | CD14+ monocytes | Chen et al., |
| rs3865444 | 19 | ENSG00000268849.1 | SIGLEC22P | 1.72E-07 | 5.23E-05 | CD14+ monocytes | Chen et al., |
| rs10838725 | 11 | ENSG00000134571.6 | MYBPC3 | 1.79E-07 | 2.18E-07 | CD14+ monocytes | Chen et al., |
| rs10838725 | 11 | ENSG00000110514.14 | MADD | 2.40E-06 | 5.68E-06 | CD14+ monocytes | Chen et al., |
| rs701006 | 12 | ENSG00000123297.11 | TSFM | 2.48E-06 | 8.05E-05 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000221838.5 | AP4M1 | 9.93E-06 | 5.90E-05 | CD14+ monocytes | Chen et al., |
| rs76904798 | 12 | ENSG00000229899.1 | AC084290.2 | 1.04E-05 | 7.97E-04 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000078319.7 | PMS2P1 | 5.41E-05 | 8.76E-04 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000166997.3 | CNPY4 | 8.00E-05 | 2.58E-03 | CD14+ monocytes | Chen et al., |
| rs701006 | 12 | ENSG00000135506.11 | OS9 | 2.50E-04 | 2.63E-02 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000087077.7 | TRIP6 | 4.80E-04 | 7.17E-03 | CD14+ monocytes | Chen et al., |
| rs76904798 | 12 | ENSG00000225342.1 | AC079630.4 | 5.25E-04 | 2.49E-02 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000196411.5 | EPHB4 | 8.54E-04 | 5.47E-03 | CD14+ monocytes | Chen et al., |
| rs701006 | 12 | ENSG00000135439.6 | AGAP2 | 8.61E-04 | 1.03E-02 | CD14+ monocytes | Chen et al., |
| rs3865444 | 19 | ENSG00000142512.10 | SIGLEC10 | 9.02E-04 | 1.47E-01 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000066923.12 | STAG3 | 3.09E-03 | 1.10E-02 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000085514.10 | PILRA | 3.10E-03 | 7.72E-02 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000106290.10 | TAF6 | 6.23E-03 | 3.42E-02 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000106330.7 | MOSPD3 | 8.13E-03 | 3.15E-02 | CD14+ monocytes | Chen et al., |
| rs1476679 | 7 | ENSG00000146830.8 | GIGYF1 | 9.97E-03 | 4.92E-02 | CD14+ monocytes | Chen et al., |
| rs701006 | 12 | 870056 | METTL21B | 4.79E-63 | 1.95E-59 | CD14+ monocytes | Fairfax et al., |
| rs1476679 | 7 | 2190541 | GATS | 1.14E-20 | 5.73E-18 | CD14+ monocytes | Fairfax et al., |
| rs701006 | 12 | 5670577 | TSPAN31 | 3.25E-15 | 1.04E-12 | CD14+ monocytes | Fairfax et al., |
| rs11175655 | 12 | 3520291 | LRRK2 | 6.66E-12 | 1.48E-09 | CD14+ monocytes | Fairfax et al., |
| rs11175655 | 12 | 2490315 | HS.306876 | 3.06E-10 | 5.51E-08 | CD14+ monocytes | Fairfax et al., |
| rs701006 | 12 | 5290239 | XRCC6BP1 | 2.68E-09 | 4.22E-07 | CD14+ monocytes | Fairfax et al., |
| rs701006 | 12 | 3130102 | TSFM | 3.24E-09 | 5.04E-07 | CD14+ monocytes | Fairfax et al., |
| rs1476679 | 7 | 2810674 | TRIM4 | 3.59E-09 | 5.54E-07 | CD14+ monocytes | Fairfax et al., |
| rs1323298 | 1 | 4490176 | RGS1 | 4.91E-09 | 7.45E-07 | CD14+ monocytes | Fairfax et al., |
| rs1476679 | 7 | 1570039 | PILRB | 9.85E-08 | 1.21E-05 | CD14+ monocytes | Fairfax et al., |
| rs1476679 | 7 | 2030088 | PILRA | 1.88E-06 | 1.82E-04 | CD14+ monocytes | Fairfax et al., |
| rs11175655 | 12 | 1820725 | LRRK2 | 3.26E-05 | 2.40E-03 | CD14+ monocytes | Fairfax et al., |
| rs1476679 | 7 | 2750168 | AP4M1 | 4.30E-05 | 3.07E-03 | CD14+ monocytes | Fairfax et al., |
| rs1476679 | 7 | 5420672 | PILRA | 6.70E-03 | 1.83E-01 | CD14+ monocytes | Fairfax et al., |
| rs701006 | 12 | 870056 | METTL21B | 6.35E-51 | 3.02E-47 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs1476679 | 7 | 2190541 | GATS | 2.22E-26 | 2.39E-23 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs11175655 | 12 | 3520291 | LRRK2 | 2.70E-24 | 2.51E-21 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs701006 | 12 | 2900725 | CYP27B1 | 1.73E-19 | 1.12E-16 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs701006 | 12 | 650348 | METTL1 | 2.68E-16 | 1.26E-13 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs11175655 | 12 | 2490315 | HS.306876 | 7.94E-13 | 2.54E-10 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs1476679 | 7 | 2470577 | ZCWPW1 | 2.31E-10 | 5.36E-08 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs1323298 | 1 | 4490176 | RGS1 | 3.98E-08 | 6.43E-06 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs701006 | 12 | 5670577 | TSPAN31 | 3.96E-07 | 5.30E-05 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs1476679 | 7 | 2810674 | TRIM4 | 8.73E-07 | 1.09E-04 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs1476679 | 7 | 2030088 | PILRA | 1.29E-06 | 1.55E-04 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs701006 | 12 | 5290239 | XRCC6BP1 | 1.75E-06 | 2.05E-04 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs701006 | 12 | 3130102 | TSFM | 4.19E-05 | 3.48E-03 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs701006 | 12 | 110180 | ARHGAP9 | 7.18E-05 | 5.57E-03 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs701006 | 12 | 7380110 | CDK4 | 3.91E-04 | 2.34E-02 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs1476679 | 7 | 4850487 | UFSP1 | 6.45E-04 | 3.52E-02 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs1476679 | 7 | 5420672 | PILRA | 7.58E-04 | 4.01E-02 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs1476679 | 7 | 2750168 | AP4M1 | 2.72E-03 | 1.06E-01 | CD14+ monocytes after 24-h LPS | Fairfax et al., |
| rs701006 | 12 | 870056 | METTL21B | 1.11E-39 | 5.75E-36 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs1476679 | 7 | 2030088 | PILRA | 1.41E-10 | 5.20E-08 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs1476679 | 7 | 5420672 | PILRA | 8.96E-10 | 2.90E-07 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs701006 | 12 | 2900725 | CYP27B1 | 5.71E-09 | 1.61E-06 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs1476679 | 7 | 650040 | PILRA | 9.19E-09 | 2.49E-06 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs1476679 | 7 | 2190541 | GATS | 5.05E-08 | 1.18E-05 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs701006 | 12 | 650348 | METTL1 | 1.16E-06 | 2.04E-04 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs11175655 | 12 | 3520291 | LRRK2 | 3.67E-06 | 5.69E-04 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs701006 | 12 | 5670577 | TSPAN31 | 7.85E-05 | 8.25E-03 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs701006 | 12 | 5290239 | XRCC6BP1 | 1.69E-04 | 1.58E-02 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs1476679 | 7 | 1470195 | MCM7 | 1.91E-04 | 1.75E-02 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs1476679 | 7 | 2810674 | TRIM4 | 2.92E-04 | 2.48E-02 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs701006 | 12 | 7380110 | CDK4 | 4.49E-03 | 1.88E-01 | CD14+ monocytes after 2-h LPS | Fairfax et al., |
| rs701006 | 12 | 870056 | METTL21B | 1.30E-49 | 3.75E-46 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs11175655 | 12 | 3520291 | LRRK2 | 9.37E-23 | 5.54E-20 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs11175655 | 12 | 2490315 | HS.306876 | 4.12E-18 | 1.69E-15 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 650348 | METTL1 | 9.74E-13 | 2.33E-10 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1323298 | 1 | 4490176 | RGS1 | 1.72E-11 | 3.53E-09 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 5670577 | TSPAN31 | 3.56E-10 | 6.16E-08 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1476679 | 7 | 2030088 | PILRA | 8.04E-08 | 9.63E-06 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs11175655 | 12 | 1820725 | LRRK2 | 2.34E-07 | 2.59E-05 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 3130102 | TSFM | 8.53E-07 | 8.48E-05 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1476679 | 7 | 2810674 | TRIM4 | 2.07E-06 | 1.90E-04 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1476679 | 7 | 2190541 | GATS | 2.44E-06 | 2.20E-04 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 5290239 | XRCC6BP1 | 1.19E-05 | 9.19E-04 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 4610066 | TMEM194A | 6.51E-05 | 4.16E-03 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1476679 | 7 | 5420672 | PILRA | 1.29E-04 | 7.55E-03 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1476679 | 7 | 6200743 | CNPY4 | 5.44E-04 | 2.55E-02 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1476679 | 7 | 3940750 | LOC100134648 | 6.57E-04 | 2.98E-02 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 1030593 | OS9 | 7.07E-04 | 3.17E-02 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 1580431 | PIP4K2C | 1.09E-03 | 4.49E-02 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1476679 | 7 | 2470577 | ZCWPW1 | 1.23E-03 | 4.95E-02 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 2900725 | CYP27B1 | 1.29E-03 | 5.13E-02 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs701006 | 12 | 110180 | ARHGAP9 | 4.40E-03 | 1.29E-01 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
| rs1476679 | 7 | 650040 | PILRA | 6.69E-03 | 1.72E-01 | CD14+ monocytes after exposure to IFN-γ | Fairfax et al., |
eQTL analysis of six neurodegenerative disease variants in CD14+ monocytes.
The associations identified by Ryan and colleagues are bolded. The rs11175655 variant is in high LD with rs76904798 (r2 = 0.83, and D′ = 1), and rs1323298 is in high LD with rs1323292 (r2 = 0.92, and D′ = 0.97). Chr, chromosome; IFN-γ, interferon-γ; LPS, lipopolysaccharide.
In the second dataset including 414 samples in the naive state, two genetic variants rs3865444 and rs10838725, as well as their proxy variants (r2 ≥ 0.8), are not available. The other four genetic variants or their proxy variants were significantly associated with the expression of nearby genes including PILRB, LRRK2, RGS1, and METTL21B. In brief, rs11175655 is in high LD with rs76904798 (r2 = 0.83, and D′ = 1), and rs1323298 is in high LD with rs1323292 (r2 = 0.92, and D′ = 0.97). Meanwhile, there also some other genes, as provided in Table 2. In addition to the naive state, some findings were also observed in other three states including 367 individuals after exposure to IFN-γ, 322 individuals after 24-h LPS, and 261 individuals after 2-h LPS, as provided in Table 3.
QTLs Analysis in Human Peripheral Blood
In the two eQTL datasets in human peripheral blood, we again found significant association of all these six genetic variants (rs3865444, rs1476679, rs10838725, rs76904798, rs1323292, and rs701006) with the expression of nearby genes including CD33, PILRB, NUP160, LRRK2, RGS1, and METTL21B, as well as other genes (Table 4).
Table 4
| SNP | Position (hg19) | Probe ID | P-value | Sample size | Reference |
|---|---|---|---|---|---|
| rs1323292 | chr1:192541021 | RGS1 | 7.36E-05 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | METTL21B | 3.54E-224 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | TSFM;RP11-571M6.15 | 1.16E-40 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | AVIL | 1.80E-32 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | TSFM | 1.38E-25 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | XRCC6BP1 | 3.26E-21 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | OS9;RP11-571M6.7 | 4.95E-15 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | MARCH9 | 2.78E-13 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | AVIL;U6 | 4.08E-13 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | TSPAN31 | 4.79E-13 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | OS9 | 6.58E-12 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | RP11-571M6.18 | 1.26E-10 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | METTL1;METTL21B | 4.04E-10 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | CDK4;TSPAN31 | 1.90E-09 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | METTL21B;RP11-571M6.15 | 2.35E-09 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | METTL1;RP11-571M6.13 | 4.52E-05 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | AGAP2 | 7.06E-05 | 2116 | Zhernakova et al., |
| rs701006 | chr12:58106836 | METTL1 | 7.42E-05 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | PILRB | 9.88E-127 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | STAG3;PVRIG | 3.96E-58 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | STAG3;GATS | 4.26E-52 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | PILRA;PILRB | 3.91E-37 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | STAG3;GATS;GATS | 1.09E-32 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | STAG3;PVRIG;AC005071.1 | 2.05E-24 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | PILRB;CTB-161A2.4 | 4.97E-22 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | GAL3ST4 | 1.21E-18 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | ZCWPW1 | 1.38E-11 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | MOSPD3 | 2.82E-09 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | GPC2 | 2.92E-08 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | PILRB;CTB-161A2.3 | 5.97E-08 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | TFR2 | 1.98E-06 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | PILRA | 6.63E-06 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | C7orf61 | 4.12E-05 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | TSC22D4;C7orf61 | 4.34E-05 | 2116 | Zhernakova et al., |
| rs1476679 | chr7:100004446 | PPP1R35;RP11-758P17.2 | 5.99E-05 | 2116 | Zhernakova et al., |
| rs10838725 | chr11:47557871 | MYBPC3 | 3.05E-47 | 2116 | Zhernakova et al., |
| rs10838725 | chr11:47557871 | C1QTNF4 | 6.18E-44 | 2116 | Zhernakova et al., |
| rs10838725 | chr11:47557871 | MADD | 8.32E-25 | 2116 | Zhernakova et al., |
| rs10838725 | chr11:47557871 | FNBP4 | 2.20E-18 | 2116 | Zhernakova et al., |
| rs10838725 | chr11:47557871 | FNBP4;Y_RNA | 2.92E-15 | 2116 | Zhernakova et al., |
| rs10838725 | chr11:47557871 | RP11-750H9.5 | 5.00E-13 | 2116 | Zhernakova et al., |
| rs10838725 | chr11:47557871 | SLC39A13 | 1.16E-09 | 2116 | Zhernakova et al., |
| rs10838725 | chr11:47557871 | RAPSN | 1.33E-07 | 2116 | Zhernakova et al., |
| rs76904798 | chr12:40614434 | LRRK2 | 1.47E-52 | 2116 | Zhernakova et al., |
| rs76904798 | chr12:40614434 | RP11-476D10.1 | 1.93E-45 | 2116 | Zhernakova et al., |
| rs3865444 | chr19:51727962 | CD33 | 1.58E-44 | 2116 | Zhernakova et al., |
| rs3865444 | chr19:51727962 | SIGLEC22P | 5.55E-24 | 2116 | Zhernakova et al., |
| rs3865444 | chr19:51727962 | VSIG10L | 5.74E-06 | 2116 | Zhernakova et al., |
| rs1323292 | chr1:192541021 | SLC41A1 | 3.44E-05 | 5257 | Joehanes et al., |
| rs1323292 | chr1:192541021 | P4HTM | 5.15E-05 | 5257 | Joehanes et al., |
| rs701006 | chr12:58106836 | AVIL;PP12719; LOC100653271 | 2.45E-38 | 5257 | Joehanes et al., |
| rs701006 | chr12:58106836 | METTL21B | 9.31E-31 | 5257 | Joehanes et al., |
| rs701006 | chr12:58106836 | ARHGAP9 | 2.71E-09 | 5257 | Joehanes et al., |
| rs701006 | chr12:58106836 | TSPAN31 | 7.01E-07 | 5257 | Joehanes et al., |
| rs701006 | chr12:58106836 | OS9 | 7.77E-07 | 5257 | Joehanes et al., |
| rs701006 | chr12:58106836 | XRCC6BP1 | 2.87E-06 | 5257 | Joehanes et al., |
| rs701006 | chr12:58106836 | GLI1 | 3.96E-06 | 5257 | Joehanes et al., |
| rs701006 | chr12:58106836 | TSFM | 7.43E-05 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | PILRB | 6.16E-69 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | PVRIG;STAG3 | 1.23E-44 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | GATS | 1.56E-26 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | EPHB4 | 4.30E-12 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | TRIM4 | 1.87E-07 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | PILRA | 1.82E-06 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | CNPY4 | 1.00E-05 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | MAN2B1 | 1.38E-05 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | ZKSCAN1 | 1.84E-05 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | OR2AE1 | 2.45E-05 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | C3orf30;IGSF11 | 4.87E-05 | 5257 | Joehanes et al., |
| rs1476679 | chr7:100004446 | AQP4 | 5.83E-05 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | MADD | 7.68E-32 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | PTPRJ | 1.06E-22 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | MTCH2 | 1.64E-12 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | MYBPC3 | 6.88E-11 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | SLC39A13 | 5.96E-06 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | SBF2 | 1.90E-05 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | NDUFS3 | 2.86E-05 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | DKK4 | 3.64E-05 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | NUP160 | 6.14E-05 | 5257 | Joehanes et al., |
| rs10838725 | chr11:47557871 | FCER2 | 7.50E-05 | 5257 | Joehanes et al., |
| rs76904798 | chr12:40614434 | TANK | 3.77E-05 | 5257 | Joehanes et al., |
| rs3865444 | chr19:51727962 | ETFB | 2.21E-07 | 5257 | Joehanes et al., |
| rs3865444 | chr19:51727962 | PITPNM2 | 5.24E-05 | 5257 | Joehanes et al., |
eQTL analysis of six neurodegenerative disease variants in whole blood.
The associations identified by Ryan and colleagues are bolded.
Discussion
In recent years, large-scale GWAS datasets have identified 94 genes associated with Alzheimer's disease, Parkinson's disease, or multiple sclerosis (Ryan et al., ). However, it is still unclear how these variants functionally affect the underlying neurodegenerative disease pathogenesis. Growing evidence shows that genetic variants may affect disease risk by regulating gene expression (Bao et al., ; Liu et al., , , ; Hu et al., ). Ryan and colleagues applied a human MDMi cellular model, and conducted an eQTL analysis to evaluate the effects of these neurodegenerative disease variants (Ryan et al., ).
In summary, Ryan et al. identified that six neurodegenerative disease variants were associated with disease susceptibility, and could alter the expression of six nearby genes. Both rs1476679 and rs76904798 variants could only regulate the expression of PILRB and LRRK2 in the MDMi cells, but not in human peripheral blood monocytes (Ryan et al., ). Ryan et al. concluded that the differentiation of monocytes into microglia-like cells could cause the acquisition of a cellular state, which could reveal the functional consequences of certain genetic variants (Ryan et al., ). Ryan et al. provided an in vitro translational tool to generate microglia-like cells quickly and easily from adult blood, and could be useful for exploring microglia function and dysfunction (Ryan et al., ). However, Ryan et al. just selected single eQTL dataset in human peripheral blood monocytes without any replication.
Here, we first performed an enhancer analysis of these six neurodegenerative disease variants. The results showed that these six genetic variants were predicted to be mainly located in enhancers of blood cell types, especially in CD14+ monocytes. The promoter analysis showed that rs3865444 and rs1323292 variants were predicted to be mainly located in promoter of blood cell types, especially in immune cell types. Hence, we further evaluated these six genetic variants using multiple eQTL datasets in human peripheral blood immune cell CD14+ monocytes. The results that showed that rs1476679 and rs76904798 variants or their proxy variants could significantly regulate the expression of PILRB and LRRK2 in immune cell CD14+ monocytes and human peripheral blood.
In summary, we believe that these findings provide important supplementary information about the regulatory mechanisms by which both variants influence PILRB and LRRK2 gene expression and neurodegenerative disease risk.
Statements
Author contributions
BS and HY conceived and initiated the project. JS and YH analyzed the data, and wrote the first draft of the manuscript. YZ, LW, and LL contributed to the interpretation of the results and critical revision of the manuscript for important intellectual content and approved the final version of the manuscript.
Acknowledgments
We thank the International Human Epigenome Consortium (IHEC), Fairfax and colleagues for the eQTL datasets. This work was supported in part by funds from the National Natural Science Foundation of China (Grant No. 81870938) and Fund of Taishan scholar project (to BS).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
BaoX.LiuG.JiangY.JiangQ.LiaoM.FengR.et al. (2015). Cell adhesion molecule pathway genes are regulated by cis-regulatory SNPs and show significantly altered expression in Alzheimer's disease brains. Neurobiol. Aging36, 2904. e1–7. 10.1016/j.neurobiolaging.2015.06.006
2
BeechamA. H.PatsopoulosN. A.XifaraD. K.DavisM. F.KemppinenA.CotsapasC.et al. (2013). Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis. Nat. Genet.45, 1353–1360. 10.1038/ng.2770
3
ChangD.NallsM. A.HallgrimsdottirI. B.HunkapillerJ.Van Der BrugM.CaiF.et al. (2017). A meta-analysis of genome-wide association studies identifies 17 new Parkinson's disease risk loci. Nat. Genet.49, 1511–1516. 10.1038/ng.3955
4
ChenL.GeB.CasaleF. P.VasquezL.KwanT.Garrido-MartinD.et al. (2016). Genetic drivers of epigenetic and transcriptional variation in human immune cells. Cell167, 1398–1414e24. 10.1016/j.cell.2016.10.026
5
FairfaxB. P.HumburgP.MakinoS.NaranbhaiV.WongD.LauE.et al. (2014). Innate immune activity conditions the effect of regulatory variants upon monocyte gene expression. Science343:1246949. 10.1126/science.1246949
6
HuY.JinS.ChengL.LiuG.JiangQ. (2017). Autoimmune disease variants regulate GSDMB gene expression in human immune cells and whole blood. Proc. Natl. Acad. Sci. U.S.A.114, E7860–E7862. 10.1073/pnas.1712127114
7
JiangQ.JinS.JiangY.LiaoM.FengR.ZhangL.et al. (2017). Alzheimer's disease variants with the genome-wide significance are significantly enriched in immune pathways and active in immune cells. Mol. Neurobiol.54, 594–600. 10.1007/s12035-015-9670-8
8
JoehanesR.ZhangX.HuanT.YaoC.YingS. X.NguyenQ. T.et al. (2017). Integrated genome-wide analysis of expression quantitative trait loci aids interpretation of genomic association studies. Genome Biol.18:16. 10.1186/s13059-016-1142-6
9
JunG. R.ChungJ.MezJ.BarberR.BeechamG. W.BennettD. A.et al. (2017). Transethnic genome-wide scan identifies novel Alzheimer's disease loci. Alzheimers Dement.13, 727–738. 10.1016/j.jalz.2016.12.012
10
LambertJ. C.Ibrahim-VerbaasC. A.HaroldD.NajA. C.SimsR.BellenguezC.et al. (2013). Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer's disease. Nat. Genet.45, 1452–1458. 10.1038/ng.2802
11
LiY.SongD.JiangY.WangJ.FengR.ZhangL.et al. (2016). CR1 rs3818361 polymorphism contributes to Alzheimer's disease susceptibility in Chinese population. Mol. Neurobiol.53, 4054–4059. 10.1007/s12035-015-9343-7
12
LiuG.BaoX.WangR. (2015). Expression quantitative trait loci regulate HNF4A and PTBP1 expression in human brains. Proc. Natl. Acad. Sci. U.S.A.112:E3975. 10.1073/pnas.1509048112
13
LiuG.HuY.JinS.JiangQ. (2017a). Genetic variant rs763361 regulates multiple sclerosis CD226 gene expression. Proc. Natl. Acad. Sci. U.S.A.114, E906–E907. 10.1073/pnas.1618520114
14
LiuG.JiangY.WangP.FengR.JiangN.ChenX.et al. (2012). Cell adhesion molecules contribute to Alzheimer's disease: multiple pathway analyses of two genome-wide association studies. J. Neurochem.120, 190–198. 10.1111/j.1471-4159.2011.07547.x
15
LiuG.LiuY.JiangQ.JiangY.FengR.ZhangL.et al. (2016). Convergent genetic and expression datasets highlight TREM2 in Parkinson's disease susceptibility. Mol. Neurobiol.53, 4931–4938. 10.1007/s12035-015-9416-7
16
LiuG.SunJ. Y.XuM.YangX. Y.SunB. L. (2017b). SORL1 variants show different association with early-onset and late-onset Alzheimer's disease risk. J. Alzheimers. Dis.58, 1121–1128. 10.3233/JAD-170005
17
LiuG.WangT.TianR.HuY.HanZ.WangP.et al. (2018a). Alzheimer's disease risk variant rs2373115 regulates GAB2 and NARS2 expression in human brain tissues. J. Mol. Neurosci.66, 37–43. 10.1007/s12031-018-1144-9
18
LiuG.XuY.JiangY.ZhangL.FengR.JiangQ. (2017c). PICALM rs3851179 variant confers susceptibility to Alzheimer's disease in Chinese population. Mol. Neurobiol.54, 3131–3136. 10.1007/s12035-016-9886-2
19
LiuG.ZhangF.HuY.JiangY.GongZ.LiuS.et al. (2017d). Genetic variants and multiple sclerosis risk gene SLC9A9 expression in distinct human brain regions. Mol. Neurobiol.54, 6820–6826. 10.1007/s12035-016-0208-5
20
LiuG.ZhangF.JiangY.HuY.GongZ.LiuS.et al. (2017e). Integrating genome-wide association studies and gene expression data highlights dysregulated multiple sclerosis risk pathways. Mult. Scler.23, 205–212. 10.1177/1352458516649038
21
LiuG.ZhangY.WangL.XuJ.ChenX.BaoY.et al. (2018b). Alzheimer's disease rs11767557 variant regulates EPHA1 gene expression specifically in human whole blood. J. Alzheimers Dis.61, 1077–1088. 10.3233/JAD-170468
22
LiuG.ZhaoY.JinS.HuY.WangT.TianR.et al. (2018c). Circulating vitamin E levels and Alzheimer's disease: a Mendelian randomization study. Neurobiol Aging72, 189.e1–189.e9. 10.1016/j.neurobiolaging.2018.08.008
23
LiuG.ZhaoY.SunJ. Y.SunB. L. (2018d). Parkinson's disease risk variant rs1109303 regulates the expression of INPP5K and CRK in human brain. Neurosci. Bull. [Epub ahead of print]. 10.1007/s12264-018-0289-4
24
NallsM. A.PankratzN.LillC. M.DoC. B.HernandezD. G.SaadM.et al. (2014). Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson's disease. Nat. Genet.46, 989–993. 10.1038/ng.3043
25
OngC. T.CorcesV. G. (2011). Enhancer function: new insights into the regulation of tissue-specific gene expression. Nat. Rev. Genet.12, 283–293. 10.1038/nrg2957
26
RyanK. J.WhiteC. C.PatelK.XuJ.OlahM.ReplogleJ. M.et al. (2017). A human microglia-like cellular model for assessing the effects of neurodegenerative disease gene variants. Sci. Transl. Med.9:eaai7635. 10.1126/scitranslmed.aai7635
27
ShlyuevaD.StampfelG.StarkA. (2014). Transcriptional enhancers: from properties to genome-wide predictions. Nat. Rev. Genet.15, 272–286. 10.1038/nrg3682
28
WardL. D.KellisM. (2012). HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res.40, D930–D934. 10.1093/nar/gkr917
29
WardL. D.KellisM. (2016). HaploReg v4: systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res.44, D877–D881. 10.1093/nar/gkv1340
30
ZhangS.ZhangD.JiangY.WuL.ShangH.LiuJ.et al. (2015). CLU rs2279590 polymorphism contributes to Alzheimer's disease susceptibility in Caucasian and Asian populations. J. Neural. Transm.122, 433–439. 10.1007/s00702-014-1260-9
31
ZhangY.WangL.JiaH.LiaoM.ChenX.XuJ.et al. (2018). Genetic variants regulate NR1H3 expression and contribute to multiple sclerosis risk. J. Neurol. Sci.390, 162–165. 10.1016/j.jns.2018.04.037
32
ZhernakovaD. V.DeelenP.VermaatM.Van ItersonM.Van GalenM.ArindrartoW.et al. (2017). Identification of context-dependent expression quantitative trait loci in whole blood. Nat. Genet.49, 139–145. 10.1038/ng.3737
Summary
Keywords
genome-wide association study, neurodegenerative disease, eQTLs, CD14+ monocytes, Alzheimer's disease
Citation
Sun J, Hou Y, Zhang Y, Wang L, Liu L, Sun B and Yuan H (2018) Genetic Variants Associated With Neurodegenerative Diseases Regulate Gene Expression in Immune Cell CD14+ Monocytes. Front. Genet. 9:666. doi: 10.3389/fgene.2018.00666
Received
06 October 2018
Accepted
04 December 2018
Published
18 December 2018
Volume
9 - 2018
Edited by
Yan Huang, Harvard Medical School, United States
Reviewed by
Feng ZHU, Zhejiang University, China; Weiwei Xue, Chongqing University, China
Updates
Copyright
© 2018 Sun, Hou, Zhang, Wang, Liu, Sun and Yuan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bao-liang Sun blsun88@163.comHui Yuan yuanhui0314@163.com
This article was submitted to Neurogenomics, a section of the journal Frontiers in Genetics
†These authors share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.