Abstract Leishmania parasites cause cutaneous leishmaniasis (CL), a pathologic disease characterized by disfiguring, ulcerative skin lesions. Both parasite and host gene expression following infection with various Leishmania species has been investigated in vitro, but global transcriptional analysis following L. major infection in vivo is lacking. Thus, we conducted a comprehensive transcriptomic profiling study combining bulk RNA sequencing (RNA-Seq) and single-cell RNA sequencing (scRNA-Seq) to identify global changes in gene expression in vivo following L. major infection. Bulk RNA-Seq analysis revealed that host immune response pathways like the antigen processing and presentation pathway were significantly enriched amongst differentially expressed genes (DEGs) upon infection, while ribosomal pathways were significantly downregulated in infected mice compared to naive controls. scRNA-Seq analyses revealed cellular heterogeneity including distinct resident and recruited cell types in the skin following murine L. major infection. Within the individual immune cell types, several DEGs indicative of many interferon induced GTPases and antigen presentation molecules were significantly enhanced in the infected ears including macrophages ( Gbp2, H2-K1, H2-Aa, H2-Ab1 ), resident macrophages ( H2-K1, H2-D1, Gbp4, Gbp8, Gbp2 ), and inflammatory monocytes ( Gbp2, Gbp5, Gbp7, Gbp3 ). Ingenuity Pathway Analysis of scRNA-Seq data indicated the antigen presentation pathway was increased with infection, while EIF2 signaling is the top downregulated pathway followed by eIF4/p70S6k and mTOR signaling in multiple cell types including macrophages, BECs, and LECs. Altogether, this transcriptomic profile highlights known recruitment of myeloid cells to lesions and recognizes a previously undefined role for EIF2 signaling in murine L. major infection in vivo. Author summary Leishmania major cause cutaneous leishmaniasis, which is characterized by disfiguring, ulcerative skin lesions. Here, we show murine L. major -directed reprogramming of the host transcriptome in vivo. Our bulk RNA-Seq analyses revealed upregulation of antigen processing and presentation pathway, while the host ribosomal pathway was downregulated following L. major infection. Similarly, scRNA-Seq analyses revealed the upregulation of transcripts responsible for antigen presentation and host defense proteins like guanylate binding proteins (GBPs) alongside the downregulation of EIF2 signalling at the site of L. major infection. Overall, our transcriptomic dataset not only provides the comprehensive list of gene expression at the single-cell resolution, and highlights a previously undefined role for EIF2 signalling during murine L. major infection in vivo.

In vivo reprogramming of murine host immune response genes following Leishmania major infection