Skip to main content

Immunological and genomic characterization of Ibizan Hound dogs in an endemic Leishmania infantum region

Abstract

Background

The Ibizan Hound is a canine breed native to the Mediterranean region, where leishmaniosis is an endemic zoonosis. Several studies indicate low prevalence of this disease in these dogs but the underlying molecular mechanism remains unknown.

Methods

In this study, qualitative immunological and genomic profiles of this breed have been analyzed.

Results

Our analysis shows relevant differences between the cytokine serum profile of Ibizan Hound dogs and previously published data from other canine strains. Additionally, several genetic risk variants related to the immune response, regulation of the immune system, and genes encoding cytokines and their receptors have been studied. The most relevant genes that presented such fixed polymorphisms were IFNG and IL6R. Other variants with frequencies ≥ 0.7 were found in the genes ARHGAP18, DAPK1, GNAI2, MITF, IL12RB1, LTBP1, SCL28A3, SCL35D2, PTPN22, CIITA, THEMIS, and CD180. Epigenetic regulatory genes such as HEY2 and L3MBTL3 showed also intronic polymorphisms.

Conclusions

Our analysis and results indicate that the regulation of immune responses is different in Ibizan Hounds compared to other breeds. Future studies are needed to elucidate whether these differences are related to the low prevalence of L. infantum infection in the Ibizan Hound.

Graphical Abstract

Background

The Ibizan Hound, also named “Podenco Ibicenco” or “Ca Eivissenc,” is the oldest of the 32 Spanish dog breeds officially recognized by the Spanish Kennel Club and the “Federation Cynologique Internationale.” The Federation classified the Ibizan Hound breed with the standard number 89, in group 5-Sect. 7 (Spitz and primitive types-primitive type hunting dogs), being accepted, without working trial in Spain since 1954, with an official valid standard since 1982 and official genealogical book since 2007 [1]. Although the Ibizan Hound breed is native to the Balearic Islands, it is also quite abundant in other Mediterranean regions in Spain such as Catalonia or Valencia, and in France, such as Roussillon or Provence. According to The American Kennel Club, the main physical characteristics of the Ibizan Hound are: (i) height ranging between 23.5 and 27.5 inches in males and 22.5 and 26 inches in females; (ii) weight ranging between 45 and 50 pounds in both males and females; (iii) mean life expectancy ranging between 11 and 14 years in both sexes [1].

Hunting dogs have an increased risk of infectious disease transmission due to (i) direct or indirect contacts with wildlife reservoirs, (ii) exposure to vector insects and ticks, and (iii) feeding habits like entrail ingestion in areas where downed game animals are field dressed [2]. Leishmaniosis is considered one of the most relevant vector-borne diseases in the Mediterranean region, where it is zoonotic, and its seroprevalence is estimated between 63 and 80% of the domestic dog population [3, 4], causing around 20,000–40,000 deaths per year in humans [5]. This disease is caused by different species of the genus Leishmania, including L. infantum, the most prevalent causal agent of leishmaniosis in the Mediterranean area, which is transmitted by Phlebotomus spp. [6].

The global seroprevalence in domestic dogs in Spain is estimated in 10.1%, one of the highest incidedences in the Mediterranean region [7]. In a recent study from Llobat’s group exploring 26 canine breeds, not including Ibizan Hounds, estimated the prevalence of L. infantum infection in the Ibizan Islands to be 4.7% (positive animals with or without clinical signs), being the most prevalent in the Doberman Pinscher and Boxer breeds, with seroprevalences of 42.9% and 39.1%, respectively [8].

Several studies have shown differences in immune response against L. infantum infection in Ibizan Hounds. Burhnam et al. [9] found higher concentrations of anti-Phlebotomus perniciosus antibodies in saliva and concomitant increase of interferon-γ (IFN-γ) serum levels in Ibizan Hounds compared to other canine breeds. Therefore, it has been proposed that such high IFN-γ values have a protective effect, since they induce the cellular immune response [10]. Moreover, Solano-Gallego et al. (2000) demonstrated a different specific immune response against L. infantum in Ibizan Hound specimens compared to other canine breeds, showing that 81% of Ibizan Hound dogs analyzed presented a cellular immune response vs. only 48% presented in the other dogs [11].

In addition, some genomic analyses suggest that several polymorphisms are involved in the molecular mechanisms associated with the resistance or susceptibility to canine leishmaniosis. Two polymorphisms present in the canine Slc11a1 gene, one intronic single-nucleotide polymorphism (SNP) and one silent SNP in exon 8, have been associated with an increased risk of canine visceral leishmaniosis [12]. Four SNPs in the canine β-defensin-1 (CBD1) gene have been associated with L. infantum infection in dogs from Italy [13].

As far as we know, no genomic studies in Ibizan Hound have been perfomed focusing on the immune response against L. infantum. Here, we present an immunological and genetic characterization of the Ibizan Hound breed to shed light on the underlying mechanisms leading to different immune responses.

Methods

Animal and epidemiological data

Twenty-eight Ibizan Hounds were analyzed. For all animals, epidemiological data were collected: sex, age, general and leishmaniosis vaccination status, external deworming status and type, living conditions (indoors or outdoors), type of parasitic feeding burden, and clinical signs.

Sample collection and cytokine levels

Ten milliliters of whole blood was taken by cephalic venipuncture using Vacutainer tubes without anticoagulant. Samples were maintained at room temperature to obtain serum aliquots, which were stored at − 20˚C until processing. Serological testing for L. infantum detection of specific antibodies was performed using the indirect immunofluorescent antibody test (IFAT) for anti-Leishmania-specific immunoglobulin G (IgG) antibodies (MegaFLUO® LEISH, Megacor Diagnostik GmbH, Hörbranz, Austria). Samples were considered seropositive with IFAT titer ≥ 1/80, following the manufacturer’s instructions. The whole blood samples were used for DNA extraction before 24 h to recovery.

