Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to visceral leishmaniasis

IFN-λ3 is induced by Leishmania donovani and can inhibit parasite growth in cell line models but not in the mouse model, while it shows a significant association with leishmaniasis in humans

The intracellular protozoan parasite Leishmania donovani causes debilitating human diseases that involve visceral and dermal manifestations. Type 3 interferons (IFNs), also referred to as lambda IFNs (IFNL, IFN-L, or IFN-λ), are known to play protective roles against intracellular pathogens at the epithelial surfaces. Herein, we show that L. donovani induces IFN-λ3 in human as well as mouse cell line-derived macrophages. Interestingly, IFN-λ3 treatment significantly decreased parasite load in infected cells, mainly by increasing reactive oxygen species production. Microscopic examination showed that IFN-λ3 inhibited uptake but not replication, while the phagocytic ability of the cells was not affected. This was confirmed by experiments that showed that IFN-λ3 could decrease parasite load only when added to the medium at earlier time points, either during or soon after parasite uptake, but had no effect on parasite load when added at 24 h post-infection, suggesting that an early event during parasite uptake was targeted. Furthermore, the parasites could overcome the inhibitory effect of IFN-λ3, which was added at earlier time points, within 2-3 days post-infection. BALB/c mice treated with IFN-λ3 before infection led to a significant increase in expression of IL-4 and ARG1 post-infection in the spleen and liver, respectively, and to different pathological changes, especially in the liver, but not to changes in parasite load. Treatment with IFN-λ3 during infection did not decrease the parasite load in the spleen either. However, IFN-λ3 was significantly increased in the sera of visceral leishmaniasis patients, and the IFNL genetic variant rs12979860 was significantly associated with susceptibility to leishmaniasis.

The Association of Human Leucocyte Antigen (HLA) Class I and II Genes with Cutaneous and Visceral Leishmaniasis in Iranian Patients: A Preliminary Case-Control Study

Background

Leishmaniasis is currently considered a re-emerging or emerging infection based on the geographic region. The outcome of leishmaniasis vastly depends on Leishmania-host interaction. This preliminary study aimed to show the association of human leukocyte antigen (HLA) class I and II genes with healed and non-healed cutaneous leishmaniasis (CL), and symptomatic and asymptomatic visceral leishmaniasis (VL) compared with control groups in Iran.

Methods

Ninety-five people, including 31 patients versus 64 individuals in the control group, were enrolled. Among them, 20 patients had confirmed CL based on amastigote observation, 10 had improved CL and 10 non-healed CL. Eleven patients were suffering from confirmed VL based on direct agglutination test (Five asymptomatic and six symptomatic VL cases). Besides, they were residents in an endemic area of VL in the northwest of Iran. To select a control group, it was ensured that they had no history of leishmaniasis. Peripheral blood samples were collected from each patient. After DNA extraction, HLA typing was conducted using polymerase chain reaction - sequence-specific priming (PCR-SSP). Subsequently, data were statistically analyzed by SPSS.

Results

There was a statistical relationship between the presence of HLA-A26 and CL, healed CL and the existence of the B38 allele, C1 allele and symptomatic VL, as well as B1.4 allele and asymptomatic VL (P<0.05).

Conclusion

This primary finding indicates that several HLA genes have a potential role in the susceptibility of Iranian people to CL and VL.

Malaria and leishmaniasis: Updates on co-infection

Malaria and leishmaniasis are endemic parasitic diseases in tropical and subtropical countries. Although the overlap of these diseases in the same host is frequently described, co-infection remains a neglected issue in the medical and scientific community. The complex relationship of concomitant infections with Plasmodium spp. and Leishmania spp. is highlighted in studies of natural and experimental co-infections, showing how this "dual" infection can exacerbate or suppress an effective immune response to these protozoa. Thus, a Plasmodium infection preceding or following Leishmania infection can impact the clinical course, accurate diagnosis, and management of leishmaniasis, and vice versa. The concept that in nature we are affected by concomitant infections reinforces the need to address the theme and ensure its due importance. In this review we explore and describe the studies available in the literature on Plasmodium spp. and Leishmania spp. co-infection, the scenarios, and the factors that may influence the course of these diseases.

Human leukocyte antigen system associations in Malassezia-related skin diseases

BACKGROUND

The human leukocyte antigen system (HLA) is divided into two classes involved in antigen presentation: class I presenting intracellular antigens and class II presenting extracellular antigens. While susceptibility to infections is correlated with the HLA system, data on associations between HLA genotypes and Malassezia-related skin diseases (MRSD) are lacking. Thus, the objective of this study was to investigate associations between HLA alleles and MRSD.

MATERIALS AND METHODS

Participants in The Danish Blood Donor Study (2010-2018) provided questionnaire data on life style, anthropometric measures, and registry data on filled prescriptions. Genotyping was done using Illumina Infinium Global Screening Array, and HLA alleles were imputed using the HIBAG algorithm. Cases and controls were defined using filled prescriptions on topical ketoconazole 2% as a proxy of MRSD. Logistic regressions assessed associations between HLA alleles and MRSD adjusted for confounders and Bonferroni corrected for multiple tests.

RESULTS

A total of 9455 participants were considered MRSD cases and 24,144 participants as controls. We identified four risk alleles B*57:01, OR 1.19 (95% CI: 1.09-1.31), C*01:02, OR 1.19 (95% CI: 1.08-1.32), C*06:02, OR 1.14 (95% CI: 1.08-1.22), and DRB1*01:01, OR 1.10 (95% CI: 1.04-1.17), and two protective alleles, DQB1*02:01, OR 0.89 (95% CI: 0.85-0.94), and DRB1*03:01, OR 0.89 (95% CI: 0.85-0.94).

CONCLUSION

Five novel associations between HLA alleles and MRSD were identified in our cohort, and one previous association was confirmed. Future studies should assess the correlation between Malassezia antigens and antigen-binding properties of the associated HLA alleles.

Human genetic polymorphism and Leishmaniasis

Leishmaniasis is a disease of the subtropical and tropical spheres of the earth and has various clinical manifestations. The different form of leishmaniasis includes cutaneous leishmaniasis, mucocutaneous leishmaniasis, most lethal visceral leishmaniasis and PKDL form. These different forms depend on many factors such as parasite and vector species, geographical, environmental conditions and population ethnicity. Host genetic factors have been widely investigated for their role in developing the disease in various infections. There are several reports on associations or resistance between candidate gene polymorphisms and the risk and outcome of Leishmania infection. Polymorphism in genes involved in both innate and adaptive immune systems, as well as genes of metabolic processes contributes to disease manifestation. The wide availability and advancement of molecular techniques permits to exploration of hereditary factors related to leishmaniasis. Many candidate gene studies were conducted on family-based and population to identify novel biomarkers for understanding disease pathogenesis pathways and possible drug targets. This comprehensive review presents an update on various human genes polymorphism that influence the outcome of different forms of Leishmania infection in endemic regions of the world. Various electronic databases were searched systematically for relevant publications and thoroughly analyzed. Most of the candidate gene studies were found with discrepancies in findings. Genetic and functional studies with adequate power are needed to validate the contribution of host genes in susceptibility or resistance towards Leishmania infection and understanding pathogenesis.

Comparison of mixed model based approaches for correcting for population substructure with application to extreme phenotype sampling

Background

Mixed models are used to correct for confounding due to population stratification and hidden relatedness in genome-wide association studies. This class of models includes linear mixed models and generalized linear mixed models. Existing mixed model approaches to correct for population substructure have been previously investigated with both continuous and case-control response variables. However, they have not been investigated in the context of extreme phenotype sampling (EPS), where genetic covariates are only collected on samples having extreme response variable values. In this work, we compare the performance of existing binary trait mixed model approaches (GMMAT, LEAP and CARAT) on EPS data. Since linear mixed models are commonly used even with binary traits, we also evaluate the performance of a popular linear mixed model implementation (GEMMA).

Results

We used simulation studies to estimate the type I error rate and power of all approaches assuming a population with substructure. Our simulation results show that for a common candidate variant, both LEAP and GMMAT control the type I error rate while CARAT’s rate remains inflated. We applied all methods to a real dataset from a Québec, Canada, case-control study that is known to have population substructure. We observe similar type I error control with the analysis on the Québec dataset. For rare variants, the false positive rate remains inflated even after correction with mixed model approaches. For methods that control the type I error rate, the estimated power is comparable.

Conclusions

The methods compared in this study differ in their type I error control. Therefore, when data are from an EPS study, care should be taken to ensure that the models underlying the methodology are suitable to the sampling strategy and to the minor allele frequency of the candidate SNPs.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08297-y).

