Mapping the genes for susceptibility and response to Leishmania tropica in mouse

Autochthonous Leishmaniasis Caused by Leishmania tropica, Identified by Using Whole-Genome Sequencing, Sri Lanka

Cutaneous leishmaniasis is atypical in Sri Lanka because Leishmania donovani, which typically causes visceral disease, is the causative agent. The origins of recently described hybrids between L. donovani and other Leishmania spp. usually responsible for cutaneous leishmaniasis remain unknown. Other endemic dermotropic Leishmania spp. have not been reported in Sri Lanka. Genome analysis of 27 clinical isolates from Sri Lanka and 32 Old World Leishmania spp. strains found 8 patient isolates clustered with L. tropica and 19 with L. donovani. The L. tropica isolates from Sri Lanka shared markers with strain LtK26 reported decades ago in India, indicating they were not products of recent interspecies hybridization. Because L. tropica was isolated from patients with leishmaniasis in Sri Lanka, our findings indicate L. donovani is not the only cause of cutaneous leishmaniasis in Sri Lanka and potentially explains a haplotype that led to interspecies dermotropic L. donovani hybrids.

Functionally distinct regions of the locus Leishmania major response 15 control IgE or IFNγ level in addition to skin lesions

Leishmaniasis, a disease caused by parasites of Leishmania spp., endangers more than 1 billion people living in endemic countries and has three clinical forms: cutaneous, mucocutaneous, and visceral. Understanding of individual differences in susceptibility to infection and heterogeneity of its pathology is largely lacking. Different mouse strains show a broad and heterogeneous range of disease manifestations such as skin lesions, splenomegaly, hepatomegaly, and increased serum levels of immunoglobulin E and several cytokines. Genome-wide mapping of these strain differences detected more than 30 quantitative trait loci (QTLs) that control the response to Leishmania major. Some control different combinations of disease manifestations, but the nature of this heterogeneity is not yet clear. In this study, we analyzed the L. major response locus Lmr15 originally mapped in the strain CcS-9 which carries 12.5% of the genome of the resistant strain STS on the genetic background of the susceptible strain BALB/c. For this analysis, we used the advanced intercross line K3FV between the strains BALB/c and STS. We confirmed the previously detected loci Lmr15, Lmr18, Lmr24, and Lmr27 and performed genetic dissection of the effects of Lmr15 on chromosome 11. We prepared the interval-specific recombinant strains 6232HS1 and 6229FUD, carrying two STS-derived segments comprising the peak linkage of Lmr15 whose lengths were 6.32 and 17.4 Mbp, respectively, and analyzed their response to L. major infection. These experiments revealed at least two linked but functionally distinct chromosomal regions controlling IFNγ response and IgE response, respectively, in addition to the control of skin lesions. Bioinformatics and expression analysis identified the potential candidate gene Top3a. This finding further clarifies the genetic organization of factors relevant to understanding the differences in the individual risk of disease.

Role of host genetics and cytokines in Leishmania infection

Cytokines and chemokines are important regulators of innate and specific responses in leishmaniasis, a disease that currently affects 12 million people. We overviewed the current information about influences of genetically engineered mouse models of cytokine and chemokine on leishmaniasis. We found that genetic background of the host, parasite species and sub-strain, as well as experimental design often modify effects of genetically engineered cytokine genes. Next we analyzed genes and QTLs (quantitative trait loci) that control response to Leishmania species in mouse in order to establish relationship between genetic control of cytokine expression and organ pathology. These studies revealed a network-like complexity of the combined effects of the multiple functionally diverse QTLs and their individual specificity. Genetic control of organ pathology and systemic immune response overlap only partially. Some QTLs control both organ pathology and systemic immune response, but the effects of genes and loci with the strongest impact on disease are cytokine-independent, whereas several loci modify cytokines levels in serum without influencing organ pathology. Understanding this genetic control might be important in development of vaccines designed to stimulate certain cytokine spectrum.

Simultaneous SNP selection and adjustment for population structure in high dimensional prediction models

