A fine mapping of single nucleotide variants and haplotype analysis of IL13 gene in patients with Leishmania guyanensis-cutaneous leishmaniasis and plasma cytokines IL-4, IL-5, and IL-13

Introduction

Leishmaniasis continues to pose a substantial health burden in 97 countries worldwide. The progression and outcome of Leishmania infection are influenced by various factors, including the cytokine milieu, the skin microbiota at the infection site, the specific Leishmania species involved, the genetic background of the host, and the parasite load. In endemic regions to leishmaniasis, only a fraction of individuals infected actually develops the disease. Overexpression of IL-13 in naturally resistant C57BL/6 mice renders them susceptible to L. major infection. Haplotypes constructed from several single nucleotide variant (SNV) along a chromosome fragment may provide insight into any SNV near the fragment that may be genuinely associated with a phenotype in genetic association studies.

Methods

We investigated nine SNVs (SNV1rs1881457A>C, SNV2rs1295687C>G, SNV3rs2069744C>T, SNV4rs2069747C>T, SNV5rs20541A>G, SNV6rs1295685A>G, SNV7rs848A>C, SNV8rs2069750G >C, and SNV9rs847T>C) spanning the entire IL13 gene in patients with L. guyanensis cutaneous leishmaniasis (Lg-CL).

Results

Our analysis did not reveal any significant association between the SNVs and susceptibility/protection against Lg-CL development. However, haplotype analysis, excluding SNV4rs2069747 and SNV8rs2069750 due to low minor allele frequency, revealed that carriers of the haplotype CCCTAAC had a 93% reduced likelihood developing Lg-CL. Similarly, the haplotypes ACCCGCT (ORadj=0.02 [95% CI 0.00-0.07]; p-value, 6.0×10-19) and AGCTAAC (ORadj=0.00[95% CI 0.00-0.00]; p-value 2.7×10-12) appeared to provide protection against the development of Lg-CL. Conversely, carriers of haplotype ACCTGCC have 190% increased likelihood of developing Lg-CL (ORadj=2.9 [95%CI 1.68-5.2]; p-value, 2.5×10-6). Similarly, haplotype ACCCAAT (ORadj=2.7 [95%CI 1.5-4.7]; p-value, 3.2×10-5) and haplotype AGCCGCC are associated with susceptibility to the development of Lg-CL (ORadj=1.7[95%CI 1.04-2.8]; p-value, 0.01). In our investigation, we also found a correlation between the genotypes of rs2069744, rs20541, rs1295685, rs847, and rs848 and plasma IL-5 levels among Lg-Cl patients. Furthermore, rs20541 showed a correlation with plasma IL-13 levels among Lg-Cl patients, while rs2069744 and rs848 showed a correlation with plasma IL-4 levels among the same group.

Conclusions

Overall, our study identifies three haplotypes of IL13 associated with resistance to disease development and three haplotypes linked to susceptibility. These findings suggest the possibility of a variant outside the gene region that may contribute, in conjunction with other genes, to differences in susceptibility and partially to the pathology.

Distinct plasma chemokines and cytokines signatures in Leishmania guyanensis-infected patients with cutaneous leishmaniasis

The immunopathology associated with Leishmaniasis is a consequence of inflammation. Upon infection with Leishmania, the type of host-immune response is determinant for the clinical manifestations that can lead to either self-healing or chronic disease. Multiple pathways may determine disease severity. A comparison of systemic immune profiles in patients with cutaneous leishmaniasis caused by L. guyanensis and healthy individuals with the same socio-epidemiological characteristics coming from the same endemic areas as the patients is performed to identify particular immune profile and pathways associated with the progression of disease development. Twenty-seven plasma soluble circulating factors were evaluated between the groups by univariate and multivariate analysis. The following biomarkers pairs IL-17/IL-9 (ρ=0,829), IL-17/IL-12 (ρ=0,786), IL-6/IL-1ra (ρ=0,785), IL-6/IL-12 (ρ=0,780), IL-1β/G-CSF (ρ=0,758) and IL-17/MIP-1β (ρ=0,754) showed the highest correlation mean among the patient while only INF-γ/IL-4 (ρ=0.740), 17/MIP-1β (ρ=0,712) and IL-17/IL-9 (ρ=0,707) exhibited positive correlation among the control group. The cytokine IL-17 and IL1β presented the greater number of positive pair correlation among the patients. The linear combinations of biomarkers displayed IP-10, IL-2 and RANTES as the variables with the higher discriminatory activity in the patient group compared to PDGF, IL-1ra and eotaxin among the control subjects. IP-10, IL-2, IL-1β, RANTES and IL-17 seem to be predictive value of progression to the development of disease among the Lg-infected individuals.

IL-23R variant rs11805303 is associated with susceptibility to the development of cutaneous leishmaniasis in Leishmania guyanensis-infected individuals

BACKGROUND

Emerging evidence suggests that the IL-17/IL-23 axis may play a role in the pathogenesis of leishmaniasis. Our aim was to investigate whether the IL-23R variant rs11805303 is a risk factor for the development of cutaneous leishmaniasis (CL) in Leishmania guyanensis-infected individuals.

METHODS

We genotyped by polymerase chain reaction-restriction fragment length polymorphism the rs11805303 C/T in 828 patients with CL and 806 healthy individuals. Plasma TNF-α, IL-6, IFG-γ, IL-1β and IL-17 were measured with the Bioplex assay.

RESULTS

The distribution of the genotypes differed between patients with CL and HC with a common odds ratio (OR) of 1.78 (p = 2.2 x10 -11) for the disease-associated T allele. L. guyanensis-infected individuals homozygous for the T allele show a 200% increased risk of progressing to disease development, with a 95% confidence interval (CI) ranging from 81% to 400% (P = 9.9 x 10 -6) in comparison to individuals homozygous for the C allele. Males homozygous for the T allele have higher plasma levels of IL-17 compared with heterozygous or homozygous CC individuals.

CONCLUSION

The present association of the IL-23R variant rs11805303 with the development of cutaneous leishmaniasis suggests that the IL-17/IL-23 axis may play an important role in the pathogenesis of CL.