Elleisy N, Rohde S, Huth A, Gittel N, Glass Ä, Möller S, Lamprecht G, Schäffler H, Jaster R. Genetic association analysis of CLEC5A and CLEC7A gene single-nucleotide polymorphisms and Crohn’s disease. World J Gastroenterol 2020; 26(18): 2194-2202 [PMID: 32476786 DOI: 10.3748/wjg.v26.i18.2194]
Corresponding Author of This Article
Robert Jaster, MD, Academic Research, Professor, Senior Postdoctoral Fellow, Senior Scientist, Department of Medicine II, Division of Gastroenterology, Rostock University Medical Center, E.-Heydemann-Str. 6, Rostock 18057, Germany. firstname.lastname@example.org
Checklist of Responsibilities for the Scientific Editor of This Article
This article is an open-access article which was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Genetic association analysis of CLEC5A and CLEC7A gene single-nucleotide polymorphisms and Crohn’s disease
Nagi Elleisy, Sarah Rohde, Astrid Huth, Nicole Gittel, Änne Glass, Steffen Möller, Georg Lamprecht, Holger Schäffler, Robert Jaster
Nagi Elleisy, Sarah Rohde, Astrid Huth, Nicole Gittel, Georg Lamprecht, Holger Schäffler, Robert Jaster, Department of Medicine II, Division of Gastroenterology, Rostock University Medical Center, Rostock 18057, Germany
Änne Glass, Steffen Möller, Institute for Biostatistics and Informatics in Medicine and Ageing Research, Rostock 18057, Germany
Author contributions: Schäffler H, Rohde S and Jaster R designed the study; Huth A, Schäffler H and Lamprecht G took reponsibility for patient care and follow-up; Elleisy N, Rohde S and Jaster R performed the experiments; Gittel N, Huth A and Schäffler H collected the samples and performed the clinical characterisation of the patients; Elleisy N, Glass Ä, Möller S and Jaster R performed the biostatistics; all authors analyzed the data; Schäffler H and Jaster R wrote the manuscript and contributed equally and share senior authorship.
Supported bythe Damp-Foundation, No. 2016-04.
Institutional review board statement: The study was approved by the ethic board of the Medical Faculty of the University of Rostock (A 2017-0137). Written informed consent was obtained from each participant prior to enrollment.
Conflict-of-interest statement: The authors declare that there is no conflict of interest.
Data sharing statement: No additional data are available.
Open-Access: This article is an open-access article that was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution NonCommercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Corresponding author: Robert Jaster, MD, Academic Research, Professor, Senior Postdoctoral Fellow, Senior Scientist, Department of Medicine II, Division of Gastroenterology, Rostock University Medical Center, E.-Heydemann-Str. 6, Rostock 18057, Germany. email@example.com
Received: February 4, 2020 Peer-review started: February 4, 2020 First decision: March 21, 2020 Revised: April 13, 2020 Accepted: April 29, 2020 Article in press: April 29, 2020 Published online: May 14, 2020
Crohn’s disease (CD) is characterized by a multifactorial etiology and a significant impact of genetic traits. While NOD2 mutations represent well established risk factors of CD, the role of other genes is incompletely understood.
To challenge the hypothesis that single nucleotide polymorphisms (SNPs) in the genes CLEC5A and CLEC7A, two members of the C-type lectin domain family of pattern recognition receptors, may be associated with CD.
SNPs in CLEC5A, CLEC7A and the known CD risk gene NOD2 were studied using real time PCR-based SNP assays. Therefore, DNA samples from 175 patients and 157 healthy donors were employed. Genotyping data were correlated with clinical characteristics of the patients and the results of gene expression data analyses.
In accordance with previous studies, rs2066844 and rs2066847 in NOD2 were found to be significantly associated with CD (allelic P values = 0.0368 and 0.0474, respectively). Intriguingly, for genotype AA of rs1285933 in CLEC5A, a potential association with CD (recessive P = 0.0523; odds ratio = 1.90) was observed. There were no associations between CD and SNPs rs2078178 and rs16910631 in CLEC7A. Variants of rs1285933 had no impact on CLEC5A gene expression. In contrast, genotype-dependent differences of CXCL5 expression in peripheral blood mononuclear cells were observed. There is no statistical interaction between the tested SNPs of NOD2 and CLEC5A, suggesting of a novel pathway contributing to the disease.
Our data encourage enlarged follow-up studies to further address an association of SNP rs1285933 in CLEC5A with CD. The C-type lectin domain family member also deserves attention regarding a potential role in the pathophysiology of CD.
Core tip: The genetic traits of Crohn’s disease (CD) are incompletely understood. Here, we report a potential association of single nucleotide polymorphism (SNP) rs1285933 in CLEC5A, a member of the C-type lectin domain family of pattern recognition receptors, with CD. Variants of SNP rs1285933 had no impact on CLEC5A gene expression in peripheral blood mononuclear cells but correlated with the expression of CXCL5. The SNPs rs2078178 and rs16910631 in CLEC7A were not associated with the disease. The role of CLEC5A in the pathophysiology of CD deserves further attention.
Citation: Elleisy N, Rohde S, Huth A, Gittel N, Glass Ä, Möller S, Lamprecht G, Schäffler H, Jaster R. Genetic association analysis of CLEC5A and CLEC7A gene single-nucleotide polymorphisms and Crohn’s disease. World J Gastroenterol 2020; 26(18): 2194-2202
Together with ulcerative colitis, Crohn’s disease (CD) represents the most common and clinically relevant inflammatory bowel disease (IBD)[1,2]. While it is generally accepted that the pathogenesis of the disease is multifactorial and involves an inappropriate activation of the mucosal immune system, the precise contribution of individual environmental factors and genetic traits remains elusive[1-3]. Mutations in the NOD2 gene represent the best-characterized genetic association of CD[4-6]. Nucleotide-binding oligomerization domain 2 (NOD2) belongs to the pattern recognition receptor (PRR) family and acts as an intracellular sensor for peptidoglycan[7,8] and its fragment muramyl dipeptide[9,10]. Downstream of NOD2, the transcription factor NF-κB plays a key role in the transduction of receptor-generated signals.
C-type lectin domain (CLEC) receptors comprise a large family of carbohydrate-binding proteins. Various CLEC family receptors are considered to exert functions as PRR since they recognize pathogen-associated molecules and may induce intracellular signaling pathways that regulate inflammatory processes. CLEC proteins are crucially involved in the immune response to fungal pathogens, but have also been implicated in anti-bacterial, anti-viral and anti-parasitic defense mechanisms[13,14]. Despite their functional similarities to NOD2, CLEC proteins have not been systematically studied in the context of IBD yet. Interestingly, a single nucleotide polymorphism (SNP) in the CLEC7A (DECTIN-1) gene, rs2078178, has been reported to be strongly linked to a severe form of ulcerative colitis, and this association was even stronger for the two-marker haplotype rs2078178 to rs16910631. For another CLEC gene, CLEC5A, we recently observed a CD-associated expression pattern with higher transcript levels in patient-derived peripheral blood mononuclear cells than in corresponding controls. Furthermore, CLEC5A showed a NOD2-dependent expression profile, supporting the hypothesis that both proteins may act in a regulatory network with a pathophysiological role in CD. Given that defective bacterial clearance may contribute to the pathogenesis of CD[17,18], it is important to note that CLEC5A has also been suggested to be essentially involved in innate immunity through neutrophil trap formation and secretion of different proinflammatory cytokines after stimulation with Listeria monocytogenes. Interestingly, the SNP rs1285933 in CLEC5A is associated with dengue severity, and CLEC5A has been shown to be critical for dengue-virus-induced lethal disease.
Here, we have addressed the question if the SNPs rs2078178 and rs16910631 in CLEC7A and rs1285933 in CLEC5A are associated with CD and have analyzed effects of rs1285933 at the level of gene expression. For comparison and a positive control, the known disease-associated SNPs rs2066844 (SNP8), rs2066845 (SNP12) and rs2066847 (SNP13)[5,6] in NOD2 were included into the investigations as well.
MATERIALS AND METHODS
From October 2015 until June 2017, 175 patients (102 females and 73 males; mean age 43.1 ± 14.7 years) with CD from the Department of Gastroenterology of Rostock University Medical Center (Rostock, Germany) were included in the study. This cohort of CD patients represents an extension of a cohort that we have previously characterized regarding relationships between mutations in the NOD2 gene, the disease phenotype and anti-tumor necrosis factor-α trough levels.
The diagnosis of CD was based on clinical, endoscopic, histological and radiological findings of the patients. The following clinical data were collected: Age, sex, age at diagnosis, duration of the disease, disease location, disease behavior, disease activity (assessed by the Crohn’s disease activity index and the Harvey–Bradshaw index), disease-specific medications, and previous history of surgery (i.e., colectomy). CD was stratified via the Montreal classification. Unrelated and healthy subjects from Germany (n = 157; 101 females and 56 males; mean age 25.3 ± 5.7 years) served as controls. The study was approved by the Local Ethics Board of the University of Rostock (A-2017-0137). We obtained written informed consent from all participants prior to their enrollment.
EDTA whole-blood samples were subjected to DNA extraction employing the QIAamp DNA blood mini kit according to the instructions of the manufacturer (Qiagen, Hilden, Germany).
Genotyping was performed using TaqMan™ SNP Genotyping Allelic Discrimination Assays with VIC- and FAM-labeled probes (Thermo Fisher Scientific, Karlsruhe, Germany) for rs1285933 (CLEC5A, Assay-ID: C__9506735_10), rs2078178 (CLEC7A; Assay-ID: C__1932439_10), rs16910631 (CLEC7A; Assay-ID: C__33748498_10), rs2066844 (NOD2, SNP8, Assay-ID: C__11717468_20), rs2066845 (NOD2; SNP12, Assay-ID: C__11717466_20), and rs2066847 (NOD2, Assay-ID: SNP13 C__60383785_10). PCR was carried out in 96-well plates, employing a ViiA 7 sequence detection system (Thermo Fisher Scientific). Thermal cycling conditions were: 95 °C for 10 min, followed by 40 cycles of 15 s at 95 °C/1 min at 60 °C. After PCR, fluorescence was detected and analyzed using TaqManGenotyper software version 1.3. Alternatively, NOD2 genotypes were determined by Sanger sequencing as described before.
In vitro studies with peripheral blood mononuclear cells
In this study, previous data from our laboratory were re-evaluated with respect to the rs1285933 genotype. Briefly, peripheral blood mononuclear cells (PBMC) had been isolated from EDTA venous blood, cultured and treated with lipopolysaccharide (1 µg/mL; Sigma-Aldrich, Deisenhofen, Germany) for 6 h. Afterwards, RNA was isolated, reversely transcribed into cDNA and subjected to real-time PCR employing standard procedures and a ViiA 7 sequence detection system. The following human-specific TaqManTM gene expression assays with fluorescently labeled MGB probes were used to quantify target cDNA levels: Hs04398399_m1 (CLEC5A), Hs01099660_g1 (CXCL5), and Hs99999905_m1 (GAPDH). PCR conditions were as follows: 95 °C for 10 min, followed by 40 cycles of 15 s at 95 °C/1 min at 60 °C.
The data were stored and analyzed employing IBM SPSS Statistics 25.0 (International Business Machines Corporation, Armonk, New York, United States). Differences between patients and controls were assessed for distributions (genotype, allele and sex) using the χ2 test or Fisher’s exact test, and for means using the t-test for independent samples (age, gene expression data), respectively. Pairwise statistical interaction between SNPs in a linear model was studied employing ANOVA. The Hardy-Weinberg equilibrium was assessed using the χ2 test with 1 degree of freedom. False discovery rates were controlled by using the Benjamini-Hochberg correction. Values of P < 0.05 were considered statistically significant.
SNP genotyping was performed on DNA samples from 175 patients with CD and 157 healthy controls. Both study groups are comparable for distribution of sex (P = 0.310), while patients with CD were older than healthy volunteers who served as controls (43.1 ± 14.7 vs 25.3 ± 5.7 years; P < 0.0001). In the context of this study, this age difference was considered acceptable. For the controls, the distribution of all individual SNP genotypes was in accordance with the Hardy-Weinberg equilibrium.
To study associations of CD with SNP genotypes or allele frequencies, four genetic models (genotype, dominant, recessive, or allelic models) were employed (Table 1). As expected, significant associations with CD were found for SNPs in NOD2, specifically rs2066844 (SNP8; genotype P = 0.0498, dominant P value = 0.0219, allelic P value = 0.0368) and rs2066847 (SNP13; allelic P value = 0.0474). Intriguingly, the genotype AA of rs1285933 in CLEC5A was also potentially associated with the disease (recessive model; P = 0.0523). The corresponding odds ratios (ORs) are shown in Table 2. For NOD2, the odds of having CD might triple in the presence of the risk allele T (rs2066844: OR = 3.29), and double with allele CC (rs2066847: OR = 2.31). Increased ORs are detectable for CLEC5A, too. Genotype AA almost doubles the odds of CD (OR = 1.90). Carrying the risk allele A increases the odds of CD by 39% (OR = 1.39), whereas allele G displays a protective effect (OR = 0.72). We could not detect significant associations between CD and the two SNPs in CLEC7A (rs2078178, rs16910631) and also not for rs2066845 (SNP12) in NOD2 (Table 1). The latter finding might be explained by the rare occurrence of the risk allele C in our cohorts of small size.
Table 1 Genotype and allele frequencies of single nucleotide polymorphisms in the genes CLEC5A, CLEC7A and NOD2 in Crohn’s disease patients and controls.
1Unadjusted for multiple testing. SNP: Single nucleotide polymorphism; NA: Not applicable (missing in controls); CI: Confidence interval.
We next compared patients with different genotypes of rs1285933 in CLEC5A (AA, AG and GG, respectively) regarding their clinical characteristics, employing the following parameters: Age, age at diagnosis, duration of the disease, disease location and behavior according to Montreal classification, Crohn’s disease activity index and Harvey–Bradshaw index, history of surgical treatment and treatment with drugs (including antibodies such as tumor necrosis factor-α inhibitors). There were no statistically significant differences between the three genotypes (data not shown).
To study potential functional effects of the rs1285933 polymorphism, we re-evaluated previously published gene expression data from our laboratory. In these studies, PBMC from CD patients and controls had been employed to measure the mRNA expression of a pre-selected set of genes. Using a combined data set from 16 CD patients and 6 healthy controls, we observed no genotype-dependent differences of CLEC5A gene expression (Figure 1A). On the other hand, we found that the genotype GG, compared to AG, was associated with significantly lower mRNA levels of the proinflammatory chemokine CXCL5 (Figure 1B, please note that ΔCt values and expression levels show an inverse and logarithmic relationship that follows the function 2-(∆Ct)).
Figure 1 Effects of the rs1285933 genotype on CLEC5A and CXCL5 gene expression.
Peripheral blood mononuclear cells were isolated from individuals with genotype AA (n = 8), GG (n = 5), and AG (n = 9), cultured and treated with lipopolysaccharide (1 µg/mL) for 6 h. Subsequently, the mRNA expression of the indicated genes and the house-keeping control GAPDH was analyzed by real-time PCR. Data are presented as averaged ΔCt values ± standard error of mean. aP < 0.05 vs genotype GG.
Located on different chromosomes, the disease-associate SNP of CLEC5A is not correlating with disease-associated SNPs of NOD2 (data not shown). Furthermore, the pairwise contributions to the disease phenotype of the CLEC5A SNP and the other SNPs are independent from each other (Table 3).
Table 3 Pairwise statistical interaction between single nucleotide polymorphisms in a linear model1.
1Disease—single nucleotide polymorphism (SNP) A + SNP B + SNP A: SNP B. The uncorrected P values for the last term of an ANOVA are presented in the table for all 30 interactions. The only significance was for two chromosomally neighboring SNPs within NOD2. SNP: Single nucleotide polymorphism; NA: Not applicable.
To the best of our knowledge, the results of this study suggest for the first time a potential association of SNP rs1285933 with CD. However, our findings need to be interpreted cautiously since they are based on a relatively small number of patients from a single center.
Given that the SNP is located within the CLEC5A gene, our data implicate a PRR beyond NOD2 into the pathogenesis of the disease. The mechanisms that underlie the effect of the CLEC5A polymorphism need to be further elucidated. To this end, reinvestigating data from our past work we can report a trans effect of rs1285933 on the expression of the chemokine CXCL5 in PBMC, but not of CLEC5A itself. These data suggest that CLEC5A might be functionally affected, e.g., with respect to its ability of ligand binding or downstream signaling. In accordance with this conclusion, SNP rs1285933 has also been suggested to modulate signaling pathways after interactions between the dengue virus and CLEC5A receptors. Other disease associations of SNP rs1285933 have not been reported yet. In a population of Taiwanese children, neither rs1285933 nor other polymorphisms of CLEC5A were associated with susceptibility to Kawasaki disease, coronary artery lesion formation, and intravenous immunoglobulin treatment response.
Of note, CLEC5A is embedded into an intronic region of another gene, MGAM, and the two transcripts are known to correlate, so that the effect of SNP rs1285933 is not necessarily exclusively related to the C-type lectin domain family member. Interestingly, decreased maltase activities in the small bowel mucosa are common in children with CD, and although this is of course no evidence for a genetic association, the role of MGAM in the context of IBD may deserve further attention as well.
We also evaluated possible associations of SNP rs1285933 with different clinical characteristics of our CD patients, including disease location, disease behavior and treatment history, but did not obtain significant results. Given that such effects have been reported for NOD2 variants[29-32], the studies are nevertheless worth to be continued in larger cohorts of patients. To this end, we conclude that the principal effect of SNP rs1285933 is modulation of CD susceptibility through a different molecular pathway than NOD2.
PRRs are key regulators of innate immune responses and inflammatory processes[13,14]. For a prominent member of this family, NOD2, a role in the pathogenesis of CD is clearly established[4-6]. Our results suggest an association of a polymorphism in another PRR, rs1285933 in CLEC5A, but not of rs2078178 and rs16910631 in CLEC7A, with CD. A systematic analysis of PRR functions in the context of CD might reveal novel pathomechanistic insights and help to identify new targets for diagnostic and therapy.
Crohn’s disease (CD) is characterized by a multifactorial etiology and a significant impact of genetic traits. While NOD2 mutations represent well established risk factors of CD, the role of other genes is incompletely understood.
A better knowledge of the molecular basis of CD is considered as an essential prerequisite for a further improvement of diagnostics and therapy.
Previous studies from our laboratory have pointed to a possible link between CD and the expression of pattern recognition receptors of the C-type lectin domain family (specifically, CLEC5A) in peripheral blood mononuclear cells (PBMC). This observation prompted us to ask if single nucleotide polymorphisms in the genes CLEC5A and CLEC7A might be associated with the disease.
DNA samples from patients with CD and healthy donors were subjected to the analysis of single nucleotide polymorphisms in the genes CLEC5A, CLEC7A and NOD2. For studies on gene expression, PBMC from subgroups of both cohorts were employed. Molecular findings were correlated with clinical characteristics of the patients.
For genotype AA of rs1285933 in CLEC5A, a potential association with CD and an increased odds ratio were detected. As expected, risk variants of NOD2 were associated with an increased occurrence of CD as well. Polymorphisms of rs1285933 correlated with CXCL5 gene expression but had no effect on CLEC5A expression in PBMC.
SNP rs1285933 in CLEC5A may represent a novel genetic association of CD. The finding, however, needs to be reproduced in multicenter studies with larger numbers of CD patients.
Pattern recognition receptors of the C-type lectin domain family deserve further attention regarding their potential role in the pathogenesis of CD and their relevance as diagnostic markers and therapeutic targets.
We thank Mrs. Katja Bergmann for expert technical assistance.
Hugot JP, Chamaillard M, Zouali H, Lesage S, Cézard JP, Belaiche J, Almer S, Tysk C, O'Morain CA, Gassull M, Binder V, Finkel Y, Cortot A, Modigliani R, Laurent-Puig P, Gower-Rousseau C, Macry J, Colombel JF, Sahbatou M, Thomas G. Association of NOD2 leucine-rich repeat variants with susceptibility to Crohn's disease.Nature. 2001;411:599-603.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 3772][Cited by in F6Publishing: 3523][Article Influence: 179.6][Reference Citation Analysis (0)]
Xavier-Carvalho C, Cezar RDDS, Freire NM, Vasconcelos CMM, Solorzano VEF, de Toledo-Pinto TG, Fialho LG, do Carmo RF, Vasconcelos LRS, Cordeiro MT, Baptista P, de Azeredo EL, da Cunha RV, de Souza LJ, Pacheco AG, Kubelka CF, Moura PMMF, Moraes MO. Association of rs1285933 single nucleotide polymorphism in CLEC5A gene with dengue severity and its functional effects.Hum Immunol. 2017;78:649-656.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 13][Cited by in F6Publishing: 11][Article Influence: 2.6][Reference Citation Analysis (0)]
Silverberg MS, Satsangi J, Ahmad T, Arnott ID, Bernstein CN, Brant SR, Caprilli R, Colombel JF, Gasche C, Geboes K, Jewell DP, Karban A, Loftus EV, Peña AS, Riddell RH, Sachar DB, Schreiber S, Steinhart AH, Targan SR, Vermeire S, Warren BF. Toward an integrated clinical, molecular and serological classification of inflammatory bowel disease: report of a Working Party of the 2005 Montreal World Congress of Gastroenterology.Can J Gastroenterol. 2005;19 Suppl A:5A-36A.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 1957][Cited by in F6Publishing: 1808][Article Influence: 244.6][Reference Citation Analysis (0)]
Schnitzler F, Friedrich M, Wolf C, Stallhofer J, Angelberger M, Diegelmann J, Olszak T, Tillack C, Beigel F, Göke B, Glas J, Lohse P, Brand S. The NOD2 Single Nucleotide Polymorphism rs72796353 (IVS4+10 A>C) Is a Predictor for Perianal Fistulas in Patients with Crohn's Disease in the Absence of Other NOD2 Mutations.PLoS One. 2015;10:e0116044.
[PubMed] [DOI][Cited in This Article: ][Cited by in Crossref: 13][Cited by in F6Publishing: 13][Article Influence: 1.9][Reference Citation Analysis (0)]