The IL-2, IL-6, IL-8, IFN-γ (Canine IL-2 ELISA kit, Canine IL-6 ELISA kit, Canine IL-8 ELISA kit, and Canine IFN-γ ELISA kit, respectively; Invitrogen, Waltham, MA, USA,), and IL-18 (Canine IL-18 ELISA kit, Mybiosource, San Diego, CA, USA) levels were measured in serum samples by commercial ELISA kit method following the manufacturers' recommendations. In brief, 50 µl serum was used for the sandwich ELISA analysis. The microplate had been pre-coated with an antibody specific to cytokines. The samples were added to the microplate wells and combined with the specific antibody. Then, a biotinylated detection antibody specific for each cytokine and avidin-horseradish peroxidase (HRP) conjugate were added successively to each microplate well and incubated. Free components were washed away. The substrate solution was added to each well. The enzyme-substrate reaction was determined by the optical density (OD) and measured spectrophotometrically at a 450-nm wavelength in a Victor-X3™ plate reader (Perkin Elmer®). The concentration of each cytokine was calculated by comparing the OD of the samples to the standard curve.

DNA extraction and whole genome analysis

Genomic DNA (gDNA) from 28 Ibizan Hound samples was isolated using a QIAamp DNA Blood Kit following the manufacturer’s protocol (QIAamp; Qiagen, Hilden, Germany). DNA was quantified using a Nanodrop spectrophotometer (Thermo Fisher Scientific), and only samples with A260/A280 ratio > 1.8 were used. gDNA concentrations for all samples were a minimum of 50 ng/µl. DNA samples were whole genome amplified for 20–24 h at 37 ºC, fragmented, precipitated and resuspended in an appropriate hybridization buffer.

Samples were genotyped using the CanineHD DNA Analysis BeadChip WG-440–1001 (Illumina, Inc., San Diego, CA, USA) and hybridized on the prepared BeadChips for 16–24 h at 48 ºC. Following the hybridization, nonspecifically hybridized samples were removed by washing, while the remaining specifically hybridized loci were processed for the single-base extension reaction, stained, and imaged on an Illumina iScan Reader. GenomeStudio 2.0.5 (Illumina Inc., San Diego, CA, USA) was used to process data generated from the iScan system for subsequent analysis, according to manufacturer guidelines. Intensity data were loaded into the Genotyping Module to primary data analysis, including raw data normalization, clustering, and genotype calling.

SNPs on sexual chromosomes and with a call rate < 95% were discarded using PLINK v1.90b6.22. The final data set included 165,480 mapped positions in 24 Ibizan Hound samples with a mean genotyping rate of 0.988.

Analysis of polymorphisms related to immune response

The allelic status of previously published polymorphisms and genomic regions described as related to the immune response to L. infantum infection in dogs [12, 14,15,16,17,18,19] were interrogated in the analyzed samples.

PLINK v1.90b6.22 was used to extract variants from the selected genome regions, encompassing 38 genes, according to the mapping information of the Canis lupus familiaris genome assembly CanFam3.1. and to calculate allele frequencies. Polymorphisms were considered fixed in the Ibizan Hound breed when presenting frequencies > 0.7. Detected polymorphisms were annotated using NCBI refseqs, release 105. The rsID information was downloaded and annotated from the European Variation Archive EVA release 3 files corresponding to the CanFam3.1 assembly.

Results

First, we studied relevant epidemiological parameters of the Ibizan Hound such as gender, age, and diet; 60.71% of the analyzed Ibizan Hound specimens were males, categorized as adults (32.14%) or elderly (39.29%). Most (82.14%) were fed commercial food (Table 1). According to their lifestyle, all animals lived outdoors in contact with other dogs, used external deworming, and had the required vaccinations. None of the subjects were vaccinated against Leishmania. Likewise, none of the animals included in this study presented a positive result for antibodies against Leishmania according to the IFAT technique. Serum levels of the analyzed cytokines are summarized in Table 2.

Table 1 Epidemiological data of the 28 Ibizan Hound samples analyzed
Table 2 Concentration values of the analyzed cytokines

Second, we investigated the incidence of genomic variance in our cohort. In Tables 3, 4 and Additional file 1, information on 217 polymorphisms has been analyzed in 38 immune response-related genes. Table 3 summarizes the cytokine genes variants found in the genomic positions and their frequencies in the analyzed samples. In the IFN-γ gene, which is encoded by a gene located at chromosome 10, one polymorphism (rs22078594) was found with a frequency of 0.86 of the alternative alleles (A). In the IL2 gene, located at chromosome 2, a different variant was identified (rs22707631) in 3’ UTR with a 0.63 frequency of the alternative allele. No variants were found in the IL6 gene in chromosome 14; however, the gene that encodes its receptor, located at chromosome 7, showed two variants. The first variant rs24423644 in the 19:17,756,202 position was the most relevant (frequency of the alternative allele = 1) and the second rs24416159 in 7:42,902,063 position (frequency of the alternative allele = 0.78).

Table 3 Genomic variants found in genes encoding analyzed cytokines and their receptors
Table 4 Genomic variants found in genes other than cytokines analyzed with frequency of the alternative allele > 0.7

