IL23R polymorphisms influence phenotype and response to therapy in patients with ulcerative colitis

Objective - We aimed to identify the clinical and genetic [IL23 receptor (IL23R) single nucleotide polymorphisms (SNPs)] predictors of response to therapy in patients with ulcerative colitis. Patients and methods - A total of 174 patients with ulcerative colitis, 99 women and 75 men, were included. The mean age of the patients was 47±15 years and the mean disease duration was 11±9 years. The number of patients classified as responders (R) or nonresponders (NR) to several therapies was as follows: 110 R and 53 NR to mesalazine (5-ASA), 28 R and 20 NR to azathioprine (AZT), 18 R and 7 NR to infliximab. Clinical and demographic variables were recorded. A total of four SNPs were studied: IL23R G1142A, C2370A, G43045A, and G9T. Genotyping was performed by real-time PCR using Taqman probes. Results - Older patients were more prone to respond to 5-ASA (P=0.004), whereas those with pancolitis were less likely to respond to such therapies (P=0.002). Patients with extraintestinal manifestations (EIMs) were less likely to respond to 5-ASA (P=0.001), AZT (P=0.03), and corticosteroids (P=0.06). Carriers of the mutant allele for IL23R SNPs had a significantly higher probability of developing EIMs (P<0.05), a higher probability of being refractory to 5-ASA (P<0.03), but a higher likelihood of responding to AZT (P=0.05). A significant synergism was observed between IL23R C2370A and EIMs with respect to nonresponse to 5-ASA (P=0.03). Conclusion - Besides extent of disease and age at disease onset, the presence of EIMs may be a marker of refractoriness to 5-ASA, corticosteroids, and AZT. IL23R SNPs are associated both with EIMs and with nonresponse to 5-ASA and corticosteroids.

No evidence of association between common autoimmunity STAT4 and IL23R risk polymorphisms and non-anterior uveitis.

OBJECTIVE: STAT4 and IL23R loci represent common susceptibility genetic factors in autoimmunity. We decided to investigate for the first time the possible role of different STAT4/IL23R autoimmune disease-associated polymorphisms on the susceptibility to develop non-anterior uveitis and its main clinical phenotypes. METHODS Four functional polymorphisms (rs3821236, rs7574865, rs7574070, and rs897200) located within STAT4 gene as well as three independent polymorphisms (rs7517847, rs11209026, and rs1495965) located within IL23R were genotyped using TaqMan® allelic discrimination in a total of 206 patients with non-anterior uveitis and 1553 healthy controls from Spain. RESULTS No statistically significant differences were found when allele and genotype distributions were compared between non-anterior uveitis patients and controls for any STAT4 (rs3821236: P=0.39, OR=1.12, CI 95%=0.87-1.43; rs7574865: P=0.59 OR=1.07, CI 95%=0.84-1.37; rs7574070: P=0.26, OR=0.89, CI 95%=0.72-1.10; rs897200: P=0.22, OR=0.88, CI 95%=0.71-1.08;) or IL23R polymorphisms (rs7517847: P=0.49, OR=1.08, CI 95%=0.87-1.33; rs11209026: P=0.26, OR=0.78, CI 95%=0.51-1.21; rs1495965: P=0.51, OR=0.93, CI 95%=0.76-1.15). CONCLUSION Our results do not support a relevant role, similar to that described for other autoimmune diseases, of IL23R and STAT4 polymorphisms in the non-anterior uveitis genetic predisposition. Further studies are needed to discard a possible weak effect of the studied variant.

rs1004819 is the main disease-associated IL23R variant in German Crohn's disease patients: combined analysis of IL23R, CARD15, and OCTN1/2 variants

