Epigenome-wide analysis of sperm cells identifies IL22 as a possible germ line risk locus for psoriatic arthritis

Remy A. Pollock; Laila Zaman; Vinod Chandran; Dafna D. Gladman

doi:10.1371/journal.pone.0212043

Abstract

Psoriasis and its associated inflammatory arthritis, psoriatic arthritis (PsA), have a clear heritable component, but a large proportion of the heritable risk remains unexplained by gene sequence variation. This study aimed to determine if epigenetic factors contribute to the missing heritability in psoriatic disease. DNA methylation profiling was performed on sperm cells from 23 probands with psoriasis without PsA (PsC), 13 PsA probands, and 18 unaffected controls. Differentially methylated CpGs and regions (DMRs) were identified and validated by pyrosequencing. Underlying AluY and copy number variation (CNV) in the HCG26 and IL22 genes, respectively, were assessed by genotyping. Array, subject’s age, age of psoriasis onset, psoriasis severity, and medication usage were found to influence methylation at many genes and were included as covariates in the analysis. Between PsC probands vs. controls, 169 DMRs were found; 754 DMRs were found between PsA probands vs. controls, and 86 between PsA and PsC probands (adjusted p<0.05). Differences in methylation across DMRs were generally subtle (<10%) but correlated well with pyrosequencing. Biological inference prioritized notable DMRs associated with skin disease (SIGLEC14, JAM3, PCOLCE, RXRB), skin and/or joint disease (MBP, OSBPL5, SNORD115, HCG26), and joint disease (IL22, ELF5, PPP2R2D, PTPRN2, HCG26). Hypermethylation of the DMR within the first exon of arthritis-associated IL22 showed significant correlation (rho = 0.34, 95% CI 0.06–0.57, p = 0.01) between paired sperm and blood samples, independent of a CNV within the same region. Further studies are needed to rule out underlying genetic causes and determine if these represent heritable, constitutional epimutations, or are the result of exposure of germ cells to endogenous or exogenous environmental factors.

Citation: Pollock RA, Zaman L, Chandran V, Gladman DD (2019) Epigenome-wide analysis of sperm cells identifies IL22 as a possible germ line risk locus for psoriatic arthritis. PLoS ONE 14(2): e0212043. https://doi.org/10.1371/journal.pone.0212043

Editor: Osman El-Maarri, University of Bonn, Institute of Experimental Hematology and Transfusion Medicine, GERMANY

Received: August 22, 2018; Accepted: January 20, 2019; Published: February 19, 2019

Copyright: © 2019 Pollock et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All datasets generated and analysed in this study have been uploaded to the Gene Expression Omnibus (GEO) repository (accession number GSE126017).

Funding: RAP was supported by a CIHR Banting and Best Doctoral Research Award and Krembil Research Institute Postdoctoral Fellowship, and LZ was supported by CIHR-IMHA and CIHR-Janssen studentships. The study was funded by the Arthritis Research Foundation, the National Psoriasis Foundation, and the Krembil Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interest exist.

Background

Psoriasis is a common inflammatory skin disease associated with significant morbidity, mortality, and poor quality of life that affects approximately 3% of Caucasians [1, 2]. Approximately 30% of psoriasis patients develop psoriatic arthritis (PsA), an inflammatory arthritis characterized by peripheral and/or axial arthritis, skin and nail disease, dactylitis, and enthesitis. The high recurrence risk ratios among first-degree relatives of psoriasis and PsA patients [3, 4], and higher disease concordance among monozygotic (62–70%) than dizygotic twins (21–23%) [5–8] suggest that both have a strong heritable component. Numerous susceptibility loci for psoriasis and PsA have been identified [9]; however together they account for only ~40% of the heritability of psoriasis and PsA [10].

The paradigm of DNA sequence variation as the sole substrate of heritability is challenged by animal studies demonstrating the ability of epigenetic marks to be inherited [11], as in the examples of agouti and axin mice, which carry DNA methylation variations within retrotransposons inserted into the agouti viable yellow (A^vy) and axin fused (Axin^FU) alleles, respectively, and result in variable expressivity of agouti and axin loci leading to variable coat colour [12] and kinked-tail [13] phenotypes. Moreover, methylation status of the Axin^FU allele can be inherited through both maternal and paternal transmissions; however, the penetrance of the kinked-tail phenotype is higher following paternal transmission [13]. The methylation status of the Axin^FU allele in somatic tissues is also reflected in sperm cells. These studies suggest that epigenetic marks can be inherited due to a failure to reset methylation of A^vy and Axin^FU in the germ line, and this resistance to resetting can result in a parental transmission bias or ‘parent-of-origin’ effect.

Parent-of-origin effects (POE) have been identified in humans in epidemiological analyses of multiple complex diseases including psoriasis and PsA. Several independent investigations, including studies of large psoriasis and PsA cohorts from the Faroe Islands, Scotland, and Canada [14–17], have consistently demonstrated a greater prevalence of psoriasis among the offspring of psoriatic fathers compared to psoriatic mothers, and a significantly greater tendency for psoriasis and PsA probands to report an affected father compared to an affected mother. Paternal transmission is accompanied by a significant reduction in age of psoriasis onset, and a tendency to manifest as the more severe PsA phenotype in subsequent generations.

Evidence of POE in psoriasis and PsA suggests that epigenetic mechanisms may contribute to the missing heritability in psoriatic disease. POE in humans may be mediated by epigenetic defects with or without an underlying genetic cause (known as secondary and primary epimutations, respectively [18, 19]). The aim of this study was to test the hypothesis that epigenetic marks contribute to the heritability of psoriatic disease by performing an epigenome-wide analysis of sperm cells from probands with PsA, cutaneous psoriasis without PsA (PsC) and controls to identify putative heritable DNA methylation variants. The second aim was to determine if these marks are present in somatic cells (whole blood) from the same patients and controls, and investigate the impact of underlying genetic variation on aberrant methylation marks.

Materials and methods

Samples

Male PsC and PsA probands were recruited from the University of Toronto Psoriatic Disease Program. All PsC patients were diagnosed by a dermatologist and examined by a rheumatologist to verify the absence of PsA. All PsA patients were diagnosed by a rheumatologist and satisfied the CASPAR criteria [20]. Unaffected male controls with no family history of PsC or PsA were recruited from the general population. All participants provided written informed consent and the study was conducted with approval from the University Health Network Research Ethics Board. Participants provided semen samples from which the motile fraction of mature spermatozoa was isolated by two-layer density gradient centrifugation using ISolate reagent (Irvine Scientific, Santa Ana, CA, USA) according to the manufacturer’s instructions.

DNA extraction and bisulfite conversion

Somatic cells were removed using a solution of 0.5% Triton X-100 and 0.1% SDS. Sperm cells were lysed with 400ul of 100mM Tris-Cl (pH 8), 10mM EDTA, 500nM NaCl, 1% SDS, and 2% B-mercaptoethanol and 100ul of Proteinase K (20mg/ml) at 55°C and 900rpm. An additional 50ul of Proteinase K was added after 2 hours, and again after 18 hours. After 20 hours of incubation, 20μl of RNase A/T1 and 10μl RNAse H were added and the samples incubated at 37°C for 30 minutes. DNA was extracted by phenol-chloroform extraction and quantity and purity were assessed by NanoDrop spectrophotometry. DNA (500ng) was treated with sodium bisulfite using the EZ DNA Methylation Kit (Zymo Research, Irvine, CA, USA) according to the manufacturer’s instructions for Illumina Infinium Methylation Assays.