Polymorphisms in genes other than cytokines were also found in the analyzed samples. More precisely, intronic variants have been found in Rho GTPase Activating Protein 18 (ARHGAP18), Death Associated Protein Kinase 1 (DAPK1), G Protein subunit Alpha 2 (GNAI2), and Melanocyte Inducing Transcription Factor (MITF) genes, some of them with alterative allele frequency of 1. Other intronic and downstream polymorphisms have been found in different genes which encode for other interleukin receptors (IL12RB1), interleukin regulators such as Latent Transforming Growth Factor-β binding protein (LTBP1), two members of Solute Carrier Family (SCL28A3 and SCL35D2), and Protein Tyrosine Phosphatase Non-Receptor Type 22 (PTPN22), several genes related to expression gene regulation by epigenetic mechanisms (HEY2, L3MBTL3, and HABP4), membrane transport (ECHDC1 and RASEF), and immune system regulators (CIITA, THEMIS, and CD180), among others (Table 4). The results presented in Table 4 show those fixed variants in the population studied. The remaining SNPs were found with lower frequencies (additional file 1: Table S1).

Discussion

Our results show the immunological and genomic profile of the Ibizan Hound, which is the oldest of the 32 Spanish dog breeds officially recognized by the Spanish Kennel Club. In our analysis, several cytokine levels were compared to those previously published indicating a differential immunological profile. For the first time, serum levels of IL-2, IL-8, and IL-18 are studied in this canine breed.

Previous studies have indicated that the Ibizan Hound has a natural genetic resistance against leishmaniosis due to the early activation of the cellular immune response related to specific cytokines [11, 17, 20]. In fact, our results showed high levels of IFN-γ, IL-2, and IL-18, whereas the serum levels of IL-8 in this breed were lower than in other breeds [21]. The activation of macrophages by IFN-γ has been associated with effective T helper 1 (Th1) cellular immunity and the elimination of intracellular amastigotes and the control of the disease [11]. In Ibizan Hound dogs [9] a higher expression level of IFN-γ was observed in animals positive to the leishmanin skin test (LST) than in dogs with healthy skin in the nonendemic area [22].

Moreover, the upregulation of IFN-γ has been correlated with low disease severity [23]. Furthermore, stimulation with Leishmania infantum soluble antigen (LSA) in Ibizan Hound dogs provokes higher IFN-γ serum levels than in other healthy or sick dog breed [24]. The IFNG gene is located at chromosome 10 and presents four exons that codify a unique protein with 166 amino acids. Genomic analyses of IFNG show a fixed variant 10:10408724G > A (rs22078594) in an intronic region. This polymorphism produces a transcript named IFNG-201 with no associated biological functions [25]. However, since another intronic variant in IFNG has been correlated with both low levels of IFN-γ in humans and susceptibility to cutaneous leishmaniasis [26], more studies on the rs22078594 variant will be needed to elucidate its potential role in the regulation of the immune response.

The host ability to control Leishmania infection requires a strong cellular immune response, which is associated with the activation of T helper (Th)-1 cells producing IFN-γ and IL-2 [27, 28]. Additional evidence suggest that both factors have a key protective role in response to Leishmania infection [22, 29, 30]. Aslan et al. [34] found that IL-2 expression was negatively correlated with splenic parasite loads in experimentally infected dogs [31]. Our results are in agreement with previous data, indicating a natural resistance against leishmaniosis in Ibizan Hounds compared to other canine breeds [11, 24, 32]. The gene encoding for IL-2 is located at chromosome 19 and presents four exons. The variant found in IL2 gene 19, 17756202 T > C (rs22707631), overlaps one transcript of this gene [33]. This polymorphism is in the 3'-UTR region of the gene, which could have consequences in its post-transcriptional gene expression. In humans, several polymorphisms have been found in IL2, but it does not seem to predict susceptibility or protection against the development of disease to Leishmania guyanensis [34]. The biological changes associated with this variant are not reported in dogs or humans; thus, further studies will be necessary to shed light on the relationship between this 3'-UTR polymorphism and differences in serum levels of cytokines.

Additionally, IL-18 is a cytokine that plays an important role in both types of immunity, innate and acquired, and mediates Th1 response by inducing IFN-γ production in T cells and natural killer cells [35, 36]. Kumar et al. [40] found a significant relationship between susceptibility to visceral leishmaniasis in humans and the cytokine serum levels [37]. In dogs, it has been reported that IL-18 has no prominent role in the infection outcome of the disease [31, 38,39,40,41]. There is no evidence about the expression levels of IL-18 in different canine breeds, but similarly to IFN-g and IL-2, it is possible that IL-18 plays an important role in the natural resistance against leishmaniosis in Ibizan Hound dogs. Previous results from our group showed higher levels of IL-18 in Ibizan Hound than in the Boxer breed (data unpublished), which would be in agreement with these studies supporting the theory that the Ibizan Hound has a higher immune response against Leishmania infection than other canine breeds. IL18 gene is located at chromosome 5 and has seven exons. The array used for this study does not include any position in this gene, but it would be interesting to further analyze this gene and its genomic alterations to dissect whether they could regulate the elevated serum levels.

Low IL-8 serum level values were found in Ibizan Hound dogs. These results are consistent with the scarce studies focused on the role of IL-8 in Leishmania infection. Notably, IL-8 is an important cytokine for the activation of immune cells, especially neutrophils. Wardini et al. [45] found lower basal levels of IL-8 production in neutrophils in healthy dogs compared with infected animals [42]. A decrease in IL-8 levels has been reported in dogs infected both naturally and experimentally [43]. However, contrarily, it has been stated that amastigotes accumulate and survive in neutrophils in animals with high IL-8 levels [44], although these results are not conclusive. In humans, IL-8 polymorphism has been considered a risk factor for the development of visceral leishmaniasis [45].