BACKGROUND: The IL23R gene has been identified as a susceptibility gene for inflammatory bowel disease (IBD) in the North American population. The aim of our study was to test this association in a large German IBD cohort and to elucidate potential interactions with other IBD genes as well as phenotypic consequences of IL23R variants. METHODS: Genomic DNA from 2670 Caucasian individuals including 833 patients with Crohn's disease (CD), 456 patients with ulcerative colitis (UC), and 1381 healthy unrelated controls was analyzed for 10 IL23R SNPs. Genotyping included the NOD2 variants p.Arg702Trp, p.Gly908Arg, and p.Leu1007fsX1008 and polymorphisms in SLC22A4/OCTN1 (1672 C-->T) and SLC22A5/OCTN2 (-207 G-->C). RESULTS: All IL23R gene variants analyzed displayed highly significant associations with CD. The strongest association was found for the SNP rs1004819 [P = 1.92x10(-11); OR 1.56; 95 % CI (1.37-1.78)]. 93.2% of the rs1004819 TT homozygous carriers as compared to 78% of CC wildtype carriers had ileal involvement [P = 0.004; OR 4.24; CI (1.46-12.34)]. The coding SNP rs11209026 (p.Arg381Gln) was protective for CD [P = 8.04x10(-8); OR 0.43; CI (0.31-0.59)]. Similar, but weaker associations were found in UC. There was no evidence for epistasis between the IL23R gene and the CD susceptibility genes CARD15 and SLC22A4/5. CONCLUSION: IL23R is an IBD susceptibility gene, but has no epistatic interaction with CARD15 and SLC22A4/5. rs1004819 is the major IL23R variant associated with CD in the German population, while the p.Arg381Gln IL23R variant is a protective marker for CD and UC.

Novel genetic risk markers for ulcerative colitis in the IL2/IL21 region are in epistasis with IL23R and suggest a common genetic background for ulcerative colitis and celiac disease

