Structural and mechanistic insights into the function of Leishmania ribosome lacking a single pseudouridine modification.

Leishmania is the causative agent of cutaneous and visceral diseases affecting millions of individuals worldwide. Pseudouridine (Ψ), the most abundant modification on rRNA, changes during the parasite life cycle. Alterations in the level of a specific Ψ in helix 69 (H69) affected ribosome function. To decipher the molecular mechanism of this phenotype, we determine the structure of ribosomes lacking the single Ψ and its parental strain at ∼2.4-3 Å resolution using cryo-EM. Our findings demonstrate the significance of a single Ψ on H69 to its structure and the importance for its interactions with helix 44 and specific tRNAs. Our study suggests that rRNA modification affects translation of mRNAs carrying codon bias due to selective accommodation of tRNAs by the ribosome. Based on the high-resolution structures, we propose a mechanism explaining how the ribosome selects specific tRNAs.

Intracellular persistence of Leishmania tarentolae in primary canine macrophage cells.

Leishmania tarentolae is a non-pathogenic species first isolated from geckoes in the Mediterranean basin. The finding that dogs test positive against both Leishmania infantum and L. tarentolae raises questions regarding the ability of the latter species to persist and adapt to new hosts. This study aimed to evaluate in vitro the capability of L. tarentolae to colonize, survive and persist in canine primary monocyte-derived mononuclear cells. Monocytes were isolated from dog whole blood samples and placed in 24-well plates for differentiation into macrophages and for incubation with L. tarentolae field-isolated strains (RI-325 and SF-178) and laboratory (LEM-124) strain; the parasite burden was assessed at different time points post-infection. The L. infantum laboratory strain (MON-1) was used as control. Infection parameters were evaluated by microscopy, counting the number of amastigotes/200 infected cells, and by duplex real-time PCR from supernatants and detached cells. Similar to L. infantum, L. tarentolae strains developed into round-shaped amastigote-like forms, with higher infection rates detected at 4 h followed by an overall decrease until 48 h. RI-325 presented also a higher infection rate at 72 h. Data showed that L. tarentolae strains infect and persist inside in vitro primary canine mononuclear cells, opening new perspectives for further laboratory studies.

Amphotericin B resistance correlates with increased fitness in vitro and in vivo in Leishmania (Mundinia) martiniquensis.

Amphotericin B (AmpB) deoxycholate is the available first-line drug used to treat visceral leishmaniasis caused by Leishmania (Mundinia) martiniquensis, however, some cases of AmpB treatment failure have been reported in Thailand. Resistance to drugs is known to affect parasite fitness with a potential impact on parasite transmission but still little is known about the effect of resistance to drugs on L. martiniquensis. Here we aimed to gain insight into the fitness changes occurring after treatment failure or in vitro-induced resistance to AmpB. L. martiniquensis parasites isolated from a patient before (LSCM1) and after relapse (LSCM1-6) were compared for in vitro and in vivo fitness changes together with an in vitro induced AmpB-resistant parasite generated from LSCM1 parasites (AmpBRP2i). Results revealed increased metacyclogenesis of the AmpBPR2i and LSCM1-6 strains (AmpB-resistant strains) compared to the LSCM1 strain and increased fitness with respect to growth and infectivity. The LSCM1-6 and AmpBRP2i strains were present in mice for longer periods compared to the LSCM1 strain, but no clinical signs of the disease were observed. These results suggest that the AmpB-resistant parasites could be more efficiently transmitted to humans and maintained in asymptomatic hosts longer than the susceptible strain. The asymptomatic hosts therefore may represent "reservoirs" for the resistant parasites enhancing transmission. The results in this study advocate an urgent need to search and monitor for AmpB-resistant L. martiniquensis in patients with relapsing leishmaniasis and in asymptomatic patients, especially, in HIV/Leishmania coinfected patients.

Leishmania allelic selection during experimental sand fly infection correlates with mutational signatures of oxidative DNA damage.