Levels of the inflammatory cytokine IL-6 were similar to previously published data for other canine breeds [21]. Several authors found a higher expression of IL-6 in healthy than in experimentally infected sick dogs [20, 39]. Abbehusen et al. [46] demonstrated that serum levels of this cytokine increased in animals after natural and experimental infection [43]. De Lima et al. [49] pointed out that this cytokine is a good marker in active disease [46]. Other authors did not find a relationship between these cytokine levels and the development of leishmaniosis [31, 38]. The IL6 gene is located at chromosome 14 and presents five exons. No variants of this gene are included in the array. However, several SNPs have been found in the coding region of the IL-6 receptor (IL6R) and IL-6 cytokine family signal transducer (IL6ST) gene. For example, one SNP in 3'-UTR and one intronic polymorphism of IL6R are located at chromosome 7 and are fixed in the Ibizan Hound population studied. The 3’UTR variant results in a transcript named IL6R-201, which encodes a protein with 468 amino acids. The intronic polymorphism overlaps the same transcript as 3’UTR. Although these two SNPs are fixed in the Ibizan Hound population, their respective biological functions are largely unknown. Given that the results obtained in the serum levels of this interleukin were similar to those found in other canine breeds [21], it is possible that these variants do not influence them. IL-6 production increases in blood from dogs living in endemic areas of canine leishmaniosis cultured with several Toll-like receptors (TLR), in particular with TLR-4a and TLR-7a, but not with TLR-3a, whereas the IFN-γ production is not affected [47]. Other TLRs are related to L. infantum infection, so TLR-2 expression increases in mononucleate cells of the dermis with LST in Ibizan Hound dogs [48]. In human macrophages, the stimulation with agonist of TLR-2 increases the IFN-γ production through the PI3K-mTOR pathway [49]. In fact, in vitro studies have shown the ability of TLR-2 to reduce anti-leishmanial response mediated by IFN-γ [50], and TLR-2 is expressed in inflammatory cells such as biopsy samples of clinically affected skin [48]. The pathways of TLR and cytokines, and the relationship between them and their expression in macrophages of infected and non-infected dogs, and the polymorphisms of genes encoding these TLRs could be studied in future experiments to elucidate the mechanisms of resistance to L. infantum infection.

Different polymorphisms of several genes are fixed in the studied cohort of Ibizan Hound. Some are related to immune regulation. More precisely, six intronic SNPs with unknown biological function have been found in ARHGAP18 gene. This gene is related to regulation of lymphocyte-specific transcription factor interferon regulatory factor 4 (IRF4) [51]. Recently, Sundararaj and Casarotto [54] published a review indicating different functions of IRF4 in immune regulation, including B cell development [52]. Seven intronic variants are fixed in LAMA2 gene. Although LAMA2 is mainly related to muscular diseases (see editorial 70), some authors have indicated functions in immune response. For example, Lama2-deficient mouse model shows elevated expression of proinflammatory cytokines, concretely, TNF-α and IL1β [53], and Li et al. [57] demonstrated that LAMA2 is a coreceptor to activate T cells upon infection by S. aureus [54]. Thirteen intronic polymorphisms have been found in LTBP1 gene, a known gene involved in T regulatory lymphocyte differentiation [14]. In fact, Batista et al. [14] performed a genome-wide association study of cell-mediate response in dogs infected by L. infantum, and their results showed several genetic markers, which explain the phenotypic variance in cytokines expression [14]. These genetic markers are in some genes where fixed variants have been identified in this study, For example, in LTBP1 and in IL12RB1, which is associated with TNF-α and involved in cytokine and chemokine signaling, with one intronic and downstream polymorphism fixed in our dog population. Another SNP fixed in the Ibizan Hound population and related to immune response is DAPK1 gene variant. DAPK1 protein is related to mTORC1 activation and CD8 + cell function, by a recently discovered mechanism and is necessary to antiviral activity [55, 56] and most likely antiparasitic immune response. Additionally, there are polymorphisms in known genes that regulate the immune response through other pathways, e.g. GNAI2, which regulates B cell motility and entrance into lymph nodes [57, 58], MITF, inhibitor of IL-1, IL-6, CCL-2, and TNF-α [59], PTPN22, a negative regulator of T cell receptor and B cell receptor signaling [60, 61], CIITA, a regulator of Il-6 and HLA-II expression [62, 63], and THEMIS, which encodes a protein necessary to maintenance of CD8 + T cell response and regulation of cytokine signal [64]. Two family members of Solute Carrier Family (SLC28A3 and SLC35D2) show intronic variants. Polymorphisms in other members of this family have been associated with susceptibility or resistance against visceral leishmaniosis in different dog populations [12, 65,66,67,68], probably by regulation of survival in host macrophages to parasite [69], which could be performed by other members of the same family. Finally, genes encoding regulatory proteins involved in epigenetic regulation also show fixed polymorphisms, such as HEY2, L3MBTL3, and HABP4. HEY2 is targeted by Notch1 in humans [70] and is involved in immune response activation [71]. Moreover, HEY2 expression is regulated by the same epigenetic pathway as NR4A1 [72], a key mediator of T cell dysfunction [73]. Notch signaling is also regulated by L3MBTL3 as epigenetic regulator [74]. Overall, our results indicate that Ibizan Hound is a canine breed with a singular immune response against L. infantum infection and probably against other pathogens. Therefore, further studies will be necessary to elucidate the molecular mechanism underlaying the immune response and its relationship with the low prevalence of leishmaniosis in this breed.

Conclusions

