Gene expression patterns associated with Leishmania panamensis infection in macrophages from BALB/c and C57BL/6 mice

Unraveling molecular mechanistic disparities in pathogenic visceral Leishmania resistance between reptiles and mammals through comparative transcriptomic analyses

Leishmaniasis is one of the most important neglected tropical parasitic diseases, manifesting various clinical forms depending on the parasite species and the genetic background of the host. In order to elucidate the underlying mechanisms of reptilian defense against pathogenic Leishmania species and to delineate the global gene expression profile alterations during host-pathogen interaction, we established experimental animal and cell models using both heterothermic lizards (Phrynocephalus przewalskii) and homothermic mammals (BALB/c mice) infected with pathogenic Leishmania infantum (high virulence HCZ strain) and Leishmania donovani (low virulence 801 strain). Overall, the lizards didn't show any obvious clinical symptoms or immune responses in vivo. Using RNA-seq methodology, differentially expressed genes identified in the HCZ and 801-comparison groups of P. przewalskii were primarily associated with arginine biosynthesis, the MAPK signaling pathway and the PI3K-Akt signaling pathway. In contrast, higher parasite loads, exacerbated hepatic inflammatory lesions and enhanced immune responses were observed in BALB/c mice, with DEGs predominantly associated with immunological diseases, innate and adaptive immune responses. By integrating transcriptional data from reptile and mammalian hosts, we elucidated the pivotal role of amino acid metabolism and lipid metabolism in parasite control. In contrast to findings from animal experiments, Leishmania parasites effectively infected peritoneal macrophages of lizards in vitro, demonstrating a high infection rate. Furthermore, we used RT-qPCR to detect changes in cytokine expression in macrophages and found that Th1-type cytokines were significantly upregulated in lizards, facilitating the clearance of the HCZ strain 24 hours post-infection. Conversely, cytokine expression was generally suppressed in BALB/c mice, allowing immune evasion by the parasites.

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.

Genetic Iron Overload Hampers Development of Cutaneous Leishmaniasis in Mice

The survival, growth, and virulence of Leishmania spp., a group of protozoan parasites, depends on the proper access and regulation of iron. Macrophages, Leishmania's host cell, may divert iron traffic by reducing uptake or by increasing the efflux of iron via the exporter ferroportin. This parasite has adapted by inhibiting the synthesis and inducing the degradation of ferroportin. To study the role of iron in leishmaniasis, we employed Hjv^-/- mice, a model of hemochromatosis. The disruption of hemojuvelin (Hjv) abrogates the expression of the iron hormone hepcidin. This allows unrestricted iron entry into the plasma from ferroportin-expressing intestinal epithelial cells and tissue macrophages, resulting in systemic iron overload. Mice were injected with Leishmania major in hind footpads or intraperitoneally. Compared with wild-type controls, Hjv^-/- mice displayed transient delayed growth of L. major in hind footpads, with a significant difference in parasite burden 4 weeks post-infection. Following acute intraperitoneal exposure to L. major, Hjv^-/- peritoneal cells manifested increased expression of inflammatory cytokines and chemokines (Il1b, Tnfa, Cxcl2, and Ccl2). In response to infection with L. infantum, the causative agent of visceral leishmaniasis, Hjv^-/- and control mice developed similar liver and splenic parasite burden despite vastly different tissue iron content and ferroportin expression. Thus, genetic iron overload due to hemojuvelin deficiency appears to mitigate the early development of only cutaneous leishmaniasis.

Host M-CSF induced gene expression drives changes in susceptible and resistant mice-derived BMdMs upon Leishmania major infection

Leishmaniases are a group of diseases with different clinical manifestations. Macrophage-Leishmania interactions are central to the course of the infection. The outcome of the disease depends not only on the pathogenicity and virulence of the parasite, but also on the activation state, the genetic background, and the underlying complex interaction networks operative in the host macrophages. Mouse models, with mice strains having contrasting behavior in response to parasite infection, have been very helpful in exploring the mechanisms underlying differences in disease progression. We here analyzed previously generated dynamic transcriptome data obtained from Leishmania major (L. major) infected bone marrow derived macrophages (BMdMs) from resistant and susceptible mouse. We first identified differentially expressed genes (DEGs) between the M-CSF differentiated macrophages derived from the two hosts, and found a differential basal transcriptome profile independent of Leishmania infection. These host signatures, in which 75% of the genes are directly or indirectly related to the immune system, may account for the differences in the immune response to infection between the two strains. To gain further insights into the underlying biological processes induced by L. major infection driven by the M-CSF DEGs, we mapped the time-resolved expression profiles onto a large protein-protein interaction (PPI) network and performed network propagation to identify modules of interacting proteins that agglomerate infection response signals for each strain. This analysis revealed profound differences in the resulting responses networks related to immune signaling and metabolism that were validated by qRT-PCR time series experiments leading to plausible and provable hypotheses for the differences in disease pathophysiology. In summary, we demonstrate that the host's gene expression background determines to a large degree its response to L. major infection, and that the gene expression analysis combined with network propagation is an effective approach to help identifying dynamically altered mouse strain-specific networks that hold mechanistic information about these contrasting responses to infection.