Complex traits are known to be influenced by a combination of environmental factors and rare and common genetic variants. However, detection of such multivariate associations can be compromised by low statistical power and confounding by population structure. Linear mixed effects models (LMM) can account for correlations due to relatedness but have not been applicable in high-dimensional (HD) settings where the number of fixed effect predictors greatly exceeds the number of samples. False positives or false negatives can result from two-stage approaches, where the residuals estimated from a null model adjusted for the subjects' relationship structure are subsequently used as the response in a standard penalized regression model. To overcome these challenges, we develop a general penalized LMM with a single random effect called ggmix for simultaneous SNP selection and adjustment for population structure in high dimensional prediction models. We develop a blockwise coordinate descent algorithm with automatic tuning parameter selection which is highly scalable, computationally efficient and has theoretical guarantees of convergence. Through simulations and three real data examples, we show that ggmix leads to more parsimonious models compared to the two-stage approach or principal component adjustment with better prediction accuracy. Our method performs well even in the presence of highly correlated markers, and when the causal SNPs are included in the kinship matrix. ggmix can be used to construct polygenic risk scores and select instrumental variables in Mendelian randomization studies. Our algorithms are available in an R package available on CRAN (https://cran.r-project.org/package=ggmix).

Novel Loci Controlling Parasite Load in Organs of Mice Infected With Leishmania major, Their Interactions and Sex Influence

Leishmaniasis is a serious health problem in many countries, and continues expanding to new geographic areas including Europe and USA. This disease, caused by parasites of Leishmania spp. and transmitted by phlebotomine sand flies, causes up to 1.3 million new cases each year and despite efforts toward its functional dissection and treatment it causes 20-50 thousands deaths annually. Dependence of susceptibility to leishmaniasis on sex and host's genes was observed in humans and in mouse models. Several laboratories defined in mice a number of Lmr (Leishmania major response) genetic loci that control functional and pathological components of the response to and outcome of L. major infection. However, the development of its most aggressive form, visceral leishmaniasis, which is lethal if untreated, is not yet understood. Visceral leishmaniasis is caused by infection and inflammation of internal organs. Therefore, we analyzed the genetics of parasite load, spread to internal organs, and ensuing visceral pathology. Using a new PCR-based method of quantification of parasites in tissues we describe a network-like set of interacting genetic loci that control parasite load in different organs. Quantification of Leishmania parasites in lymph nodes, spleen and liver from infected F2 hybrids between BALB/c and recombinant congenic strains CcS-9 and CcS-16 allowed us to map two novel parasite load controlling Leishmania major response loci, Lmr24 and Lmr27. We also detected parasite-controlling role of the previously described loci Lmr4, Lmr11, Lmr13, Lmr14, Lmr15, and Lmr25, and describe 8 genetic interactions between them. Lmr14, Lmr15, Lmr25, and Lmr27 controlled parasite load in liver and lymph nodes. In addition, Leishmania burden in lymph nodes but not liver was influenced by Lmr4 and Lmr24. In spleen, parasite load was controlled by Lmr11 and Lmr13. We detected a strong effect of sex on some of these genes. We also mapped additional genes controlling splenomegaly and hepatomegaly. This resulted in a systematized insight into genetic control of spread and load of Leishmania parasites and visceral pathology in the mammalian organism.

Transforming Growth Factor Beta (TGFβ1) and Epidermal Growth Factor (EGF) as Biomarkers of Leishmania (V) braziliensis Infection and Early Therapeutic Response in Cutaneous Leishmaniasis: Studies in Hamsters

Introduction: In cutaneous leishmaniasis, the host immune response is responsible for the development of skin injuries but also for resolution of the disease especially after antileishmanial therapy. The immune factors that participate in the regulation of inflammation, remodeling of the extracellular matrix, cell proliferation and differentiation may constitute biomarkers of diseases or response to treatment. In this work, we analyzed the production of the growth factors EGF, TGFβ1, PDGF, and FGF during the infection by Leishmania parasites, the development of the injuries and the early response to treatment. Methodology: Golden hamsters were infected with L. (V) braziliensis. The growth factors were detected in skin scrapings and biopsies every 2 weeks after infected and then at day 7 of treatment with different drug candidates by RT-qPCR. The parasitic load was also quantified by RT-qPCR in skin biopsies sampled at the end of the study. Results: The infection by L. (V) braziliensis induced the expression of all the growth factors at day 15 of infection. One month after infection, EGF and TGFβ1 were expressed in all hamsters with inverse ratio. While the EGF and FGF levels decreased between day 15 and 30 of infection, the TGFβ1 increased and the PGDF levels did not change. The relative expression of EGF and TGFβ1 increased notably after treatment. However, the increase of EGF was associated with clinical cure while the increase of TGFβ1 was associated with failure to treatment. The amount of parasites in the cutaneous lesion at the end of the study decreased according to the clinical outcome, being lower in the group of cured hamsters and higher in the group of hamsters that had a failure to the treatment. Conclusions: A differential profile of growth factor expression occurred during the infection and response to treatment. Higher induction of TGFβ1 was associated with active disease while the higher levels of EGF are associated with adequate response to treatment. The inversely EGF/TGFβ1 ratio may be an effective biomarker to identify establishment of Leishmania infection and early therapeutic response, respectively. However, further studies are needed to validate the utility of the proposed biomarkers in field conditions.