Epigenome-wide methylation analysis

Bisulfite converted DNA samples were interrogated on Infinium HumanMethylation 450k v1 BeadChips (Illumina, San Diego, CA, USA) according to the manufacturer’s protocol. Arrays were scanned on the iScan system (Illumina, San Diego, CA, USA). Fluorescence intensities were quantified, and quality control was performed in GenomeStudio Version 2011.1 (Illumina, San Diego, CA, USA) using the HumanMethylation450_15017482_v.1.2 annotation file. Data were normalized against controls and background subtracted. The datasets generated and analysed during the current study are available in the Gene Expression Omnibus (GEO) repository (accession number GSE126017).

Bioinformatics and statistical analyses

Data were imported into the lumi package for colour balance adjustment and quantile normalization [21–23]. Type I vs II probe bias was corrected using the BMIQ method implemented through ChAMP [24]. A total of 485,577 CpG sites were initially assessed on the arrays. Probes were filtered based on detection p value (>0.01), cross-hybridization to multiple genomic locations, and presence of SNPs (MAF>5%) at the CpG or single base extension site. Furthermore, the least variable 25% of probes based on interquartile range (IQR) were removed [25, 26]. After filtering, 331,193 CpG sites were carried forward for analysis. The DMRcate package was used to identify contiguous genomic regions that differ between groups, or differentially methylated regions (DMRs), by kernel smoothing [27]. Comparisons planned a priori were: PsC probands vs. controls, PsA probands vs. controls, and PsA vs. PsC probands. DMRs were determined by linear regression with adjustment for clinical and technical covariates which significantly affected methylation. WebGestalt [28] was used to perform Overrepresentation Enrichment Analysis of gene ontologies, pathways, diseases, phenotypes, and chromosomal regions. Lists of DMRs from each comparison were tested for enrichment against a reference list of all of the CpGs analyzed. Enriched chromosomal regions were visualized using Integrated Genomics Viewer v2.3.

Pyrosequencing

For technical validation of array results, TPPP (cg22936884), CSMD1 (cg12027248), PTDSS2 (cg22993527), and HCG26 (cg16166826) methylation was assessed in bisulfite converted sperm DNA samples on a PyroMark Q24 system. For whole blood validation, whole blood DNA was extracted using a modified salting out technique (Gentra PureGene Blood kit), bisulfite converted, and methylation was measured at PTPRN2, ELF5, IL22, OSBPL5, H19, MBP, SNORD115, CARS2, and JAM3 as described above.

Transposable element and CNV genotyping

Genotyping of the 298bp HCG26 AluY insertion was performed by amplifying the surrounding region using primers flanking the insertion site (S1 Table). PCR was performed using 1X PCR buffer (Quanta Biosciences Inc.), 1.0mM MgCl₂, 200uM dNTPs, 500nM forward and reverse primers, 1U AccuStart II Taq polymerase (Quanta), and 100ng of DNA. Products were amplified on a 2720 thermal cycler (Applied Biosystems) at 94°C for 1 minute, 30 cycles of 94°C for 30 seconds, 60°C for 30 seconds, and 72°C for 90 seconds. Alleles containing the insertion yielded a 1697bp product, while alleles without the insertion yielded a 1399bp product, which were resolved on a 1.5% agarose gel. Genotyping of the IL22 CNV was performed using an inventoried Taqman Copy Number Assay (Hs00146600_cn), Taqman RNAseP Copy Number Reference Assay, and 1X Taqman Genotyping Master Mix (Life Technologies) following the manufacturer’s instructions and cycling conditions. Assays were performed in quadruplicate on a 7900HT system and analyzed using CopyCaller v2.1 software (Life Technologies).

Results

Investigation of clinical and demographic factors that affect sperm methylation

Following quality control, 54 subjects (23 PsC, 13 PsA probands, and 18 unaffected controls) were included in the analysis (Table 1). Clinical, demographic, and technical characteristics affecting methylation across all CpG sites were explored by linear regression. Methylation levels at many genes were significantly associated with array chip, subject’s age, age of psoriasis onset, Psoriasis Area and Severity Index (PASI- a clinical measure of psoriasis severity), and medication usage (non-steroidal anti-inflammatory drugs, disease modifying anti-rheumatic drugs, or biologics) (S1 Fig). For this reason, in the planned comparisons of PsC patients vs. controls and PsA patients vs. controls, age and slide were included as covariates, while in the comparison of PsA vs. PsC, age, slide, medications, age of psoriasis onset, and PASI were included as covariates.

Download:

Table 1. Demographic and clinical characteristics of the study subjects.

https://doi.org/10.1371/journal.pone.0212043.t001

Genome-wide trends in the psoriatic disease sperm methylome

The 331,193 CpGs assessed showed a bimodal distribution of methylation across all sperm samples analyzed, with the majority of sites showing either high methylation or low methylation levels (beta [β]>85% or <15%), consistent with previous reports [29]. Between PsC patients vs. controls there were 574 differentially methylated (DM) CpGs, 2,467 between PsA patients vs. controls, and 342 between PsA and PsC patients (mean FDR-adjusted p value [q]<0.05). The spatial distribution of these CpGs relative to genes and CpG islands is shown in Fig 1A and 1B.

Download:

Fig 1. The spatial distribution of DM CpGs by gene-centric and CpG island-centric annotations.

(A) No specific enrichment of DM CpG sites among gene-centric annotations was found. (B) In contrast, relative to CpG islands, DM CpG sites in PsC vs. controls were enriched in both north and south CpG island shores but depleted in north and south shelves and CpG islands (χ² = 17.51, df = 4, p = 0.0015). DM CpG sites in PsA vs. controls and PsA vs. PsC were not enriched in any specific CpG island-centric annotation (χ² = 5.58, df = 4, p = 0.233 and χ² = 9.41, df = 4, p = 0.052, respectively).

https://doi.org/10.1371/journal.pone.0212043.g001

Significant DM CpG sites were assembled into regions of contiguous differential methylation wherever the distance to the next significant probe was less than 1kb to prioritize functionally relevant regions of methylation change. Between PsC probands vs. controls the DM CpGs grouped into 169 differentially methylated regions (DMRs) of variable size, 754 DMRs in PsA probands vs. controls, and 86 DMRs between PsA and PsC probands (mean adjusted p<0.05 of constituent CpG sites). In all 3 comparisons, the majority of DMRs showed subtle changes in methylation (from -10% hypo to 10% hypermethylated) (Fig 2A–2C).

Download:

Fig 2. Distribution of the methylation changes observed within each set of significant DMRs.

(A) PsC vs. Controls (B) PsA vs. Controls (C) PsA vs. PsC. Data is plotted in bins of 0.1 (10%) change in methylation, representing the max. beta fold change among its constituent CpG sites, with each bin indicated by its middle value (i.e. changes ranging from 0–0.1 are indicated by the bin labelled 0.05). In all 3 comparisons, the majority of DMRs represented subtle changes in methylation (from -10% hypomethylated to 10% hypermethylated). Relatively few DMRs were beyond this range in PsC probands vs. controls. Differences larger than 10% were observed in PsA probands vs. controls, which were skewed towards hypermethylation, and in PsA vs. PsC probands which were normally distributed about 0 and ranged from approximately -45% hypomethylated to 35% hypermethylated.

