Genetic association analysis of CLEC5A and CLEC7A gene singlenucleotide 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

Abstract

Key words: Crohn's disease; Single nucleotide polymorphisms; NOD2; CLEC5A; Gene expression; CXCL5

INTRODUCTION

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 theNOD2gene 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[11].

C-type lectin domain (CLEC) receptors comprise a large family of carbohydratebinding proteins[12]. 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 theCLEC7A(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[15]. For anotherCLECgene,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,CLEC5Ashowed aNOD2-dependent expression profile, supporting the hypothesis that both proteins may act in a regulatory network with a pathophysiological role in CD[16]. Given that defective bacterial clearance may contribute to the pathogenesis of CD[17,18], it is important to note thatCLEC5Ahas also been suggested to be essentially involved in innate immunity through neutrophil trap formation and secretion of different proinflammatory cytokines after stimulation withListeria monocytogenes[19].Interestingly, the SNP rs1285933 inCLEC5Ais associated with dengue severity[20], andCLEC5Ahas been shown to be critical for dengue-virus-induced lethal disease[21].

Here, we have addressed the question if the SNPs rs2078178 and rs16910631 inCLEC7Aand rs1285933 inCLEC5Aare 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]inNOD2were included into the investigations as well.

MATERIALS AND METHODS

Patients

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 theNOD2gene, the disease phenotype and anti-tumor necrosis factor-α trough levels[22].

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[23]and the Harvey-Bradshaw index[24]), disease-specific medications, and previous history of surgery (i.e., colectomy). CD was stratifiedviathe Montreal classification[25]. 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.

DNA extraction

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

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,NOD2genotypes were determined by Sanger sequencing as described before[22].

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[16]. 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 humanspecific TaqManTMgene 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.

Statistical analysis

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χ2test or Fisher's exact test, and for means using thet-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χ2test with 1 degree of freedom.False discovery rates were controlled by using the Benjamini-Hochberg correction.Values ofP＜ 0.05 were considered statistically significant.

RESULTS

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.7vs25.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 inNOD2,specifically rs2066844 (SNP8; genotypeP= 0.0498, dominantPvalue = 0.0219, allelicPvalue = 0.0368) and rs2066847 (SNP13; allelicPvalue = 0.0474). Intriguingly, the genotype AA of rs1285933 inCLEC5Awas also potentially associated with the disease(recessive model;P= 0.0523). The corresponding odds ratios (ORs) are shown in Table 2. ForNOD2, 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 forCLEC5A, 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 inCLEC7A(rs2078178,rs16910631) and also not for rs2066845 (SNP12) inNOD2(Table 1). The latter finding might be explained by the rare occurrence of the risk allele C in our cohorts of small size.

We next compared patients with different genotypes of rs1285933 inCLEC5A(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 reevaluated 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[16]. Using a combined data set from 16 CD patients and 6 healthy controls, we observed no genotype-dependent differences ofCLEC5Agene 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 chemokineCXCL5(Figure 1B, please note that ΔCt values and expression levels show an inverse and logarithmic relationship that follows the function 2-(ΔCt)).

Located on different chromosomes, the disease-associate SNP ofCLEC5Ais not correlating with disease-associated SNPs ofNOD2(data not shown). Furthermore, the pairwise contributions to the disease phenotype of theCLEC5ASNP and the other SNPs are independent from each other (Table 3).

DISCUSSION

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 theCLEC5Agene, our data implicate a PRR beyondNOD2into the pathogenesis of the disease. The mechanisms that underlie the effect of theCLEC5Apolymorphism need to be further elucidated. To this end,reinvestigating data from our past work[16]we can report atranseffect of rs1285933 on the expression of the chemokineCXCL5in PBMC, but not ofCLEC5Aitself. 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[20]. Other disease associations of SNP rs1285933 have not been reported yet. In a population of Taiwanese children, neither rs1285933 nor other polymorphisms ofCLEC5Awereassociated with susceptibility to Kawasaki disease, coronary artery lesion formation,and intravenous immunoglobulin treatment response[26].

Table 1 Genotype and allele frequencies of single nucleotide polymorphisms in the genes CLEC5A, CLEC7A and NOD2 in Crohn's disease patients and controls

Of note,CLEC5Ais embedded into an intronic region of another gene,MGAM, and the two transcripts are known to correlate[27], 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[28], and although this is of course no evidence for a genetic association, the role ofMGAMin 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 forNOD2variants[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 thanNOD2.

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 inCLEC5A, but not of rs2078178 and rs16910631 inCLEC7A, 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.

Table 2 Odds ratios of genotypes and alleles of single nucleotide polymorphisms in the genes CLEC5A and NOD2

Table 3 Pairwise statistical interaction between single nucleotide polymorphisms in a linear model1

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.

ARTICLE HIGHLIGHTS

Research background

Crohn's disease (CD) is characterized by a multifactorial etiology and a significant impact of genetic traits. WhileNOD2mutations represent well established risk factors of CD, the role of other genes is incompletely understood.

Research motivation

A better knowledge of the molecular basis of CD is considered as an essential prerequisite for a further improvement of diagnostics and therapy.

Research objectives

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 genesCLEC5AandCLEC7Amight be associated with the disease.

Research methods

DNA samples from patients with CD and healthy donors were subjected to the analysis of single nucleotide polymorphisms in the genesCLEC5A,CLEC7AandNOD2. For studies on gene expression, PBMC from subgroups of both cohorts were employed. Molecular findings were correlated with clinical characteristics of the patients.

Research results

For genotype AA of rs1285933 inCLEC5A, a potential association with CD and an increased odds ratio were detected. As expected, risk variants ofNOD2were associated with an increased occurrence of CD as well. Polymorphisms of rs1285933 correlated withCXCL5gene expression but had no effect onCLEC5Aexpression in PBMC.

Research conclusions

SNP rs1285933 inCLEC5Amay represent a novel genetic association of CD. The finding,however, needs to be reproduced in multicenter studies with larger numbers of CD patients.

Research perspectives

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.

ACKNOWLEDGEMENTS

We thank Mrs. Katja Bergmann for expert technical assistance.