A novel locus on mouse chromosome 7 that influences survival after infection with tick-borne encephalitis virus

Background

Tick-borne encephalitis (TBE) is the main tick-borne viral infection in Eurasia. Its manifestations range from inapparent infections and fevers with complete recovery to debilitating or fatal encephalitis. The basis of this heterogeneity is largely unknown, but part of this variation is likely due to host genetic. We have previously found that BALB/c mice exhibit intermediate susceptibility to the infection of TBE virus (TBEV), STS mice are highly resistant, whereas the recombinant congenic strain CcS-11, carrying 12.5% of the STS genome on the background of the BALB/c genome is even more susceptible than BALB/c. Importantly, mouse orthologs of human TBE controlling genes Oas1b, Cd209, Tlr3, Ccr5, Ifnl3 and Il10, are in CcS-11 localized on segments derived from the strain BALB/c, so they are identical in BALB/c and CcS-11. As they cannot be responsible for the phenotypic difference of the two strains, we searched for the responsible STS-derived gene-locus. Of course the STS-derived genes in CcS-11 may operate through regulating or epigenetically modifying these non-polymorphic genes of BALB/c origin.

Methods

To determine the location of the STS genes responsible for susceptibility of CcS-11, we analyzed survival of TBEV-infected F₂ hybrids between BALB/c and CcS-11. CcS-11 carries STS-derived segments on eight chromosomes. These were genotyped in the F₂ hybrid mice and their linkage with survival was tested by binary trait interval mapping. We have sequenced genomes of BALB/c and STS using next generation sequencing and performed bioinformatics analysis of the chromosomal segment exhibiting linkage with TBEV survival.

Results

Linkage analysis revealed a novel suggestive survival-controlling locus on chromosome 7 linked to marker D7Nds5 (44.2 Mb). Analysis of this locus for polymorphisms between BALB/c and STS that change RNA stability and genes’ functions led to detection of 9 potential candidate genes: Cd33, Klk1b22, Siglece, Klk1b16, Fut2, Grwd1, Abcc6, Otog, and Mkrn3. One of them, Cd33, carried a nonsense mutation in the STS strain.

Conclusions

The robust genetic system of recombinant congenic strains of mice enabled detection of a novel suggestive locus on chromosome 7. This locus contains 9 candidate genes, which will be focus of future studies not only in mice but also in humans.

Genetic Regulation of Guanylate-Binding Proteins 2b and 5 during Leishmaniasis in Mice

Interferon-induced GTPases [guanylate-binding proteins (GBPs)] play an important role in inflammasome activation and mediate innate resistance to many intracellular pathogens, but little is known about their role in leishmaniasis. We therefore studied expression of Gbp2b/Gbp1 and Gbp5 mRNA in skin, inguinal lymph nodes, spleen, and liver after Leishmania major infection and in uninfected controls. We used two different groups of related mouse strains: BALB/c, STS, and CcS-5, CcS-16, and CcS-20 that carry different combinations of BALB/c and STS genomes, and strains O20, C57BL/10 (B10) and B10.O20, OcB-9, and OcB-43 carrying different combinations of O20 and B10 genomes. The strains were classified on the basis of size and number of infection-induced skin lesions as highly susceptible (BALB/c, CcS-16), susceptible (B10.O20), intermediate (CcS-20), and resistant (STS, O20, B10, OcB-9, OcB-43). Some uninfected strains differed in expression of Gbp2b/Gbp1 and Gbp5, especially of Gbp2b/Gbp1 in skin. Uninfected BALB/c and STS did not differ in their expression, but in CcS-5, CcS-16, and CcS-20, which all carry BALB/c-derived Gbp gene-cluster, expression of Gbp2b/Gbp1 exceeds that of both parents. These data indicate trans-regulation of Gbps. Infection resulted in approximately 10× upregulation of Gbp2b/Gbp1 and Gbp5 mRNAs in organs of both susceptible and resistant strains, which was most pronounced in skin. CcS-20 expressed higher level of Gbp2b/Gbp1 than both parental strains in skin, whereas CcS-16 expressed higher level of Gbp2b/Gbp1 than both parental strains in skin and liver. This indicates a trans-regulation present in infected mice CcS-16 and CcS-20. Immunostaining of skin of five strains revealed in resistant and intermediate strains STS, CcS-5, O20, and CcS-20 tight co-localization of Gbp2b/Gbp1 protein with most L. major parasites, whereas in the highly susceptible strain, BALB/c most parasites did not associate with Gbp2b/Gbp1. In conclusion, expression of Gbp2b/Gbp1 and Gbp5 was increased even in organs of clinically asymptomatic resistant mice. It suggests a hidden inflammation, which might contribute to control of persisting parasites. This is supported by the co-localization of Gbpb2/Gbp1 protein and L. major parasites in skin of resistant and intermediate but not highly susceptible mice.