The issue beyond resistance: Methicillin-resistant Staphylococcus epidermidis biofilm formation is induced by subinhibitory concentrations of cloxacillin, cefazolin, and clindamycin

Staphylococcus epidermis is one of the most frequent causes of device-associated infections due to biofilm formation. Current reports noted that subinhibitory concentrations of antibiotics induce biofilm production in some bacteria. Accordingly, we evaluated the effect of exposure of different subinhibitory concentrations of cloxacillin, cefazolin, clindamycin, and vancomycin on the biofilm formation of methicillin-resistant S. epidermidis (MRSE). Antimicrobial susceptibility testing and minimum inhibitory/bactericidal concentration of antimicrobial agents were determined. MRSE isolates were selected, and their biofilm formation ability was evaluated. The effect of subinhibitory concentrations of cloxacillin, cefazolin, clindamycin, and vancomycin, antibiotics selected among common choices in the clinic, on MRSE biofilm formation was determined by the microtitre method. Besides, the effect of subinhibitory concentrations of cloxacillin, cefazolin, clindamycin, and vancomycin on the expression of the biofilm-associated genes icaA and atlE was evaluated by Reverse-transcription quantitative real-time polymerase chain reaction (RT-qPCR). Antimicrobial susceptibility patterns of MRSE strains showed a high level of resistance as follows: 80%, 53.3%, 33.3%, 33.3%, and 26.6%, for erythromycin, trimethoprim-sulfamethoxazole, tetracycline, clindamycin, and gentamicin, respectively. Besides, 73.3% of S. epidermidis strains were Multidrug-resistant (MDR). Minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC) values were in the range of 0.5 to512 μg/mL and 1 to1024 μg/mL for cloxacillin, 0.125 to256 μg/mL and 1 to512 μg/mL for cefazolin, 0.125 to64 μg/mL and 4 to>1024 μg/mL for clindamycin, and 2 to32 μg/mL and 4 to32 μg/mL for vancomycin, respectively. The findings showed that subinhibitory concentrations of cloxacillin, cefazolin, and clindamycin induce biofilm production in MRSE strains. In particular, the OD values of strains were in the range of 0.09–0.95, 0.05–0.86, and 0.06–1 toward cloxacillin, cefazolin, and clindamycin, respectively. On the other hand, exposure to subinhibitory vancomycin concentrations did not increase the biofilm formation in MRSE strains. The findings also demonstrated that sub-MIC of antibiotics up-regulated biofilm-associated genes. In particular, atlE and icaA were up-regulated 0.062 to 1.16 and 0.078 to 1.48 folds, respectively, for cloxacillin, 0.11 to 0.8, and 0.1 to 1.3 folds for cefazolin, 0.18 to 0.98, and 0.19 to 1.4 folds, respectively, for clindamycin. In contrast, the results showed that sub-MIC of vancomycin did not increase the biofilm-associated genes. These findings overall show that exposure to sub-MIC of traditional antibiotics can cause biofilm induction in MRSE, thereby increasing the survival and persistence on various surfaces that worsen the condition of comorbid infections.

Dissection of the macrophage response towards infection by the Leishmania-viral endosymbiont duo and dynamics of the type I interferon response

Leishmania RNA virus 1 (LRV1) is a double-stranded RNA virus found in some strains of the human protozoan parasite Leishmania, the causative agent of leishmaniasis, a neglected tropical disease. Interestingly, the presence of LRV1 inside Leishmania constitutes an important virulence factor that worsens the leishmaniasis outcome in a type I interferon (IFN)–dependent manner and contributes to treatment failure. Understanding how macrophages respond toward Leishmania alone or in combination with LRV1 as well as the role that type I IFNs may play during infection is fundamental to oversee new therapeutic strategies. To dissect the macrophage response toward infection, RNA sequencing was performed on murine wild-type and Ifnar-deficient bone marrow–derived macrophages infected with Leishmania guyanensis (Lgy) devoid or not of LRV1. Additionally, macrophages were treated with poly I:C (mimetic virus) or with type I IFNs. By implementing a weighted gene correlation network analysis, the groups of genes (modules) with similar expression patterns, for example, functionally related, coregulated, or the members of the same functional pathway, were identified. These modules followed patterns dependent on Leishmania, LRV1, or Leishmania exacerbated by the presence of LRV1. Not only the visualization of how individual genes were embedded to form modules but also how different modules were related to each other were observed. Thus, in the context of the observed hyperinflammatory phenotype associated to the presence of LRV1, it was noted that the biomarkers tumor-necrosis factor α (TNF-α) and the interleukin 6 (IL-6) belonged to different modules and that their regulating specific Src-family kinases were segregated oppositely. In addition, this network approach revealed the strong and sustained effect of LRV1 on the macrophage response and genes that had an early, late, or sustained impact during infection, uncovering the dynamics of the IFN response. Overall, this study contributed to shed light and dissect the intricate macrophage response toward infection by the Leishmania-LRV1 duo and revealed the crosstalk between modules made of coregulated genes and provided a new resource that can be further explored to study the impact of Leishmania on the macrophage response.