https://doi.org/10.1371/journal.pone.0212043.g002

Differentially methylated regions in psoriatic disease are enriched in the major histocompatibility complex

Overrepresentation enrichment analysis yielded no significant enrichment of gene ontologies, pathways, networks, diseases and phenotypes among the lists of DMRs (q<0.05). However, DMRs in PsA probands vs. controls and PsC probands vs. controls were significantly enriched within the major histocompatibility complex on 6p21.3, a region which contains several loci associated with PsA and PsC, including HLA-B, HLA-C, MICA, HCP5, and TNFA [30]. Thirty-one DMRs between PsA probands vs. controls mapped this region (Hg19 coordinates chr6:26033821–35659141), 11 of which mapped to annotated genes (q = 4.31x10^-4). Eighteen DMRs between PsC probands vs. controls mapped to a smaller region contained therein (hg19 coordinates chr6:28889358–33677751), 10 of which mapped to annotated genes (q = 1.36x10^-4) (S2 Table). DMRs from both comparisons tended to occur outside the main psoriatic disease susceptibility region containing HLA-C and HLA-B. In PsA probands vs. controls, the closest DMRs to this region were found in HCG26 and LTB centromeric to HLA-B, and near HCG22 telomeric to HLA-C, while in PsC vs. controls, the closest DMRs to this region were found in HCG26, and just upstream of PSORS1C1 telomeric to HLA-C (S2 Fig).

Top differentially methylated regions in sperm cells contain biologically relevant genes

We prioritized DMRs for downstream analysis by limiting our analysis to those showing ≥10% change in methylation (hypo or hypermethylation). Of these, biological inferences were used to further prioritize biologically relevant sperm DMRs (Table 2). Several genes relevant to psoriatic disease pathogenesis, such as those involved in innate immunity, and skin and bone biology were identified. The top hit associated with skin disease (PsC probands vs. controls) was within the promoter of the immune-activating receptor SIGLEC14 (sialic acid binding Ig like lectin 14), which was hypermethylated by 10% (q = 4.32x10^-6) [31]. DMRs within the 3’UTR of junctional adhesion molecule 3 (JAM3), which regulates neutrophil transepithelial migration [32], and in the promoter region of procollagen-C endopeptidase enhancer (PCOLCE), a glycoprotein which binds and promotes the cleavage of type I procollagen to yield mature fibrillar collagen type I, a major component of the dermis, were 16% hypomethylated and 11% hypermethylated, respectively. Furthermore, a DMR within the body of retinoid X receptor beta (RXRB), which contains a 3'+140A polymorphism associated with chronic plaque psoriasis patients with a positive family history of disease [33], was 11% hypermethylated.

Download:

Table 2. Biologically relevant DMRs with ≥10% change in methylation identified in each of the 3 comparisons.

https://doi.org/10.1371/journal.pone.0212043.t002

Top hits representing skin and/or joint disease (PsA probands vs. controls) included a 14% hypermethylated DMR spanning the promoter and body of myelin basic protein (MBP), an autoantigen in multiple sclerosis. The promoter of oxysterol binding protein like 5 (OSBPL5), a paternally imprinted gene in the placenta that encodes an intracellular lipid receptor which plays a role in maintaining cholesterol balance, was also hypermethylated by 11%. The small nucleolar RNA 115 (SNORD115), another paternally imprinted gene located in the Prader-Willi Syndrome (PWS) region of chromosome 15, was hypermethylated by 16%. Notably, the promoter and body of HLA complex group 26 (HCG26) was hypomethylated by 25%. The function of HCG26 is unknown; however, it is located within the MHC between putative PsA risk loci MICA and MICB and lies in an intron of the adjacent locus HCP5, which has also been associated with PsA in GWAS.

DMRs associated with joint disease (PsA vs. PsC probands) included the promoter and body of interleukin-22 (IL22), a cytokine produced by IL-23-driven Th17 cells, activated γδ T cells, CD8+ T cells, and monocytes [34], which was 10% hypermethylated in PsA probands. The promoter of E74 like ETS transcription factor 5 (ELF5), a transcription factor expressed solely in epithelium which regulates terminal differentiation of keratinocytes, was 14% hypermethylated in PsA probands. The body of protein phosphatase 2 regulatory subunit B delta (PPP2R2D) was 14% hypomethylated in PsA probands. Knockout of this gene in cancer cells has been found to inhibit apoptosis and enhance proliferation and cytokine production in T cells [35]. The body of protein tyrosine phosphatase, receptor type N2 (PTPRN2), which contains a SNP previously associated with PsA in GWAS [36], was 17% hypomethylated in PsA probands, while the promoter and body of HCG26 was also significantly hypomethylated by 31% in PsA vs. PsC probands.

Pyrosequencing validates array results for several genes except HCG26 due to a transposable element insertion

