Comparative transcriptome analysis of Peromyscus leucopus and C3H mice infected with the Lyme disease pathogen

Comparative reservoir competence of Peromyscus leucopus, C57BL/6J, and C3H/HeN for Borrelia burgdorferi B31

Borrelia burgdorferi, a Lyme disease spirochete, causes a range of acute and chronic maladies in humans. However, a primary vertebrate reservoir in the United States, the white-footed deermouse Peromyscus leucopus, is reported not to have reduced fitness following infection. Although laboratory strains of Mus musculus mice have successfully been leveraged to model acute human Lyme disease, the ability of these rodents to model B. burgdorferi-P. leucopus interactions remains understudied. Here, we compared infection of P. leucopus with B. burgdorferi B31 with infection of the traditional B. burgdorferi murine models-C57BL/6J and C3H/HeN Mus musculus, which develop signs of inflammation akin to human disease. We find that B. burgdorferi was able to reach much higher burdens (10- to 30-times higher) in multiple M. musculus skin sites and that the overall dynamics of infection differed between the two rodent species. We also found that P. leucopus remained transmissive to larval Ixodes scapularis for a far shorter period than either M. musculus strain. In line with these observations, we found that P. leucopus does launch a modest but sustained inflammatory response against B. burgdorferi in the skin, which we hypothesize leads to reduced bacterial viability and rodent-to-tick transmission in these hosts. Similarly, we also observe evidence of inflammation in infected P. leucopus hearts. These observations provide new insight into reservoir species and the B. burgdorferi enzootic cycle.IMPORTANCEA Lyme disease-causing bacteria, Borrelia burgdorferi, must alternate between infecting a vertebrate host-usually rodents or birds-and ticks. In order to be successful in that endeavor, the bacteria must avoid being killed by the vertebrate host before it can infect a new larval tick. In this work, we examine how B. burgdorferi and one of its primary vertebrate reservoirs, Peromyscus leucopus, interact during an experimental infection. We find that B. burgdorferi appears to colonize its natural host less successfully than conventional laboratory mouse models, which aligns with a sustained seemingly anti-bacterial response by P. leucopus against the microbe. These data enhance our understanding of P. leucopus host-pathogen interactions and could potentially serve as a foundation to uncover ways to disrupt the spread of B. burgdorferi in nature.

Estimating Enzyme Expression and Metabolic Pathway Activity in Borreliella-Infected and Uninfected Mice

Evaluating changes in metabolic pathway activity is essential for studying disease mechanisms and developing new treatments, with significant benefits extending to human health. Here, we propose EMPathways2, a maximum likelihood pipeline that is based on the expectation-maximization algorithm, which is capable of evaluating enzyme expression and metabolic pathway activity level. We first estimate enzyme expression from RNA-seq data that is used for simultaneous estimation of pathway activity levels using enzyme participation levels in each pathway. We implement the novel pipeline to RNA-seq data from several groups of mice, which provides a deeper look at the biochemical changes occurring as a result of bacterial infection, disease, and immune response. Our results show that estimated enzyme expression, pathway activity levels, and enzyme participation levels in each pathway are robust and stable across all samples. Estimated activity levels of a significant number of metabolic pathways strongly correlate with the infected and uninfected status of the respective rodent types.

Cellular and transcriptome signatures unveiled by single-cell RNA-Seq following ex vivo infection of murine splenocytes with Borrelia burgdorferi

Introduction

Lyme disease, the most common tick-borne infectious disease in the US, is caused by a spirochetal pathogen Borrelia burgdorferi (Bb). Distinct host responses are observed in susceptible and resistant strains of inbred of mice following infection with Bb reflecting a subset of inflammatory responses observed in human Lyme disease. The advent of post-genomic methodologies and genomic data sets enables dissecting the host responses to advance therapeutic options for limiting the pathogen transmission and/or treatment of Lyme disease.

Methods

In this study, we used single-cell RNA-Seq analysis in conjunction with mouse genomics exploiting GFP-expressing Bb to sort GFP+ splenocytes and GFP- bystander cells to uncover novel molecular and cellular signatures that contribute to early stages of immune responses against Bb.

Results

These data decoded the heterogeneity of splenic neutrophils, macrophages, NK cells, B cells, and T cells in C3H/HeN mice in response to Bb infection. Increased mRNA abundance of apoptosis-related genes was observed in neutrophils and macrophages clustered from GFP+ splenocytes. Moreover, complement-mediated phagocytosis-related genes such as C1q and Ficolin were elevated in an inflammatory macrophage subset, suggesting upregulation of these genes during the interaction of macrophages with Bb-infected neutrophils. In addition, the role of DUSP1 in regulating the expression of Casp3 and pro-inflammatory cytokines Cxcl1, Cxcl2, Il1b, and Ccl5 in Bb-infected neutrophils were identified.

Discussion

These findings serve as a growing catalog of cell phenotypes/biomarkers among murine splenocytes that can be exploited for limiting spirochetal burden to limit the transmission of the agent of Lyme disease to humans via reservoir hosts.