Precision Medicine in Control of Visceral Leishmaniasis Caused by L. donovani

Precision medicine and precision global health in visceral leishmaniasis (VL) have not yet been described and could take into account how all known determinants improve diagnostics and treatment for the individual patient. Precision public health would lead to the right intervention in each VL endemic population for control, based on relevant population-based data, vector exposures, reservoirs, socio-economic factors and other determinants. In anthroponotic VL caused by L. donovani, precision may currently be targeted to the regional level in nosogeographic entities that are defined by the interplay of the circulating parasite, the reservoir and the sand fly vector. From this 5 major priorities arise: diagnosis, treatment, PKDL, asymptomatic infection and transmission. These 5 priorities share the immune responses of infection with L. donovani as an important final common pathway, for which innovative new genomic and non-genomic tools in various disciplines have become available that provide new insights in clinical management and in control. From this, further precision may be defined for groups (e.g. children, women, pregnancy, HIV-VL co-infection), and eventually targeted to the individual level.

Family-Based Genome-Wide Association Study of Autism Spectrum Disorder in Middle Eastern Families

Autism spectrum disorder (ASD) is a neurodevelopmental disease characterized by abnormalities in language and social communication with substantial clinical heterogeneity. Genetic factors play an important role in ASD with heritability estimated between 70% to 80%. Genome-wide association studies (GWAS) have identified multiple loci associated with ASD. However, most studies were performed on European populations and little is known about the genetic architecture of ASD in Middle Eastern populations. Here, we report the first GWAS of ASD in the Middle eastern population of Qatar. We analyzed 171 families with ASD, using linear mixed models adjusting for relatedness and other confounders. Results showed that common single nucleotide polymorphisms (SNP) in seven loci are associated with ASD (p < 1 × 10-5). Although the identified loci did not reach genome-wide significance, many of the top associated SNPs are located within or near genes that have been implicated in ASD or related neurodevelopmental disorders. These include GORASP2, GABBR2, ANKS6, THSD4, ERCC6L, ARHGEF6, and HDAC8. Additionally, three of the top associated SNPs were significantly associated with gene expression. We also found evidence of association signals in two previously reported ASD-susceptibility loci (rs10099100 and rs4299400). Our results warrant further functional studies and replication to provide further insights into the genetic architecture of ASD.

Statistical Learning Methods Applicable to Genome-Wide Association Studies on Unbalanced Case-Control Disease Data

Despite the fact that imbalance between case and control groups is prevalent in genome-wide association studies (GWAS), it is often overlooked. This imbalance is getting more significant and urgent as the rapid growth of biobanks and electronic health records have enabled the collection of thousands of phenotypes from large cohorts, in particular for diseases with low prevalence. The unbalanced binary traits pose serious challenges to traditional statistical methods in terms of both genomic selection and disease prediction. For example, the well-established linear mixed models (LMM) yield inflated type I error rates in the presence of unbalanced case-control ratios. In this article, we review multiple statistical approaches that have been developed to overcome the inaccuracy caused by the unbalanced case-control ratio, with the advantages and limitations of each approach commented. In addition, we also explore the potential for applying several powerful and popular state-of-the-art machine-learning approaches, which have not been applied to the GWAS field yet. This review paves the way for better analysis and understanding of the unbalanced case-control disease data in GWAS.

HLA-DRB1 Alleles Associated with Lower Leishmaniasis Susceptibility Share Common Amino Acid Polymorphisms and Epitope Binding Repertoires

Susceptibility for leishmaniasis is largely dependent on host genetic and immune factors. Despite the previously described association of human leukocyte antigen (HLA) gene cluster variants as genetic susceptibility factors for leishmaniasis, little is known regarding the mechanisms that underpin these associations. To better understand this underlying functionality, we first collected all known leishmaniasis-associated HLA variants in a thorough literature review. Next, we aligned and compared the protection- and risk-associated HLA-DRB1 allele sequences. This identified several amino acid polymorphisms that distinguish protection- from risk-associated HLA-DRB1 alleles. Subsequently, T cell epitope binding predictions were carried out across these alleles to map the impact of these polymorphisms on the epitope binding repertoires. For these predictions, we used epitopes derived from entire proteomes of multiple Leishmania species. Epitopes binding to protection-associated HLA-DRB1 alleles shared common binding core motifs, mapping to the identified HLA-DRB1 amino acid polymorphisms. These results strongly suggest that HLA polymorphism, resulting in differential antigen presentation, affects the association between HLA and leishmaniasis disease development. Finally, we established a valuable open-access resource of putative epitopes. A set of 14 HLA-unrestricted strong-binding epitopes, conserved across species, was prioritized for further epitope discovery in the search for novel subunit-based vaccines.

Genetics, Transcriptomics and Meta-Taxonomics in Visceral Leishmaniasis

Visceral leishmaniasis (VL) caused by parasites of the Leishmania donovani complex can be fatal in susceptible individuals. Understanding the interactions between host and pathogen is one way to obtain leads to develop better drugs and for vaccine development. In recent years multiple omics-based approaches have assisted researchers to gain a more global picture of this interaction in leishmaniasis. Here we review results from studies using three omics-based approaches to study VL caused by L. donovani in India: (i) chip-based analysis of single nucleotide variants in the first genome-wide association study of host genetic risk factors for VL, followed by analysis of epitope binding to HLA DRB1 risk versus protective alleles; (ii) transcriptional profiling demonstrating pathways important in Amphotericin B treated compared to active VL cases, including demonstration that anti-interleukin-10 unleashes a storm of chemokines and cytokines in whole blood responses to soluble leishmania antigen in active cases; and (iii) a meta-taxonomic approach based on sequencing amplicons derived from regions of 16S ribosomal RNA (16S rRNA) and 18S rRNA genes that allowed us to determine composition of both prokaryotic and eukaryotic gut microflora in VL cases compared to endemic controls. Overall, our omics-based approaches demonstrate that global analyses of genetic risk factors, host responses to infection, and the interaction between host, parasite and the microbiome can point to the most critical factors that determine the outcome of infection.

Post kala-azar dermal leishmaniasis: A threat to elimination program

Leishmaniasis remains a public health concern around the world that primarily affects poor folks of the developing world spanning across 98 countries with mortality of 0.2 million to 0.4 million annually. Post kala-azar dermal leishmaniasis (PKDL) is the late skin manifestation of visceral leishmaniasis (VL). It has been reported that about 2.5% to 20% of patients recovered from VL develop PKDL having stilted macular or nodular lesions with parasites. In the Indian subcontinent (ISC), it manifests a few months after recovery from VL, though in Africa it can occur simultaneously with VL or a little later. New cases of PKDL are also observed without prior VL in the ISC. These individuals with PKDL represent an important but largely neglected reservoir of infection that perpetuates anthroponotic Leishmania donovani transmission in the ISC and can jeopardize the VL elimination program as these cases can infect the sand flies and spread the endemic. Therefore, it becomes imperative to eradicate PKDL as a part of the VL elimination program. With the limited treatment options besides little knowledge on PKDL, this review stands out in focusing on different aspects that should be dealt for sustained VL elimination.

Fatal progression of experimental visceral leishmaniasis is associated with intestinal parasitism and secondary infection by commensal bacteria, and is delayed by antibiotic prophylaxis

Leishmania donovani causes visceral leishmaniasis (VL), which is typically fatal without treatment. There is substantial variation between individuals in rates of disease progression, response to treatment and incidence of post-treatment sequelae, specifically post-kala-azar dermal leishmaniasis (PKDL). Nevertheless, the majority of infected people are asymptomatic carriers. Hamsters and mice are commonly used as models of fatal and non-fatal VL, respectively. Host and parasite genetics are likely to be important factors, but in general the reasons for heterogeneous disease presentation in humans and animal models are poorly understood. Host microbiota has become established as a factor in cutaneous forms of leishmaniasis but this has not been studied in VL. We induced intestinal dysbiosis in mice and hamsters by long-term treatment with broad-spectrum antibiotics in their drinking water. There were no significant differences in disease presentation in dysbiotic mice. In contrast, dysbiotic hamsters infected with L. donovani had delayed onset and progression of weight loss. Half of control hamsters had a rapid progression phenotype compared with none of the ABX-treated animals and the nine-month survival rate was significantly improved compared to untreated controls (40% vs. 10%). Antibiotic-treated hamsters also had significantly less severe hepatosplenomegaly, which was accompanied by a distinct cytokine gene expression profile. The protective effect was not explained by differences in parasite loads or haematological profiles. We further found evidence that the gut-liver axis is a key aspect of fatal VL progression in hamsters, including intestinal parasitism, bacterial translocation to the liver, malakoplakia and iron sequestration, none of which occurred in non-progressing murine VL. Diverse bacterial genera were cultured from VL affected livers, of which Rodentibacter was specifically absent from ABX-treated hamsters, indicating this pathobiont may play a role in promoting disease progression. The results provide experimental support for antibiotic prophylaxis against secondary bacterial infections as an adjunct therapy in human VL patients.

