Modulation of Host Immune Response during Leishmania infantum Natural Infection: A Whole-Transcriptome Analysis of the Popliteal Lymph Nodes in Dogs

Exploring the impact of epidemiological and clinical factors on the progression of canine leishmaniosis by statistical and whole genome analyses: from breed predisposition to comorbidities

Canine leishmaniosis (CanL), caused by Leishmania infantum, is a complex disease of growing importance in Europe. Clinical manifestations result from the down-modulation of the host immune response through multiple host-parasite interactions. Although several factors might influence CanL progression, this is the first known study evaluating risk factors for its different clinical stages in a large referral hospital population (n = 35.669) from an endemic area, over a 20 year period. Genome-wide scans for selection signatures were also conducted to explore the genomic component of clinical susceptibility to L. infantum infection. The prevalence of CanL was 3.2% (16.7% stage I; 43.6% stage II; 32.1% stage III; 7.6% stage IV). Dog breed (crossbreed), bodyweight (<10 kg), living conditions (indoors), regular deworming treatment, and being vaccinated against Leishmania significantly decreased the transmission risk and the risk for developing severe clinical forms. Conversely, the detection of comorbidities was associated with advanced clinical forms, particularly chronic kidney disease, neoplasia, cryptorchidism, infectious tracheobronchitis and urate urolithiasis, although those did not impact the clinical outcome. Significant associations between an increased risk of severe clinical stages and findings in the anamnesis (renal or skin-related manifestations) and physical examination (ocular findings) were also detected, highlighting their diagnostic value in referred cases of CanL. Sixteen breeds were found to be significantly more susceptible to developing severe stages of leishmaniosis (e.g. Great Dane, Rottweiler, English Springer Spaniel, Boxer, American Staffordshire Terrier, Golden Retriever), while 20 breeds displayed a clinical resistantance phenotype and, thus, are more likely to mount an efficient immune response against L. infantum (e.g. Pointer, Samoyed, Spanish Mastiff, Spanish Greyhound, English Setter, Siberian Husky). Genomic analyses of these breeds retrieved 12 regions under selection, 63 candidate genes and pinpointed multiple biological pathways such as the IRE1 branch of the unfolded protein response, which could play a critical role in clinical susceptibility to L. infantum infection.

Transcriptional profiling of canine osteosarcoma identifies prognostic gene expression signatures with translational value for humans

Canine osteosarcoma is increasingly recognized as an informative model for human osteosarcoma. Here we show in one of the largest clinically annotated canine osteosarcoma transcriptional datasets that two previously reported, as well as de novo gene signatures devised through single sample Gene Set Enrichment Analysis (ssGSEA), have prognostic utility in both human and canine patients. Shared molecular pathway alterations are seen in immune cell signaling and activation including TH1 and TH2 signaling, interferon signaling, and inflammatory responses. Virtual cell sorting to estimate immune cell populations within canine and human tumors showed similar trends, predominantly for macrophages and CD8+ T cells. Immunohistochemical staining verified the increased presence of immune cells in tumors exhibiting immune gene enrichment. Collectively these findings further validate naturally occurring osteosarcoma of the pet dog as a translationally relevant patient model for humans and improve our understanding of the immunologic and genomic landscape of the disease in both species.

The role of atypical MAP kinase 4 in the host interaction with Cryptosporidium parvum

Cryptosporidium parvum is an apicomplexan parasite that causes severe zoonotic diarrhea in humans and calves. Since there are no effective treatments or vaccines for infants or immunocompromised patients, it is important to understand the molecular mechanisms of the parasite-host interaction for novel drug discovery. Mitogen-activated protein kinase (MAP kinase) is a key host factor in interactions between host and various pathogens, including parasites. Although the function of conventional MAP kinases against parasite infection has been investigated, that of atypical MAP kinases remains largely unknown. Therefore, we focused on one of the atypical MAP kinases, MAPK4, and its effect on C. parvum infection in human intestinal cells. Here, we report that MAPK4-deficient intestinal cells showed a significant reduction in C. parvum infection. We also show that host MAPK4 has a role in host cell survival from C. parvum infection. In addition, we show that C. parvum requires host MAPK4 for its successful invasion and asexual reproduction. Taken together, our data suggest that MAPK4 is an important host factor contributing to C. parvum infection in human intestinal cells.

Comparative transcriptomic analysis of long noncoding RNAs in Leishmania-infected human macrophages

It is well established that infection with Leishmania alters the host cell's transcriptome. Since mammalian cells have multiple mechanisms to control gene expression, different molecules, such as noncoding RNAs, can be involved in this process. MicroRNAs have been extensively studied upon Leishmania infection, but whether long noncoding RNAs (lncRNAs) are also altered in macrophages is still unexplored. We performed RNA-seq from THP-1-derived macrophages infected with Leishmania amazonensis (La), L. braziliensis (Lb), and L. infantum (Li), investigating a previously unappreciated fraction of macrophage transcriptome. We found that more than 24% of the total annotated transcripts and 30% of differentially expressed (DE) RNAs in Leishmania-infected macrophage correspond to lncRNAs. LncRNAs and protein coding RNAs with altered expression are similar among macrophages infected with the Leishmania species. Still, some species-specific alterations could occur due to distinct pathophysiology in which Li infection led to a more significant number of exclusively DE RNAs. The most represented classes among DE lncRNAs were intergenic and antisense lncRNAs. We also found enrichment for immune response-related pathways in the DE protein coding RNAs, as well as putative targets of the lncRNAs. We performed a coexpression analysis to explore potential cis regulation of coding and antisense noncoding transcripts. We identified that antisense lncRNAs are similarly regulated as its neighbor protein coding genes, such as the BAALC/BAALC-AS1, BAALC/BAALC-AS2, HIF1A/HIF1A-AS1, HIF1A/HIF1A-AS3 and IRF1/IRF1-AS1 pairs, which can occur as a species-specific modulation. These findings are a novelty in the field because, to date, no study has focused on analyzing lncRNAs in Leishmania-infected macrophage. Our results suggest that lncRNAs may account for a novel mechanism by which Leishmania can control macrophage function. Further research must validate putative lncRNA targets and provide additional prospects in lncRNA function during Leishmania infection.

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.