Trypanosomatid pathogens are transmitted by blood-feeding insects, causing devastating human infections. These parasites show important phenotypic shifts that often impact parasite pathogenicity, tissue tropism, or drug susceptibility. The evolutionary mechanisms that allow for the selection of such adaptive phenotypes remain only poorly investigated. Here, we use Leishmania donovani as a trypanosomatid model pathogen to assess parasite evolutionary adaptation during experimental sand fly infection. Comparing the genome of the parasites before and after sand fly infection revealed a strong population bottleneck effect as judged by allele frequency analysis. Apart from random genetic drift caused by the bottleneck effect, our analyses revealed haplotype and allelic changes during sand fly infection that seem under natural selection given their convergence between independent biological replicates. Our analyses further uncovered signature mutations of oxidative DNA damage in the parasite genomes after sand fly infection, suggesting that Leishmania suffers from oxidative stress inside the insect digestive tract. Our results propose a model of Leishmania genomic adaptation during sand fly infection, with oxidative DNA damage and DNA repair processes likely driving haplotype and allelic selection. The experimental and computational framework presented here provides a useful blueprint to assess evolutionary adaptation of other eukaryotic pathogens inside their insect vectors, such as Plasmodium spp, Trypanosoma brucei, and Trypanosoma cruzi.

Trifloxystrobin blocks the growth of Theileria parasites and is a promising drug to treat Buparvaquone resistance.

Theileria parasites are responsible for devastating cattle diseases, causing major economic losses across Africa and Asia. Theileria spp. stand apart from other apicomplexa parasites by their ability to transform host leukocytes into immortalized, hyperproliferating, invasive cells that rapidly kill infected animals. The emergence of resistance to the theilericidal drug Buparvaquone raises the need for new anti-Theileria drugs. We developed a microscopy-based screen to reposition drugs from the open-access Medicines for Malaria Venture (MMV) Pathogen Box. We show that Trifloxystrobin (MMV688754) selectively kills lymphocytes or macrophages infected with Theileria annulata or Theileria parva parasites. Trifloxystrobin treatment reduced parasite load in vitro as effectively as Buparvaquone, with similar effects on host gene expression, cell proliferation and cell cycle. Trifloxystrobin also inhibited parasite differentiation to merozoites (merogony). Trifloxystrobin inhibition of parasite survival is independent of the parasite TaPin1 prolyl isomerase pathway. Furthermore, modeling studies predicted that Trifloxystrobin and Buparvaquone could interact distinctly with parasite Cytochrome B and we show that Trifloxystrobin was still effective against Buparvaquone-resistant cells harboring TaCytB mutations. Our study suggests that Trifloxystrobin could provide an effective alternative to Buparvaquone treatment and represents a promising candidate for future drug development against Theileria spp.

Retrospective Analysis of Leishmaniasis in Sicily (Italy) from 2013 to 2021: One-Health Impact and Future Control Strategies.

Leishmaniasis is an important vector-borne disease that represents a serious public health problem, including in Sicily (Italy), which is considered an endemic area. We collected canine, feline and human data from 2013 to 2021 in Sicily, while entomological surveys were conducted only in 2013 and 2021. Overall, 23,794/74,349 (34.4%) of dogs and 274/4774 (11.8%) of cats were positive in one or more diagnostic tests. A total of 467 cases of human Leishmaniasis were reported, with 71% showing cutaneous and 29% visceral involvement. The provinces with the largest number of patients were Agrigento (45.4%) and Palermo (37%). In 2013, Phlebotomus perfiliewi was the dominant sandfly species in Sicily (68.7%), followed by Phlebotomus perniciosus (17.2%) and Sergentomya minuta (14%). In 2021, Phlebotomusperfiliewi was confirmed as the most common species (61.6%), followed by Phlebotomusperniciosus (33.1%) and Sergentomyaminuta (4.7%). Of particular interest was the identification of Phlebotomus papatasi (0.41%) in Agrigento. Our retrospective study can inform health authorities for the development of appropriate screening, treatment and control strategies to reduce Leishmania incidence rate. This study examined the present state of Leishmaniasis control, surveillance, and prevention in Sicily, but also highlighted deficiencies that could be addressed through the application of One-Health principles.

Effects of Structurally Different HDAC Inhibitors against Trypanosoma cruzi, Leishmania, and Schistosoma mansoni.

