- Open Access
A functional SNP associated with atopic dermatitis controls cell type-specific methylation of the VSTM1 gene locus
- Dilip Kumar†1,
- Kia Joo Puan†1,
- Anand Kumar Andiappan†1,
- Bernett Lee1,
- Geertje H. A. Westerlaken2,
- Doreen Haase1,
- Rossella Melchiotti1,
- Zhuang Li1,
- Nurhashikin Yusof1,
- Josephine Lum1,
- Geraldine Koh1,
- Shihui Foo1,
- Joe Yeong1, 16,
- Alexessander Couto Alves3,
- Juha Pekkanen4,
- Liang Dan Sun5,
- Astrid Irwanto6,
- Benjamin P. Fairfax7, 8,
- Vivek Naranbhai7, 8,
- John E. A. Common9,
- Mark Tang12,
- Chin Keh Chuang17, 18,
- Marjo-Riitta Jarvelin3, 13, 14, 15,
- Julian C. Knight7,
- Xuejun Zhang5,
- Fook Tim Chew19,
- Shyam Prabhakar6,
- Liu Jianjun6,
- De Yun Wang10, 11,
- Francesca Zolezzi1,
- Michael Poidinger1,
- E. Birgitte Lane9,
- Linde Meyaard†2Email author and
- Olaf Rötzschke†1Email author
© The Author(s). 2017
- Received: 10 August 2016
- Accepted: 11 January 2017
- Published: 20 February 2017
Expression quantitative trait loci (eQTL) databases represent a valuable resource to link disease-associated SNPs to specific candidate genes whose gene expression is significantly modulated by the SNP under investigation. We previously identified signal inhibitory receptor on leukocytes-1 (SIRL-1) as a powerful regulator of human innate immune cell function. While it is constitutively high expressed on neutrophils, on monocytes the SIRL-1 surface expression varies strongly between individuals. The underlying mechanism of regulation, its genetic control as well as potential clinical implications had not been explored yet.
Whole blood eQTL data of a Chinese cohort was used to identify SNPs regulating the expression of VSTM1, the gene encoding SIRL-1. The genotype effect was validated by flow cytometry (cell surface expression), correlated with electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP) and bisulfite sequencing (C-methylation) and its functional impact studied the inhibition of reactive oxygen species (ROS).
We found a significant association of a single CpG-SNP, rs612529T/C, located in the promoter of VSTM1. Through flow cytometry analysis we confirmed that primarily in the monocytes the protein level of SIRL-1 is strongly associated with genotype of this SNP. In monocytes, the T allele of this SNP facilitates binding of the transcription factors YY1 and PU.1, of which the latter has been recently shown to act as docking site for modifiers of DNA methylation. In line with this notion rs612529T associates with a complete demethylation of the VSTM1 promoter correlating with the allele-specific upregulation of SIRL-1 expression. In monocytes, this upregulation strongly impacts the IgA-induced production of ROS by these cells. Through targeted association analysis we found a significant Meta P value of 1.14 × 10–6 for rs612529 for association to atopic dermatitis (AD).
Low expression of SIRL-1 on monocytes is associated with an increased risk for the manifestation of an inflammatory skin disease. It thus underlines the role of both the cell subset and this inhibitory immune receptor in maintaining immune homeostasis in the skin. Notably, the genetic regulation is achieved by a single CpG-SNP, which controls the overall methylation state of the promoter gene segment.
- Signal inhibitory receptor on leukocytes-1 (SIRL-1)
- Expression quantitative trait loci (eQTL)
- Atopic dermatitis
- Reactive oxygen species (ROS)
Genome-wide association studies (GWAS) in the recent past have identified thousands of disease-related genetic variants and trait-associated polymorphisms [1–3]. The National Genome Research Institute (NHGRI) GWAS catalog alone comprises data of more than 2000 different studies organized with association P values from hundreds of disease phenotypes and complex traits . However, translating these associated variants to causality and potential mechanisms of disease have been difficult. Recent advances in high throughput genotyping and gene expression analysis have allowed for an extensive mapping of expressed quantitative trait loci (eQTL) [5, 6]. Although many of these eQTL are cell-type–specific and also vary between the ethnic groups, they represent a valuable resource for candidate proteins, whose messenger RNA (mRNA) expression is, at least in part, genetically controlled [5–7]. These eQTL datasets provide an excellent opportunity to address and understand the mechanisms underlying GWAS hits. Zhu et al. have recently demonstrated this at a whole genome level by providing a list of priority candidate genes for further functional evaluation of the five complex human traits . However, once the priority candidate genes have been identified, it needs to be validated using well designed in vitro, in vivo, or ex vivo models using techniques, which go beyond assaying the genome and transcriptome [9, 10]. Additionally, this also requires knowledge on the protein molecule under investigation or pathways identified, cell/tissue specificity, and functional consequence on the immune system. Hence, there is a need for targeted but yet comprehensive candidate gene-based studies looking at various aspects of the potential disease gene from DNA to protein, also in pure populations of immune cells and/or tissues with appropriate stimuli and functional readouts. In this study, we present an in-depth, well-characterized study of the SIRL-1 (signal inhibitory receptor on leukocytes-1) gene using various complementary techniques.
The “immune inhibitory receptors” represent a heterogeneous family of transmembrane receptors that function as negative regulators of immune cells [11, 12]. Upon ligation, the inhibitory signal is transduced through immunoreceptor tyrosine-based inhibitory motifs (ITIMs) located in the cytoplasmic tail of the receptor . One member of this family is SIRL-1. While the natural ligand is not known yet, cross-linking of the receptor with agonistic antibodies results in the tyrosine phosphorylation of the ITIMs by the Src family kinases (SFK). The phosphorylated ITIMs act as a docking site for the recruitment of the Src homology 2 (SH2) domain-containing tyrosine phosphatases SHP-1 and/or SHP-2. These phosphatases antagonize the effects of activating SFK and Syk family kinases, triggered for instance by the engagement of IgE-specific Fc-epsilon receptors . Therefore, these immune inhibitory receptors have been implicated to play a role in regulating antibody-mediated immune responses. SIRL-1 is known to primarily express on human myeloid cells, where it inhibits crucial pro-inflammatory functions, such as Fc receptor (FcR) mediated respiratory burst and neutrophil extracellular trap (NET) formation [12, 14, 16, 17]. In a prior study, we made an observation reporting striking differences in the expression level of SIRL-1 on various immune subsets. Granulocyte subsets neutrophils and eosinophils constitutively expressed high levels SIRL-1 on most individuals, while the expression level of the receptor on monocytes varied strongly between individuals . The cause of these variations as well as the potential clinical implications was not known.
In this study, we used mRNA and genotype data from two independent cohorts (Chinese, Caucasian) to link population variations of the molecule on monocytes to a single single nucleotide polymorphism (SNP), rs612529T/C, a CpG-SNP located in the promoter region of the SIRL-1 encoding gene VSTM1. By employing flow cytometry (FACS), functional assays, electrophoretic mobility shift assay (EMSA), bisulfite DNA sequencing, and chromatin immune precipitation (ChIP) analysis, we could further show that SIRL-1 upregulation correlates with the allele dependent binding of PU.1 and YY1 to the CpG-SNP, which in turn is associated with a demethylation of the entire promoter region. Importantly, genetic association analysis suggested that the cell-type–specific fluctuations of SIRL-1 indeed might have pathologic consequences. We observed that in three independent cohorts, the C allele (indicative of low expression) is consistently associated with an increased risk for atopic dermatitis (AD), an inflammatory skin condition.
Human blood samples and cohorts
The cohorts used in this study were outlined in Additional file 1. eQTL analysis was carried out with ethnic Chinese cohorts collected in Singapore (“blood eQTL cohort” and “monocyte & B cell eQTL cohort”) of 202 (unpublished) and 15 individuals , respectively. In addition, data were used from a “Monocyte & B cell cohort” of 279 individuals of Caucasian ethnicity that had been reported before by Fairfax et al.  With the exception of the “functional cohort” (see below), details on all cohorts had been published previously. Data of three additional cohorts (Singapore Chinese population , Chinese Han Population  and Caucasian population ) that had been published previously were used for the genetic association of the VSTM1-SNP rs612529 with AD.
Another genotyped Singaporean cohort consisting of ethnic Chinese individuals was used for the functional characterization of VSTM1 by FACS (n = 44), reactive oxygen species (ROS) assay (n = 30), DNA-methylation (n = 8), and ChIP assay (n = 9) (“functional cohort”). For all individuals of the “functional cohort,” the rs612529T/C genotype was determined by using high resolution melting (HRM) analysis. The functional cohort was used for the rs612529 genotype/phenotype association, measurement of NADPH oxidase mediated ROS production assay, ChIP assay and bisulfite treatment, and methylation analysis as described below.
DNA was isolated from whole blood or peripheral blood mononuclear cells (PBMCs) using DNeasy Blood & Tissue Kits (Qiagen) according to the manufacturer’s instructions. For targeted genotyping, the HRM analysis was performed by real-time polymerase chain reaction (PCR) using a CFX96 Real-Time Detection System (Bio-Rad). The primers used for rs612529 genotyping were 5′-TCTTGGCAGAACTTCAGATAAGGT-3′ (forward) and 5′-ACAAGAAGCCGTCGATGATAACT-3′ (reverse). Amplification was carried out using the following protocol: 3 min at 95 °C, 40 cycles of 5 s at 95 °C, 5 s at 50 °C, and final extension for 10 s at 95 °C. A melting curve was generated in the range of 65–95 °C (in 0.2 °C increments) with 10 s/step. Heterozygosity and homozygosity of the allelic state were deduced using Precision Melting Analysis software (Bio-Rad).
RNA extraction and global gene expression for whole blood samples
Total RNA was extracted from human whole blood collected in Tempus Blood RNA tubes (ThermoFisher) using the MagMAX™ for Stabilized Blood Tubes RNA Isolation Kit (ThermoFisher). Total RNA was subjected to globin mRNA depletion using the Ambion GLOBINclear kit (ThermoFisher). Agilent Bioanalyzer was used to evaluate total RNA integrity and the RNA Integrity Number (RIN) was calculated. All samples had RIN in the range of 4.7–8.8; however, the median value was 7. For genome-wide gene expression profiling, fluorescent-labeled PCR products were prepared according to Illumina Human Whole-Genome Gene Expression DASL Assay Guide. The labeled products were hybridized onto the Illumina HumanHT-12-v4 Expression Bead ChIP for 16 h at 58 °C. These arrays were washed, coated based on the Assay Guide, and scanned using Bead Array Scanner 500GX at BSF Microarray Facility.
The whole blood cohort used for eQTL analysis has not been published yet. The median age is 21 years and 55% male. Whole blood gene expression data were processed using the Bioconductor lumi package in R 3.1.2. The gene expression data were quantile normalized and log2 transformed prior to further analysis. Whole blood gene expression data for genes located 200kbp upstream and downstream of the SNP rs612529 were extracted and combined with the SNPs located within this same region. One-way ANOVAs were then done to determine if any of the SNP genotypes associated with the gene expression levels. The one-way ANOVA was conducted in R 3.1.2 and all resulting P values were corrected for multiple testing using the method of Benjamini and Hochberg. The eQTL data for the B cells and monocytes for the 15 samples for the Chinese samples from Singapore were previously described in Melchiotti et al. . Data for the B cells and monocytes for the Caucasian cohort were extracted from Fairfax et al. .
Association analysis for rs612529
The association analysis for the SNP rs612529 was performed using statistics from the various cohorts using Stouffer z-trend method using the Meta P program (http://igm.cumc.columbia.edu/MetaP/metap.php). The Stouffer’s z trend method considers P values, sample sizes, and effect directions. For more information on the program, refer to Whitlock et al. .
Flow cytometry analysis
Blood samples were collected in either BD K2EDTA or citrate vacutainer tubes. FACS staining was performed on either using whole blood sample or isolated PBMCs. Erythrocytes were removed from whole blood by incubation in a 14-mL erythrocyte lysis buffer (155 mM NH4Cl, 10 mM KHCO3, 0.1 mM EDTA) for 10 min at room temperature. After erythrocytes lysis, the cells were centrifuged at 860 × g for 3 min. The cell pellets were re-suspended in phosphate-buffered saline (PBS) and washed once with PBS by centrifugation at 860 × g for 3 min at room temperature. For PBMC isolation, whole blood was layered on Ficoll-Paque (GE Healthcare) and centrifuged for 500 × g for 30 min without a break. PBMCs were harvested at the interface between Ficoll and plasma layers. To discriminate the cells live from the dead cells, lyzed whole blood samples or PBMC cell pellets were incubated with 100 μL of LIVE/DEAD Fixable Aqua Dead kit (Life Technologies) in PBS for 10 min at room temperature. The cells were washed once with 100 μL MACS buffer (0.5% BSA, 2 mM EDTA in PBS) and transferred into 96 V-bottom plates for centrifugation at 1000 × g for 3 min. For the SIRL-1 phenotyping, the cell pellets were incubated with for 15 min at 4 °C anti-SIRL-1 FITC (1A5; as previously described ), anti-FcεRI PE (AER37; eBioscience), anti-CD123 PerCP Cy5.5 (6H6; eBioscience), anti-CD1c BV421 (L161; Biolegend), anti-CD14 PE-CF594 (MΦP9; BD Biosciences), anti-CCR3 Alexa Fluor647 (5E8; Biolegend) and anti-HLA-DR APC-H7 (L243; BD Biosciences). For the correlation of SIRL-1 with Fc receptor expression, the staining was carried out with anti-FcεR1 (AER37), anti-CD16 (3G8; Biolegend), and anti-CD89 (A59; BD Bioscience) monoclonal antibodies. After incubation, the cell pellets were washed once with MACS buffer, re-suspended in the 250 μL of the same buffer, and SIRL-1 expression level on different cell types was determined using a LSRII or a Fortessa flow cytometer (BD Biosciences).
Monocytes used for nuclear protein extraction were isolated from apheresis blood and for ChIP and bisulfite sequencing monocytes were isolated from PBMCs derived from functional cohort by positive selection using MACS Human CD14 Microbeads (Miltenyi Biotec) according to the instructions recommended by manufacturer. Neutrophils used for the bisulfite sequencing were isolated by Ficoll-Paque density gradient centrifugation. The neutrophils were isolated from the layer immediately above erythrocyte layer and lyzing the erythrocyte with erythrocyte lysis buffer (155 mM NH4Cl, 10 mM KHCO3, 0.1 mM EDTA) for 10 min at room temperature on a roller mixer. The neutrophils were washed once with cold RMPI1640 supplemented with 10% fetal bovine serum.
Measurement of NADPH oxidase mediated ROS production
ROS production was measured in primary monocytes isolated from frozen PBMC samples of genotyped donors of the functional cohort. All experiments were performed using technical duplicates. NADPH oxidase activity was, unless otherwise indicated, assessed as H2O2 formation and determined by Amplex Red (Molecular Probes), which reacts with H2O2 to produce fluorescent resorufin in the presence of HRP (Sigma-Aldrich). The assay was performed in microfluor white plates (Thermo Scientific). Isolated monocytes were re-suspended in HEPES buffer (20 mM HEPES, 132 mM NaCl, 6 mM KCl, 1 mM MgSO4, 1.2 mM KH2PO4, pH 7.4 supplemented with 0.5% human serum albumin (Sigma), 1 mM Ca2+, and 5 mM D-glucose) containing 50 μM Amplex Red and 1 U/mL HRP and stimulated (1 × 106 cells/mL) with 10 μg/mL plate-bound human IgA (ChromPure Human IgA, Jackson ImmunoResearch Laboratories) together with 10 μg/mL plate-bound anti-SIRL-1 (1A5) or an isotype-matched control mAb (functional grade mouse IgG clone:P3, eBioscience). Formation of H2O2 was measured directly from the start of stimulation, when cells were added to coated wells, and continued for 1 h. Fluorescence was measured at 1-min intervals (λEx/λEm = 545/590 nm) using a Tecan Infinite M200 PRO fluorescence reader (Tecan). Measurements were corrected for spontaneous H2O2 production, as determined by incubation of cells in non-coated wells. The percentage of SIRL-1-mediated inhibition of ROS production was calculated by sub-traction of background H2O2 production and subsequent calculation of the area under the curve of all samples. The AUC of IgA co-coated with an isotype-matched control mAb was set at 100% and compared with the AUC of IgA co-coated with anti-SIRL-1 mAb.
Nuclear protein isolation
Freshly isolated monocytes from leukapheresis blood were washed with ice cold PBS, then buffer A (20 mM HEPES, pH 7.9, 20% Glycerol, 10 mM NaCl, 0.2 mM EDTA (pH 8), 1 mM DTT, 0.1% Triton X-100) supplemented with Pierce Proteases and Phosphatase Inhibitor Mini Tablets. The tablet contained aprotinin, bestatin, E-64, EDTA, leupeptin, sodium fluoride, sodium orthovanadate, sodium pyrophosphate, beta-glycerophosphate, and metalloproteinase inhibitor. After 15 min incubation on ice, homogenates were centrifuged at 2000 rpm at 4 °C for 15 min, and the resultant “nuclear pellets” were re-suspended in buffer A containing 500 mM NaCl. The nuclear proteins were incubated for 1 h on ice with intermittent tapping followed by centrifugation at 13,000 rpm (4 °C) for 15 min. The supernatants were then aliquoted and snap-frozen at –80 °C until used. Quantitation of nuclear proteins was performed using Bio-Rad protein quantitation kit (Bio-Rad).
Electrophoretic mobility shift assay
The double-strand oligonucleotides had a length of 25 bp (the sequences are listed in Additional file 2). Allele-specific oligonucleotides corresponding to VSTM-1 promoter region –34 to –58 were designed with the C/T polymorphism of rs612529 in the center (termed as “C probe” and “T probe”). C and T probes, as well as MZF-1 consensus/mutant oligonucleotides, were synthesized by 1st BASE, Singapore; all remaining oligonucleotides were purchased from Santa Cruz Biotechnology. 32P labeling was introduced using T4 Polynucleotide Kinase. For EMSA reactions, 5 μg of the monocyte nuclear protein extract was incubated with 40,000 cpm of the 32P-labeled oligonucleotide probe in a total volume of 20 μL 5 mM HEPES (pH 7.9), 5% glycerol, 0.5 mM EDTA, 1 μg dIdC, 1 mM DTT, and 80 mM NaCl at 4 °C for 45 min. For “cold” competition experiments, 12.5–50 ng of the unlabeled oligonucleotide competitor 15 min were added before adding the respective radiolabeled probe, representing a 12.5–50-fold excess of the unlabeled competitor. For supershift assays, 2 uL of either rabbit polyclonal anti-PU.1 IgG (H-135 X, Santa Cruz Biotechnology) or rabbit polyclonal anti-YY1 IgG (C-20, Santa Cruz Biotechnology) or rabbit IgG (Abcam) were added 15 min before adding the respective radiolabeled probe. After incubation, the DNA–protein complexes were separated for 90 min with a 5% non-denaturing polyacrylamide gel at 240 V on a protein II gel-apparatus (Bio-Rad) using 0.5X Tris/Borate/EDTA (TBE). After electrophoresis, the gel was transferred to a Whatman paper and dried on a vacuum gel dryer (Gel Dryer 583, Bio-Rad) for 2 h at 80 °C. The dried gel was then exposed on an imaging plate and scanned using a FLA-5000 phospho-image scanner (Fujifilm Life Science).
ChIP assays were performed on monocytes isolated from four CC-donors and five TT-donors (rs612529) of the functional cohort. ChIP was carried out essentially as previously described . Briefly, DNA/protein cross-linking was achieved by incubating the cells for 8 min at 37 °C in 1% formaldehyde. Ultrasound sonication of the chromatin was performed in lysis buffer (50 mM Tris-Cl pH 8.0, 10 mM EDTA, 1% SDS) with 3-min cycles (30 s “ON”, 30 s “OFF”) for a total of ten cycles using a Bioruptor UCD-300 sonicator (Diagenode). After sonication, samples were diluted in fivefold dilution buffer (0.01% SDS,1.1% Triton X-100, 1.1 mM EDTA, 20 mM Tris-Cl pH 8.0) followed by overnight incubation at 4 °C with 10 μg of either rabbit polyclonal anti-PU.1 antibody IgG (H-135 X, Santa Cruz Biotechnology), or rabbit polyclonal IgG (Abcam). Real-time PCRs of genomic regions containing the putative PU.1-binding site were performed in triplicate by using iTaq SYBR green supermix (Bio-Rad) with SIRL-1 specific primers (forward primer: 5′-TCTTGGCAGAACTTCAGATAAGGT-3′; reverse primer: 5′-ACAAGAAGCCGTCGATGATAACT-3′). The relative occupancy of the immunoprecipitated factor at SIRL-1 locus is estimated by using the comparative threshold method 2ˆ(Ctmock–Ctspecific), where Ctmock and Ctspecific are mean threshold cycles of PCR done in triplicate on DNA samples from mock and specific immunoprecipitation .
Bisulfite treatment and DNA methylation analysis
Primers for the bisulfite genomic sequencing (BGS) were designed to cover a region 361 bp upstream and 335 bp downstream to the transcription start site of VSTM1 (Additional file 3). DNA bisulfite treatment of the genomic DNA isolated from freshly isolated monocytes and neutrophils of respectively four donors of the CC and TT genotype (functional cohort) was performed according to the protocol described in EZ DNA Methylation-Direct™ Kit (Zymo Research). For BGS, bisulfite-treated DNA was amplified using the SIRL-1 promoter specific primers
F2: 5′- ACTCTAACACTCTAATTCTCTACCCCAC-3′,;
B1:5′ AGAGTGGGGTAGAGAATTAGAGTGTTA and
B2: 5′- ACCTAAACTACCACTCCCACATAAAT-3′) with optimized conditions. The PCR products were cloned into the TA cloning vector (Invitrogen, Life technologies). For both monocytes and neutrophils three clones from each donor and sequenced using M13 forward primer.
CpG islands (CGI) prediction
A Hidden Markov Model (HMM) based CGI prediction for the VSTM1 gene locus was obtained from UCSC genome browser as custom tracks (http://genome.ucsc.edu/cgi-bin/hgTracks) . CpG islands (CGI) from DNA sequences and custom tracks were uploaded from Rafael Irizarry lab (http://rafalab.jhsph.edu/CGI/).
VSTM1 promoter-associated SNP rs612529 is an eQTL for SIRL-1
Cell-type–specific regulation of SIRL-1 protein expression in monocytes
Impact of rs612529 on the inhibition of ROS
rs612529 regulates the recruitment of PU.1 and YY1 to the VSTM1 promoter
In silico analysis by JASPAR (http://jaspar.binf.ku.dk/) and TRANSFAC (http://www.gene-regulation.com/) predicted binding by PU.1, MZF-1, and YY1 (Additional file 6). EMSA competition of the radiolabeled T probe with unlabeled oligonucleotides representing the consensus binding sequences of these transcription factors suggested band 1 to be formed by PU.1, while band 2 derived from YY1 (Fig. 4b; Additional files 2 and 7). The identity of the two factors was formally confirmed by super-shift assays, in which the PU.1-specific and YY1-specific antibodies super-shifted and/or blocked the complex formation with the T probe (Fig. 4c).
The allele-specific interaction of PU.1 with the VSTM1 promoter region in monocyte-derived nuclear samples was also confirmed in ChIP assays (Fig. 4d). DNA fragments precipitated from four CC and five TT donors with PU.1-specific antibodies were analyzed by quantitative PCR with a primer pair covering a region of the VSTM1 promoter including the rs612529 SNP (Additional file 3). In line with the results from the EMSA experiments, a significant increase in the amount of VSTM1 DNA was detected in the PU.1 precipitate from TT donors compared to CC donors (Fig. 4d). Thus, in monocytes PU.1 is strongly associated with the analyzed VSTM1 promoter segment when rs612529 consists of the T allele.
The allelic state of rs612529 is correlated with the methylation of VSTM1 in monocytes
In order to study a potential role of the rs612529 CpG-SNP in the control of the monocyte-specific SIRL-1 expression, we analyzed cytosine methylation by bisulfite sequencing. A total of eight donors (four CC genotype and four TT genotype) were used to determine the influence of the allelic state in both eQTL relevant monocytes and allele-independent neutrophils (Fig. 5b and c). While the CpG-SNP itself would appear to be a potential site for an allele-specific methylation, it was apparently not a target for modification, as it was found to be located in a 54-bp core region containing six additional CpG pairs, which were constitutively demethylated in all samples analyzed. While the CpG-SNP was not directly affected by epigenetic modifications, it had a strong and cell-type–specific effect on the general methylation state in the analyzed gene region. In line with the bright but allele-independent surface expression detected by FACS (compare Fig. 2) neutrophils exhibit demethylated CpG elements both in the CC and the TT donors (Fig. 5b). In monocytes, however, only TT donors had a demethylated VSTM1 promoter, while the gene segment was almost completely methylated in CC donors (Fig. 5c). Thus, in monocytes the formation of a CpG pair by the C allele of rs612529 seems to promote the methylation of other CpG pairs in the region results in a silencing of the locus (compare Fig. 2).
Genetic association of rs612529 to AD phenotype
Association of rs612529 with AD
In this study, we significantly associate the previously reported cell-type–specific variations in SIRL-1 expression to a genetic cause. A single SNP located in the promoter region mediates the epigenetic silencing of the entire gene locus, a process presumably mediated by the allele-specific recruitment of PU.1/YY1. The allelic effect is particularly evident in monocytes, where the C allele of rs612529 is associated with a near complete lack of SIRL-1 expression. This in turn results in a diminished control of pro-inflammatory reactions, evidenced by the failure of SIRL-1 agonists to dampen the production of Fc-receptor induced ROS in these cells. Genetic association results in independent cohorts suggest that this is associated consistently with an increased risk of AD.
We then demonstrated allele-specific binding of PU.1 to the promoter region of VSTM1 using follow-up functional experiments such EMSA and ChIP assays. This transcription factor PU.1 has been shown to be a master regulator of myeloid cell development, which could explain the observed specificity of the SIRL-1 expression in cells of this lineage [31–33]. PU.1 forms complexes with various regulatory factors chromatin remodelers, including DNA methyltransferases and demethylases to activate or repress transcription . A recent study has shown that in non-differentiated monocytes, PU.1 can act as docking site for the demethylase, Tet2 . In these cells the recruitment of Tet2 by PU.1 leads to the targeted demethylation of the respective gene locus. In line with this observation we therefore hypothesize that the allele-dependent methylation of the VSTM1 gene reported here can be explained by the conditional recruitment of Tet2 or other demethylases as a result of the allele-specific binding of PU.1 to rs612529 (Additional file 8). In monocytes rs612529 would thus represent a genetic master switch for the epigenetic control of SIRL-1 expression by CpG demethylation.
Mechanistically, the methylation of CpG promotes the compaction of chromatin by providing interaction sites for methyl-CpG-binding proteins (MBPs) [24, 36, 37]. Hence, CpG-SNPs, which allow methylation in an allele-dependent way, are often associated with function [38–40]. In this context, it is interesting to note that rs612529 is a CpG-SNP where the C allele forming the CpG pair is associated with the epigenetic silencing of the VSTM1 gene. However, in contrast to conventional CpG-SNPs, it is not a direct target for methylation. Instead of acting as a simple “switch” to close the chromatin, it rather functions as a genetic/epigenetic regulator that modulates the DNA methylation state by the recruitment of CpG modifying enzymes. In this way, the effect of the CpG-SNP becomes potentiated, as the methylation affects not just a single CpG but also many other CpG sites in its vicinity. Subtle changes associated with the allelic switching between C and T therefore translate into extensive epigenetic alterations of the gene, which makes rs612529 a master genetic regulator of the VSTM1 gene in monocytes.
We then evaluated if the causal SNP we identified, rs612529, has potential pathological implications, given the striking variations in SIRL-1 surface expression observed on monocytes. By using the allelic states of the SNP as proxy for the expression levels in these cells we could establish AD as being associated with a low SIRL-1 expression in these cells. AD is a chronic inflammatory skin disease affecting 60 million people in Europe and the USA . The pathophysiology of AD is a result from a complex interplay of host genetic factors, environmental triggers, skin integrity, immune balance, and the skin microbiome . Thus, a tight control of the immune response is needed to ensure effective pathogen defense without collateral damage to the protective skin layers. IgE is one of the determinants predisposing factor for AD. However, earlier reports have shown that in early childhood, elevated levels of IgE do not seem to correlate with disease manifestations or severity of AD . Moreover, some evidence pointing towards a possible role of IgA derives from an earlier report indicating that patients with selective IgA deficiency have an increased incidence of atopic dermatitis .
SIRL-1 is an effective inhibitor of monocyte-derived pro-inflammatory reactions such as FcR-induced ROS production. In our study, we show that blood monocytes represent a heterogeneous population characterized by different Fc receptors. In Fig. 3a and Additional file 5, FcαR+ monocytes co-segregate with a high expression of SIRl-1 expression. This is not observed in FcγRIII+ and FcεR+ monocytes, the latter being expressed on a very small subset of monocytes. As such, we examined if co-activating of FcαR and SIRL-1 with human IgA and anti-SIRL-1 mAb would influence the capacity of monocytes to generate ROS. We show the capacity of SIRL-1 to modulate ROS production in monocytes is a function of VSTM1 genotype. The association of the low expressing C allele with AD thus leads us to hypothesize that inhibitory regulation by SIRL-1 of monocytes or monocyte-derived cells plays a key role in skin homeostasis.
In summary, our study introduces with rs612529 an important genetic factor influencing the risk for skin inflammation. The identification of this link was made possible by defining an eQTL associated with these fluctuations from the whole blood data of our Chinese cohort and using the genotypes to mine pre-existing genetic association databases. A growing number of eQTL and GWAS data have been published already but only a fraction of the GWAS datasets are publicly accessible so far. Given the power the approach has to uncover new and unexpected links between genetically controlled phenotypes and disease manifestation, it appears to be mandatory that more of the “large data,” including eQTL and GWAS, has to be compiled and curated in easily accessible repositories which are open for the general scientific community. In the case of SIRL-1 it not only helped to understand the role of the receptor in controlling immune cells, it also identified a novel candidate target for the treatment of AD. The identification of the as-yet unknown ligand(s) of SIRL-1 may thus open up new avenues to control or suppress the manifestation of the inflammatory skin disease.
The authors thank all apheresis blood donors from the Health Sciences Authority of Singapore for consenting to give blood for this study. We thank the volunteers and family members for participating in this study, as well as members of the Functional Genomics Laboratory at Department of Biological Sciences, National University of Singapore for their help with sample collection. We also thank Ivy Low, Nurhidaya Binte Shadan, and Seri Munriah Mustafah from the Flow Cytometry Facility. Esther Mok and Joni Chong from the Immunomonitoring platform at the Singapore Immunology Network for helping with the sample processing.
The study was supported by the Singapore Ministry of Education Academic Research Fund (R-154-000-404-112, R-154-000-553-112, R-154-000-565-112, R-154-000-630-112), National Medical Research Council (NMRC; Singapore, NMRC/1150/2008), and Singapore Immunology Network, Agency for Science, Technology and Research (A*STAR) Singapore (SIgN-06-006, SIgN-08-020, and SIgN-10-029).
Availability of data and materials
For the data relevant for the Caucasian eQTL study, the gene expression data are already available at www.ebi.ac.uk (E-MTAB-945), while the genotyping data are accessible at European Genome-Phenome Archive (EGA; EGAS00000000109) and available on request.
DK, KJP, and AKA performed conception and design; DK, KJP, GHAW, DH, ZL, and NY carried out the functional experiments while AKA, BL, RM, MP, DK, and KJP carried out analysis and interpretation; JL, GK, and FS, were involved in large-scale data generation while JY, CKC, SP, and DK were involved in the CHIP assay and interpretation of the epigenetics. ACA, JP, LDS, AI, JEAC, EBL, WDY, LJ, MT, MRJ, XZ, FTC, and FZ were involved in the association analysis and BPF, VN, and JK performed eQTL analysis in Caucasian cohorts. LM and OR were the lead investigators who led the study and facilitated the study design and execution. All authors were involved in drafting the manuscript for submission and provided feedback for subsequent revisions. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All Singaporean cohorts are part of the Singapore Chinese Cohort Study, which was approved by the Institutional Review Board at the National University Hospital (IRB nos. NUS 07-023, NUS 10-445, and NUS 09-256) and are in compliance with the Helsinki declaration. Written informed consent was obtained from all volunteers prior to the collection of blood samples. All human apheresis cone blood used for the monocyte isolation was obtained from the Health Sciences Authority of Singapore and was approved by the National University of Singapore (IRB no. NUS-IRB 10-250) and is in compliance with the Helsinki declaration. Written informed consent was obtained for these donors and the samples were anonymized before they were accessible for research.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Visscher PM, Brown MA, McCarthy MI, Yang J. Five years of GWAS discovery. Am J Hum Genet. 2012;90:7–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Manolio TA. Genomewide association studies and assessment of the risk of disease. N Engl J Med. 2010;363:166–76.View ArticlePubMedGoogle Scholar
- Chan SL, Jin S, Loh M, Brunham LR. Progress in understanding the genomic basis for adverse drug reactions: a comprehensive review and focus on the role of ethnicity. Pharmacogenomics. 2015;16:1161–78.View ArticlePubMedGoogle Scholar
- Welter D, MacArthur J, Morales J, Burdett T, Hall P, Junkins H, et al. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014;42:D1001–6.View ArticlePubMedGoogle Scholar
- Lappalainen T. Functional genomics bridges the gap between quantitative genetics and molecular biology. Genome Res. 2015;25:1427–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Pickrell JK, Marioni JC, Pai AA, Degner JF, Engelhardt BE, Nkadori E, et al. Understanding mechanisms underlying human gene expression variation with RNA sequencing. Nature. 2010;464:768–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Nica AC, Dermitzakis ET. Expression quantitative trait loci: present and future. Philos Trans R Soc Lond B Biol Sci. 2013;368:20120362.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhu Z, Zhang F, Hu H, Bakshi A, Robinson MR, Powell JE, et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat Genet. 2016;48:481–7.View ArticlePubMedGoogle Scholar
- Edwards SL, Beesley J, French JD, Dunning AM. Beyond GWASs: illuminating the dark road from association to function. Am J Hum Genet. 2013;93:779–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang Q. Genetic study of complex diseases in the post-GWAS era. J Genet Genomics. 2015;42:87–98.View ArticlePubMedGoogle Scholar
- Steevels TAM, Meyaard L. Immune inhibitory receptors: Essential regulators of phagocyte function. Eur J Immunol. 2011;41:575–87.View ArticlePubMedGoogle Scholar
- van Avondt K, Fritsch-Stork R, Derksen RHWM, Meyaard L. Ligation of signal inhibitory receptor on leukocytes-1 suppresses the release of neutrophil extracellular traps in systemic lupus erythematosus. PLoS ONE. 2012;8, e78459.View ArticleGoogle Scholar
- Vivier E, Daëron M. Immunoreceptor tyrosine-based inhibition motifs. Immunol Today. 1997;18:286–91.View ArticlePubMedGoogle Scholar
- Steevels TAM, Lebbink RJ, Westerlaken GHA, Coffer PJ, Meyaard L. Signal inhibitory receptor on leukocytes-1 is a novel functional inhibitory immune receptor expressed on human phagocytes. J Immunol. 2010;184:4741–8.View ArticlePubMedGoogle Scholar
- Ravetch JV, Lanier LL. Immune inhibitory receptors. Science. 2000;290:84–9.View ArticlePubMedGoogle Scholar
- Steevels TAM, van Avondt K, Westerlaken GHA, Stalpers F, Walk J, Bont L, et al. Signal inhibitory receptor on leukocytes-1 (SIRL-1) negatively regulates the oxidative burst in human phagocytes. Eur J Immunol. 2013;43:1297–308.View ArticlePubMedGoogle Scholar
- van Avondt K, van der Linden M, Naccache PH, Egan DA, Meyaard L. Signal inhibitory receptor on leukocytes-1 limits the formation of neutrophil extracellular traps, but preserves intracellular bacterial killing. J Immunol. 2016;196:3686–94.View ArticlePubMedGoogle Scholar
- Melchiotti R, Puan KJ, Andiappan AK, Poh TY, Starke M, Zhuang L, et al. Genetic analysis of an allergic rhinitis cohort reveals an intercellular epistasis between FAM134B and CD39. BMC Med Genet. 2014;15:73.View ArticlePubMedPubMed CentralGoogle Scholar
- Fairfax BP, Makino S, Radhakrishnan J, Plant K, Leslie S, Dilthey A, et al. Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles. Nat Genet. 2012;44:502–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen H, Common JEA, Haines RL, Balakrishnan A, Brown SJ, Goh CSM, et al. Wide spectrum of filaggrin-null mutations in atopic dermatitis highlights differences between Singaporean Chinese and European populations. Br J Dermatol. 2011;165:106.View ArticlePubMedGoogle Scholar
- Sun L-D, Xiao F-L, Li Y, Zhou W-M, Tang H-Y, Tang X-F, et al. Genome-wide association study identifies two new susceptibility loci for atopic dermatitis in the Chinese Han population. Nat Genet. 2011;43:690–4.View ArticlePubMedGoogle Scholar
- Rantakallio P. Groups at risk in low birth weight infants and perinatal mortality. Acta Paediatr Scand. 1969;193:1+.Google Scholar
- Whitlock MC. Combining probability from independent tests: the weighted Z-method is superior to Fisher’s approach. J Evol Biol. 2005;18:1368–73.View ArticlePubMedGoogle Scholar
- Lunyak VV, Burgess R, Prefontaine GG, Nelson C, Sze S-H, Chenoweth J, et al. Corepressor-dependent silencing of chromosomal regions encoding neuronal genes. Science. 2002;298:1747–52.View ArticlePubMedGoogle Scholar
- Wu H, Caffo B, Jaffee HA, Irizarry RA, Feinberg AP. Redefining CpG islands using hidden Markov models. Biostatistics. 2010;11:499–514.View ArticlePubMedPubMed CentralGoogle Scholar
- Jaenisch R, Bird A. Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003;33(Suppl):245–54.View ArticlePubMedGoogle Scholar
- Li T, Guo X, Wang W, Mo X, Wang P, Han W. V-set and transmembrane domain-containing 1 is silenced in human hematopoietic malignancy cell lines with promoter methylation and has inhibitory effects on cell growth. Mol Med Rep. 2014;11:1344–51.PubMedGoogle Scholar
- Nakai K, Yoneda K, Kubota Y. Oxidative stress in allergic and irritant dermatitis: from basic research to clinical management. Recent Pat Inflamm Allergy Drug Discov. 2012;6:202–9.View ArticlePubMedGoogle Scholar
- Okayama Y. Oxidative stress in allergic and inflammatory skin diseases. Curr Drug Targets Inflamm Allergy. 2005;4:517–9.View ArticlePubMedGoogle Scholar
- Andiappan AK, Wang DY, Anantharaman R, Parate PN, Suri BK, Low HQ, et al. Genome-wide association study for atopy and allergic rhinitis in a Singapore Chinese population. PLoS ONE. 2011;6:e19719.View ArticlePubMedPubMed CentralGoogle Scholar
- Scott EW, Simon MC, Anastasi J, Singh H. Requirement of transcription factor PU.1 in the development of multiple hematopoietic lineages. Science. 1994;265:1573–7.View ArticlePubMedGoogle Scholar
- Vangala RK, Heiss-Neumann MS, Rangatia JS, Singh SM, Schoch C, Tenen DG, et al. The myeloid master regulator transcription factor PU.1 is inactivated by AML1-ETO in t(8;21) myeloid leukemia. Blood. 2003;101:7.View ArticleGoogle Scholar
- Carotta S, Wu L, Nutt SL. Surprising new roles for PU.1 in the adaptive immune response. Immunol Rev. 2010;238:63–75.View ArticlePubMedGoogle Scholar
- van Riel B, Rosenbauer F. Epigenetic control of hematopoiesis: the PU.1 chromatin connection. Biol Chem. 2014;395:1265–74.PubMedGoogle Scholar
- la Rica De L, Rodríguez-Ubreva J, García M, Islam ABMMK, Urquiza JM, Hernando H, et al. PU.1 target genes undergo Tet2-coupled demethylation and DNMT3b-mediated methylation in monocyte-to-osteoclast differentiation. Genome Biol. 2012;14:R99–9.View ArticleGoogle Scholar
- Deaton AM, Bird A. CpG islands and the regulation of transcription. Genes Dev. 2011;25:1010–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Klose RJ, Bird AP. Genomic DNA methylation: the mark and its mediators. Trends Biochem Sci. 2006;31:89–97.View ArticlePubMedGoogle Scholar
- Zhi D, Aslibekyan S, Irvin MR, Claas SA, Borecki IB, Ordovas JM, et al. SNPs located at CpG sites modulate genome-epigenome interaction. Epigenetics. 2013;8:802–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Ziller MJ, Gu H, Müller F, Donaghey J, Tsai LT-Y, Kohlbacher O, et al. Charting a dynamic DNA methylation landscape of the human genome. Nature. 2013;500:477–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Shoemaker R, Deng J, Wang W, Zhang K. Allele-specific methylation is prevalent and is contributed by CpG-SNPs in the human genome. Genome Res. 2010;20:883–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Guttman-Yassky E, Nograles KE, Krueger JG. Contrasting pathogenesis of atopic dermatitis and psoriasis--part I: clinical and pathologic concepts. J Allergy Clin Immunol. 2011;127:1110–8.View ArticlePubMedGoogle Scholar
- Leung DYM. New insights into atopic dermatitis: role of skin barrier and immune dysregulation. Allergol Int. 2013;62:151–61.View ArticlePubMedGoogle Scholar
- Illi S, von Mutius E, Lau S, Nickel R, Grüber C, Niggemann B, et al. The natural course of atopic dermatitis from birth to age 7 years and the association with asthma. J Allergy Clin Immunol. 2004;113:925–31.View ArticlePubMedGoogle Scholar
- Aghamohammadi A, Cheraghi T, Gharagozlou M, Movahedi M, Rezaei N, Yeganeh M, et al. IgA deficiency: correlation between clinical and immunological phenotypes. J Clin Immunol. 2008;29:130–6.View ArticlePubMedGoogle Scholar