Canine susceptibility to visceral leishmaniasis: A systematic review upon genetic aspects, considering breed factors and immunological concepts

Dogs have different susceptibility degrees to leishmaniasis; however, genetic research on this theme is scarce, manly on visceral form. The aims of this systematic review were to describe and discuss the existing scientific findings on genetic susceptibility to canine leishmaniasis, as well as to show the gaps of the existing knowledge. Twelve articles were selected, including breed immunological studies, genome wide associations or other gene polymorphism or gene sequencing studies, and transcription approaches. As main results of literature, there was a suggestion of genetic clinical resistance background for Ibizan Hound dogs, and alleles associated with protection or susceptibility to visceral leishmaniasis in Boxer dogs. Genetic markers can explain phenotypic variance in both pro- and anti-inflammatory cytokines and in cellular immune responses, including antigen presentation. Many gene segments are involved in canine visceral leishmaniasis phenotype, with Natural Resistance Associated Macrophage Protein 1 (NRAMP1) as the most studied. This was related to both protection and susceptibility. In comparison with murine and human genetic approaches, lack of knowledge in dogs is notorious, with many possibilities for new studies, revealing a wide field to be assessed on canine leishmaniasis susceptibility research.

Gene-specific sex effects on eosinophil infiltration in leishmaniasis

Background

Sex influences susceptibility to many infectious diseases, including some manifestations of leishmaniasis. The disease is caused by parasites that enter to the skin and can spread to the lymph nodes, spleen, liver, bone marrow, and sometimes lungs. Parasites induce host defenses including cell infiltration, leading to protective or ineffective inflammation. These responses are often influenced by host genotype and sex. We analyzed the role of sex in the impact of specific gene loci on eosinophil infiltration and its functional relevance.

Methods

We studied the genetic control of infiltration of eosinophils into the inguinal lymph nodes after 8 weeks of Leishmania major infection using mouse strains BALB/c, STS, and recombinant congenic strains CcS-1,-3,-4,-5,-7,-9,-11,-12,-15,-16,-18, and -20, each of which contains a different random set of 12.5% genes from the parental “donor” strain STS and 87.5% genes from the “background” strain BALB/c. Numbers of eosinophils were counted in hematoxylin-eosin-stained sections of the inguinal lymph nodes under a light microscope. Parasite load was determined using PCR-ELISA.

Results

The lymph nodes of resistant STS and susceptible BALB/c mice contained very low and intermediate numbers of eosinophils, respectively. Unexpectedly, eosinophil infiltration in strain CcS-9 exceeded that in BALB/c and STS and was higher in males than in females. We searched for genes controlling high eosinophil infiltration in CcS-9 mice by linkage analysis in F₂ hybrids between BALB/c and CcS-9 and detected four loci controlling eosinophil numbers. Lmr14 (chromosome 2) and Lmr25 (chromosome 5) operate independently from other genes (main effects). Lmr14 functions only in males, the effect of Lmr25 is sex independent. Lmr15 (chromosome 11) and Lmr26 (chromosome 9) operate in cooperation (non-additive interaction) with each other. This interaction was significant in males only, but sex-marker interaction was not significant. Eosinophil infiltration was positively correlated with parasite load in lymph nodes of F₂ hybrids in males, but not in females.

Conclusions

We demonstrated a strong influence of sex on numbers of eosinophils in the lymph nodes after L. major infection and present the first identification of sex-dependent autosomal loci controlling eosinophilic infiltration. The positive correlation between eosinophil infiltration and parasite load in males suggests that this sex-dependent eosinophilic infiltration reflects ineffective inflammation.

Mast cells have no impact on cutaneous leishmaniasis severity and related Th2 differentiation in resistant and susceptible mice