Neglected tropical diseases (NTDs), including trypanosomiasis, leishmaniasis, and schistosomiasis, result in a significant burden in terms of morbidity and mortality worldwide every year. Current antiparasitic drugs suffer from several limitations such as toxicity, no efficacy toward all of the forms of the parasites' life cycle, and/or induction of resistance. Histone-modifying enzymes play a crucial role in parasite growth and survival; thus, the use of epigenetic drugs has been suggested as a strategy for the treatment of NTDs. We tested structurally different HDACi 1-9, chosen from our in-house library or newly synthesized, against Trypanosoma cruzi, Leishmania spp, and Schistosoma mansoni. Among them, 4 emerged as the most potent against all of the tested parasites, but it was too toxic against host cells, hampering further studies. The retinoic 2'-aminoanilide 8 was less potent than 4 in all parasitic assays, but as its toxicity is considerably lower, it could be the starting structure for further development. In T. cruzi, compound 3 exhibited a single-digit micromolar inhibition of parasite growth combined with moderate toxicity. In S. mansoni, 4's close analogs 17-20 were tested in new transformed schistosomula (NTS) and adult worms displaying high death induction against both parasite forms. Among them, 17 and 19 exhibited very low toxicity in human retinal pigment epithelial (RPE) cells, thus being promising compounds for further optimization.

Experimental evolution links post-transcriptional regulation to Leishmania fitness gain.

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

Going ballistic: Leishmania nuclear subversion of host cell plasticity.

Intracellular parasites have evolved intricate strategies to subvert host cell functions for their own survival. These strategies are particularly damaging to the host if the infection involves immune cells, as illustrated by protozoan parasites of the genus Leishmania that thrive inside mononuclear phagocytic cells, causing devastating immunopathologies. While the impact of Leishmania infection on host cell phenotype and functions has been well documented, the regulatory mechanisms underlying host cell subversion were only recently investigated. Here we summarize the current knowledge on how Leishmania infection affects host nuclear activities and propose thought-provoking new concepts on the reciprocal relationship between epigenetic and transcriptional regulation in host cell phenotypic plasticity, its potential subversion by the intracellular parasite, and its relevance for host-directed therapy.

GIP: an open-source computational pipeline for mapping genomic instability from protists to cancer cells.

Genome instability has been recognized as a key driver for microbial and cancer adaptation and thus plays a central role in many diseases. Genome instability encompasses different types of genomic alterations, yet most available genome analysis software are limited to just one type of mutation. To overcome this limitation and better understand the role of genetic changes in enhancing pathogenicity we established GIP, a novel, powerful bioinformatic pipeline for comparative genome analysis. Here, we show its application to whole genome sequencing datasets of Leishmania, Plasmodium, Candida and cancer. Applying GIP on available data sets validated our pipeline and demonstrated the power of our tool to drive biological discovery. Applied to Plasmodium vivax genomes, our pipeline uncovered the convergent amplification of erythrocyte binding proteins and identified a nullisomic strain. Re-analyzing genomes of drug adapted Candida albicans strains revealed correlated copy number variations of functionally related genes, strongly supporting a mechanism of epistatic adaptation through interacting gene-dosage changes. Our results illustrate how GIP can be used for the identification of aneuploidy, gene copy number variations, changes in nucleic acid sequences, and chromosomal rearrangements. Altogether, GIP can shed light on the genetic bases of cell adaptation and drive disease biomarker discovery.

Genome instability drives epistatic adaptation in the human pathogen Leishmania.

How genome instability is harnessed for fitness gain despite its potential deleterious effects is largely elusive. An ideal system to address this important open question is provided by the protozoan pathogen Leishmania, which exploits frequent variations in chromosome and gene copy number to regulate expression levels. Using ecological genomics and experimental evolution approaches, we provide evidence that Leishmania adaptation relies on epistatic interactions between functionally associated gene copy number variations in pathways driving fitness gain in a given environment. We further uncover posttranscriptional regulation as a key mechanism that compensates for deleterious gene dosage effects and provides phenotypic robustness to genetically heterogenous parasite populations. Finally, we correlate dynamic variations in small nucleolar RNA (snoRNA) gene dosage with changes in ribosomal RNA 2'-O-methylation and pseudouridylation, suggesting translational control as an additional layer of parasite adaptation. Leishmania genome instability is thus harnessed for fitness gain by genome-dependent variations in gene expression and genome-independent compensatory mechanisms. This allows for polyclonal adaptation and maintenance of genetic heterogeneity despite strong selective pressure. The epistatic adaptation described here needs to be considered in Leishmania epidemiology and biomarker discovery and may be relevant to other fast-evolving eukaryotic cells that exploit genome instability for adaptation, such as fungal pathogens or cancer.