Human genetics of leishmania infections

Identifying genetic risk factors for parasitic infections such as the leishmaniases could provide important leads for improved therapies and vaccines. Until recently most genetic studies of human leishmaniasis were underpowered and/or not replicated. Here, we focus on recent genome-wide association studies of visceral leishmaniasis (VL) and cutaneous leishmaniasis (CL). For VL, analysis across 2287 cases and 2692 controls from three cohorts identified a single major peak of genome-wide significance (P_combined = 2.76 × 10^–17) at HLA-DRB1–HLA-DQA1. HLA-DRB1*1501 and DRB1*1404/DRB1*1301 were the most significant protective versus risk alleles, respectively, with specific residues at amino acid positions 11 and 13 unique to protective alleles. Epitope-binding studies showed higher frequency of basic AAs in DRB1*1404-/*1301-specific epitopes compared to hydrophobic and polar AAs in DRB1*1501-specific epitopes at anchor residues P4 and P6 which interact with residues at DRB1 positions 11 and 13. For CL, genome-wide significance was not achieved in combined analysis of 2066 cases and 2046 controls across 2 cohorts. Rather, multiple top hits at P < 5 × 10^–5 were observed, amongst which IFNG-AS1 was of specific interest as a non-coding anti-sense RNA known to influence responses to pathogens by increasing IFN-γ secretion. Association at LAMP3 encoding dendritic cell lysosomal associated membrane protein 3 was also interesting. LAMP3 increases markedly upon activation of dendritic cells, localizing to the MHC Class II compartment immediately prior to translocation of Class II to the cell surface. Together these GWAS results provide firm confirmation for the importance of antigen presentation and the regulation of IFNγ in determining the outcome of Leishmania infections.

Leishmania amazonensis resistance in murine macrophages: Analysis of possible mechanisms

Leishmaniasis encompass a group of infectious parasitic diseases occurring in 97 endemic countries where over one billion people live in areas at risk of infection. It is in the World Health Organization list of neglected diseases and it is considered a serious public health problem, with more than 20,000 deaths a year and high morbidity. Infection by protozoa from the genus Leishmania can cause several forms of the disease, which may vary from a self-healing ulcer to fatal visceral infection. Leishmania species, as well as host immune response and genetics can modulate the course of the disease. Leishmania sp are obligatory intracellular parasites that have macrophages as their main host cell. Depending on the activation phenotype, these cells may have distinct roles in disease development, acting in parasite control or proliferation. Therefore, the purpose of this work was to analyze Leishmania amazonensis infection in primary macrophage cells obtained from mice with two distinct genetic backgrounds, ie. different susceptibility to the infection; evaluating the cause for that difference. After infection, peritoneal macrophages from the resistant C3H/He strain presented lower parasite load when compared to susceptible BALB/c macrophages. The same was also true when cells received a Th2 stimulus after infection, but the difference was abrogated under Th1 stimulus. Nitric oxide production and arginase activity was different between the strains under Th1 or Th2 stimulus, respectively, but iNOS inhibition was unable to suppress C3H/He resistance. Hydrogen peroxide production was also higher in C3H/He than BALB/c under Th1 stimulus, but it could not account for differences in susceptibility. These results led us to conclude that, although they have an important role in parasite control, neither NO nor H2O2 production can explain C3H/He resistance to infection. Other studies are needed to uncover different mechanisms of resistance/susceptibility to L. amazonensis.

HLA diversity in Brazil

Brazil is the fifth largest country in the world in area and the fifth most populous. The Brazilian voluntary Bone Marrow Donor Registry is the third largest in terms of number of donors in the world, being a valuable source of HLA genetics to characterize the donor population of Brazil as well. The genetic background of the Brazilian population is quite heterogeneous, resulting from 5 centuries of admixture among Native Americans, Europeans and Africans, making the Brazilian population unique in terms of genetic ancestry. The unique characteristics of populations in different Brazilian regions make them an exciting focus for genetic diversity studies. Studies on HLA genetic diversity of Brazilian populations have been conducted since the late 1980s and, in this review, we highlight the main findings from studies carried out in Brazil based on classical HLA. In addition, we calculated the genetic distance from the molecular data of the studies included in this review in order to have a broader view of the HLA diversity in Brazilian populations. We emphasize that characterization of HLA diversity is not only important for transplantation programs, but can shed a light on ancestry, history and other demographic patterns with or without association with autoimmune disease.

Genomic medicine in Chagas disease

Genetic approaches have been proposed for improving the understanding of the causes of differential susceptibility to Trypanosoma cruzi infection and Chagas disease outcome. Polymorphisms in genes involved in the immune/inflammatory response are being studied in order to clarify their possible role in the occurrence or severity of the cardiac and/or gastrointestinal complications. However still today, the number of significant associated genes is limited and the pathophysiological mechanisms underlying this condition are unknown. This article review the information currently available from the published scientific literature regarding the genetic variants of molecules of the immune system and other variants that can contribute to the clinical presentation of the disease. Genomic medicine will improve our knowledge about the molecular basis of Chagas disease, will open new avenues for developing biomarkers of disease progression, new therapeutic strategies to suit the requirements of individual patients, and will contribute to the control of one of the infections with the greatest socio-economic impact in the Americas.

Susceptibility to leishmaniasis is affected by host SLC11A1 gene polymorphisms: a systematic review and meta-analysis

Leishmaniases are cutaneous, mucocutaneous, and visceral diseases affecting humans and domesticated animals mostly in the tropical and subtropical areas of the planet. Host genetics have been widely investigated for their role in developing various infectious diseases. The SLC11A1 gene has been reported to play a role in neutrophil function and is associated with susceptibility to infectious and inflammatory diseases such as tuberculosis or rheumatoid arthritis. In the present meta-analysis, we investigate the genetic association of SLC11A1 polymorphisms with susceptibility to leishmaniasis. Genotypes and other risk-related data were collected from 13 case-control and family-based studies (after literature search). Conventional random-effects meta-analysis was performed using STATA 13. To pool case-control and family-based data, the weighted Stouffer's method was also applied. Eight polymorphisms were investigated: rs2276631, rs3731865, rs3731864, rs17221959, rs201565523, rs2279015, rs17235409, and rs17235416. We found that rs17235409 (D543N) and rs17235416 (1729 + 55del4) are significantly associated with a risk for cutaneous leishmaniasis (CL), whereas rs17221959, rs2279015, and rs17235409 are associated with visceral leishmaniasis (VL). Our results suggest that polymorphisms in SLC11A1 affect susceptibility to CL and VL. These findings open new pathways in understanding macrophage response to Leishmania infection and the genetic factors predisposing to symptomatic CL or VL that can lead to the usage of predictive biomarkers in populations at risk.

Determinants for progression from asymptomatic infection to symptomatic visceral leishmaniasis: A cohort study

BACKGROUND

Asymptomatic Leishmania donovani infections outnumber clinical presentations, however the predictors for development of active disease are not well known. We aimed to identify serological, immunological and genetic markers for progression from L. donovani infection to clinical Visceral Leishmaniasis (VL).

METHODS

We enrolled all residents >2 years of age in 27 VL endemic villages in Bihar (India). Blood samples collected on filter paper on two occasions 6-12 months apart, were tested for antibodies against L. donovani with rK39-ELISA and DAT. Sero converters, (negative for both tests in the first round but positive on either of the two during the second round) and controls (negative on both tests on both occasions) were followed for three years. At the start of follow-up venous blood was collected for the following tests: DAT, rK39- ELISA, Quantiferon assay, SNP/HLA genotyping and L.donovani specific quantitative PCR.

RESULTS

Among 1,606 subjects enrolled,17 (8/476 seroconverters and 9/1,130 controls) developed VL (OR 3.1; 95% CI 1.1-8.3). High DAT and rK39 ELISA antibody titers as well as positive qPCR were strongly and significantly associated with progression from seroconversion to VL with odds ratios of 19.1, 30.3 and 20.9 respectively. Most VL cases arose early (median 5 months) during follow-up.

CONCLUSION

We confirmed the strong association between high DAT and/or rK39 titers and progression to disease among asymptomatic subjects and identified qPCR as an additional predictor. Low predictive values do not warrant prophylactic treatment but as most progressed to VL early during follow-up, careful oberservation of these subjects for at least 6 months is indicated.

Asymptomatic Leishmania infection in blood donors from the Southern of Spain

OBJECTIVES

To investigate the proportion of asymptomatic infection among blood donors in a region endemic for Leishmania; and to ascertain epidemiological and genetic factors associated with this condition.

METHODS

We studied 1260 blood donors in the Province of Granada in the Southern Spain. After obtaining informed consent in each participant, a poll about habits, housing and contact with animals were carried out. Blood samples were obtained for determining antileishmanial antibodies and a PCR assay. HLA typing was performed in a randomly sample among the donors with positive serology.

RESULTS

We have found that L. infantum antibodies were present in 7.9% of blood donors and DNA in blood was detected in 2.5% of donors. There was no concordance between both determinations, except in one patient. Taking into consideration both techniques, 129 participants were considered to have asymptomatic Leishmania infection. No participant in this study developed clinical leishmaniasis during a follow-up period of 2 years. HLA were typed in 51 donors. Asymptomatic Leishmania infection might be associated with certain HLA antigens. A multivariate analysis was done with the variables obtained through the participants' interview. The contact with livestock (goats, pigs, and sheep), but not dogs, either at home or in the environment, was significantly and independently associated with asymptomatic leishmania infection.

CONCLUSIONS

Asymptomatic leishmanial infection among blood donors is frequent in the Granada Province, south of Spain. The presence of livestock in this region is related to this infection, perhaps influencing vector density of this disease. Some HLA genes might be associated with asymptomatic leishmanial state.

Immunogenicity and Protective Efficacy of T-Cell Epitopes Derived From Potential Th1 Stimulatory Proteins of Leishmania (Leishmania) donovani

Development of a suitable vaccine against visceral leishmaniasis (VL), a fatal parasitic disease, is considered to be vital for maintaining the success of kala-azar control programs. The fact that Leishmania-infected individuals generate life-long immunity offers a viable proposition in this direction. Our prior studies demonstrated that T-helper1 (Th1) type of cellular response was generated by six potential recombinant proteins viz. elongation factor-2 (elF-2), enolase, aldolase, triose phosphate isomerase (TPI), protein disulfide isomerase (PDI) and p45, derived from a soluble antigenic fraction (89.9–97.1 kDa) of Leishmania (Leishmania) donovani promastigote, in treated Leishmania patients and golden hamsters and showed significant prophylactic potential against experimental VL. Moreover, since, it is well-known that our immune system, in general, triggers production of specific protective immunity in response to a small number of amino acids (peptide), this led to the identification of antigenic epitopes of the above-stated proteins utilizing immunoinformatics. Out of thirty-six, three peptides-P-10 (enolase), P-14, and P-15 (TPI) elicited common significant lymphoproliferative as well as Th1-biased cytokine responses both in golden hamsters and human subjects. Further, immunization with these peptides plus BCG offered 75% prophylactic efficacy with boosted cellular immune response in golden hamsters against Leishmania challenge which is indicative of their candidature as potential vaccine candidates.

Variations in macrophage migration inhibitory factor gene are not associated with visceral leishmaniasis in India

BACKGROUND

The host genetic factors play important role in determining the outcome of visceral leishmaniasis (VL). Macrophage migration inhibitory factor (MIF) is an important host cytokine, which is a key regulator of innate immune system. Genetic variants in MIF gene have been found to be associated with several inflammatory and infectious diseases. Role of MIF is well documented in leishmaniasis diseases, including Indian visceral leishmaniasis, where elevated level of serum MIF has been associated with VL phenotypes. However, there was no genetic study to correlate MIF variants in VL, therefore, we aimed to study the possible association of three reported MIF gene variants -794 CATT, -173G > C and non-coding RNA gene LOC284889 in Indian VL phenotype.

METHODS

Study subjects comprised of 214 VL patients along with ethnically and demographically matched 220 controls from VL endemic regions of Bihar state in India.

RESULTS

We found no significant difference between cases and controls in allelic, genotypic and haplotype frequency of the markers analysed [-794 CATT repeats (χ2=0.86; p=0.35; OR=0.85; 95% CI=0.61-1.19); -173 G>C polymorphism (χ2=1.11; p=0.29; OR=0.83; 95% CI=0.59-1.16); and LOC284889 (χ2=0.78; p=0.37; OR=0.86; 95% CI=0.61-1.20)].

CONCLUSION

Since we did not find any significant differences between case and control groups, we conclude that sequencing of complete MIF gene and extensive study on innate and adaptive immunity genes may help in identifying genetic variations that are associated with VL susceptibility/resistance among Indians.

HLA-DR Class II expression on myeloid and lymphoid cells in relation to HLA-DRB1 as a genetic risk factor for visceral leishmaniasis

Genetic variation at HLA-DRB1 is a risk factor for visceral leishmaniasis (VL) caused by Leishmania donovani. The single nucleotide polymorphism rs9271252 upstream of the DRB1 gene provides a perfect tag for protective versus risk HLA-DRB1 four-digit alleles. In addition to the traditional role of the membrane-distal region of HLA class II molecules in antigen presentation and CD4 T-cell activation, the membrane-proximal region mediates 'non-traditional' multi-functional activation, differentiation, or death signals, including in DR-expressing T cells. To understand how HLA-DR contributes to disease pathogenesis, we examined expression at the protein level in circulating myeloid (CD14⁺ , CD16⁺ ) and lymphoid (CD4⁺ , CD8⁺ , CD19⁺ ) cells of VL patients (pre- and post-treatment) compared with endemic healthy controls (EHC). Although DR expression is reduced in circulating myeloid cells in active disease relative to EHC and post-treatment groups, expression is enhanced on CD4⁺  DR⁺ and CD8⁺  DR⁺ T cells consistent with T-cell activation. Cells of all myeloid and lymphoid populations from active cases were refractory to stimulation of DR expression with interferon-γ (IFN-γ). In contrast, all populations except CD19⁺ B cells from healthy blood bank controls showed enhanced DR expression following IFN-γ stimulation. The rs9271252 genotype did not impact significantly on IFN-γ-activated DR expression in myeloid, B or CD8⁺ T cells, but CD4⁺ T cells from healthy individuals homozygous for the risk allele were particularly refractory to activated DR expression. Further analysis of DR expression on subsets of CD4⁺ T cells regulating VL disease could uncover additional ways in which pleiotropy at HLA DRB1 contributes to disease pathogenesis.

Coping with family structure in genome-wide association studies: a comparative evaluation

In this paper, a fully statistical investigation of the control of family structure as random effects is analyzed and discussed, using both the genome-wide association studies (GWAS) data and simulated data. Three modeling strategies are proposed and the analysis results suggest the hybrid use of results from all possible models should be combined in practice.

An HLA class II locus, previously identified by a genome-wide association study, is also associated with susceptibility to pulmonary tuberculosis in a Chinese population

OBJECTIVE

Genome-wide association study (GWAS) in Icelanders identified HLA class II sequence variants on chromosome 6p21 as tuberculosis (TB) susceptibility loci. To evaluate the role of these loci in other populations with different ancestry, we conducted a case-control study in Chinese population.

METHODS

We genotyped two genetic variants (rs9272461 and rs9271300) on the reported chromosome 6p21 in 739 pulmonary tuberculosis (PTB) cases and 749 healthy controls from Chinese Han population using TaqMan allelic discrimination assay. Logistic regression was applied to evaluate the association between genetic variants and PTB risk and to estimate corresponding odds ratios (ORs) and 95% confidence intervals (95%CIs).

RESULTS

We found that rs9272461 was significantly associated with the risk of PTB in various genetic models (dominant OR = 0.75, 95%CI: 0.61-0.92; recessive OR = 0.64, 95%CI: 0.46-0.90, and additive OR = 0.78, 95%CI: 0.67-0.90). Moreover, in the stratified analysis in additive model, the association was also significant in the old (age ≥ 48 years) (OR = 0.76, 95%CI: 0.62-0.93; P = .008), men (OR = 0.71, 95%CI: 0.59-0.85; P < .001), and new PTB cases (OR = 0.76, 95%CI: 0.65-0.90; P = .001). The association results were similar between the microbiologically negative (OR = 0.78, 95%CI: 0.64-0.94; P = .008) and positive cases (OR = 0.77, 95%CI: 0.64-0.93; P = .008). We did not observe significant association for rs9271300 neither in the overall analysis (additive model: OR = 0.98, 95%CI: 0.85-1.13; P = .776) nor in the stratified analysis.

CONCLUSIONS

Our findings indicate that the HLA class II locus also affects the susceptibility to PTB in Chinese population. Further validation studies and function experiments are required to confirm the roles of the discovered variant.

Canine leishmaniasis: Genome-wide analysis and antibody response to Lutzomyia longipalpis saliva

The anti-inflammatory properties of sand fly saliva favor the establishment of the Leishmania infantum infection. In contrast, an antibody response against Lutzomyia longipalpis saliva is often associated with a protective cell-mediated response against canine visceral leishmaniasis. Genetic studies may demonstrate to what extent the ability to secrete anti-saliva antibodies depends on genetic or environmental factors. However, the genetic basis of canine antibody response against sand fly saliva has not been assessed. The aim of this study was to identify chromosomal regions associated with the anti-Lu. longipalpis salivary IgG response in 189 dogs resident in endemic areas in order to provide information for prophylactic strategies. Dogs were classified into five groups based on serological and parasitological diagnosis and clinical evaluation. Anti-salivary gland homogenate (SGH) IgG levels were assessed by Enzyme-Linked Immunosorbent Assay (ELISA). Genomic DNA was isolated from blood samples and genotyped using a SNP chip with 173,662 single nucleotide polymorphism (SNP) markers. The following linear regression model was fitted: IgG level = mean + origin + sex + age + use of a repellent collar, and the residuals were assumed as pseudo-phenotypes for the association test between phenotypes and genotypes (GWA). A component of variance model that takes into account polygenic and sample structure effects (EMMAX) was employed for GWA. Phenotypic findings indicated that anti-SGH IgG levels remained higher in exposed and subclinically infected dogs than in severely diseased dogs even in regression model residuals. Five associated markers were identified on chromosomes 2, 20 and 31. The mapped genes included CD180 (RP105) and MITF related to the rapid activation of B lymphocytes and differentiation into antibody-secreting plasma cells. The findings pointed to chromosomal segments useful for functional confirmation studies and a search for adjuvant molecules of the anti-saliva response.

Transformation of Summary Statistics from Linear Mixed Model Association on All-or-None Traits to Odds Ratio

Genome-wide association studies (GWAS) have identified thousands of loci that are robustly associated with complex diseases. The use of linear mixed model (LMM) methodology for GWAS is becoming more prevalent due to its ability to control for population structure and cryptic relatedness and to increase power. The odds ratio (OR) is a common measure of the association of a disease with an exposure (e.g., a genetic variant) and is readably available from logistic regression. However, when the LMM is applied to all-or-none traits it provides estimates of genetic effects on the observed 0-1 scale, a different scale to that in logistic regression. This limits the comparability of results across studies, for example in a meta-analysis, and makes the interpretation of the magnitude of an effect from an LMM GWAS difficult. In this study, we derived transformations from the genetic effects estimated under the LMM to the OR that only rely on summary statistics. To test the proposed transformations, we used real genotypes from two large, publicly available data sets to simulate all-or-none phenotypes for a set of scenarios that differ in underlying model, disease prevalence, and heritability. Furthermore, we applied these transformations to GWAS summary statistics for type 2 diabetes generated from 108,042 individuals in the UK Biobank. In both simulation and real-data application, we observed very high concordance between the transformed OR from the LMM and either the simulated truth or estimates from logistic regression. The transformations derived and validated in this study improve the comparability of results from prospective and already performed LMM GWAS on complex diseases by providing a reliable transformation to a common comparative scale for the genetic effects.

Epitope-Binding Characteristics for Risk versus Protective DRB1 Alleles for Visceral Leishmaniasis

HLA-DRB1 is the major genetic risk factor for visceral leishmaniasis (VL). We used SNP2HLA to impute HLA-DRB1 alleles and SNPTEST to carry out association analyses in 889 human cases and 977 controls from India. NetMHCIIpan 2.1 was used to map epitopes and binding affinities across 49 Leishmania vaccine candidates, as well as across peptide epitopes captured from dendritic cells treated with crude Leishmania Ag and identified using mass spectrometry and alignment to amino acid sequences of a reference Leishmania genome. Cytokines were measured in peptide-stimulated whole blood from 26 cured VL cases and eight endemic healthy controls. HLA-DRB1*1501 and DRB1*1404/DRB1*1301 were the most significant protective and risk alleles, respectively, with specific residues at aa positions 11 and 13 unique to protective alleles. We observed greater peptide promiscuity in sequence motifs for 9-mer core epitopes predicted to bind to risk (*1404/*1301) compared with protective (*1501) DRB1 alleles. There was a higher frequency of basic amino acids in DRB1*1404/*1301-specific epitopes compared with hydrophobic and polar amino acids in DRB1*1501-specific epitopes at anchor residues pocket 4 and pocket 6, which interact with residues at DRB1 positions 11 and 13. Cured VL patients made variable, but robust, IFN-γ, TNF, and IL-10 responses to 20-mer peptides based on captured epitopes, with peptides based on DRB1*1501-captured epitopes resulting in a higher proportion (odds ratio 2.23, 95% confidence interval 1.17-4.25, p = 0.017) of patients with IFN-γ/IL-10 ratios > 2-fold compared with peptides based on DRB1*1301-captured epitopes. Our data provide insight into the molecular mechanisms underpinning the association of HLA-DRB1 alleles with risk versus protection in VL in humans.

Detection of Leishmania donovani in peripheral blood of asymptomatic individuals in contact with patients with visceral leishmaniasis

BACKGROUND

The majority of individuals infected with Leishmania donovani complex remain asymptomatic. They may act as transmission reservoirs for visceral leishmaniasis (VL). We investigated sero-prevalence of L. donovani complex amongst those closely associated with patients with VL and whether these sero-reactive individuals had Leishmania parasites in their peripheral blood. Other risk factors were also investigated.

METHODS

A total of 257 individuals in contact with patients with VL were tested for anti-Leishmania antibodies by rK39 immunochromatographic test (rK39 ICT), ELISA using promastigote antigen (p-ELISA) and indirect fluorescent antibody test (IFAT). Buffy coats of rK39 ICT positive individuals were cultured; sero-reactive buffy coats were tested for Leishmania DNA by ITS1 PCR. DNA obtained from culture was sequenced to confirm Leishmania species. Risk factors were evaluated for each sero-positive sample.

RESULTS

The results showed 29.2% (75/257) prevalence by serological tests: 14.4% (37/257) were positive by rK39 ICT, 25.3% (65/257) by p-ELISA, 18.3% (47/257) by IFAT and 10.9% (28/257) by all three serological methods. Ten percent (3/30) of cultures were positive for Leishmania promastigotes. Only 3% (2/74) sero-reactive buffy coats were positive for DNA; sequence analysis revealed L. donovani species. Significant risk factors were age, working as farmers, domestic animals in household and proximity to animal shelters.

CONCLUSIONS

Asymptomatic family members of patients with VL can carry live L. donovani in peripheral blood and may act as potential reservoirs.

GENBANK ACCESSION NUMBER

BankIt1863680 Leishmania KT921417 (DNA sequences of the ribosomal ITS1 region of L. donovani).

Genome-wide linkage and association study implicates the 10q26 region as a major genetic contributor to primary nonsyndromic vesicoureteric reflux

Vesicoureteric reflux (VUR) is the commonest urological anomaly in children. Despite treatment improvements, associated renal lesions – congenital dysplasia, acquired scarring or both – are a common cause of childhood hypertension and renal failure. Primary VUR is familial, with transmission rate and sibling risk both approaching 50%, and appears highly genetically heterogeneous. It is often associated with other developmental anomalies of the urinary tract, emphasising its etiology as a disorder of urogenital tract development. We conducted a genome-wide linkage and association study in three European populations to search for loci predisposing to VUR. Family-based association analysis of 1098 parent-affected-child trios and case/control association analysis of 1147 cases and 3789 controls did not reveal any compelling associations, but parametric linkage analysis of 460 families (1062 affected individuals) under a dominant model identified a single region, on 10q26, that showed strong linkage (HLOD = 4.90; ZLRLOD = 4.39) to VUR. The ~9Mb region contains 69 genes, including some good biological candidates. Resequencing this region in selected individuals did not clearly implicate any gene but FOXI2, FANK1 and GLRX3 remain candidates for further investigation. This, the largest genetic study of VUR to date, highlights the 10q26 region as a major genetic contributor to VUR in European populations.

Combining epidemiology with basic biology of sand flies, parasites, and hosts to inform leishmaniasis transmission dynamics and control

Quantitation of the nonlinear heterogeneities in Leishmania parasites, sand fly vectors, and mammalian host relationships provides insights to better understand leishmanial transmission epidemiology towards improving its control. The parasite manipulates the sand fly via production of promastigote secretory gel (PSG), leading to the “blocked sand fly” phenotype, persistent feeding attempts, and feeding on multiple hosts. PSG is injected into the mammalian host with the parasite and promotes the establishment of infection. Animal models demonstrate that sand flies with the highest parasite loads and percent metacyclic promastigotes transmit more parasites with greater frequency, resulting in higher load infections that are more likely to be both symptomatic and efficient reservoirs. The existence of mammalian and sand fly “super-spreaders” provides a biological basis for the spatial and temporal clustering of clinical leishmanial disease. Sand fly blood-feeding behavior will determine the efficacies of indoor residual spraying, topical insecticides, and bed nets. Interventions need to have sufficient coverage to include transmission hot spots, especially in the absence of field tools to assess infectiousness. Interventions that reduce sand fly densities in the absence of elimination could have negative consequences, for example, by interfering with partial immunity conferred by exposure to sand fly saliva. A deeper understanding of both sand fly and host biology and behavior is essential to ensuring effectiveness of vector interventions.

We review recent research that sheds light on the quantitative biology of leishmanial transmission between sand flies and mammalian hosts and use these insights to better understand transmission, the observed epidemiology of the disease, and their implications in choice of control strategy. Using animal models, we show how the parasite-induced processes manipulate sand fly blood-feeding behavior and the infectious metacyclic dose to promote host infection and to differentially regulate the onward transmission potential of individual vectors and hosts. The existence of subpopulations of mammalian and sand fly “super-spreaders” provides a biological basis for the spatial and temporal clustering of clinical leishmanial disease. While tools are unavailable to distinguish these individuals in mixed populations, blanket interventions will be necessary to ensure inclusion of transmission hot spots. Interventions that reduce sand fly densities without elimination could interfere with vector—host dynamics and conferred partial immunity to host populations.

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.

Genetic susceptibility to infectious diseases: Current status and future perspectives from genome-wide approaches

Genome-wide association studies (GWASs) have been widely applied to identify genetic factors that affect complex diseases or traits. Presently, the GWAS Catalog includes > 2800 human studies. Of these, only a minority have investigated the susceptibility to infectious diseases or the response to therapies for the treatment or prevention of infections. Despite their limited application in the field, GWASs have provided valuable insights by pinpointing associations to both innate and adaptive immune response loci, as well as novel unexpected risk factors for infection susceptibility. Herein, we discuss some issues and caveats of GWASs for infectious diseases, we review the most recent findings ensuing from these studies, and we provide a brief summary of selected GWASs for infections in non-human mammals. We conclude that, although the general trend in the field of complex traits is to shift from GWAS to next-generation sequencing, important knowledge on infectious disease-related traits can be still gained by GWASs, especially for those conditions that have never been investigated using this approach. We suggest that future studies will benefit from the leveraging of information from the host's and pathogen's genomes, as well as from the exploration of models that incorporate heterogeneity across populations and phenotypes. Interactions within HLA genes or among HLA variants and polymorphisms located outside the major histocompatibility complex may also play an important role in shaping the susceptibility and response to invading pathogens.

•

Relatively few GWASs for infectious diseases were performed.

•

Phenotype heterogeneity and case/control misclassification can affect GWAS power.

•

Adaptive and innate immunity loci were identified in several infectious disease GWASs.

•

Unexpected loci (e.g., lncRNAs) were also associated with infection susceptibility.

•

GWASs should integrate host and pathogen diversity and use complex association models.

Exome Sequencing Identifies Two Variants of the Alkylglycerol Monooxygenase Gene as a Cause of Relapses in Visceral Leishmaniasis in Children, in Sudan

Background

Visceral leishmaniasis (kala-azar, KA) is the most severe form of leishmaniasis, characterized by fever, weight loss, hepatosplenomegaly, and lymphadenopathy. During an outbreak of KA in Babar El Fugara (Sudan), 5.7% of cured patients displayed relapses, with familial clustering in half the cases.

Methods

We performed whole-exome sequencing on 10 relapsing individuals and 11 controls from 5 nuclear families.

Results

Rare homozygous and compound-heterozygous nonsense (c.1213C > T, rs139309795, p.Arg405*) and missense (c.701A > G, rs143439626, p.Lys234Arg) mutations of the alkylglycerol monooxygenase (AGMO) gene were associated with KA relapse in 3 families. Sequencing in additional family members confirmed the segregation of these mutations with relapse and revealed an autosomal dominant mode of transmission. These mutations were detected heterozygous in 2 subjects among 100 unrelated individuals with KA who never relapsed after cure, suggesting incomplete penetrance of AGMO deficiency. AGMO is expressed in hematopoietic cells, and is strongly expressed in the liver. AGMO modulates PAF production by mouse macrophages, suggesting that it may act through the PAF/PAF receptor pathway previously shown to have anti-Leishmania activity.

Conclusions

This is the first demonstration that relapses after a first episode of KA are due to differences in human genetic susceptibility and not to modifications of parasite pathogenicity.

Expression profiling of Sudanese visceral leishmaniasis patients pre- and post-treatment with sodium stibogluconate

Visceral leishmaniasis (VL) in Sudan caused by Leishmania donovani is fatal in susceptible individuals if untreated. Treatment with sodium stibogluconate (SSG) leads to post-kala-azar dermal leishmaniasis (PKDL) in 58% of patients. Here, Affymetrix microarrays were used to identify genes differentially expressed in lymph nodes (N=9 paired samples) pre- and post-treatment with SSG. Using the Bioconductor package limma, 438 genes from 28 869 post-quality-control probe sets were differentially expressed (P_nominal ≤.02) post- vs pretreatment. Canonical pathway analysis using Ingenuity Pathway Analysis™ identified "role of nuclear factor of activated T-cell in regulation of immune response" (P_nominal =1.35×10^-5 ; P_BH-adjusted =4.79×10^-3 ), "B-cell development" (P_nominal =2.04×10^-4 ; P_BH-adjusted =.024), "Fcγ receptor-mediated phagocytosis in macrophages and monocytes" (P_nominal =2.04×10^-4 ; P_BH-adjusted =.024) and "OX40 signalling" (P_nominal =2.82×10^-4 ; P_BH-adjusted =.025) as pathways differentially regulated post- vs pretreatment. Major network hub genes included TP53, FN1, MYC, BCL2, JUN, SYK, RUNX2, MMP1 and ACTA2. Top endogenous upstream regulators included IL-7 (P=2.28×10^-6 ), TNF (P=4.26×10^-6 ), Amyloid Precursor Protein (P=4.23×10^-5 ) and SPI1/PI.1 (P=1.17×10^-7 ). Top predicted chemical drug regulators included the flavonoid genistein (P=4.56×10^-7 ) and the quinoline alkaloid camptothecin (P=5.14×10^-5 ). These results contribute to our understanding of immunopathology associated with VL and response to SSG treatment. Further replication could identify novel therapeutic strategies that improve on SSG treatment and reduce the likelihood of progression to PKDL.

The MHC locus and genetic susceptibility to autoimmune and infectious diseases

In the past 50 years, variants in the major histocompatibility complex (MHC) locus, also known as the human leukocyte antigen (HLA), have been reported as major risk factors for complex diseases. Recent advances, including large genetic screens, imputation, and analyses of non-additive and epistatic effects, have contributed to a better understanding of the shared and specific roles of MHC variants in different diseases. We review these advances and discuss the relationships between MHC variants involved in autoimmune and infectious diseases. Further work in this area will help to distinguish between alternative hypotheses for the role of pathogens in autoimmune disease development.

The MHC locus and genetic susceptibility to autoimmune and infectious diseases

In the past 50 years, variants in the major histocompatibility complex (MHC) locus, also known as the human leukocyte antigen (HLA), have been reported as major risk factors for complex diseases. Recent advances, including large genetic screens, imputation, and analyses of non-additive and epistatic effects, have contributed to a better understanding of the shared and specific roles of MHC variants in different diseases. We review these advances and discuss the relationships between MHC variants involved in autoimmune and infectious diseases. Further work in this area will help to distinguish between alternative hypotheses for the role of pathogens in autoimmune disease development.

Comprehensive candidate gene analysis for symptomatic or asymptomatic outcomes of Leishmania infantum infection in Brazil

Genetic risk factors contribute to asymptomatic versus symptomatic visceral leishmaniasis (VL) outcomes following infection with Leishmania infantum. We therefore carried out a family-based (n = 918 post-quality control fully genotyped and phenotyped individuals) candidate gene study for symptomatic VL or asymptomatic delayed-type hypersensitivity (DTH) skin test phenotypes in highly endemic neighborhoods of northeast Brazil. A total of 248 SNPs were genotyped in 42 genes selected as candidates on the basis of prior genetic, immunological, and transcriptional profiling studies. The most significant association with the VL phenotype was with SNP rs6785358 (P = 5.7e-04; p_corrected = 0.026) 3.8 kb upstream of TGFBR2, the gene encoding the type 2 receptor for transforming growth factor beta (TGFβ). A second inhibitory member of the TGBβ superfamily signaling pathway, SMAD7, was associated with the DTH phenotype (SNP rs7238442: P = 0.001; p_corrected = 0.051). The most significant association for the DTH phenotype was with SNP rs10800309 (P = -8.4e-06; p_corrected = 3.9e-04) situated 3.1 kb upstream of FCGR2A, the gene encoding the low-affinity IIa receptor for the Fc fragment of IgG. Overall, our results imply a role for IgG-mediated inflammation in determining DTH associated with asymptomatic infection and contribute to growing evidence that the TGFβ pathway is important in the immunopathogenesis of VL.

Blood Transcriptional Profiling Reveals Immunological Signatures of Distinct States of Infection of Humans with Leishmania infantum

Visceral leishmaniasis (VL) can be lethal if untreated; however, the majority of human infections with the etiological agents are asymptomatic. Using Illumina Bead Chip microarray technology, we investigated the patterns of gene expression in blood of active VL patients, asymptomatic infected individuals, patients under remission of VL and controls. Computational analyses based on differential gene expression, gene set enrichment, weighted gene co-expression networks and cell deconvolution generated data demonstrating discriminative transcriptional signatures. VL patients exhibited transcriptional profiles associated with pathways and gene modules reflecting activation of T lymphocytes via MHC class I and type I interferon signaling, as well as an overall down regulation of pathways and gene modules related to myeloid cells, mainly due to differences in the relative proportions of monocytes and neutrophils. Patients under remission of VL presented heterogeneous transcriptional profiles associated with activation of T lymphocytes via MHC class I, type I interferon signaling and cell cycle and, importantly, transcriptional activity correlated with activation of Notch signaling pathway and gene modules that reflected increased proportions of B cells after treatment of disease. Asymptomatic and uninfected individuals presented similar gene expression profiles, nevertheless, asymptomatic individuals exhibited particularities which suggest an efficient regulation of lymphocyte activation and a strong association with a type I interferon response. Of note, we validated a set of target genes by RT-qPCR and demonstrate the robustness of expression data acquired by microarray analysis. In conclusion, this study profiles the immune response during distinct states of infection of humans with Leishmania infantum with a novel strategy that indicates the molecular pathways that contribute to the progression of the disease, while also providing insights into transcriptional activity that can drive protective mechanisms.

Infections of humans with the protozoan parasites L. donvani and L. infantum can lead to the development of the disease visceral leishmaniasis, but also to an asymptomatic status. However, the mechanisms that result in these clinical outcomes after infection are poorly understood. In this study, we applied a data-driven approach to obtain insights into the immunological processes linked to the progression of the disease or to protective mechanisms. For this purpose, we evaluated the patterns of expression for genes that code proteins from the entire human genome in the peripheral blood from patients with visceral leishmaniasis, from individuals who remained asymptomatic after infections with L. infantum, from patients who were recovering from disease after treatment and from uninfected individuals. By employing computational analysis to evaluate the blood transcriptional activity of each group, we identified transcriptional signatures that correlate with previous findings obtained through different analytical methods. Moreover, our analyses uncovered hitherto unidentified molecular pathways and gene networks associated with the transcriptional profiles of individuals recovering from disease or that did not develop symptoms after infection. This suggests that activation of protective responses can be useful targets for the development of new therapies for visceral leishmaniasis.

Homology modelling of frequent HLA class-II alleles: A perspective to improve prediction of HLA binding peptide and understand the HLA associated disease susceptibility

Human leukocyte antigen (HLA) plays significant role via the regulation of immune system and contribute in the progression and protection of many diseases. HLA molecules bind and present peptides to T- cell receptors which generate the immune response. HLA peptide interaction and molecular function of HLA molecule is the key to predict peptide binding and understanding its role in different diseases. The availability of accurate three dimensional (3D) structures is the initial step towards this direction. In the present work, homology modelling of important and frequent HLA-DRB1 alleles (07:01, 11:01 and 09:01) was done and acceptable models were generated. These modelled alleles were further refined and cross validated by using several methods including Ramachandran plot, Z-score, ERRAT analysis and root mean square deviation (RMSD) calculations. It is known that numbers of allelic variants are related to the susceptibility or protection of various infectious diseases. Difference in amino acid sequences and structures of alleles were also studied to understand the association of HLA with disease susceptibility and protection. Susceptible alleles showed more amino acid variations than protective alleles in three selected diseases caused by different pathogens. Amino acid variations at binding site were found to be more than other part of alleles. RMSD values were also higher at variable positions within binding site. Higher RMSD values indicate that mutations occurring at peptide binding site alter protein structure more than rest of the protein. Hence, these findings and modelled structures can be used to design HLA-DRB1 binding peptides to overcome low prediction accuracy of HLA class II binding peptides. Furthermore, it may help to understand the allele specific molecular mechanisms involved in susceptibility/resistance against pathogenic diseases.

Association of HLA class I and II genes with cutaneous leishmaniasis: a case control study from Sri Lanka and a systematic review

Background

The outcome of leishmaniasis is an interplay between Leishamania and the host. Identifying contributory host genetic factors is complicated by the variability in phenotype, ethnicity and parasite species. Leishmaniasis is caused exclusively by L. donovani in Sri Lanka with localized cutaneous leishmaniasis (LCL) being the predominant form. We report here an association study of human leucocyte antigen (HLA) class I and II genes with LCL in Sri Lanka, the first on HLA associations in cutaneous leishmaniasis in a South Asian population.

Methods

An existing DNA repository of 200 each of patients and controls was typed for HLA-DQ by PCR-SSP. Next generation sequencing-based typing for HLA-A, HLA-B and HLA-DRB1 alleles was done in a subset of 280 samples. Association tests were performed on 28,489 genotyped and imputed SNPs spanning a region of 1.4 Mb across the HLA genes. To compare our results with similar studies, we carried out a systematic review to document all HLA associations reported to-date for cutaneous and muco-cutaneous leishmaniasis.

Results

DRB1*04 DQB1*02 (P = 0.03; Pc = 0.09), DRB1*07 DQB1*02 (P = 0.03; Pc = 0.09) haplotypes were absent in patients. B*07 (P = 0.007; Pc = 0.13; OR = 0.36; 95 % CI = 0.17–0.77) allele and DRB1*15 DQB1*06 (P = 0.00; Pc < 0.01; OR = 0.3; 95 % CI = 0.2–.0.6) haplotype were over represented in controls and DRB1*15 (P = 0.002; Pc = 0.01) allele was over represented in patients. Two SNPs (rs281864595/rs1050517) in the antigen recognition region of HLA-B, comprised a haplotype more frequent in controls (P = 0.04). The alleles identified by the systematic review to predispose or to protect from cutaneous/mucocutaneous leishmaniasis remained highly heterogeneous in different populations studied.

Conclusions

Our preliminary findings suggest a role for some class I and class II HLA genes in determining predisposition to LCL in this population which should be corroborated with further studies. The systematic review reiterates this need, as the purported susceptibility or protection gained by certain HLA alleles or haplotypes has rarely been independently verified.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-1626-8) contains supplementary material, which is available to authorized users.

The Potential Use of Forensic DNA Methods Applied to Sand Fly Blood Meal Analysis to Identify the Infection Reservoirs of Anthroponotic Visceral Leishmaniasis

Background

In the Indian sub-continent, visceral leishmaniasis (VL), also known as kala azar, is a fatal form of leishmaniasis caused by the kinetoplastid parasite Leishmania donovani and transmitted by the sand fly Phlebotomus argentipes. VL is prevalent in northeast India where it is believed to have an exclusive anthroponotic transmission cycle. There are four distinct cohorts of L. donovani exposed individuals who can potentially serve as infection reservoirs: patients with active disease, cured VL cases, patients with post kala azar dermal leishmaniasis (PKDL), and asymptomatic individuals. The relative contribution of each group to sustaining the transmission cycle of VL is not known.

Methodology/Principal Findings

To answer this critical epidemiological question, we have addressed the feasibility of an approach that would use forensic DNA methods to recover human DNA profiles from the blood meals of infected sand flies that would then be matched to reference DNA sampled from individuals living or working in the vicinity of the sand fly collections. We found that the ability to obtain readable human DNA fingerprints from sand flies depended entirely on the size of the blood meal and the kinetics of its digestion. Useable profiles were obtained from most flies within the first 24 hours post blood meal (PBM), with a sharp decline at 48 hours and no readable profiles at 72 hours. This early time frame necessitated development of a sensitive, nested-PCR method compatible with detecting L. donovani within a fresh, 24 hours blood meal in flies fed on infected hamsters.

Conclusion/Significance

Our findings establish the feasibility of the forensic DNA method to directly trace the human source of an infected blood meal, with constraints imposed by the requirement that the flies be recovered for analysis within 24 hours of their infective feed.

Visceral leishmaniasis, also known as kala-azar, is a fatal form of leishmaniasis that is caused by the protozoan parasite Leishmania donovani. In the Indian sub-continent, the parasite is transmitted between people by the sand fly Phlebotomus argentipes. There are four different groups of infected people that can carry the parasite: active cases, cured cases, cases of post-kala azar dermal leishmanisis, and people with asymptomatic infections. The relative contribution of these different groups to sustaining the transmission cycle is not known. In the current work, we explore the feasibility of forensic DNA methods that could be applied to directly trace the human source of an infected blood meal in flies collected from high transmission areas. We report that a readable human DNA profile can be generated from the majority of sand flies so long as the flies are sampled within 24 hours after blood feeding. We developed a highly sensitive nested-PCR method that was able to detect parasites present in a fresh, 24 hour blood meal after feeding on an infected host.

Fine mapping under linkage peaks for symptomatic or asymptomatic outcomes of Leishmania infantum infection in Brazil

Infection with the protozoan Leishmania infantum can lead to asymptomatic infection and protective immunity, or to the progressive and potentially fatal disease visceral leishmaniasis (VL). Published studies show host genetic background determines in part whether infected individuals will develop a symptomatic or asymptomatic outcome. The purpose of the current study was to fine map chromosome regions previously linked with risk for symptomatic (chromosome 9) or asymptomatic (chromosomes 15 and 19) manifestations of L. infantum infection. We conducted a family-based genetic study of VL and asymptomatic infection (detected by a DTH skin test) with a final post quality control sample of 961 individuals with full genotype and phenotype information from highly endemic neighborhoods of northeast Brazil. A total of 5485 SNPs under the linkage peaks on chromosomes 9, 15 and 19 were genotyped. No strong SNP associations were observed for the DTH phenotype. The most significant associations with the VL phenotype were with SNP rs1470217 (p=5.9e-05; pcorrected=0.057) on chromosome 9, and with SNP rs8107014 (p=1.4e-05; pcorrected=0.013) on chromosome 19. SNP rs1470217 is situated in a 180kb intergenic region between TMEM215 (Transmembrane protein 215) and APTX (Aprataxin). SNP rs8107014 lies in the intron between exons 26 and 27 of a 34 exon transcript (ENST00000204005) of LTBP4, (Latent transforming growth factor-beta-binding protein 4a). The latter supports growing evidence that the transforming growth factor-beta pathway is important in the immunopathogenesis of VL.

Control for Population Structure and Relatedness for Binary Traits in Genetic Association Studies via Logistic Mixed Models

Linear mixed models (LMMs) are widely used in genome-wide association studies (GWASs) to account for population structure and relatedness, for both continuous and binary traits. Motivated by the failure of LMMs to control type I errors in a GWAS of asthma, a binary trait, we show that LMMs are generally inappropriate for analyzing binary traits when population stratification leads to violation of the LMM's constant-residual variance assumption. To overcome this problem, we develop a computationally efficient logistic mixed model approach for genome-wide analysis of binary traits, the generalized linear mixed model association test (GMMAT). This approach fits a logistic mixed model once per GWAS and performs score tests under the null hypothesis of no association between a binary trait and individual genetic variants. We show in simulation studies and real data analysis that GMMAT effectively controls for population structure and relatedness when analyzing binary traits in a wide variety of study designs.

Multiple testing correction in linear mixed models

BACKGROUND

Multiple hypothesis testing is a major issue in genome-wide association studies (GWAS), which often analyze millions of markers. The permutation test is considered to be the gold standard in multiple testing correction as it accurately takes into account the correlation structure of the genome. Recently, the linear mixed model (LMM) has become the standard practice in GWAS, addressing issues of population structure and insufficient power. However, none of the current multiple testing approaches are applicable to LMM.

RESULTS

We were able to estimate per-marker thresholds as accurately as the gold standard approach in real and simulated datasets, while reducing the time required from months to hours. We applied our approach to mouse, yeast, and human datasets to demonstrate the accuracy and efficiency of our approach.

CONCLUSIONS

We provide an efficient and accurate multiple testing correction approach for linear mixed models. We further provide an intuition about the relationships between per-marker threshold, genetic relatedness, and heritability, based on our observations in real data.

Successful Treatment of Human Visceral Leishmaniasis Restores Antigen-Specific IFN-γ, but not IL-10 Production

One of the key immunological characteristics of active visceral leishmaniasis (VL) is a profound immunosuppression and impaired production of Interferon-γ (IFN-γ). However, recent studies from Bihar in India showed using a whole blood assay, that whole blood cells have maintained the capacity to produce IFN-γ. Here we tested the hypothesis that a population of low-density granulocytes (LDG) might contribute to T cell responses hyporesponsiveness via the release of arginase. Our results show that this population is affected by the anticoagulant used to collect blood: the frequency of LDGs is significantly lower when the blood is collected with heparin as compared to EDTA; however, the anticoagulant does not impact on the levels of arginase released. Next, we assessed the capacity of whole blood cells from patients with active VL to produce IFN-γ and IL-10 in response to antigen-specific and polyclonal activation. Our results show that whole blood cells produce low or levels below detection limit of IFN-γ and IL-10, however, after successful treatment of VL patients, these cells gradually regain their capacity to produce IFN-γ, but not IL-10, in response to activation. These results suggest that in contrast to VL patients from Bihar, India, whole blood cells from VL patients from Gondar, Ethiopia, have lost their ability to produce IFN-γ during active VL and that active disease is not associated with sustained levels of IL-10 production following stimulation.

The leishmaniases, a group of diseases caused by Leishmania parasites, belong to the most neglected tropical diseases: they are mainly found in low-income countries and affect the poorest populations. These parasites infect cells of the immune system called macrophages, which can kill the intracellular parasites in response to soluble mediators they receive from other cells of the immune system, the lymphocytes. Visceral leishmaniasis is the most severe form of the leishmaniases and is characterized by enlarged liver and spleen, fever, weight-loss and anaemia and represents a major public health problem in Ethiopia. Currently there is no vaccine available, the existing treatment has many severe side effects and drug-resistance is increasing. In the present study, we worked with patients suffering from visceral leishmaniasis. This form of the disease is fatal if the patients are not treated. We studied the ability of lymphocytes isolated from their blood to produce soluble mediators before and at different times after the end of treatment. Our results show that the lymphocytes have an impaired capacity to produce the soluble mediator required to instruct infected cells to kill the intracellular parasites, but that this lack of response is gradually restored with time after successful treatment.

Progress and promise in understanding the genetic basis of common diseases

Susceptibility to common human diseases is influenced by both genetic and environmental factors. The explosive growth of genetic data, and the knowledge that it is generating, are transforming our biological understanding of these diseases. In this review, we describe the technological and analytical advances that have enabled genome-wide association studies to be successful in identifying a large number of genetic variants robustly associated with common disease. We examine the biological insights that these genetic associations are beginning to produce, from functional mechanisms involving individual genes to biological pathways linking associated genes, and the identification of functional annotations, some of which are cell-type-specific, enriched in disease associations. Although most efforts have focused on identifying and interpreting genetic variants that are irrefutably associated with disease, it is increasingly clear that--even at large sample sizes--these represent only the tip of the iceberg of genetic signal, motivating polygenic analyses that consider the effects of genetic variants throughout the genome, including modest effects that are not individually statistically significant. As data from an increasingly large number of diseases and traits are analysed, pleiotropic effects (defined as genetic loci affecting multiple phenotypes) can help integrate our biological understanding. Looking forward, the next generation of population-scale data resources, linking genomic information with health outcomes, will lead to another step-change in our ability to understand, and treat, common diseases.

Tolerance to Trypanosomatids: A Threat, or a Key for Disease Elimination?

So far, research on trypanosomatid infections has been driven by 'disease by disease' approaches, leading to different concepts and control strategies. It is, however, increasingly clear that they share common features such as the ability to generate long-lasting asymptomatic infections in their mammalian hosts. Trypanotolerance, long integrated in animal African trypanosomiasis control, historically refers to the ability of cattle breeds to limit Trypanosoma infection and pathology, but has only recently been recognized in humans. Whilst trypanotolerance is absent from the vocabulary on leishmaniasis and Chagas disease, asymptomatic infections also occur. We review the concept of trypanotolerance across the trypanosomatids and discuss the importance of asymptomatic carriage in the current context of elimination.