Several studies indicate that the Ibizan Hound seems to present resistance against L. infantum infection, with lower prevalence of leishmaniosis than other canine breeds. This qualitative study shows high levels of IFN-γ, IL-2, and IL-18 in Ibizan Hound dogs, whereas the serum levels of IL-8 in this breed are lower than in other breeds. The genomic analysis reveals many variants which are fixed in this dog population and related to expression and regulation of the immune system. These results indicate that this breed presents a specific immune response, which could explain this low prevalence of leishmaniosis and would be related to resistance against other infection pathogens. Future studies will be necessary to discover the specific biological function of the genetic SNPs found here, more precisely in genes related to expression and regulation of cytokines levels.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

  1. Ibizan Hound Dog Breed Information. American Kennel Club. https://www.akc.org/dog-breeds/ibizan-hound/. Accessed on 30 May 2022.

  2. Ridgway M. Hunting dogs. Vet Clin North Am Small Anim Pract. 2021;51:877–90.

    Article  PubMed  Google Scholar 

  3. Leontides LS, Saridomichelakis MN, Billinis C, Kontos V, Koutinas AF, Galatos AD, et al. A cross-sectional study of Leishmania spp. infection in clinically healthy dogs with polymerase chain reaction and serology in Greece. Vet Parasitol. 2002;109:19–27.

    Article  PubMed  Google Scholar 

  4. Mohebali M, Malmasi A, Khodabakhsh M, Zarei Z, Akhoundi B, Hajjaran H, et al. Feline leishmaniosis due to Leishmania infantum in Northwest Iran: the role of cats in endemic areas of visceral leishmaniosis. Vet Parasitol Reg Stud Rep. 2017;9:13–6.

    Google Scholar 

  5. Alvar J, Vélez ID, Bern C, Herrero M, Desjeux P, Cano J, et al. Leishmaniasis worldwide and global estimates of its incidence. PLoS ONE. 2012;7:e35671.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. González E, Molina R, Iriso A, Ruiz S, Aldea I, Tello A, et al. Opportunistic feeding behaviour and Leishmania infantum detection in Phlebotomus perniciosus females collected in the human leishmaniasis focus of Madrid, Spain (2012–2018). PLoS Negl Trop Dis. 2021;15:e0009240.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Gálvez R, Montoya A, Cruz I, Fernández C, Martín O, Checa R, et al. Latest trends in Leishmania infantum infection in dogs in Spain, part I: mapped seroprevalence and sand fly distributions. Parasit Vectors. 2020;13:204.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Edo M, Marín-García PJ, Llobat L. Is the prevalence of Leishmania infantum linked to breeds in dogs? Characterization of seropositive dogs in Ibiza. Animals (Basel). 2021;11:2579.

    Article  PubMed  Google Scholar 

  9. Burnham AC, Ordeix L, Alcover MM, Martínez-Orellana P, Montserrat-Sangrà S, Willen L, et al. Exploring the relationship between susceptibility to canine leishmaniosis and anti-Phlebotomus perniciosus saliva antibodies in Ibizan hounds and dogs of other breeds in Mallorca. Spain Parasit Vectors. 2020;13:129.

    Article  CAS  PubMed  Google Scholar 

  10. Abdeladhim M, Kamhawi S, Valenzuela JG. What’s behind a sand fly bite? The profound effect of sand fly saliva on host hemostasis, inflammation and immunity. Infect Genet Evol. 2014;28:691–703.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Solano-Gallego L, Llull J, Ramos G, Riera C, Arboix M, Alberola J, et al. The Ibizian hound presents a predominantly cellular immune response against natural Leishmania infection. Vet Parasitol. 2000;90:37–45.

    Article  CAS  PubMed  Google Scholar 

  12. Sanchez-Robert E, Altet L, Utzet-Sadurni M, Giger U, Sanchez A, Francino O. Slc11a1 (formerly Nramp1) and susceptibility to canine visceral leishmaniasis. Vet Res. 2008;39:36.

    Article  PubMed  Google Scholar 

  13. da Silva LG, Costa-Júnior CRL, Figueiredo-Júnior CAS, Leal-Balbino TC, Crovella S, Otranto D, et al. Canine β-defensin-1 (CBD1) gene as a possible marker for Leishmania infantum infection in dogs. Parasit Vectors. 2017;10:199.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Batista LFS, Utsunomiya YT, Silva TBF, Dias RA, Tomokane TY, Pacheco AD, et al. Genome-wide association study of cell-mediated response in dogs naturally infected by Leishmania infantum. Infect Immun. 2016;84:3629–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Batista LFS, Torrecilha RBP, Silva RB, Utsunomiya YT, Silva TBF, Tomokane TY, et al. Chromosomal segments may explain the antibody response cooperation for canine leishmaniasis pathogenesis. Vet Parasitol. 2020;288:109276.

    Article  CAS  PubMed  Google Scholar 

  16. Batista LFS, Utsunomiya YT, Silva TBF, Carneiro MM, Paiva JSF, Silva RB, et al. Canine leishmaniasis: genome-wide analysis and antibody response to Lutzomyia longipalpis saliva. PLoS ONE. 2018;13:e0197215.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Soutter F, Solano-Gallego L, Attipa C, Gradoni L, Fiorentino E, Foglia Manzillo V, et al. An investigation of polymorphisms in innate and adaptive immune response genes in canine leishmaniosis. Vet Parasitol. 2019;269:34–41.

    Article  CAS  PubMed  Google Scholar 

  18. Utsunomiya YT, Ribeiro ÉS, Quintal APN, Sangalli JR, Gazola VR, Paula HB, et al. Genome-wide scan for visceral leishmaniasis in mixed-breed dogs identifies candidate genes involved in T helper cells and macrophage signaling. PLoS ONE. 2015;10:e0136749.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Altet L, Francino O, Solano-Gallego L, Renier C, Sánchez A. Mapping and sequencing of the canine NRAMP1 gene and identification of mutations in leishmaniasis-susceptible dogs. Infect Immun. 2002;70:2763–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Martínez-Orellana P, Marí-Martorell D, Montserrat-Sangrà S, Ordeix L, Baneth G, Solano-Gallego L. Leishmania infantum-specific IFN-γ production in stimulated blood from dogs with clinical leishmaniosis at diagnosis and during treatment. Vet Parasitol. 2017;248:39–47.

    Article  PubMed  Google Scholar 

  21. Marín-García PJ, Llobat L. Canine cytokines profile in an endemic region of L. infantum: related factors. Vet Sci. 2022. https://doi.org/10.3390/vetsci9060305.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Ordeix L, Montserrat-Sangrà S, Martínez-Orellana P, Solano-Gallego L. Toll-Like Receptors 2, 4, and 7, Interferon-gamma, interleukin 10, and programmed death ligand 1 transcripts in leishmanin skin test-positive reactions of Ibizan Hound dogs. J Immunol Res. 2020;2020:9602576.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Esteve LO. The spectrum of cutaneous manifestations in canine leishmaniosis: insights into diagnosis and immune responses. http://purl.org/dc/dcmitype/Text. Universitat Autonomade Barcelona. 2018

  24. Martínez-Orellana P, González N, Baldassarre A, Álvarez-Fernández A, Ordeix L, Paradies P, et al. Humoral responses and ex vivo IFN-γ production after canine whole blood stimulation with Leishmania infantum antigen or KMP11 recombinant protein. Vet Sci. 2022;9:116.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Transcript: ENSCAFT00845024391.1 (IFNG-201) - Variants - Canis_lupus_familiaris - Ensembl genome browser 107. http://www.ensembl.org/Canis_lupus_familiaris/Transcript/ProtVariations?db=core;g=ENSCAFG00845013680;r=10:10530421-10535258;t=ENSCAFT00845024391;vdb=variation;vf=10918. Accessed on 13 Jul 2022.

  26. da Silva GAV, Mesquita TG, Souza VC, Junior J, do ES, Gomes de Souza ML, Talhari AC, et al. A Single haplotype of IFNG correlating with low circulating levels of interferon-γ is associated with susceptibility to cutaneous leishmaniasis caused by Leishmania guyanensis. Clin Infect Dis. 2020;71:274–81.

    Article  PubMed  Google Scholar 

  27. Maia C, Campino L. Biomarkers associated with Leishmania infantum exposure, infection, and disease in dogs. Front Cell Infect Microbiol. 2018;8:302.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Maia C, Campino L. Cytokine and phenotypic cell profiles of Leishmania infantum infection in the dog. J Trop Med. 2012;2012:541571.

    Article  PubMed  Google Scholar 

  29. Ordeix L, Montserrat-Sangrà S, Martínez-Orellana P, Baxarias M, Solano-Gallego L. Toll-like receptors 2, 4 and 7, interferon-gamma and interleukin 10, and programmed death ligand 1 transcripts in skin from dogs of different clinical stages of leishmaniosis. Parasit Vectors. 2019;12:575.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Carrillo E, Moreno J. Cytokine profiles in canine visceral leishmaniasis. Vet Immunol Immunopathol. 2009;128:67–70.

    Article  CAS  PubMed  Google Scholar 

  31. Aslan H, Oliveira F, Meneses C, Castrovinci P, Gomes R, Teixeira C, et al. New insights into the transmissibility of Leishmania infantum from dogs to sand flies: experimental vector-transmission reveals persistent parasite depots at bite sites. J Infect Dis. 2016;213:1752–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Martínez-Orellana P, Quirola-Amores P, Montserrat-Sangrà S, Ordeix L, Llull J, Álvarez-Fernández A, et al. The inflammatory cytokine effect of Pam3CSK4 TLR2 agonist alone or in combination with Leishmania infantum antigen on ex-vivo whole blood from sick and resistant dogs. Parasit Vectors. 2017;10:123.

    Article  PubMed  PubMed Central  Google Scholar 

  33. rs22707631 (SNP) - Explore this variant - Canis_lupus_familiaris - Ensembl genome browser 107. http://www.ensembl.org/Canis_lupus_familiaris/Variation/Explore?db=core;g=ENSCAFG00845020168;r=19:17869103-17874712;t=ENSCAFT00845035554;v=rs22707631;vdb=variation;vf=1375817. Accessed on 13 Jul 2022.

  34. de Araujo FJ, da Silva LDO, Mesquita TG, Pinheiro SK, Vital WdS, Chrusciak-Talhari A, et al. Polymorphisms in the TOLLIP gene influence susceptibility to cutaneous leishmaniasis caused by Leishmania guyanensis in the Amazonas State of Brazil. PLoS Negl Trop Dis. 2015. https://doi.org/10.1371/journal.pntd.0003875.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Okamura H, Kashiwamura S, Tsutsui H, Yoshimoto T, Nakanishi K. Regulation of interferon-gamma production by IL-12 and IL-18. Curr Opin Immunol. 1998;10:259–64.

    Article  CAS  PubMed  Google Scholar 

  36. Takeda K, Tsutsui H, Yoshimoto T, Adachi O, Yoshida N, Kishimoto T, et al. Defective NK cell activity and Th1 response in IL-18-deficient mice. Immunity. 1998;8:383–90.

    Article  CAS  PubMed  Google Scholar 

  37. Kumar D, Tiwary P, Chakravarty J, Sundar S. Association of interleukin-18 gene polymorphism with susceptibility to visceral leishmaniasis in endemic area of Bihar, an Indian population. ScientificWorldJournal. 2014;2014:852104.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Pinelli E, Killick-Kendrick R, Wagenaar J, Bernadina W, del Real G, Ruitenberg J. Cellular and humoral immune responses in dogs experimentally and naturally infected with Leishmania infantum. Infect Immun. 1994;62:229–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Chamizo C, Moreno J, Alvar J. Semi-quantitative analysis of cytokine expression in asymptomatic canine leishmaniasis. Vet Immunol Immunopathol. 2005;103:67–75.

    Article  CAS  PubMed  Google Scholar 

  40. Manna L, Reale S, Viola E, Vitale F, Foglia Manzillo V, Pavone LM, et al. Leishmania DNA load and cytokine expression levels in asymptomatic naturally infected dogs. Vet Parasitol. 2006;142:271–80.

    Article  CAS  PubMed  Google Scholar 

  41. Carrillo E, Ahmed S, Goldsmith-Pestana K, Nieto J, Osorio Y, Travi B, et al. Immunogenicity of the P-8 amastigote antigen in the experimental model of canine visceral leishmaniasis. Vaccine. 2007;25:1534–43.

    Article  CAS  PubMed  Google Scholar 

  42. Wardini AB, Pinto-da-Silva LH, Nadaes NR, Nascimento MT, Roatt BM, Reis AB, et al. Neutrophil properties in healthy and Leishmania infantum-naturally infected dogs. Sci Rep. 2019;9:6247.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Abbehusen MMC, Almeida VDA, Solcà MdS, Pereira LdS, Costa DJ, Gil-Santana L, et al. Clinical and immunopathological findings during long term follow-up in Leishmania infantum experimentally infected dogs. Sci Rep. 2017;7:15914.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Sanz CR, Miró G, Sevane N, Reyes-Palomares A, Dunner S. Modulation of Host Immune Response during Leishmania infantum natural infection: a whole-transcriptome analysis of the popliteal lymph nodes in dogs. Front Immunol. 2021;12:794627.

    Article  CAS  PubMed  Google Scholar 

  45. Hajilooi M, Abasi M, Bazmani A, Ahmadi A, Matini M, Solgi G, et al. Evaluation of interleukin-8 -251 t/a polymorphisms in visceral leishmaniasis. J Res Health Sci. 2015;15:59–61.

    PubMed  Google Scholar 

  46. de Lima VMF, Peiro JR, de Oliveira VR. IL-6 and TNF-alpha production during active canine visceral leishmaniasis. Vet Immunol Immunopathol. 2007;115:189–93.

    Article  PubMed  Google Scholar 

  47. Martínez-Orellana P, Montserrat-Sangrà S, Quirola-Amores P, González N, Solano-Gallego L. Cytokine effect of TLR3, TLR4, and TLR7 agonists alone or associated with Leishmania infantum antigen on blood from dogs. Biomed Res Int. 2018;2018:5693736.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Ordeix L, Silva JEDS, Llull J, Quirola P, Montserrat-Sangrà S, Martínez-Orellana P, et al. Histological and immunological description of the leishmanin skin test in Ibizan Hounds. J Comp Pathol. 2018;158:56–65.

    Article  CAS  PubMed  Google Scholar 

  49. Gajanayaka N, Dong SXM, Ali H, Iqbal S, Mookerjee A, Lawton DA, et al. TLR-4 agonist induces IFN-γ production selectively in proinflammatory human M1 macrophages through the PI3K-mTOR- and JNK-MAPK-activated p70S6K pathway. J Immunol. 2021;207:2310–24.

    Article  CAS  PubMed  Google Scholar 

  50. Srivastava S, Pandey SP, Jha MK, Chandel HS, Saha B. Leishmania expressed lipophosphoglycan interacts with toll-like receptor (TLR)-2 to decrease TLR-9 expression and reduce anti-leishmanial responses. Clin Exp Immunol. 2013;172:403–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Wang L, Yao ZQ, Moorman JP, Xu Y, Ning S. Gene expression profiling identifies IRF4-associated molecular signatures in hematological malignancies. PLoS ONE. 2014;9:e106788.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Sundararaj S, Casarotto MG. Molecular interactions of IRF4 in B cell development and malignancies. Biophys Rev. 2021;13:1219–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Jeudy S, Wardrop KE, Alessi A, Dominov JA. Bcl-2 inhibits the innate immune response during early pathogenesis of murine congenital muscular dystrophy. PLoS ONE. 2011;6:e22369.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Li Z, Zeppa JJ, Hancock MA, McCormick JK, Doherty TM, Hendy GN, et al. Staphylococcal superantigens use LAMA2 as a coreceptor to activate T cells. J Immunol. 2018;200:1471–9.

    Article  CAS  PubMed  Google Scholar 

  55. Wei Z, Du Q, Li P, Liu H, Xia M, Chen Y, et al. Death-associated protein kinase 1 (DAPK1) controls CD8+ T cell activation, trafficking, and antitumor activity. FASEB J. 2021;35:e21138.

    Article  CAS  PubMed  Google Scholar 

  56. Wei Z, Li P, He R, Liu H, Liu N, Xia Y, et al. DAPK1 (death associated protein kinase 1) mediates mTORC1 activation and antiviral activities in CD8+ T cells. Cell Mol Immunol. 2021;18:138–49.

    Article  CAS  PubMed  Google Scholar 

  57. Hwang I-Y, Park C, Kehrl JH. Impaired trafficking of Gnai2+/- and Gnai2-/- T lymphocytes: implications for T cell movement within lymph nodes. J Immunol. 2007;179:439–48.

    Article  CAS  PubMed  Google Scholar 

  58. Han S-B, Moratz C, Huang N-N, Kelsall B, Cho H, Shi C-S, et al. Rgs1 and Gnai2 regulate the entrance of B lymphocytes into lymph nodes and B cell motility within lymph node follicles. Immunity. 2005;22:343–54.

    Article  CAS  PubMed  Google Scholar 

  59. Ballotti R, Cheli Y, Bertolotto C. The complex relationship between MITF and the immune system: a melanoma immunotherapy (response) factor? Mol Cancer. 2020;19:170.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Arechiga AF, Habib T, He Y, Zhang X, Zhang Z-Y, Funk A, et al. Cutting edge: the PTPN22 allelic variant associated with autoimmunity impairs B cell signaling. J Immunol. 2009;182:3343–7.

    Article  CAS  PubMed  Google Scholar 

  61. Gjörloff-Wingren A, Saxena M, Williams S, Hammi D, Mustelin T. Characterization of TCR-induced receptor-proximal signaling events negatively regulated by the protein tyrosine phosphatase PEP. Eur J Immunol. 1999;29:3845–54.

    Article  PubMed  Google Scholar 

  62. Ramia E, Chiaravalli AM, Bou Nasser Eddine F, Tedeschi A, Sessa F, Accolla RS, et al. CIITA-related block of HLA class II expression, upregulation of HLA class I, and heterogeneous expression of immune checkpoints in hepatocarcinomas: implications for new therapeutic approaches. Oncoimmunology. 2019;8:1548243.

    Article  PubMed  Google Scholar 

  63. Adhikari A, Cobb B, Eddington S, Becerra N, Kohli P, Pond A, et al. IFN-γ and CIITA modulate IL-6 expression in skeletal muscle. Cytokine X. 2020;2:100023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Brzostek J, Gautam N, Zhao X, Chen EW, Mehta M, Tung DWH, et al. T cell receptor and cytokine signal integration in CD8+ T cells is mediated by the protein Themis. Nat Immunol. 2020;21:186–98.

    Article  CAS  PubMed  Google Scholar 

  65. Sanchez-Robert E, Altet L, Sanchez A, Francino O. Polymorphism of Slc11a1 (Nramp1) gene and canine leishmaniasis in a case-control study. J Hered. 2005;96:755–8.

    Article  CAS  PubMed  Google Scholar 

  66. Braliou GG, Kontou PI, Boleti H, Bagos PG. Susceptibility to leishmaniasis is affected by host SLC11A1 gene polymorphisms: a systematic review and meta-analysis. Parasitol Res. 2019;118:2329–42.

    Article  PubMed  Google Scholar 

  67. Ejghal R, Hida M, Idrissi ML, Hessni AE, Lemrani M. SLC11A1 polymorphisms and susceptibility to visceral leishmaniasis in Moroccan patients. Acta Trop. 2014;140:130–6.

    Article  CAS  PubMed  Google Scholar 

  68. Bueno R, Carvalho Neta AV, Xavier MN, Oliveira RG, Diniz SA, Melo MN, et al. cDNA sequencing and expression of Nramp1 (Slc11a1) in dogs phenotypically resistant or susceptible to visceral leishmaniasis. Vet Immunol Immunopathol. 2009;127:332–9.

    Article  CAS  PubMed  Google Scholar 

  69. Singh N, Gedda MR, Tiwari N, Singh SP, Bajpai S, Singh RK. Solute carrier protein family 11 member 1 (Slc11a1) activation efficiently inhibits Leishmania donovani survival in host macrophages. J Parasit Dis. 2017;41:671–7.

    Article  PubMed  Google Scholar 

  70. Bourdon M, Santulli P, Doridot L, Jeljeli M, Chêne C, Chouzenoux S, et al. Immune cells and Notch1 signaling appear to drive the epithelial to mesenchymal transition in the development of adenomyosis in mice. Mol Hum Reprod. 2021. https://doi.org/10.1093/molehr/gaab053.

    Article  PubMed  Google Scholar 

  71. Castro RC, Zambuzi FA, Fontanari C, de Morais FR, Bollela VR, Kunkel SL, et al. NOTCH1 and DLL4 are involved in the human tuberculosis progression and immune response activation. Tuberculosis (Edinb). 2020;124:101980.

    Article  CAS  PubMed  Google Scholar 

  72. Paauw ND, Lely AT, Joles JA, Franx A, Nikkels PG, Mokry M, et al. H3K27 acetylation and gene expression analysis reveals differences in placental chromatin activity in fetal growth restriction. Clin Epigenetics. 2018;10:85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Liu X, Wang Y, Lu H, Li J, Yan X, Xiao M, et al. Genome-wide analysis identifies NR4A1 as a key mediator of T cell dysfunction. Nature. 2019;567:525–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Xu T, Park S-S, Giaimo BD, Hall D, Ferrante F, Ho DM, et al. RBPJ/CBF1 interacts with L3MBTL3/MBT1 to promote repression of Notch signaling via histone demethylase KDM1A/LSD1. EMBO J. 2017;36:3232–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We sincerely thank the “Podenco Ibicenco Association” for allowing their data and samples of Ibizan Hound dogs. We are grateful to Darío Gabanes for the graphics, to Precision Medicine Unit of INCLIVA, and to the Veterinary Medicine Faculty of Universidad Cardenal Herrera-CEU.