Nitric oxide controls proliferation of Leishmania major by inhibiting the recruitment of permissive host cells.

Nitric oxide (NO) is an important antimicrobial effector but also prevents unnecessary tissue damage by shutting down the recruitment of monocyte-derived phagocytes. Intracellular pathogens such as Leishmania major can hijack these cells as a niche for replication. Thus, NO might exert containment by restricting the availability of the cellular niche required for efficient pathogen proliferation. However, such indirect modes of action remain to be established. By combining mathematical modeling with intravital 2-photon biosensors of pathogen viability and proliferation, we show that low L. major proliferation results not from direct NO impact on the pathogen but from reduced availability of proliferation-permissive host cells. Although inhibiting NO production increases recruitment of these cells, and thus pathogen proliferation, blocking cell recruitment uncouples the NO effect from pathogen proliferation. Therefore, NO fulfills two distinct functions for L. major containment: permitting direct killing and restricting the supply of proliferation-permissive host cells.

Cutaneous Leishmaniasis in Algeria; Highlight on the Focus of M'Sila.

Algeria ranks second after Afghanistan for the incidence of cutaneous leishmaniasis (CL) worldwide. Here, we report a 34-years retrospective analysis of CL in Algeria and focused on the most affected region, the M'Sila province. All 66 cutaneous isolates corresponded to Leishmania (L.) major. Our study of the sandfly and rodent fauna further highlighted the high density of Phlebotomus papatasi and additional phlebotomine species of medical importance, not previously identified in M'Sila. Wild rodents belonging to nine species were trapped in M'Sila, and Psammomys obesus and Meriones shawi were found infected by L. major. In addition, Leishmania infantum was isolated from two visceral leishmaniasis cases, one dog and its proven vectors (P. perniciosus, P. longicuspis, and P. perfiliewi) inventoried during the survey. The high incidence of CL in the M'Sila province is likely a consequence of the increase in minimum temperatures recorded that constitutes suitable conditions for establishing a high endemicity and leads to an explosive rise in leishmaniases cases in this region. A thorough investigation of the underlying risk factors is urgently needed to detect new cases earlier. All these would improve the preparedness to fight the disease.

Dangerous Duplicity: The Dual Functions of Casein Kinase 1 in Parasite Biology and Host Subversion.

Casein Kinase 1 (CK1) family members are serine/threonine protein kinases that are involved in many biological processes and highly conserved in eukaryotes from protozoan to humans. Even though pathogens exploit host CK1 signaling pathways to survive, the role of CK1 in infectious diseases and host/pathogen interaction is less well characterized compared to other diseases, such as cancer or neurodegenerative diseases. Here we present the current knowledge on CK1 in protozoan parasites highlighting their essential role for parasite survival and their importance for host-pathogen interactions. We also discuss how the dual requirement of CK1 family members for parasite biological processes and host subversion could be exploited to identify novel antimicrobial interventions.

Colonization and genetic diversification processes of Leishmania infantum in the Americas.

Leishmania infantum causes visceral leishmaniasis, a deadly vector-borne disease introduced to the Americas during the colonial era. This non-native trypanosomatid parasite has since established widespread transmission cycles using alternative vectors, and human infection has become a significant concern to public health, especially in Brazil. A multi-kilobase deletion was recently detected in Brazilian L. infantum genomes and is suggested to reduce susceptibility to the anti-leishmanial drug miltefosine. We show that deletion-carrying strains occur in at least 15 Brazilian states and describe diversity patterns suggesting that these derive from common ancestral mutants rather than from recurrent independent mutation events. We also show that the deleted locus and associated enzymatic activity is restored by hybridization with non-deletion type strains. Genetic exchange appears common in areas of secondary contact but also among closely related parasites. We examine demographic and ecological scenarios underlying this complex L. infantum population structure and discuss implications for disease control.

Identification of the Host Substratome of Leishmania-Secreted Casein Kinase 1 Using a SILAC-Based Quantitative Mass Spectrometry Assay.