Array data were validated by pyrosequencing of the same sperm samples initially analyzed. Four DMRs representing a range of methylation differences were selected for validation (CSMD1, PTDSS2, TPPP, and HCG26). Array and pyrosequencing data were strongly correlated (Pearson’s r = 0.97) for three out of the four DMRs (CSMD1, PTDSS2, and TPPP, but not HCG26 [r = 0.53] (S3 Fig). HCG26 was found to contain a 298bp AluY transposable element insertion which is polymorphic, being present in some HCG26 transcripts but not others [37]. Pyrosequencing of the HCG26 DMR revealed 100% methylation of all CpG sites in all subjects, while genotyping of the AluY insertion showed a perfect correlation between zygosity and methylation measured by Infinium arrays (S3 Table), suggesting that the insertion interfered with the binding of the Infinium probe yielding inaccurate measurements on the arrays.

Methylation of IL22 in whole blood correlates with sperm and is independent of copy number variation

Two rounds of epigenome-wide reprogramming occur between generations, first in the blastocyst after fertilization, during which gamete-specific marks are erased to ensure totipotency, and second in early embryogenesis, during which imprinted marks are erased and re-set in the primordial germ cells depending on the sex of the embryo [38]. For a proband to inherit an epigenetic mark from a previous generation, it must become ‘imprinted-like’ and resist the first round of reprogramming [39], and for the same proband to transmit it to the next generation, it must also resist the second round. DMRs present in sperm may be explained by the fact that they were inherited from a previous generation and have resisted the first and second rounds of epigenetic reprogramming.

To investigate whether sperm DMRs were inherited, we measured them in whole blood, a somatic tissue, collected from the same individuals who provided semen samples. Thirteen sperm DMRs were chosen based on adjusted p-value, β difference, biological or functional relevance, and the density of significant probes within the DMR. Successful pyrosequencing assays were developed for 9 of the 13 DMRs. Due to fact that they are distinct tissues, there were large and significant differences (p<0.0001) in methylation at all DMRs in paired blood compared to sperm samples. Most DMRs, apart from H19, showed nearly 100% methylation (80–100%) in blood, while showing a range in methylation in sperm samples from 0–40% methylation (IL22, ELF5, and PTPRN2), to 40–80% methylation (JAM3, CARS2, SNORD115, OBSPL5, H19, and MBP) (Fig 3). H19, which is paternally imprinted and normally heavily methylated on the paternal allele, was found to have higher methylation levels in sperm cells (60–80%) compared to blood cells (~50%).

Download:

Fig 3. Comparison of methylation levels in blood and sperm samples.

Methylation was measured in paired blood samples collected from the same individuals who provided semen samples by pyrosequencing. Differences between tissues were analyzed by paired t-test. ****, p<0.0001.

https://doi.org/10.1371/journal.pone.0212043.g003

Despite these tissue-specific differences, there was a significant positive correlation between sperm and blood methylation for IL22 (ρ = 0.34 [95% CI 0.06–0.57], p = 0.015, Fig 4), but not the other loci (S4 Table). Exon 1 of IL22 contains a CNV previously associated with psoriasis vulgaris [40]. This CNV encompasses the binding site of an Infinium probe within the IL22 DMR. Due to the potential of CNVs to bias methylation measurements, the IL22 exon 1 CNV was genotyped in all sperm and blood samples. No variation in copy number was found (all samples contained two copies, S5 Table), indicating that IL22 methylation cannot be attributed to the underlying CNV and reflect true methylation levels at this locus.

Download:

Fig 4. Correlation between sperm and blood methylation of IL22.

Methylation of exon 1 of IL22 was significantly correlated across sperm and blood samples, as measured by Infinium arrays (sperm) and pyrosequencing (blood).

https://doi.org/10.1371/journal.pone.0212043.g004

Discussion

Despite several large-scale efforts to identify genetic underpinnings of psoriatic disease, much of its heritability remains unexplained. Here, we explored whether epigenetic factors contribute to its heritability by taking the novel approach of genome-wide methylation profiling of sperm cells. An advantage of beginning our investigation with sperm cells is that they were easily attainable and could be purified into a homogeneous cell preparation, allowing us to avoid the confounding effects of cellular heterogeneity.

As sperm cells of psoriasis and PsA patients have never been studied, we were interested in large genomic and small-scale differences between patients and controls. One study of human sperm cells found that the largest degree of epigenetic variation occurs at functionally important promoter-associated CpG islands [31]. When associated with transcriptional start sites, CpG island methylation is associated with long-term gene silencing such as imprinting, or with genes expressed predominantly in germ cells and some tissue-specific genes [33]. We similarly found that the DM CpG sites in PsC vs. controls were enriched in such regions, as well as in north and south CpG island shores—areas whose functional relevance is still unclear [33].

DM CpG sites were assembled into larger DMRs to aid in identifying functionally relevant regions of methylation change [27]. As this effectively averages methylation across large regions, DMRs were overall quite subtle between groups, with most changes <10% in each comparison. Other studies comparing methylation at CpG site level have also noted that although there is considerable epigenetic variation in human sperm cells [41], these variations are quite subtle, or are present in very low frequencies of cells (<1%) [41, 42], consistent with our observations. Although few DMRs exceeded a 10% change in methylation, it is noteworthy that many of those which did play roles in the pathogenesis of skin disease (SIGLEC14, JAM3, PCOLCE, RXRB, ELF5, IL22), joint disease (MBP, HCG26, IL22, PPP2R2D, PTPRN2) or are known to be imprinted (OSBPL5, SNORD115) (Table 2). DMRs were also found to be enriched within the MHC, but were generally not in close proximity to the Class I psoriatic disease susceptibility loci HLA-C and HLA-B. Whether these DMRs represent distal regulatory elements for known psoriatic disease susceptibility loci remains to be determined.

We investigated candidate DMRs in blood to provide further evidence that they were inherited, but also to determine if observations in sperm can be translated to blood, a more pathogenically relevant tissue. The finding of IL22 as significantly hypermethylated in sperm of PsA compared to PsC probands and its significant correlation with methylation levels in the blood is particularly interesting given its role in inflammatory skin and joint disease. IL-22 is an IL-10 family cytokine that is secreted by a diverse repertoire of cells including Th1, Th17, Th22, and NK cells in response to IL-23. It is recognized by a heterodimeric receptor consisting of IL-10R2 and IL-22R1, with the latter being expressed specifically in lung, gastrointestinal, and skin epithelia but not on immune cells, making IL-22 a terminal effector cytokine [43]. IL-22 is suspected to play an important role in psoriatic skin disease by inducing hyperplasia and migration of epidermal keratinocytes, and expression of antimicrobial peptides, suggesting a role in wound healing and antimicrobial responses [43, 44]. It is highly overexpressed in psoriatic skin lesions and blood of psoriasis patients, and expression correlates with disease severity [45]. In mice, IL-22 is highly expressed in the entheses, and has been found to promote entheseal and periosteal bone formation [46]. Furthermore, IL-22 is expressed at higher levels in PsA compared to OA synovial fluid, and activated T cells from the synovium of PsA patients produce more IL-22 than OA patients, suggesting a role for IL-22 in PsA [47]. In this context, our finding of IL22 hypermethylation seems paradoxical given that it is usually associated with transcriptional repression. However, the DMR identified encompassed the gene body, which can be associated with moderate expression [48]. Unfortunately, we could not assess expression of IL22 in whole blood of the probands tested due to its low expression level in unstimulated ex vivo cells. Future studies will aim to validate IL22 hypermethylation and correlate with expression in stimulated cells to demonstrate the functional consequences of the identified epigenetic variants on transcription.

Although a DMR within the first exon of IL22 showed a significant correlation between blood and sperm cells, which is independent of an overlapping CNV in IL22, at present we cannot rule out the possibility that other genetic variants in IL-22 affect methylation of this DMR. Sequencing studies are necessary to demonstrate that methylation is independent of cis or trans-acting genetic mutations. If methylation is indeed independent of genetic variation in IL-22, it may be considered a putative primary epimutation warranting further study. For example, it will be necessary to know whether these germ line variants are present in normal somatic tissues derived from the three germ layers (ectoderm, mesoderm, and endoderm), as this would suggest that it is a ‘constitutional’ epimutation derived from the germ line and present in all tissues [18], lending further support for inheritance from the previous generation. Similarly, it would also be helpful to demonstrate the presence of these variants in the germ line of the parents, or the somatic tissues of the offspring of the patients studied here.

Conclusions

This study provides preliminary evidence of epigenetic variations in human sperm cells that are associated with PsC and PsA. These variations are generally subtle but occur near or within several genes that have potential pathogenic relevance. We identified a hypermethylated DMR within the first exon of IL22 which was correlated between sperm and whole blood as a putative epimutation associated with joint disease, but further studies are needed to rule out an underlying genetic cause and determine if this is widely distributed in somatic cells, which would suggest a heritable, constitutional epimutation. We also identified several other DMRs between psoriasis patients, PsA patients, and controls which may be the result of exposure of germ cells to endogenous or exogenous environmental factors.

Supporting information

S1 Table. PCR primers for amplifying the 298bp AluY insertion in HCG26.

https://doi.org/10.1371/journal.pone.0212043.s001

(PDF)

S2 Table. Enriched DMRs between PsA probands vs. controls and PsC probands vs. controls in the MHC that map to annotated genes.

https://doi.org/10.1371/journal.pone.0212043.s002

(PDF)

S3 Table. Genotyping results of the AluY insertion in HCG26 and corresponding Infinium array results.

0 copies of the insertion results in >90% methylation at the HCG26 locus, 1 copy results in a readout of ~55% methylation, while 2 copies results in a readout of ~20% methylation.

https://doi.org/10.1371/journal.pone.0212043.s003

(PDF)

S4 Table. Correlation between sperm and whole blood methylation across candidate DMRs.

https://doi.org/10.1371/journal.pone.0212043.s004

(PDF)

S5 Table. IL22 CNV genotyping results in blood and sperm samples.

https://doi.org/10.1371/journal.pone.0212043.s005

(PDF)

S1 Fig. Frequency of p values for the association of clinical, demographic, and technical characteristics with sperm CpG methylation.

https://doi.org/10.1371/journal.pone.0212043.s006

(PDF)

S2 Fig. Genetic map of DMRs in PsA and PsC probands vs. controls enriched in the MHC on 6p21.3.

Top-bottom, high-level to detailed view of chromosome 6 showing enrichment within the MHC. Dashes indicate location of DMRs between PsA probands vs. controls (red) and PsC probands vs. controls (blue). Bar graphs depict beta fold change and–log10 of the q values for each DMR. Green box indicates psoriatic disease susceptibility region including HLA-C and HLA-B.

https://doi.org/10.1371/journal.pone.0212043.s007

(PDF)

S3 Fig. Technical validation of arrays in sperm samples.

Correlation between Infinium % methylation and gold-standard bisulfite pyrosequencing (x = y included for comparison).

https://doi.org/10.1371/journal.pone.0212043.s008

(PDF)

Acknowledgments

We would like to thank Daniel Pereira and Sutha Shanmugarajah from the Psoriatic Arthritis Program for their help in subject recruitment.

References

1. Chandran V, Raychaudhuri SP. Geoepidemiology and environmental factors of psoriasis and psoriatic arthritis. J Autoimmun. 2010; 34:J314–321. pmid:20034760
- View Article
- PubMed/NCBI
- Google Scholar
2. Gladman D. Psoriatic Arthritis. In: Moderate to Severe Psoriasis. 2013, Informa Health Care: New York. p. 315–334.
3. Chandran V, Schentag CT, Brockbank JE, Pellett FJ, Shanmugarajah S, Toloza SM, et al. Familial aggregation of psoriatic arthritis. Ann Rheum Dis. 2009; 68:664–667. pmid:18524791
- View Article
- PubMed/NCBI
- Google Scholar
4. Rahman P,Elder JT. Genetic epidemiology of psoriasis and psoriatic arthritis. Ann Rheum Dis. 2005; 64 Suppl 2:ii37–39; discussion ii40-31.
- View Article
- Google Scholar
5. Pedersen OB, Svendsen AJ, Ejstrup L, Skytthe A,Junker P. On the heritability of psoriatic arthritis. Disease concordance among monozygotic and dizygotic twins. Ann Rheum Dis. 2008; 67:1417–1421. pmid:18218666
- View Article
- PubMed/NCBI
- Google Scholar
6. Bowcock AM,Krueger JG. Getting under the skin: the immunogenetics of psoriasis. Nat Rev Immunol. 2005; 5:699–711. pmid:16138103
- View Article
- PubMed/NCBI
- Google Scholar
7. Lonnberg AS, Skov L, Skytthe A, Kyvik KO, Pedersen OB,Thomsen SF. Heritability of psoriasis in a large twin sample. Br J Dermatol. 2013; 169:412–416. pmid:23574549
- View Article
- PubMed/NCBI
- Google Scholar
8. Duffy DL, Spelman LS,Martin NG. Psoriasis in Australian twins. J Am Acad Dermatol. 1993; 29:428–434. pmid:8349859
- View Article
- PubMed/NCBI
- Google Scholar
9. Eder L, Chandran V,Gladman DD. What have we learned about genetic susceptibility in psoriasis and psoriatic arthritis? Curr Opin Rheumatol. 2015; 27:91–98. pmid:25415529
- View Article
- PubMed/NCBI
- Google Scholar
10. Ellinghaus D, Jostins L, Spain SL, Cortes A, Bethune J, Han B, et al. Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nat Genet. 2016; 48:510–518. pmid:26974007
- View Article
- PubMed/NCBI
- Google Scholar
11. Petronis A. Epigenetics as a unifying principle in the aetiology of complex traits and diseases. Nature. 2010; 465:721–727. pmid:20535201
- View Article
- PubMed/NCBI
- Google Scholar
12. Morgan HD, Sutherland HG, Martin DI,Whitelaw E. Epigenetic inheritance at the agouti locus in the mouse. Nat Genet. 1999; 23:314–318. pmid:10545949
- View Article
- PubMed/NCBI
- Google Scholar
13. Rakyan VK, Chong S, Champ ME, Cuthbert PC, Morgan HD, Luu KV, et al. Transgenerational inheritance of epigenetic states at the murine Axin(Fu) allele occurs after maternal and paternal transmission. Proc Natl Acad Sci U S A. 2003; 100:2538–2543. pmid:12601169
- View Article
- PubMed/NCBI
- Google Scholar
14. Traupe H, Van Gurp PJM, Happle R, Boezeman J,Van de Kerkhof PCM. Psoriasis vulgaris, fetal growth, and genomic imprinting. American Journal of Medical Genetics. 1992; 42:649–654. pmid:1632432
- View Article
- PubMed/NCBI
- Google Scholar
15. Burden AD, Javed S, Bailey M, Hodgins M, Connor M,Tillman D. Genetics of psoriasis: paternal inheritance and a locus on chromosome 6p. J Invest Dermatol. 1998; 110:958–960. pmid:9620305
- View Article
- PubMed/NCBI
- Google Scholar
16. Rahman P, Gladman DD, Schentag CT,Petronis A. Excessive paternal transmission in psoriatic arthritis. Arthritis Rheum. 1999; 42:1228–1231. pmid:10366116
- View Article
- PubMed/NCBI
- Google Scholar
17. Pollock RA TA, Pellett F, Chandran V, Petronis A, Rahman P, Gladman DD. Further Evidence Supporting a Parent-of-Origin Effect in Psoriatic Disease. Arthritis Care Res (Hoboken). 2015; [Submitted, Under Review]:
- View Article
- Google Scholar
18. Oey H,Whitelaw E. On the meaning of the word 'epimutation'. Trends Genet. 2014; 30:519–520. pmid:25301328
- View Article
- PubMed/NCBI
- Google Scholar
19. Horsthemke B. Epimutations in human disease. Curr Top Microbiol Immunol. 2006; 310:45–59. pmid:16909906
- View Article
- PubMed/NCBI
- Google Scholar
20. Taylor W, Gladman D, Helliwell P, Marchesoni A, Mease P, Mielants H, et al. Classification criteria for psoriatic arthritis: development of new criteria from a large international study. Arthritis Rheum. 2006; 54:2665–2673. pmid:16871531
- View Article
- PubMed/NCBI
- Google Scholar
21. Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010; 11:587. pmid:21118553
- View Article
- PubMed/NCBI
- Google Scholar
22. Du P, Kibbe WA,Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008; 24:1547–1548. pmid:18467348
- View Article
- PubMed/NCBI
- Google Scholar
23. Pan Du GF, Spencer Huang, Warren , Simon Lin. Analyze Illumina Infinium methylation microarray data. R package version 2201. 2014;
- View Article
- Google Scholar
24. Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, et al. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics. 2014; 30:428–430. pmid:24336642
- View Article
- PubMed/NCBI
- Google Scholar
25. Gentleman R CV, Huber W and Hahne F. genefilter: genefilter: methods for filtering genes from high-throughput experiments. R package version 1500.
26. Chen YA, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013; 8:203–209. pmid:23314698
- View Article
- PubMed/NCBI
- Google Scholar
27. Peters TJ, Buckley MJ, Statham AL, Pidsley R, Samaras K, R VL, et al. De novo identification of differentially methylated regions in the human genome. Epigenetics Chromatin. 2015; 8:6. pmid:25972926
- View Article
- PubMed/NCBI
- Google Scholar
28. Wang J, Vasaikar S, Shi Z, Greer M,Zhang B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 2017; 45:W130–W137. pmid:28472511
- View Article
- PubMed/NCBI
- Google Scholar
29. Krausz C, Sandoval J, Sayols S, Chianese C, Giachini C, Heyn H, et al. Novel insights into DNA methylation features in spermatozoa: stability and peculiarities. PLoS One. 2012; 7:e44479. pmid:23071498
- View Article
- PubMed/NCBI
- Google Scholar
30. Pollock R, Chandran V. Genetics of psoriasis and psoriatic arthritis. In: Psoriasis and Psoriatic Arthritis: Pathophysiology, Therapeutic Intervention, and Complementary Medicine. Boca Raton: CRC Press; 2017. pp. 9–34.
31. Ali SR, Fong JJ, Carlin AF, Busch TD, Linden R, Angata T, et al. Siglec-5 and Siglec-14 are polymorphic paired receptors that modulate neutrophil and amnion signaling responses to group B Streptococcus. J Exp Med. 2014; 211:1231–1242. pmid:24799499
- View Article
- PubMed/NCBI
- Google Scholar
32. Chavakis T, Keiper T, Matz-Westphal R, Hersemeyer K, Sachs UJ, Nawroth PP, et al. The junctional adhesion molecule-C promotes neutrophil transendothelial migration in vitro and in vivo. J Biol Chem. 2004; 279:55602–55608. pmid:15485832
- View Article
- PubMed/NCBI
- Google Scholar
33. Vasku V, Bienertova Vasku J, Pavkova Goldbergova M,Vasku A. Three retinoid X receptor gene polymorphisms in plaque psoriasis and psoriasis guttata. Dermatology. 2007; 214:118–124. pmid:17341859
- View Article
- PubMed/NCBI
- Google Scholar
34. Zheng Y, Danilenko DM, Valdez P, Kasman I, Eastham-Anderson J, Wu J, et al. Interleukin-22, a T(H)17 cytokine, mediates IL-23-induced dermal inflammation and acanthosis. Nature. 2007; 445:648–651. pmid:17187052
- View Article
- PubMed/NCBI
- Google Scholar
35. Zhou P, Shaffer DR, Alvarez Arias DA, Nakazaki Y, Pos W, Torres AJ, et al. In vivo discovery of immunotherapy targets in the tumour microenvironment. Nature. 2014; 506:52–57. pmid:24476824
- View Article
- PubMed/NCBI
- Google Scholar
36. Stuart PE, Nair RP, Tsoi LC, Tejasvi T, Das S, Kang HM, et al. Genome-wide Association Analysis of Psoriatic Arthritis and Cutaneous Psoriasis Reveals Differences in Their Genetic Architecture. Am J Hum Genet. 2015; 97:816–836. pmid:26626624
- View Article
- PubMed/NCBI
- Google Scholar
37. Xing J, Zhang Y, Han K, Salem AH, Sen SK, Huff CD, et al. Mobile elements create structural variation: analysis of a complete human genome. Genome Res. 2009; 19:1516–1526. pmid:19439515
- View Article
- PubMed/NCBI
- Google Scholar
38. Kota SK,Feil R. Epigenetic transitions in germ cell development and meiosis. Dev Cell. 2010; 19:675–686. pmid:21074718
- View Article
- PubMed/NCBI
- Google Scholar
39. Martos SN, Tang WY,Wang Z. Elusive inheritance: Transgenerational effects and epigenetic inheritance in human environmental disease. Prog Biophys Mol Biol. 2015; 118:44–54. pmid:25792089
- View Article
- PubMed/NCBI
- Google Scholar
40. Prans E, Kingo K, Traks T, Silm H, Vasar E,Koks S. Copy number variations in IL22 gene are associated with Psoriasis vulgaris. Hum Immunol. 2013; 74:792–795. pmid:23395647
- View Article
- PubMed/NCBI
- Google Scholar
41. Flanagan JM, Popendikyte V, Pozdniakovaite N, Sobolev M, Assadzadeh A, Schumacher A, et al. Intra- and interindividual epigenetic variation in human germ cells. Am J Hum Genet. 2006; 79:67–84. pmid:16773567
- View Article
- PubMed/NCBI
- Google Scholar
42. Suter CM, Martin DI,Ward RL. Germline epimutation of MLH1 in individuals with multiple cancers. Nat Genet. 2004; 36:497–501. pmid:15064764
- View Article
- PubMed/NCBI
- Google Scholar
43. Johnston A,Gudjonsson JE. 22 again: IL-22 as a risk gene and important mediator in psoriasis. J Invest Dermatol. 2014; 134:1501–1503. pmid:24825058
- View Article
- PubMed/NCBI
- Google Scholar
44. Boniface K, Bernard FX, Garcia M, Gurney AL, Lecron JC,Morel F. IL-22 inhibits epidermal differentiation and induces proinflammatory gene expression and migration of human keratinocytes. J Immunol. 2005; 174:3695–3702. pmid:15749908
- View Article
- PubMed/NCBI
- Google Scholar
45. Boniface K, Guignouard E, Pedretti N, Garcia M, Delwail A, Bernard FX, et al. A role for T cell-derived interleukin 22 in psoriatic skin inflammation. Clin Exp Immunol. 2007; 150:407–415. pmid:17900301
- View Article
- PubMed/NCBI
- Google Scholar
46. Sherlock JP, Joyce-Shaikh B, Turner SP, Chao CC, Sathe M, Grein J, et al. IL-23 induces spondyloarthropathy by acting on ROR-gammat+ CD3+CD4-CD8- entheseal resident T cells. Nat Med. 2012; 18:1069–1076. pmid:22772566
- View Article
- PubMed/NCBI
- Google Scholar
47. Kotenko SV, Izotova LS, Mirochnitchenko OV, Esterova E, Dickensheets H, Donnelly RP, et al. Identification of the functional interleukin-22 (IL-22) receptor complex: the IL-10R2 chain (IL-10Rbeta) is a common chain of both the IL-10 and IL-22 (IL-10-related T cell-derived inducible factor, IL-TIF) receptor complexes. J Biol Chem. 2001; 276:2725–2732. pmid:11035029
- View Article
- PubMed/NCBI
- Google Scholar
48. Suzuki MM,Bird A. DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008; 9:465–476. pmid:18463664
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Chandran V, Raychaudhuri SP. Geoepidemiology and environmental factors of psoriasis and psoriatic arthritis. J Autoimmun. 2010; 34:J314–321. pmid:20034760
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Gladman D. Psoriatic Arthritis. In: Moderate to Severe Psoriasis. 2013, Informa Health Care: New York. p. 315–334.

[ref3] 3. Chandran V, Schentag CT, Brockbank JE, Pellett FJ, Shanmugarajah S, Toloza SM, et al. Familial aggregation of psoriatic arthritis. Ann Rheum Dis. 2009; 68:664–667. pmid:18524791
View Article
PubMed/NCBI
Google Scholar

[7] View Article

[8] PubMed/NCBI

[9] Google Scholar

[ref4] 4. Rahman P,Elder JT. Genetic epidemiology of psoriasis and psoriatic arthritis. Ann Rheum Dis. 2005; 64 Suppl 2:ii37–39; discussion ii40-31.
View Article
Google Scholar

[11] View Article

[12] Google Scholar

[ref5] 5. Pedersen OB, Svendsen AJ, Ejstrup L, Skytthe A,Junker P. On the heritability of psoriatic arthritis. Disease concordance among monozygotic and dizygotic twins. Ann Rheum Dis. 2008; 67:1417–1421. pmid:18218666
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref6] 6. Bowcock AM,Krueger JG. Getting under the skin: the immunogenetics of psoriasis. Nat Rev Immunol. 2005; 5:699–711. pmid:16138103
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref7] 7. Lonnberg AS, Skov L, Skytthe A, Kyvik KO, Pedersen OB,Thomsen SF. Heritability of psoriasis in a large twin sample. Br J Dermatol. 2013; 169:412–416. pmid:23574549
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref8] 8. Duffy DL, Spelman LS,Martin NG. Psoriasis in Australian twins. J Am Acad Dermatol. 1993; 29:428–434. pmid:8349859
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref9] 9. Eder L, Chandran V,Gladman DD. What have we learned about genetic susceptibility in psoriasis and psoriatic arthritis? Curr Opin Rheumatol. 2015; 27:91–98. pmid:25415529
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref10] 10. Ellinghaus D, Jostins L, Spain SL, Cortes A, Bethune J, Han B, et al. Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nat Genet. 2016; 48:510–518. pmid:26974007
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref11] 11. Petronis A. Epigenetics as a unifying principle in the aetiology of complex traits and diseases. Nature. 2010; 465:721–727. pmid:20535201
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref12] 12. Morgan HD, Sutherland HG, Martin DI,Whitelaw E. Epigenetic inheritance at the agouti locus in the mouse. Nat Genet. 1999; 23:314–318. pmid:10545949
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref13] 13. Rakyan VK, Chong S, Champ ME, Cuthbert PC, Morgan HD, Luu KV, et al. Transgenerational inheritance of epigenetic states at the murine Axin(Fu) allele occurs after maternal and paternal transmission. Proc Natl Acad Sci U S A. 2003; 100:2538–2543. pmid:12601169
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref14] 14. Traupe H, Van Gurp PJM, Happle R, Boezeman J,Van de Kerkhof PCM. Psoriasis vulgaris, fetal growth, and genomic imprinting. American Journal of Medical Genetics. 1992; 42:649–654. pmid:1632432
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref15] 15. Burden AD, Javed S, Bailey M, Hodgins M, Connor M,Tillman D. Genetics of psoriasis: paternal inheritance and a locus on chromosome 6p. J Invest Dermatol. 1998; 110:958–960. pmid:9620305
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref16] 16. Rahman P, Gladman DD, Schentag CT,Petronis A. Excessive paternal transmission in psoriatic arthritis. Arthritis Rheum. 1999; 42:1228–1231. pmid:10366116
View Article
PubMed/NCBI
Google Scholar