Funding

This work was funded by the grants for emerging research groups of Generalitat Valenciana (GV/2021/104). The funding agency was not involved in the design of the study, data collection, analysis of data or manuscript writing.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, L.L.; methodology and formal analysis, L.A. and L.L.; data curation, P-J. M-G. and L.L.; writing—original draft preparation, L. A., P.J.M.-G., and L.L.; writing—review and editing, L.L.; supervision, L.L.; project administration, L.L.; funding acquisition, L.L. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Lola Llobat.

Ethics declarations

Ethics approval

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Animal Experimentation Ethics Committee of the Universidad Cardenal Herrera CEU (2020/VSC/PEA/0216). The owners of the dogs used in the study signed an informed consent for participation in the study.

Consent for publication

All authors agree with the The experiments involving animals were approved by the Animal Experimentation Ethics Committee of the Universidad Cardenal Herrera CEU with code 2020/VSC/PEA/0216.publication of manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Genomic variants found in non-cytokine genes with frequency of alternative allele minor than 0.7. Table shows the identify number (rsID), chromosome position, reference and alternative alleles, frequency of alternative variant in the analyzed dataset, and functional class of variant.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Álvarez, L., Marín-García, PJ. & Llobat, L. Immunological and genomic characterization of Ibizan Hound dogs in an endemic Leishmania infantum region. Parasites Vectors 15, 445 (2022). https://doi.org/10.1186/s13071-022-05504-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-022-05504-3

Keywords