Leishmaniasis is a severe public health problem, caused by the protozoan Leishmania. This parasite has two developmental forms, extracellular promastigote in the insect vector and intracellular amastigote in the mammalian host where it resides inside the phagolysosome of macrophages. Little is known about the virulence factors that regulate host-pathogen interactions and particularly host signalling subversion. All the proteomes of Leishmania extracellular vesicles identified the presence of Leishmania casein kinase 1 (L-CK1.2), a signalling kinase. L-CK1.2 is essential for parasite survival and thus might be essential for host subversion. To get insights into the functions of L-CK1.2 in the macrophage, the systematic identification of its host substrates is crucial, we thus developed an easy method to identify substrates, combining phosphatase treatment, in vitro kinase assay and Stable Isotope Labelling with Amino acids in Cell (SILAC) culture-based mass spectrometry. Implementing this approach, we identified 225 host substrates as well as a potential novel phosphorylation motif for CK1. We confirmed experimentally the enrichment of our substratome in bona fide L-CK1.2 substrates and showed they were also phosphorylated by human CK1δ. L-CK1.2 substratome is enriched in biological processes such as "viral and symbiotic interaction," "actin cytoskeleton organisation" and "apoptosis," which are consistent with the host pathways modified by Leishmania upon infection, suggesting that L-CK1.2 might be the missing link. Overall, our results generate important mechanistic insights into the signalling of host subversion by these parasites and other microbial pathogens adapted for intracellular survival.

In vitro identification of imidazo[1,2-a]pyrazine-based antileishmanial agents and evaluation of L. major casein kinase 1 inhibition.

Leishmaniasis constitutes a severe public health problem, with an estimated prevalence of 12 million cases. This potentially fatal disease has a worldwide distribution and in 2012, the fatal Visceral Leishmaniasis (VL) was declared as new emerging disease in Europe, mainly due to global warming, with expected important public health impact. The available treatments are toxic, costly or lead to parasite resistance, thus there is an urgent need for new drugs with new mechanism of action. Previously, we reported the discovery of CTN1122, a potent imidazo[1,2-a]pyrazine-based antileishmanial hit compound targeting L-CK1.2 at low micromolar ranges. Here, we described structurally related, safe and selective compounds endowed with antiparasitic properties, better than miltefosine, the reference therapy by oral route. L-CK1.2 homology model gave the first structural explanations of the role of 4-pyridyl (CTN1122) and 2-aminopyrimidin-4-yl (compound 21) moieties, at the position 3 of the central core, in the low micromolar to nanomolar L-CK1.2 inhibition, whereas N-methylpyrazole derivative 11 remained inactive against the parasite kinase.

Comparative transcriptomic analysis of antimony resistant and susceptible Leishmania infantum lines.

BACKGROUND: One of the major challenges to leishmaniasis treatment is the emergence of parasites resistant to antimony. To study differentially expressed genes associated with drug resistance, we performed a comparative transcriptomic analysis between wild-type and potassium antimonyl tartrate (SbIII)-resistant Leishmania infantum lines using high-throughput RNA sequencing. METHODS: All the cDNA libraries were constructed from promastigote forms of each line, sequenced and analyzed using STAR for mapping the reads against the reference genome (L. infantum JPCM5) and DESeq2 for differential expression statistical analyses. All the genes were functionally annotated using sequence similarity search. RESULTS: The analytical pipeline considering an adjusted p-value < 0.05 and fold change > 2.0 identified 933 transcripts differentially expressed (DE) between wild-type and SbIII-resistant L. infantum lines. Out of 933 DE transcripts, 504 presented functional annotation and 429 were assigned as hypothetical proteins. A total of 837 transcripts were upregulated and 96 were downregulated in the SbIII-resistant L. infantum line. Using this DE dataset, the proteins were further grouped in functional classes according to the gene ontology database. The functional enrichment analysis for biological processes showed that the upregulated transcripts in the SbIII-resistant line are associated with protein phosphorylation, microtubule-based movement, ubiquitination, host-parasite interaction, cellular process and other categories. The downregulated transcripts in the SbIII-resistant line are assigned in the GO categories: ribonucleoprotein complex, ribosome biogenesis, rRNA processing, nucleosome assembly and translation. CONCLUSIONS: The transcriptomic profile of L. infantum showed a robust set of genes from different metabolic pathways associated with the antimony resistance phenotype in this parasite. Our results address the complex and multifactorial antimony resistance mechanisms in Leishmania, identifying several candidate genes that may be further evaluated as molecular targets for chemotherapy of leishmaniasis.