[58] View Article

[59] PubMed/NCBI

[60] Google Scholar

[ref17] 17. Pollock RA TA, Pellett F, Chandran V, Petronis A, Rahman P, Gladman DD. Further Evidence Supporting a Parent-of-Origin Effect in Psoriatic Disease. Arthritis Care Res (Hoboken). 2015; [Submitted, Under Review]:
View Article
Google Scholar

[62] View Article

[63] Google Scholar

[ref18] 18. Oey H,Whitelaw E. On the meaning of the word 'epimutation'. Trends Genet. 2014; 30:519–520. pmid:25301328
View Article
PubMed/NCBI
Google Scholar

[65] View Article

[66] PubMed/NCBI

[67] Google Scholar

[ref19] 19. Horsthemke B. Epimutations in human disease. Curr Top Microbiol Immunol. 2006; 310:45–59. pmid:16909906
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref20] 20. Taylor W, Gladman D, Helliwell P, Marchesoni A, Mease P, Mielants H, et al. Classification criteria for psoriatic arthritis: development of new criteria from a large international study. Arthritis Rheum. 2006; 54:2665–2673. pmid:16871531
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref21] 21. Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010; 11:587. pmid:21118553
View Article
PubMed/NCBI
Google Scholar

[77] View Article

[78] PubMed/NCBI