The genus leishmania comprises different protozoan parasites which are causative agents of muco-cutaneous and systemic, potentially lethal diseases. After infection with the species Leishmania major, resistant mice expand Th1 cells which stimulate macrophages for Leishmania destruction. In contrast, susceptible mice generate Th2 cells which deactivate macrophages, leading to systemic spread of the pathogens. Th-cell differentiation is determined within the first days, and Th2 cell differentiation requires IL-4, whereby the initial IL-4 source is often unknown. Mast cells are potential sources of IL-4, and hence their role in murine leishmaniasis has previously been studied in mast cell-deficient Kit mutant mice, although these mice display immunological phenotypes beyond mast cell deficiency. We therefore readdressed this question by infecting Kit-independent mast cell-deficient mice that are Th1 (C57BL/6 Cpa(Cre) ) or Th2 (BALB/c Cpa(Cre) ) prone with L. major. Using different parasite doses and intra- or subcutaneous infection routes, the results demonstrate no role of mast cells on lesion size development, parasite load, immune cell phenotypes expanding in draining lymph nodes, and cytokine production during murine cutaneous leishmaniasis. Thus, other cell types such as ILCs or T cells have to be considered as primary source of Th2-driving IL-4.

A Review: The Current In Vivo Models for the Discovery and Utility of New Anti-leishmanial Drugs Targeting Cutaneous Leishmaniasis

The current in vivo models for the utility and discovery of new potential anti-leishmanial drugs targeting Cutaneous Leishmaniasis (CL) differ vastly in their immunological responses to the disease and clinical presentation of symptoms. Animal models that show similarities to the human form of CL after infection with Leishmania should be more representative as to the effect of the parasite within a human. Thus, these models are used to evaluate the efficacy of new anti-leishmanial compounds before human clinical trials. Current animal models aim to investigate (i) host–parasite interactions, (ii) pathogenesis, (iii) biochemical changes/pathways, (iv) in vivo maintenance of parasites, and (v) clinical evaluation of drug candidates. This review focuses on the trends of infection observed between Leishmania parasites, the predictability of different strains, and the determination of parasite load. These factors were used to investigate the overall effectiveness of the current animal models. The main aim was to assess the efficacy and limitations of the various CL models and their potential for drug discovery and evaluation. In conclusion, we found that the following models are the most suitable for the assessment of anti-leishmanial drugs: L. major–C57BL/6 mice (or–vervet monkey, or–rhesus monkeys), L. tropica–CsS-16 mice, L. amazonensis–CBA mice, L. braziliensis–golden hamster (or–rhesus monkey). We also provide in-depth guidance for which models are not suitable for these investigations.

Multilocus microsatellite typing reveals a genetic relationship but, also, genetic differences between Indian strains of Leishmania tropica causing cutaneous leishmaniasis and those causing visceral leishmaniasis

BACKGROUND

Leishmaniases are divided into cutaneous (CL) and visceral leishmaniasis (VL). In the Old World, CL is caused by Leishmania (L.) major, L. tropica and L. aethiopica. L. tropica can also visceralize and cause VL. In India, the large epidemics of VL are caused by L. donovani and cases of CL are caused by L. major and L. tropica. However, strains of L. tropica have also been isolated from Indian cases of VL.This study was done to see if Indian strains of L. tropica isolated from human cases of CL are genetically identical to or different from Indian strains of L. tropica isolated from human cases of VL and to see if any genetic differences found correlated with clinical outcome presenting as either CL or VL.

METHODS

Multilocus microsatellite typing (MLMT), employing 12 independent genetic markers specific to L. tropica, was used to characterize and identify eight strains of L. tropica isolated from human cases of CL examined in clinics in Bikaner City, Rajasthan State, north-west India. Their microsatellite profiles were compared to those of 156 previously typed strains of L. tropica from various geographical locations that were isolated from human cases of CL and VL, hyraxes and sand fly vectors.

RESULTS

Bayesian, distance-based and factorial correspondence analyses revealed two confirmed populations: India/Asia and Israel/Palestine that subdivided, respectively, into two and three subpopulations. A third population, Africa/Galilee, as proposed by Bayesian analysis was not supported by the other applied methods. The strains of L. tropica from Bikaner isolated from human cases of CL fell into one of the subpopulations in the population India/Asia together with strains from other Asian foci. Indian strains isolated from human cases of VL fell into the same sub-population but were not genetically identical to the Bikaner strains of L. tropica.

CONCLUSIONS

It seems that the genetic diversity encountered between the two groups of Indian strains is mainly owing to their geographical origins rather than their different times of isolation. Also, the genetic differences seen between the dermatotropic and viscerotropic strains might be connected with the difference in pathogenicity.