Nuclear and mitochondrial genome sequencing of North-African Leishmania infantum isolates from cured and relapsed visceral leishmaniasis patients reveals variations correlating with geography and phenotype.

Although several studies have investigated genetic diversity of Leishmania infantum in North Africa, genome-wide analyses are lacking. Here, we conducted comparative analyses of nuclear and mitochondrial genomes of seven L. infantum isolates from Tunisia with the aim to gain insight into factors that drive genomic and phenotypic adaptation. Isolates were from cured (n=4) and recurrent (n=3) visceral leishmaniasis (VL) cases, originating from northern (n=2) and central (n=5) Tunisia, where respectively stable and emerging VL foci are observed. All isolates from relapsed patients were from Kairouan governorate (Centre); one showing resistance to the anti-leishmanial drug Meglumine antimoniate. Nuclear genome diversity of the isolates was analysed by comparison to the L. infantum JPCM5 reference genome. Kinetoplast maxi and minicircle sequences (1 and 59, respectively) were extracted from unmapped reads and identified by blast analysis against public data sets. The genome variation analysis grouped together isolates from the same geographical origins. Strains from the North were very different from the reference showing more than 34 587 specific single nucleotide variants, with one isolate representing a full genetic hybrid as judged by variant frequency. Composition of minicircle classes within isolates corroborated this geographical population structure. Read depth analysis revealed several significant gene copy number variations correlating with either geographical origin (amastin and Hsp33 genes) or relapse (CLN3 gene). However, no specific gene copy number variation was found in the drug-resistant isolate. In contrast, resistance was associated with a specific minicircle pattern suggesting Leishmania mitochondrial DNA as a potential novel source for biomarker discovery.

Leishmania amazonensis Subverts the Transcription Factor Landscape in Dendritic Cells to Avoid Inflammasome Activation and Stall Maturation.

Leishmania parasites are the causative agents of human leishmaniases. They infect professional phagocytes of their mammalian hosts, including dendritic cells (DCs) that are essential for the initiation of adaptive immune responses. These immune functions strictly depend on the DC's capacity to differentiate from immature, antigen-capturing cells to mature, antigen-presenting cells-a process accompanied by profound changes in cellular phenotype and expression profile. Only little is known on how intracellular Leishmania affects this important process and DC transcriptional regulation. Here, we investigate these important open questions analyzing phenotypic, cytokine profile and transcriptomic changes in murine, immature bone marrow-derived DCs (iBMDCs) infected with antibody-opsonized and non-opsonized Leishmania amazonensis (L.am) amastigotes. DCs infected by non-opsonized amastigotes remained phenotypically immature whereas those infected by opsonized parasites displayed a semi-mature phenotype. The low frequency of infected DCs in culture led us to use DsRed2-transgenic parasites allowing for the enrichment of infected BMDCs by FACS. Sorted infected DCs were then subjected to transcriptomic analyses using Affymetrix GeneChip technology. Independent of parasite opsonization, Leishmania infection induced expression of genes related to key DC processes involved in MHC Class I-restricted antigen presentation and alternative NF-κB activation. DCs infected by non-opsonized parasites maintained an immature phenotype and showed a small but significant down-regulation of gene expression related to pro-inflammatory TLR signaling, the canonical NF-kB pathway and the NLRP3 inflammasome. This transcriptomic profile was further enhanced in DCs infected with opsonized parasites that displayed a semi-mature phenotype despite absence of inflammasome activation. This paradoxical DC phenotype represents a Leishmania-specific signature, which to our knowledge has not been observed with other opsonized infectious agents. In conclusion, systems-analyses of our transcriptomics data uncovered important and previously unappreciated changes in the DC transcription factor landscape, thus revealing a novel Leishmania immune subversion strategy directly acting on transcriptional control of gene expression. Our data raise important questions on the dynamic and reciprocal interplay between trans-acting and epigenetic regulators in establishing permissive conditions for intracellular Leishmania infection and polarization of the immune response.