[79] Google Scholar

[ref22] 22. Du P, Kibbe WA,Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008; 24:1547–1548. pmid:18467348
View Article
PubMed/NCBI
Google Scholar

[81] View Article

[82] PubMed/NCBI

[83] Google Scholar

[ref23] 23. Pan Du GF, Spencer Huang, Warren , Simon Lin. Analyze Illumina Infinium methylation microarray data. R package version 2201. 2014;
View Article
Google Scholar

[85] View Article

[86] Google Scholar

[ref24] 24. Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, et al. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics. 2014; 30:428–430. pmid:24336642
View Article
PubMed/NCBI
Google Scholar

[88] View Article

[89] PubMed/NCBI

[90] Google Scholar

[ref25] 25. Gentleman R CV, Huber W and Hahne F. genefilter: genefilter: methods for filtering genes from high-throughput experiments. R package version 1500.

[ref26] 26. Chen YA, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013; 8:203–209. pmid:23314698
View Article
PubMed/NCBI
Google Scholar

[93] View Article

[94] PubMed/NCBI

[95] Google Scholar

[ref27] 27. Peters TJ, Buckley MJ, Statham AL, Pidsley R, Samaras K, R VL, et al. De novo identification of differentially methylated regions in the human genome. Epigenetics Chromatin. 2015; 8:6. pmid:25972926
View Article
PubMed/NCBI
Google Scholar