OBJECTIVES: Recently, a genome-wide association study showed that single-nucleotide polymorphisms (SNPs) in the chromosome 4q27 region containing IL2 and IL21 are associated with celiac disease. Given the increased prevalence of inflammatory bowel disease (IBD) among celiac disease patients, we investigated the possible involvement of these SNPs in IBD. METHODS: Five SNPs strongly associated with celiac disease within the KIAA1109/TENR/IL2/IL21 linkage disequilibrium block on chromosome 4q27 and one coding SNP within the IL21 gene were analyzed in a large German IBD cohort. The study population comprised a total of 2,948 Caucasian individuals, including 1,461 IBD patients (ulcerative colitis (UC): n=514, Crohn's disease (CD): n=947) and 1,487 healthy unrelated controls. RESULTS: Three of the five celiac disease risk markers had a protective effect on UC susceptibility, and this effect remained significant after correcting for multiple testing: rs6840978: P=0.0082, P(corr)=0.049, odds ratio (OR) 0.77, 95% confidence interval (CI) 0.63-0.93; rs6822844: P=0.0028, P(corr)=0.017, OR 0.73, 95% CI 0.59-0.90; rs13119723: P=0.0058, P(corr)=0.035, OR 0.75, 95% CI 0.61-0.92. A haplotype consisting of the six SNPs tested was markedly associated with UC susceptibility (P=0.0025, P(corr)=0.015, OR 0.72, 95% CI 0.58-0.89). Moreover, in UC, epistasis was observed between the IL23R SNP rs1004819 and three SNPs in the KIAA1109/TENR/IL2/IL21 block (rs13151961, rs13119723, and rs6822844). CONCLUSIONS: Similar to other autoimmune diseases such as celiac disease, rheumatoid arthritis, type 1 diabetes, Graves' disease, and psoriatic arthritis, genetic variation in the chromosome 4q27 region predisposes to UC, suggesting a common genetic background for these diseases.

Investigation of IL23 (p19, p40) and IL23R highlights nuclear IL23p19 expression as a marker of indolent tumor features and favorable outcome in colorectal cancer patients

Objective: IL23 is involved in chronic inflammation but its role in cancer progression is not fully elucidated. Here we characterize IL23 subunits p40, p19 and IL23 receptor (IL23R) in the normal-adenoma-carcinomametastasis cascade of colorectal cancers and their relationship to clinicopathological and outcome data. Method: Immunohistochemistry for IL23R, IL12p40, IL23 and IL23p19 (monoclonal) was performed on a multi-punch tissue microarray (n=213 patients). Expression differences between normal-adenomas-cancerslymph nodes were evaluated. Correlation with clinicopathological and outcome data was undertaken. Results were validated on an independent cohort (n=341 patients). Results: An increased expression from normal-adenoma-cancer was observed (p<0.0001; all) followed by a marked reduction in lymph nodes (p<0.0001; all). Cytoplasmic and/or membranous staining of all markers was unrelated to outcome. Nuclear IL23p19 staining occurred in 23.1%and was associated with smaller tumor diameter (p=0.0333), early pT (p=0.0213), early TNM (p=0.0186), absence of vascular (p=0.0124) and lymphatic invasion (p=0.01493) and favorable survival (univariate (p=0.014) and multivariable (p=0.0321) analysis). All IL23p19 positive patients were free of distant metastasis (p=0.0146). Survival and metastasis results could be validated in Cohort 2. Conclusion: The presence of nuclear IL23p19 is related to indolent tumor features and favorable outcome supporting a more ‘protective’ role of this protein in colorectal cancer progression

Clinical and genetic factors predicting response to therapy in patients with Crohn’s disease

Aim - To identify clinical and/or genetic predictors of response to several therapies in Crohn’s disease (CD) patients. Methods - We included 242 patients with CD (133 females) aged (mean ± standard deviation) 39 ± 12 years and a disease duration of 12 ± 8 years. The single-nucleotide polymorphisms (SNPs) studied were ABCB1 C3435T and G2677T/A, IL23R G1142A, C2370A, and G9T, CASP9 C93T, Fas G670A and LgC844T, and ATG16L1 A898G. Genotyping was performed with real-time PCR with Taqman probes. Results - Older patients responded better to 5-aminosalicylic acid (5-ASA) and to azathioprine (OR 1.07, p = 0.003 and OR 1.03, p = 0.01, respectively) while younger ones responded better to biologicals (OR 0.95, p = 0.06). Previous surgery negatively influenced response to 5-ASA compounds (OR 0.25, p = 0.05), but favoured response to azathioprine (OR 2.1, p = 0.04). In respect to genetic predictors, we observed that heterozygotes for ATGL16L1 SNP had a significantly higher chance of responding to corticosteroids (OR 2.51, p = 0.04), while homozygotes for Casp9 C93T SNP had a lower chance of responding both to corticosteroids and to azathioprine (OR 0.23, p = 0.03 and OR 0.08, p = 0.02,). TT carriers of ABCB1 C3435T SNP had a higher chance of responding to azathioprine (OR 2.38, p = 0.01), while carriers of ABCB1 G2677T/A SNP, as well as responding better to azathioprine (OR 1.89, p = 0.07), had a lower chance of responding to biologicals (OR 0.31, p = 0.07), which became significant after adjusting for gender (OR 0.75, p = 0.005). Conclusions - In the present study, we were able to identify a number of clinical and genetic predictors of response to several therapies which may become of potential utility in clinical practice. These are preliminary results that need to be replicated in future pharmacogenomic studies.

Ankylosing spondylitis: genomic and functional characterization of candidate genes and their repercussion in clinical practice

RESUMO: Introdução: A espondilite anquilosante (EA) é uma doença inflamatória crónica caracterizada pela inflamação das articulações sacroilíacas e da coluna. A anquilose progressiva motiva uma deterioração gradual da função física e da qualidade de vida. O diagnóstico e o tratamento precoces podem contribuir para um melhor prognóstico. Neste contexto, a identificação de biomarcadores, assume-se como sendo muito útil para a prática clínica e representa hoje um grande desafio para a comunidade científica. Objetivos: Este estudo teve como objetivos: 1 - caracterizar a EA em Portugal; 2 - investigar possíveis associações entre genes, MHC e não-MHC, com a suscetibilidade e as características fenotípicas da EA; 3 - identificar genes candidatos associados a EA através da tecnologia de microarray. Material e Métodos: Foram recrutados doentes com EA, de acordo com os critérios modificados de Nova Iorque, nas consultas de Reumatologia dos diferentes hospitais participantes. Colecionaram-se dados demográficos, clínicos e radiológicos e colhidas amostras de sangue periférico. Selecionaram-se de forma aleatória, doentes HLA-B27 positivos, os quais foram tipados em termos de HLA classe I e II por PCR-rSSOP. Os haplótipos HLA estendidos foram estimados pelo algoritmo Expectation Maximization com recurso ao software Arlequin v3.11. As variantes alélicas dos genes IL23R, ERAP1 e ANKH foram estudadas através de ensaios de discriminação alélica TaqMan. A análise de associação foi realizada utilizando testes da Cochrane-Armitage e de regressão linear, tal como implementado pelo PLINK, para variáveis qualitativas e quantitativas, respetivamente. O estudo de expressão génica foi realizado por Illumina HT-12 Whole-Genome Expression BeadChips. Os genes candidatos foram validados usando qPCR-based TaqMan Low Density Arrays (TLDAs). Resultados: Foram incluídos 369 doentes (62,3% do sexo masculino, com idade média de 45,4 ± 13,2 anos, duração média da doença de 11,4 ± 10,5 anos). No momento da avaliação, 49,9% tinham doença axial, 2,4% periférica, 40,9% mista e 7,1% entesopática. A uveíte anterior aguda (33,6%) foi a manifestação extra-articular mais comum. Foram positivos para o HLA-B27, 80,3% dos doentes. Os haplótipo A*02/B*27/Cw*02/DRB1*01/DQB1*05 parece conferir suscetibilidade para a EA, e o A*02/B*27/Cw*01/DRB1*08/DQB1*04 parece conferir proteção em termos de atividade, repercussão funcional e radiológica da doença. Três variantes (2 para IL23R e 1 para ERAP1) mostraram significativa associação com a doença, confirmando a associação destes genes com a EA na população Portuguesa. O mesmo não se verificou com as variantes estudadas do ANKH. Não se verificou associação entre as variantes génicas não-MHC e as manifestações clínicas da EA. Foi identificado um perfil de expressão génica para a EA, tendo sido validados catorze genes - alguns têm um papel bem documentado em termos de inflamação, outros no metabolismo da cartilagem e do osso. Conclusões: Foi estabelecido um perfil demográfico e clínico dos doentes com EA em Portugal. A identificação de variantes génicas e de um perfil de expressão contribuem para uma melhor compreensão da sua fisiopatologia e podem ser úteis para estabelecer modelos com relevância em termos de diagnóstico, prognóstico e orientação terapêutica dos doentes. -----------ABSTRACT: Background: Ankylosing Spondylitis (AS) is a chronic inflammatory disorder characterized by inflammation in the spine and sacroiliac joints leading to progressive joint ankylosis and in progressive deterioration of physical function and quality of life. An early diagnosis and early therapy may contribute to a better prognosis. The identification of biomarkers would be helpful and represents a great challenge for the scientific community. Objectives: The present study had the following aims: 1- to characterize the pattern of AS in Portuguese patients; 2- to investigate MHC and non-MHC gene associations with susceptibility and phenotypic features of AS and; 3- to identify candidate genes associated with AS by means of whole-genome microarray. Material and Methods: AS was defined in accordance to the modified New York criteria and AS cases were recruited from hospital outcares patient clinics. Demographic and clinical data were recorded and blood samples collected. A random group of HLA-B27 positive patients and controls were selected and typed for HLA class I and II by PCR-rSSOP. The extended HLA haplotypes were estimated by Expectation Maximization Algorithm using Arlequin v3.11 software. Genotyping of IL23R, ERAP1 and ANKH allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests as implemented in PLINK, for dichotomous and quantitative variables, respectively. Gene expression profile was carried out using Illumina HT-12 Whole-Genome Expression BeadChips and candidate genes were validated using qPCR-based TaqMan Low Density Arrays (TLDAs). Results: A total of 369 patients (62.3% male; mean age 45.4±13.2 years; mean disease duration 11.4±10.5 years), were included. Regarding clinical disease pattern, at the time of assessment, 49.9% had axial disease, 2.4% peripheral disease, 40.9% mixed disease and 7.1% isolated enthesopathic disease. Acute anterior uveitis (33.6%) was the most common extra-articular manifestation. 80.3% of AS patients were HLA-B27 positive. The haplotype A*02/B*27/Cw*02/DRB1*01/DQB1*05 seems to confer susceptibility to AS, whereas A*02/B*27/Cw*01/DRB1*08/DQB1*04 seems to provide protection in terms of disease activity, functional and radiological repercussion. Three markers (two for IL23R and one for ERAP1) showed significant single-locus disease associations. Association of these genes with AS in the Portuguese population was confirmed, whereas ANKH markers studied did not show an association with AS. No association was seen between non-MHC genes and clinical manifestations of AS. A gene expression signature for AS was established; among the fourteen validated genes, a number of them have a well-documented inflammatory role or in modulation of cartilage and bone metabolism. Conclusions: A demographic and clinical profile of patients with AS in Portugal was established. Identification of genetic variants of target genes as well as gene expression signatures could provide a better understanding of AS pathophysiology and could be useful to establish models with relevance in terms of susceptibility, prognosis, and potential therapeutic guidance.

HLA and non-HLA genes in Behçet's disease: a multicentric study in the Spanish population.

INTRODUCTION According to genome wide association (GWA) studies as well as candidate gene approaches, Behçet's disease (BD) is associated with human leukocyte antigen (HLA)-A and HLA-B gene regions. The HLA-B51 has been consistently associated with the disease, but the role of other HLA class I molecules remains controversial. Recently, variants in non-HLA genes have also been associated with BD. The aims of this study were to further investigate the influence of the HLA region in BD and to explore the relationship with non-HLA genes recently described to be associated in other populations. METHODS This study included 304 BD patients and 313 ethnically matched controls. HLA-A and HLA-B low resolution typing was carried out by PCR-SSOP Luminex. Eleven tag single nucleotide polymorphisms (SNPs) located outside of the HLA-region, previously described associated with the disease in GWA studies and having a minor allele frequency in Caucasians greater than 0.15 were genotyped using TaqMan assays. Phenotypic and genotypic frequencies were estimated by direct counting and distributions were compared using the χ(2) test. RESULTS In addition to HLA-B*51, HLA-B*57 was found as a risk factor in BD, whereas, B*35 was found to be protective. Other HLA-A and B specificities were suggestive of association with the disease as risk (A*02 and A*24) or protective factors (A*03 and B*58). Regarding the non-HLA genes, the three SNPs located in IL23R and one of the SNPs in IL10 were found to be significantly associated with susceptibility to BD in our population. CONCLUSION Different HLA specificities are associated with Behçet's disease in addition to B*51. Other non-HLA genes, such as IL23R and IL-10, play a role in the susceptibility to the disease.

Follicular helper T cell in immunity and autoimmunity

The traditional concept that effector T helper (Th) responses are mediated by Th1/Th2 cell subtypes has been broadened by the recent demonstration of two new effector T helper cells, the IL-17 producing cells (Th17) and the follicular helper T cells (Tfh). These new subsets have many features in common, such as the ability to produce IL-21 and to express the IL-23 receptor (IL23R), the inducible co-stimulatory molecule ICOS, and the transcription factor c-Maf, all of them essential for expansion and establishment of the final pool of both subsets. Tfh cells differ from Th17 by their ability to home to B cell areas in secondary lymphoid tissue through interactions mediated by the chemokine receptor CXCR5 and its ligand CXCL13. These CXCR5+ CD4+ T cells are considered an effector T cell type specialized in B cell help, with a transcriptional profile distinct from Th1 and Th2 cells. The role of Tfh cells and its primary product, IL-21, on B-cell activation and differentiation is essential for humoral immunity against infectious agents. However, when deregulated, Tfh cells could represent an important mechanism contributing to exacerbated humoral response and autoantibody production in autoimmune diseases. This review highlights the importance of Tfh cells by focusing on their biology and differentiation processes in the context of normal immune response to infectious microorganisms and their role in the pathogenesis of autoimmune diseases.

IL-23 receptor and IL-12 receptor expression is restricted to distinct cell types in the IL-23R-GFP reporter mouse

Les maladies inflammatoires de l'intestin (MII) sont caractérisées par des réponses immunitaires incontrôlées dans l'intestin. Des études génétiques ont associé un polymorphisme dans le gène de l'IL23R à la résistance aux MII. IL23R code pour la protéine de l’IL-23r, une sous-unité du récepteur à l’IL-23 (IL-23R). Ce récepteur appartient à la famille de l’IL-12R, contenant plusieurs récepteurs hétérodimériques. D’ailleurs, IL-12R et IL-23R partagent la sous-unité IL12Rb1. Néanmoins, ces deux récepteurs favorisent des réponses immunitaires distinctes (Th1 vs Th17). Ce mémoire caractérise les dynamiques d’expression cellulaires de l’IL-23R et l’IL-12R, afin d’élucider leurs rôles dans l’inflammation. Nous avons établi qu’IL-23R et IL-12R ne sont jamais co-exprimés, malgré qu’ils partagent la sous-unité IL-12Rβ1. Parmi les cellules de rates de souris, la protéine IL-23r est trouvée dans certaines cellules T TCRγδ ou T CD4+, quelques cellules B et des cellules Lti-like. La protéine IL-12Rβ2 est exprimée par quelques cellules B. L’analyse de l’expression de l’IL-23R et l’IL-12R dans différents organes révéla que la plus grande proportion de cellules exprimant l’IL-23R se retrouve dans la lamina propria de l'intestin grêle, alors que les cellules exprimant l’IL-12Rβ2 ont été retrouvées en proportion équivalente dans tous les organes lymphoïdes. Ces observations appuient les études génétiques suggérant un rôle prédominant de l’IL23R dans les intestins. Finalement, des cultures in vitro suggèrent que l’IL-23R ou l’IL-12R avaient des réactions croisées à l’IL-12 ou l’IL-23. L’étude de l’IL-23R dans les MII devrait donc être complémentée par l’étude de l’IL-12R, car les deux récepteurs pourraient avoir des rôles complémentaires.

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.

β-glucan triggers spondylarthritis and Crohn's disease-like ileitis in SKG mice

The spondylarthritides (SpA), including ankylosing spondylitis (AS), psoriatic arthritis (PsA), reactive arthritis, and arthritis associated with inflammatory bowel disease, cause chronic inflammation of the large peripheral and axial joints, eyes, skin, ileum, and colon. Genetic studies reveal common candidate genes for AS, PsA, and Crohn's disease, including IL23R, IL12B, STAT3, and CARD9, all of which are associated with interleukin-23 (IL-23) signaling downstream of the dectin 1 β-glucan receptor. In autoimmune-prone SKG mice with mutated ZAP-70, which attenuates T cell receptor signaling and increases the autoreactivity of T cells in the peripheral repertoire, IL-17-dependent inflammatory arthritis developed after dectin 1-mediated fungal infection. This study was undertaken to determine whether SKG mice injected with 1,3-β-glucan (curdlan) develop evidence of SpA, and the relationship of innate and adaptive autoimmunity to this process.

Major histocompatibility complex and other allergy-related candidate genes associated with insect bite hypersensitivity in Icelandic horses

Insect bite hypersensitivity (IBH) is an allergic dermatitis of horses caused by bites of insects. IBH is a multifactorial disease with contribution of genetic and environmental factors. Candidate gene association analysis of IBH was performed in a group of 89 Icelandic horses all born in Iceland and imported to Europe. Horses were classified in IBH-affected and non-affected based on clinical signs and history of recurrent dermatitis, and on the results of an in vitro sulfidoleukotriene (sLT)-release assay with Culicoides nubeculosus and Simulium vittatum extract. Different genetic markers were tested for association with IBH by the Fisher's exact test. The effect of the major histocompatibility complex (MHC) gene region was studied by genotyping five microsatellites spanning the MHC region (COR112, COR113, COR114, UM011 and UMN-JH34-2), and exon 2 polymorphisms of the class II Eqca-DRA gene. Associations with Eqca-DRA and COR113 were identified (p < 0.05). In addition, a panel of 20 single nucleotide polymorphisms (SNPs) in 17 candidate allergy-related genes was tested. During the initial screen, no marker from the panel was significantly (p < 0.05) associated with IBH. Five SNPs associated with IBH at p < 0.10 were therefore used for analysis of combined genotypes. Out of them, SNPs located in the genes coding for the CD14 receptor (CD14), interleukin 23 receptor (IL23R), thymic stromal lymphopoietin (TSLP) and transforming growth factor beta 3 (TGFB3) molecules were associated with IBH as parts of complex genotypes. These results are supported by similar associations and by expression data from different horse populations and from human studies.