[97] View Article

[98] PubMed/NCBI

[99] Google Scholar

[ref28] 28. Wang J, Vasaikar S, Shi Z, Greer M,Zhang B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 2017; 45:W130–W137. pmid:28472511
View Article
PubMed/NCBI
Google Scholar

[101] View Article

[102] PubMed/NCBI

[103] Google Scholar

[ref29] 29. Krausz C, Sandoval J, Sayols S, Chianese C, Giachini C, Heyn H, et al. Novel insights into DNA methylation features in spermatozoa: stability and peculiarities. PLoS One. 2012; 7:e44479. pmid:23071498
View Article
PubMed/NCBI
Google Scholar

[105] View Article

[106] PubMed/NCBI

[107] Google Scholar

[ref30] 30. Pollock R, Chandran V. Genetics of psoriasis and psoriatic arthritis. In: Psoriasis and Psoriatic Arthritis: Pathophysiology, Therapeutic Intervention, and Complementary Medicine. Boca Raton: CRC Press; 2017. pp. 9–34.

[ref31] 31. Ali SR, Fong JJ, Carlin AF, Busch TD, Linden R, Angata T, et al. Siglec-5 and Siglec-14 are polymorphic paired receptors that modulate neutrophil and amnion signaling responses to group B Streptococcus. J Exp Med. 2014; 211:1231–1242. pmid:24799499
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref32] 32. Chavakis T, Keiper T, Matz-Westphal R, Hersemeyer K, Sachs UJ, Nawroth PP, et al. The junctional adhesion molecule-C promotes neutrophil transendothelial migration in vitro and in vivo. J Biol Chem. 2004; 279:55602–55608. pmid:15485832
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref33] 33. Vasku V, Bienertova Vasku J, Pavkova Goldbergova M,Vasku A. Three retinoid X receptor gene polymorphisms in plaque psoriasis and psoriasis guttata. Dermatology. 2007; 214:118–124. pmid:17341859
View Article
PubMed/NCBI
Google Scholar

[118] View Article

[119] PubMed/NCBI

[120] Google Scholar

[ref34] 34. Zheng Y, Danilenko DM, Valdez P, Kasman I, Eastham-Anderson J, Wu J, et al. Interleukin-22, a T(H)17 cytokine, mediates IL-23-induced dermal inflammation and acanthosis. Nature. 2007; 445:648–651. pmid:17187052
View Article
PubMed/NCBI
Google Scholar

[122] View Article

[123] PubMed/NCBI

[124] Google Scholar

[ref35] 35. Zhou P, Shaffer DR, Alvarez Arias DA, Nakazaki Y, Pos W, Torres AJ, et al. In vivo discovery of immunotherapy targets in the tumour microenvironment. Nature. 2014; 506:52–57. pmid:24476824
View Article
PubMed/NCBI
Google Scholar

[126] View Article

[127] PubMed/NCBI

[128] Google Scholar

[ref36] 36. Stuart PE, Nair RP, Tsoi LC, Tejasvi T, Das S, Kang HM, et al. Genome-wide Association Analysis of Psoriatic Arthritis and Cutaneous Psoriasis Reveals Differences in Their Genetic Architecture. Am J Hum Genet. 2015; 97:816–836. pmid:26626624
View Article
PubMed/NCBI
Google Scholar

[130] View Article

[131] PubMed/NCBI

[132] Google Scholar

[ref37] 37. Xing J, Zhang Y, Han K, Salem AH, Sen SK, Huff CD, et al. Mobile elements create structural variation: analysis of a complete human genome. Genome Res. 2009; 19:1516–1526. pmid:19439515
View Article
PubMed/NCBI
Google Scholar

[134] View Article

[135] PubMed/NCBI

[136] Google Scholar

[ref38] 38. Kota SK,Feil R. Epigenetic transitions in germ cell development and meiosis. Dev Cell. 2010; 19:675–686. pmid:21074718
View Article
PubMed/NCBI
Google Scholar

[138] View Article

[139] PubMed/NCBI

[140] Google Scholar

[ref39] 39. Martos SN, Tang WY,Wang Z. Elusive inheritance: Transgenerational effects and epigenetic inheritance in human environmental disease. Prog Biophys Mol Biol. 2015; 118:44–54. pmid:25792089
View Article
PubMed/NCBI
Google Scholar

[142] View Article

[143] PubMed/NCBI

[144] Google Scholar

[ref40] 40. Prans E, Kingo K, Traks T, Silm H, Vasar E,Koks S. Copy number variations in IL22 gene are associated with Psoriasis vulgaris. Hum Immunol. 2013; 74:792–795. pmid:23395647
View Article
PubMed/NCBI
Google Scholar

[146] View Article

[147] PubMed/NCBI

[148] Google Scholar

[ref41] 41. Flanagan JM, Popendikyte V, Pozdniakovaite N, Sobolev M, Assadzadeh A, Schumacher A, et al. Intra- and interindividual epigenetic variation in human germ cells. Am J Hum Genet. 2006; 79:67–84. pmid:16773567
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref42] 42. Suter CM, Martin DI,Ward RL. Germline epimutation of MLH1 in individuals with multiple cancers. Nat Genet. 2004; 36:497–501. pmid:15064764
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref43] 43. Johnston A,Gudjonsson JE. 22 again: IL-22 as a risk gene and important mediator in psoriasis. J Invest Dermatol. 2014; 134:1501–1503. pmid:24825058
View Article
PubMed/NCBI
Google Scholar

[158] View Article

[159] PubMed/NCBI

[160] Google Scholar

[ref44] 44. Boniface K, Bernard FX, Garcia M, Gurney AL, Lecron JC,Morel F. IL-22 inhibits epidermal differentiation and induces proinflammatory gene expression and migration of human keratinocytes. J Immunol. 2005; 174:3695–3702. pmid:15749908
View Article
PubMed/NCBI
Google Scholar

[162] View Article

[163] PubMed/NCBI

[164] Google Scholar

[ref45] 45. Boniface K, Guignouard E, Pedretti N, Garcia M, Delwail A, Bernard FX, et al. A role for T cell-derived interleukin 22 in psoriatic skin inflammation. Clin Exp Immunol. 2007; 150:407–415. pmid:17900301
View Article
PubMed/NCBI
Google Scholar

[166] View Article

[167] PubMed/NCBI

[168] Google Scholar

[ref46] 46. Sherlock JP, Joyce-Shaikh B, Turner SP, Chao CC, Sathe M, Grein J, et al. IL-23 induces spondyloarthropathy by acting on ROR-gammat+ CD3+CD4-CD8- entheseal resident T cells. Nat Med. 2012; 18:1069–1076. pmid:22772566
View Article
PubMed/NCBI
Google Scholar

[170] View Article

[171] PubMed/NCBI

[172] Google Scholar

[ref47] 47. Kotenko SV, Izotova LS, Mirochnitchenko OV, Esterova E, Dickensheets H, Donnelly RP, et al. Identification of the functional interleukin-22 (IL-22) receptor complex: the IL-10R2 chain (IL-10Rbeta) is a common chain of both the IL-10 and IL-22 (IL-10-related T cell-derived inducible factor, IL-TIF) receptor complexes. J Biol Chem. 2001; 276:2725–2732. pmid:11035029
View Article
PubMed/NCBI
Google Scholar

[174] View Article

[175] PubMed/NCBI

[176] Google Scholar

[ref48] 48. Suzuki MM,Bird A. DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008; 9:465–476. pmid:18463664
View Article
PubMed/NCBI
Google Scholar

[178] View Article

[179] PubMed/NCBI

[180] Google Scholar

Figures

Abstract

Background

Materials and methods

Samples

DNA extraction and bisulfite conversion

Epigenome-wide methylation analysis

Bioinformatics and statistical analyses

Pyrosequencing

Transposable element and CNV genotyping

Results

Investigation of clinical and demographic factors that affect sperm methylation

Genome-wide trends in the psoriatic disease sperm methylome

Differentially methylated regions in psoriatic disease are enriched in the major histocompatibility complex

Top differentially methylated regions in sperm cells contain biologically relevant genes

Pyrosequencing validates array results for several genes except HCG26 due to a transposable element insertion

Methylation of IL22 in whole blood correlates with sperm and is independent of copy number variation

Discussion

Conclusions

Supporting information

S1 Table. PCR primers for amplifying the 298bp AluY insertion in HCG26.

S2 Table. Enriched DMRs between PsA probands vs. controls and PsC probands vs. controls in the MHC that map to annotated genes.

S3 Table. Genotyping results of the AluY insertion in HCG26 and corresponding Infinium array results.

S4 Table. Correlation between sperm and whole blood methylation across candidate DMRs.

S5 Table. IL22 CNV genotyping results in blood and sperm samples.

S1 Fig. Frequency of p values for the association of clinical, demographic, and technical characteristics with sperm CpG methylation.

S2 Fig. Genetic map of DMRs in PsA and PsC probands vs. controls enriched in the MHC on 6p21.3.

S3 Fig. Technical validation of arrays in sperm samples.

Acknowledgments

References